Review Article  Open Access
Versatility of Approximating SingleParticle Electron Microscopy Density Maps Using Pseudoatoms and ApproximationAccuracy Control
Abstract
Threedimensional Gaussian functions have been shown useful in representing electron microscopy (EM) density maps for studying macromolecular structure and dynamics. Methods that require setting a desired number of Gaussian functions or a maximum number of iterations may result in suboptimal representations of the structure. An alternative is to set a desired error of approximation of the given EM map and then optimize the number of Gaussian functions to achieve this approximation error. In this article, we review different applications of such an approach that uses spherical Gaussian functions of fixed standard deviation, referred to as pseudoatoms. Some of these applications use EMmap normal mode analysis (NMA) with elastic network model (ENM) (applications such as predicting conformational changes of macromolecular complexes or exploring actual conformational changes by normalmodebased analysis of experimental data) while some other do not use NMA (denoising of EM density maps). In applications based on NMA and ENM, the advantage of using pseudoatoms in EMmap coarsegrain models is that the ENM springs are easily assigned among neighboring grains thanks to their spherical shape and uniformed size. EMmap denoising based on the map coarsegraining was so far only shown using pseudoatoms as grains.
1. Introduction
Singleparticle analysis is an electron microscopy (EM) technique that allows determining the structure at nearatomic resolutions for a large range of macromolecular complexes [1–16]. Also, it allows studying conformational variability of macromolecular complexes by determining their different conformations [17–22]. These different conformations are usually obtained by analyzing heterogeneity with methods that assume a small number of discrete conformations coexisting in the specimen [23–28], while several methods have been recently developed to help analyzing continuous conformational changes [29–33].
EMmap representations with a reduced number of points or with a set of 3D Gaussian functions have been shown useful in studying macromolecular structure and dynamics [30, 33–44]. The process of representing EM maps with a set of points or 3D Gaussian functions (grains) is sometimes referred to as coarsegraining of EM maps. A typical approach to coarsegraining is a neural network clustering approach that quantizes the given EM map so that the probability density of the grains closely resembles the probability density of the given data, which makes the coarsegrain representation retain the overall shape of the structure from the given EM map [34, 36, 38–40]. This approach is referred to as Vector Quantization (VQ). A different approach is to parametrize a Gaussian Mixture Model (GMM) of the probability density function using expectationmaximization algorithm [41, 45]. All these approaches require setting a desired (target) number of grains or a maximum number of iterations to stop the iterative procedure, which may result in suboptimal representations. Indeed, the use of a small target number of grains or a small maximum number of iterations may lead to a small final number of grains resulting in a model with overrepresented high density regions and underrepresented low density regions. Furthermore, in the case of symmetrical structures, the inadequately small final number of grains can result in representations that are overall nonuniform (asymmetrical). A difficulty is thus to choose the stopping parameter that will produce a sufficiently high number of grains to appropriately represent all density regions.
An alternative is to set a desired (target) error of approximation of the given EM map and then optimize the number of Gaussian functions, their position, and their weights to achieve the target approximation error, as in the approach that we introduced in [43]. In each iteration, this approach adds some Gaussian functions (grains) while removing some (the grains with small weights or distances will be removed). We have found that this strategy of minimizing the global representation error, involving controlled adding and removing grains, allows placing new grains where they are most needed and adapting the grains near the removed ones to better represent the local intensity in the input EM map [43], which helps overcoming the underrepresentation problem. For instance, we have found that symmetry is preserved in EMmap approximations with this strategy for typical values of the target approximation error such as 1–15% [30, 33, 42–44]. This method uses spherical Gaussian functions of fixed standard deviation that we refer to as pseudoatoms. Its versatility has been shown in applications such as predicting conformational changes of macromolecular complexes, exploring actual conformational changes, analyzing continuous conformational changes, and denoising of EM density maps [30, 33, 42–44]. Some of these applications are based on EMmap normal mode analysis (NMA) with elastic network model (ENM) [46, 47] (e.g., predicting conformational changes of macromolecular complexes or exploring actual conformational changes using normalmodebased analysis of experimental data). In some other applications, NMA is not used (e.g., denoising of EM density maps).
The advantage of using pseudoatoms in applications based on NMA and ENM, with respect to other types of grains (Table 1), is their uniformity over the molecule that is a prerequisite for a simple application of the ENM. Indeed, as pseudoatoms have spherical shape and uniformed size over the molecule, they allow an easy setting of springs among neighboring grains in the ENM. On the contrary, 3D points (the socalled codebook vectors) obtained with VQ can be regarded as Gaussian functions whose standard deviation can vary over the molecule (each codebook vector is associated with a data subregion known as Voronoi cell [34] and different subregions can have different sizes), which may make the ENM setting more complicated than with Gaussian functions of the same standard deviation. Similarly, NMA should be more complicated with ellipsoidal Gaussian distribution obtained with the GMM approach (to the best of our knowledge, such NMA has not been reported so far). Regarding EMmap denoising (not based on NMA), only pseudoatoms were so far reported as grains in that application of coarsegraining [44].

In this article, we review the mentioned applications of this EMmap approximation method while only briefly reminding the method. For algorithmic details (e.g., related to adding/removing grains), the reader is addressed to [43] that describes this method in detail.
2. Background
We start this section with a brief background on the approach for converting EM maps in sets of pseudoatoms. Then, we provide a brief reminder on NMA that is used in applications of the conversion approach to studying conformational changes of macromolecular complexes.
2.1. Approximation of EM Maps Using Pseudoatoms and Approximation Error Control
A function can be approximated using Gaussian radial basis functions (RBFs) by , where is the RBF kernel that is a Gaussian function with the standard deviation and the amplitude of 1, to which we refer as pseudoatom, is the number of pseudoatoms, is the vector of the center coordinates of the th pseudoatom, is the Euclidean distance between the vectors and , and is the weight (contribution) of the th pseudoatom. Given an EM density map , a Gaussianfunction standard deviation , and a target approximation error , our approach determines the number of pseudoatoms , their positions , and weights such that the approximation error satisfies . Here, is the effective range of values in the EM map, is the voxel location at which the given EM map is compared with its approximation, and is the total number of evaluated voxels (the evaluation can be done in a region of interest defined by a mask). To avoid getting trapped into local minima of the error , new pseudoatoms are added progressively in regions with large errors and weights and positions of the current number of pseudoatoms are determined by a gradient descent minimization of . It should be noted that pseudoatom positions do not necessarily coincide with voxel positions because the pseudoatom positions vary continuously within the EM density map. Also, it should be noted that can be expressed in angstroms, but it is here expressed in voxels. The typical values of are between 1 and 2 (voxels).
Different pseudoatom representations can be obtained with this approach by varying the Gaussianfunction standard deviation () and the target approximation error () (Figure 1). Smaller values of and result in larger numbers of pseudoatoms and vice versa. However, and should be chosen taking into account resolution and noise of the given EM map (e.g., large and may give sufficiently good approximations of low resolution maps) and the target application, as we show here.
(a)
(b)
(c)
(d)
2.2. Normal Mode Analysis
NMA models complex motions by linear combinations of harmonic oscillations around a minimumenergy conformation. NMA is often based on the standard ENM [46], which is also the case in our approach. More precisely, we perform NMA of pseudoatom representations of EM maps using the software developed by Tama et al. [47]. In this approach, the given EM map is assumed to contain the structure in the minimumenergy conformation and no energy minimization is required [46]. It has been shown that NMA of EM density maps of intermediate resolution result in a good approximation of normal modes of atomicresolution structures, in particular at low frequencies at which the motions were experimentally observed [47]. An experimentally observed motion here means a transition between two different conformations of the same complex obtained at atomic resolution (e.g., by Xray crystallography). This transition and the conformational states along it can be modeled (simulated) using lowfrequency normal modes of any of the two given conformations. For instance, the lowfrequency modes having the highest overlap with the difference between the two given atomicresolution conformations usually capture 60–70% of the conformational change [47, 48]. Normal modes of EM density maps can be used to obtain other possible conformational states of the same complex (to predict lowfrequency motions of the complex). Note here that the conformational states are discrete and can be regarded as discrete samples of a continuous trajectory of conformational transition. NMA of intermediateresolution EM density maps is especially useful when atomicresolution structures cannot be obtained [37, 38, 47], but it requires coarsegrain representations of the density maps.
Nodes of the elastic network model are 3D point particles. Each node is connected, via harmonic springs, with other nodes within a sphere of a given radius (the radius is referred to as interaction cutoff distance). In our approach, the coordinates of nodes of the ENM are the center coordinates of pseudoatoms with which the given EM density map is represented. Given nodes of the ENM, NMA requires a diagonalization of a matrix of second derivatives of the potential energy (Hessian matrix), which is performed via eigenanalysis of the Hessian matrix. Normal modes are eigenvectors of the Hessian matrix while the normal mode frequencies are the square roots of eigenvalues of the Hessian matrix.
A displacement of nodes of the ENM along normal modes modifies the given conformation, which is used in simulations of structural flexibility (Figure 2). The coordinates of normal modes are expressed in angstroms while the displacement amplitudes along normal modes have no units. Six lowestfrequency modes are related to rigidbody movements and are usually not used for the displacement.
(a)
(b)
(c)
(d)
3. Applications
3.1. Prediction of Conformational Changes
As explained in Background, NMA of intermediateresolution EM maps provides normal modes that can be used to simulate (predict) different conformational states of the same complex. Normal modes of EM maps are calculated based on the EMmap coarsegrained representations with pseudoatoms. Both EMmap coarsegraining and NMA can be performed using 3DEM Loupe web server [42].
3DEM Loupe [42] is currently the only web server that allows interactive NMA of EM density maps. Other web servers perform NMA of atomicresolution structures (e.g., ElNemo [49]) or do not allow the user to input his/her own EM density map for NMA (e.g., EMotion [38]). The workflow of 3DEM Loupe consists of the following three steps: (1) conversion of the input EM density map into a pseudoatomic structure (coarsegraining) based on the approximation error control; (2) calculation of normal modes of the pseudoatomic structure; and (3) animation of the obtained normal modes (animated displacement of pseudoatomic structures along normal modes). The pseudoatomic structure, normal modes, and animations obtained by 3DEM Loupe can be downloaded for their further analysis on a local computer. The server also contains precomputed results for several EM maps, such as 70S ribosome, GroEL, ribosomebound termination factor RF2, and connector of bacteriophage T7.
3.2. Exploring Actual Conformational Changes Including Continuous Conformational Changes
Coarsegrain representations of EM density maps and normal modes of those coarsegrain representations can be used to analyze experimental EM data. More precisely, conformations actually present in EM data can be interpreted using simulated conformations (NMAbased simulations). In the framework of conformational heterogeneity analysis, normal modes were used for image analysis first in [50], but the first procedure capable of automatically processing large series of images using as many “test” conformations as needed is HEMNMA [30, 51]. In this subsection, HEMNMA, the approach for NMAbased conformational analysis of a series of images (using normal modes of a given density map), is presented together with an approach for NMAbased conformational analysis of a series of density maps (StructMap).
3.2.1. Analysis of EM Images: HEMNMA
Given an EM density map, HEMNMA [30, 51] computes its pseudoatom representation and performs NMA of the obtained pseudoatomic structure, as does 3DEM Loupe [42]. This pseudoatomic structure and its normal modes are then used to analyze each singleparticle image by elastic projection matching. More precisely, HEMNMA deforms the pseudoatomic structure using some amplitudes of the displacement along normal modes and then compares the image with projections of the pseudoatomic structure converted into a density map. The image is assigned the orientation, translation, and normal mode amplitudes of the bestmatching projection. The normal mode amplitudes obtained for all images are finally mapped onto a lowdimensional distance space (usually, 1D, 2D, or 3D), which allows seeing the distribution of conformations. This space can be explored to detect whether the conformational change is discrete (points are grouped in a few clusters) or continuous (absence of clusters, points spread over the space). In the case of continuous conformational changes, conformational trajectories may be identified by analyzing the most populated regions in this space. Along the trajectories, the given EM map and the obtained pseudoatomic structure can be animated and 3D reconstructions can be calculated. It has been shown that this methodology can help detect the motions undetectable with methods that assume a small number of different coexisting conformations in the specimen [30] (Figure 3). In this context, the conformational heterogeneity of E. coli 70S ribosome, DNA polymerase Pol αB subunit complex of the eukaryotic primosome, and Tomato Bushy Stunt Virus was described more extensively with HEMNMA than with other methods [30].
(a)
(b)
(c)
It should be noted that HEMNMA can also analyze a series of images using an atomicresolution structure instead of a density map. It should also be noted that the conformational variability shown in Figure 3(b) is rarely seen in EM because the majority of computational methods assume specimens with few different conformations of complexes (e.g., Figures 3(a) and 4(a)). Additionally, biochemical specimen preparation protocols are usually optimized so as to reduce the number of different conformations in the specimen. To achieve atomic resolution of reconstruction, the computational methods are usually based on different rounds of 2D and 3D classifications (involving computations of average particle views and conformations), during which many particle images are removed and only particles assigned to highresolution classes (particles with most consistent views and conformations) are kept for final 3D reconstructions. Such biochemical and computational “selection” of particles may obscure full conformational variability of the complex. On the contrary, HEMNMA method estimates the conformational variability distribution (e.g., Figure 3(b)) using raw images and no classification, which facilitates deciphering the full range of conformational variability that is a quasicontinuum of conformational states (a large number of discrete samples of a continuous conformational transition trajectory). Analysis of continuous conformational transitions by EM is currently a research field in expansion [29–32, 52], which will be reviewed in a separate publication.
(a)
(b)
3.2.2. Analysis of EM Maps: StructMap
StructMap [33] automatically analyzes a set of EM density maps to map them onto a common lowdimensional distance space (usually, 1D, 2D, or 3D). It requires representing each density map with pseudoatoms and computing normal modes of each obtained pseudoatomic structure. The actual conformation in each EM density map (reference map) is estimated by elastic fitting of this map with the pseudoatomic structures of all other maps from the given set of EM maps. More precisely, pseudoatoms are displaced along normal modes and the displacement amplitudes are adjusted to maximize the crosscorrelation coefficient between the map calculated from displaced pseudoatoms and the reference map. The obtained maximum crosscorrelation coefficients are subtracted from 1 and then mapped onto a lowdimensional space, which shows the differences (distances) among given EM density maps that cannot be interpreted using normal modes. In the case of continuous conformational distribution, EM density maps can be regarded as discrete, unordered samples of continuous trajectories and StructMap can be used to get a rough idea of these trajectories (potential sequences of conformational changes) [33]. In some heterogeneity cases, such as a combined conformational and compositional heterogeneity of states along the elongation cycle of human 80S ribosome [53], this task can be particularity challenging, but StructMap can still help visualizing the differences between the states and analyzing potential trajectories [33] (Figure 4).
In the context of revealing sequences of conformational changes, StructMap can also be combined with HEMNMA. For instance, instead of using HEMNMA to analyze all images with all given EM density maps, a few EM maps can be chosen from the distance space obtained by StructMap so as to explore more finely particular regions in this space.
3.3. Denoising
CryoEM images have a low signaltonoise ratio. Thus, a large number of such images must be averaged in 3D space to reduce noise. However, even highresolution 3D reconstructed EM maps may still contain significant amounts of noise. We have shown that pseudoatom representations of EM density maps can be used for EMmap denoising [44]. We assume that the object reconstructed from images is correct enough, meaning that effects of noise dominate potential reconstruction artifacts and potential effects of heterogeneity of images that were used for this reconstruction. Smaller target approximation errors () and smaller Gaussian standard deviations () generally produce better approximations of EM maps (i.e., the density maps computed from pseudoatomic structures have more details), with a risk of reproducing the EMmap noise in the density map computed from pseudoatoms. Setting or to higher values can reduce this risk. These two parameters should be set in accordance with the resolution of the given EM map. We have shown that (voxels) and produce satisfactory denoising results in most cases of higherresolution density maps (resolutions higher than 6 Å, goldstandard FSC 0.143), whereas denoising of lowerresolution density maps is usually achieved using larger values of ( = 5%–15%) and the values of that may need to be adjusted around its default value () [44].
The method was used to denoise EM density maps of several complexes obtained at subnanometer resolutions by singleparticle analysis or subtomogram averaging (betagalactosidase, ribosome, and empty and full virus particles) [44] (Figure 5).
(a)
(b)
(c)
(d)
4. Discussion
The method for EMmap approximation using Gaussian functions of standard deviation and amplitude 1 and using a target approximation error [43] has applications in predicting conformational changes of macromolecular complexes, exploring actual (discrete or continuous) conformational changes, and denoising of EM density maps. In this article, we reviewed these applications, together with the software that uses this method and the examples of results that can be obtained with this software.
In this approach, the EMmap approximation is a weighted sum of Gaussian functions of given standard deviation and amplitude 1 (referred to as pseudoatoms) whose number, location, and weights are determined by minimizing the EMmap approximation error towards the given target error . The shape and density distribution of a given macromolecular complex are thus fully represented with the distribution of Gaussian functions. The centers of pseudoatoms are control points that can be displaced (e.g., along normal modes) to modify a given conformation of the complex. The modified conformation, represented with the pseudoatoms, can be converted into a density map by computing the weighted sum of Gaussian functions at new locations (the same weight as before the modification).
The target approximation error control allows different applications of this method. We could learn from experiments how to choose this target approximation error to suit these different applications. For instance, in EMmap denoising applications, target approximation errors smaller than are recommended for EM maps of higher resolution (higher than 6 Å according to the goldstandard FSC 0.143 criterion), whereas larger target approximation errors (5%–15%) are recommended for EM maps of lower resolution [44]. In applications such as simulations of conformational changes (NMA of EM density maps) or elastic fitting among several EM maps (elastic 3Dto3D fitting), it is recommended to use target errors of 10%–15% with EM maps that are noisy or that have different (low and high) resolutions among each other, whereas target errors of 1%–5% are recommended with clean and highresolution EM maps [33, 42]. In elastic 3Dto2D fitting applications, smaller target errors (1%–5%) are recommended to guarantee a good quality of projections of the density map from pseudoatoms (a good projection quality is required for an accurate projection matching with images) [30]. In these different applications, the standard deviation of Gaussian functions is usually in the range between 1 and 2 (voxels) and the value of 1.5 (voxels) can often be chosen as a default value, though adjustments are required around this value to achieve optimal results (quality can be evaluated by comparing target and achieved approximation errors).
Advantages of structural representations with pseudoatoms have recently been explored in the context of 3D reconstruction. For instance, a Bayesian approach to ab initio, low resolution 3D reconstruction proposed in [54] is based on estimating a pseudoatomic model using a random model for initialization and a given number of pseudoatoms. It should be noticed that the Gaussianfunction standard deviation in that approach is estimated together with other parameters of the pseudoatomic model (position and weights of pseudoatoms) as well as the parameters of image orientation and translation [54]. That approach makes calculations using class averages instead of raw images and it does not consider a possible conformational heterogeneity of the given set of images. We are currently working on a method based on pseudoatoms with a fixed Gaussianfunction standard deviation to iteratively refine a preliminary model taking into account conformational heterogeneity of the data set. This new method should result in a 3D reconstruction (model) optimally representing all given (raw) images via the determined elastic (normalmodebased) transformations between the model and the images.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This research was partially supported by the CNRS (France) and the CSIC (Spain) [Projet International de Coopération Scientifique, PICS 2011]; the French National Research Agency ANR [ANR11BSV801004]; the European Social Fund and the Ministerio de Educación y Ciencia [“Ramón y Cajal” fellowship to COSS]; the Spanish Ministry of Economy and Competitiveness [AIC–A–2011–0638 and BIO201344647R]; and the Comunidad de Madrid [CAM S2010/BMD2305]. HPC resources from GENCI (France) [Grants to project 072174 (2011–2016)] are also gratefully acknowledged.
References
 X. Yu, L. Jin, and Z. H. Zhou, “3.88 Å structure of cytoplasmic polyhedrosis virus by cryoelectron microscopy,” Nature, vol. 453, pp. 415–419, 2008. View at: Publisher Site  Google Scholar
 X. Zhang, E. Settembre, C. Xu et al., “Nearatomic resolution using electron cryomicroscopy and singleparticle reconstruction,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 6, pp. 1867–1872, 2008. View at: Publisher Site  Google Scholar
 M. Liao, E. Cao, D. Julius, and Y. Cheng, “Structure of the TRPV1 ion channel determined by electron cryomicroscopy,” Nature, vol. 504, no. 7478, pp. 107–112, 2013. View at: Publisher Site  Google Scholar
 M. Allegretti, D. J. Mills, G. McMullan, W. Kühlbrandt, and J. Vonck, “Atomic model of the F420reducing [NiFe] hydrogenase by electron cryomicroscopy using a direct electron detector,” eLife, vol. 2014, no. 3, Article ID e01963, 2014. View at: Publisher Site  Google Scholar
 A. Bartesaghi, D. Matthies, S. Banerjee, A. Merk, and S. Subramaniam, “Structure of βgalactosidase at 3.2Å resolution obtained by cryoelectron microscopy,” Proceedings of the National Academy of Sciences of the United States of America, vol. 111, no. 32, pp. 11709–11714, 2014. View at: Publisher Site  Google Scholar
 X.C. Bai, C. Yan, G. Yang et al., “An atomic structure of human γsecretase,” Nature, vol. 525, no. 7568, pp. 212–217, 2015. View at: Publisher Site  Google Scholar
 A. Bartesaghi, A. Merk, S. Banerjee et al., “2.2 å resolution cryoEM structure of βgalactosidase in complex with a cellpermeant inhibitor,” Science, vol. 348, no. 6239, pp. 1147–1151, 2015. View at: Publisher Site  Google Scholar
 N. Fischer, P. Neumann, A. L. Konevega et al., “Structure of the E. coli ribosome—EFTu complex at <3 Å resolution by C_{s}corrected cryoEM,” Nature, vol. 520, no. 7548, pp. 567–570, 2015. View at: Publisher Site  Google Scholar
 H. Khatter, A. G. Myasnikov, S. K. Natchiar, and B. P. Klaholz, “Structure of the human 80S ribosome,” Nature, vol. 520, no. 7549, pp. 640–645, 2015. View at: Publisher Site  Google Scholar
 B. Liang, Z. Li, S. Jenni et al., “Structure of the L protein of vesicular stomatitis virus from electron cryomicroscopy,” Cell, vol. 162, no. 2, pp. 314–327, 2015. View at: Publisher Site  Google Scholar
 Y. Pang, X. C. Bai, C. Yan et al., “Structure of the apoptosome: mechanistic insights into activation of an initiator caspase from Drosophil,” Genes & Development, vol. 29, pp. 277–287, 2015. View at: Google Scholar
 J. Von Der Ecken, M. Müller, W. Lehman, D. J. Manstein, P. A. Penczek, and S. Raunser, “Structure of the Factintropomyosin complex,” Nature, vol. 519, no. 7541, pp. 114–117, 2015. View at: Publisher Site  Google Scholar
 Z. Yan, X.C. Bai, C. Yan et al., “Structure of the rabbit ryanodine receptor RyR1 at nearatomic resolution,” Nature, vol. 517, no. 7532, pp. 50–55, 2015. View at: Publisher Site  Google Scholar
 M. K. Hospenthal, A. Redzej, K. Dodson et al., “Structure of a chaperoneusher pilus reveals the molecular basis of rod uncoiling,” Cell, vol. 164, no. 12, pp. 269–278, 2016. View at: Publisher Site  Google Scholar
 Z. Liu, F. Guo, F. Wang, T.C. Li, and W. Jiang, “2.9 Å resolution CryoEM 3D reconstruction of closepacked virus particles,” Structure, vol. 24, no. 2, pp. 319–328, 2015. View at: Publisher Site  Google Scholar
 A. C. Walls, M. A. Tortorici, B. J. Bosch et al., “Cryoelectron microscopy structure of a coronavirus spike glycoprotein trimer,” Nature, vol. 531, no. 7592, pp. 114–117, 2016. View at: Publisher Site  Google Scholar
 X.C. Bai, E. Rajendra, G. Yang, Y. Shi, and S. H. W. Scheres, “Sampling the conformational space of the catalytic subunit of human gsecretase,” eLife, vol. 4, no. 2015, Article ID e11182, 2015. View at: Publisher Site  Google Scholar
 A. Zhou, A. Rohou, D. G. Schep et al., “Structure and conformational states of the bovine mitochondrial ATP synthase by cryoEM,” eLife, vol. 4, no. 2015, Article ID e10180, 2015. View at: Publisher Site  Google Scholar
 P. D. Abeyrathne, C. S. Koh, T. Grant, N. Grigorieff, and A. A. Korostelev, “Ensemble cryoEM uncovers inchwormlike translocation of a viral IRES through the ribosome,” Elife, vol. 5, 2016. View at: Google Scholar
 S. Banerjee, A. Bartesaghi, A. Merk et al., “2.3 Å resolution cryoEM structure of human p97 and mechanism of allosteric inhibition,” Science, vol. 351, no. 6275, pp. 871–875, 2016. View at: Publisher Site  Google Scholar
 Y. He, C. Yan, J. Fang et al., “Nearatomic resolution visualization of human transcription promoter opening,” Nature, vol. 533, no. 7603, pp. 359–365, 2016. View at: Publisher Site  Google Scholar
 T. H. D. Nguyen, W. P. Galej, X.C. Bai et al., “CryoEM structure of the yeast U4/U6.U5 trisnRNP at 3.7 Å resolution,” Nature, vol. 530, no. 7590, pp. 298–302, 2016. View at: Publisher Site  Google Scholar
 P. A. Penczek, J. Frank, and C. M. T. Spahn, “A method of focused classification, based on the bootstrap 3D variance analysis, and its application to EFGdependent translocation,” Journal of Structural Biology, vol. 154, no. 2, pp. 184–194, 2006. View at: Publisher Site  Google Scholar
 S. H. W. Scheres, H. Gao, M. Valle et al., “Disentangling conformational states of macromolecules in 3DEM through likelihood optimization,” Nature Methods, vol. 4, no. 1, pp. 27–29, 2007. View at: Publisher Site  Google Scholar
 E. V. Orlova and H. R. Saibil, “Methods for threedimensional reconstruction of heterogeneous assemblies,” Methods in Enzymology, vol. 482, pp. 321–341, 2010. View at: Publisher Site  Google Scholar
 S. H. W. Scheres, “A bayesian view on cryoEM structure determination,” Journal of Molecular Biology, vol. 415, no. 2, pp. 406–418, 2012. View at: Publisher Site  Google Scholar
 D. Lyumkis, A. F. Brilot, D. L. Theobald, and N. Grigorieff, “Likelihoodbased classification of cryoEM images using FREALIGN,” Journal of Structural Biology, vol. 183, no. 3, pp. 377–388, 2013. View at: Publisher Site  Google Scholar
 D. Elmlund and H. Elmlund, “Cryogenic electron microscopy and singleparticle analysis,” Annual Review of Biochemistry, vol. 84, pp. 499–517, 2015. View at: Publisher Site  Google Scholar
 A. Dashti, P. Schwander, R. Langlois et al., “Trajectories of the ribosome as a Brownian nanomachine,” Proceedings of the National Academy of Sciences of the United States of America, vol. 111, no. 49, pp. 17492–17497, 2014. View at: Publisher Site  Google Scholar
 Q. Jin, C. O. S. Sorzano, J. M. de la RosaTrevín et al., “Iterative elastic 3Dto2D alignment method using normal modes for studying structural dynamics of large macromolecular complexes,” Structure, vol. 22, no. 3, pp. 496–506, 2014. View at: Publisher Site  Google Scholar
 E. Katsevich, A. Katsevich, and A. Singer, “Covariance matrix estimation for the cryoEM heterogeneity problem,” SIAM Journal on Imaging Sciences, vol. 8, no. 1, pp. 126–185, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 H. D. Tagare, A. Kucukelbir, F. J. Sigworth, H. Wang, and M. Rao, “Directly reconstructing principal components of heterogeneous particles from cryoEM images,” Journal of Structural Biology, vol. 191, no. 2, pp. 245–262, 2015. View at: Publisher Site  Google Scholar
 C. O. Sanchez Sorzano, A. L. AlvarezCabrera, M. Kazemi, J. M. Carazo, and S. Jonic, “StructMap: elastic distance analysis of electron microscopy maps for studying conformational changes,” Biophysical Journal, vol. 110, no. 8, pp. 1753–1765, 2016. View at: Publisher Site  Google Scholar
 W. Wriggers, R. A. Milligan, K. Schulten, and J. A. McCammon, “Selforganizing neural networks bridge the biomolecular resolution gap,” Journal of Molecular Biology, vol. 284, no. 5, pp. 1247–1254, 1998. View at: Publisher Site  Google Scholar
 J. García de la Torre, O. Llorca, J. L. Carrascosa, and J. M. Valpuesta, “HYDROMIC: prediction of hydrodynamic properties of rigid macromolecular structures obtained from electron microscopy images,” European Biophysics Journal, vol. 30, no. 6, pp. 457–462, 2001. View at: Publisher Site  Google Scholar
 P. A. DeAlarcón, A. PascualMontano, A. Gupta, and J. M. Carazo, “Modeling shape and topology of lowresolution density maps of biological macromolecules,” Biophysical Journal, vol. 83, no. 2, pp. 619–632, 2002. View at: Publisher Site  Google Scholar
 D. Ming, Y. Kong, M. A. Lambert, Z. Huang, and J. Ma, “How to describe protein motion without amino acid sequence and atomic coordinates,” Proceedings of the National Academy of Sciences, vol. 99, no. 13, pp. 8620–8625, 2002. View at: Publisher Site  Google Scholar
 P. Chacón, F. Tama, and W. Wriggers, “MegaDalton biomolecular motion captured from electron microscopy reconstructions,” Journal of Molecular Biology, vol. 326, no. 2, pp. 485–492, 2003. View at: Publisher Site  Google Scholar
 N. JiménezLozano, M. Chagoyen, J. CuencaAlba, and J. M. Carazo, “FEMME database: topologic and geometric information of macromolecules,” Journal of Structural Biology, vol. 144, no. 12, pp. 104–113, 2003. View at: Publisher Site  Google Scholar
 S. Birmanns and W. Wriggers, “Multiresolution anchorpoint registration of biomolecular assemblies and their components,” Journal of Structural Biology, vol. 157, no. 1, pp. 271–280, 2007. View at: Publisher Site  Google Scholar
 T. Kawabata, “Multiple subunit fitting into a lowresolution density map of a macromolecular complex using a Gaussian mixture model,” Biophysical Journal, vol. 95, no. 10, pp. 4643–4658, 2008. View at: Publisher Site  Google Scholar
 R. NogalesCadenas, S. Jonic, F. Tama et al., “3DEM Loupe: analysis of macromolecular dynamics using structures from electron microscopy,” Nucleic acids research, vol. 41, pp. W363–W367, 2013. View at: Publisher Site  Google Scholar
 S. Jonic and C. O. S. Sorzano, “Coarsegraining of volumes for modeling of structure and dynamics in electron microscopy: algorithm to automatically control accuracy of approximation,” IEEE Journal of Selected Topics in Signal Processing, vol. 10, no. 1, pp. 161–173, 2016. View at: Publisher Site  Google Scholar
 S. Jonic, J. Vargas, R. Melero, J. GomezBlanco, J. M. Carazo, and C. O. Sorzano, “Denoising of highresolution singleparticle electronmicroscopy density maps by their approximation using threedimensional Gaussian functions,” Journal of Structural Biology, vol. 194, no. 3, pp. 423–433, 2016. View at: Publisher Site  Google Scholar
 G. McLachlan and D. Peel, Finite mixture models, Wiley Series in Probability and Statistics: Applied Probability and Statistics, WileyInterscience, New York, NY, USA, 2000. View at: Publisher Site  MathSciNet
 M. M. Tirion, “Large amplitude elastic motions in proteins from a singleparameter, atomic analysis,” Physical Review Letters, vol. 77, no. 9, pp. 1905–1908, 1996. View at: Publisher Site  Google Scholar
 F. Tama, W. Wriggers, and C. L. Brooks III, “Exploring global distortions of biological macromolecules and assemblies from lowresolution structural information and elastic network theory,” Journal of Molecular Biology, vol. 321, no. 2, pp. 297–305, 2002. View at: Publisher Site  Google Scholar
 F. Tama and Y.H. Sanejouand, “Conformational change of proteins arising from normal mode calculations,” Protein Engineering, vol. 14, no. 1, pp. 1–6, 2001. View at: Publisher Site  Google Scholar
 K. Suhre and Y.H. Sanejouand, “ElNémo: a normal mode web server for protein movement analysis and the generation of templates for molecular replacement,” Nucleic Acids Research, vol. 32, pp. W610–W614, 2004. View at: Publisher Site  Google Scholar
 J. Brink, S. J. Ludtke, Y. Kong, S. J. Wakil, J. Ma, and W. Chiu, “Experimental verification of conformational variation of human fatty acid synthase as predicted by normal mode analysis,” Structure, vol. 12, no. 2, pp. 185–191, 2004. View at: Publisher Site  Google Scholar
 C. O. S. Sorzano, J. M. de la RosaTrevín, F. Tama, and S. Jonić, “Hybrid electron microscopy normal mode analysis graphical interface and protocol,” Journal of Structural Biology, vol. 188, no. 2, pp. 134–141, 2014. View at: Publisher Site  Google Scholar
 J. Frank and A. Ourmazd, “Continuous changes in structure mapped by manifold embedding of singleparticle data in cryoEM,” Methods, vol. 100, pp. 61–67, 2015. View at: Publisher Site  Google Scholar
 E. Behrmann, J. Loerke, T. V. Budkevich et al., “Structural snapshots of actively translating human ribosomes,” Cell, vol. 161, no. 4, pp. 845–857, 2015. View at: Publisher Site  Google Scholar
 P. Joubert and M. Habeck, “Bayesian inference of initial models in cryoelectron microscopy using pseudoatoms,” Biophysical Journal, vol. 108, no. 5, pp. 1165–1175, 2015. View at: Publisher Site  Google Scholar
 S. C. Harrison, A. J. Olson, C. E. Schutt, F. K. Winkler, and G. Bricogne, “Tomato bushy stunt virus at 2.9 Å resolution,” Nature, vol. 276, no. 5686, pp. 368–373, 1978. View at: Publisher Site  Google Scholar
 W. Zhang, J. A. Dunkle, and J. H. D. Cate, “Structures of the ribosome in intermediate states of ratcheting,” Science, vol. 325, no. 5943, pp. 1014–1017, 2009. View at: Publisher Site  Google Scholar
 R. Aramayo, C. Mérigoux, E. Larquet et al., “Divalent iondependent swelling of tomato bushy stunt virus: a multiapproach study,” Biochimica et Biophysica Acta—General Subjects, vol. 1724, no. 3, pp. 345–354, 2005. View at: Publisher Site  Google Scholar
 I. K. Robinson and S. C. Harrison, “Structure of the expanded state of tomato bushy stunt virus,” Nature, vol. 297, no. 5867, pp. 563–568, 1982. View at: Publisher Site  Google Scholar
 S. Klinge, R. ÑezRamírez, O. Llorca, and L. Pellegrini, “3D architecture of DNA Pol α reveals the functional core of multisubunit replicative polymerases,” The EMBO Journal, vol. 28, no. 13, pp. 1978–1987, 2009. View at: Publisher Site  Google Scholar
 G. S. Hansman, D. W. Taylor, J. S. McLellan et al., “Structural basis for broad detection of genogroup II noroviruses by a monoclonal antibody that binds to a site occluded in the viral particle,” Journal of Virology, vol. 86, no. 7, pp. 3635–3646, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Slavica Jonić and Carlos Oscar S. Sorzano. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.