Application of template matching technique to particle detection in electron micrographs

Template matching together with the comprehensive theory of image formation in electron microscope provides an optimal (in Bayesian sense) tool for solving one of the outstanding problems in single particle analysis, i.e., automatic selection of particle views from noisy micrograph fields. The method is based on the assumption that the reference three-dimensional structure is known and that the relevant parameters of the model of the image formation process can be estimated. In the first stage of the procedure, a set of possible particle views is generated using the available reference structure. The template images are constructed as linear combinations of available particle views using a clustering technique. Next, the micrograph noise characteristic is established using an automated contrast transfer function (CTF) estimation procedure. Finally, the CTF parameters calculated are used to construct a matched filter and correlation functions corresponding to the available template images are calculated. In order to alleviate the problem of the biased caused by varying image formation conditions, a decision making strategy based on the predicted distribution of correlation coefficients is proposed. It is demonstrated that due to the inclusion of CTF considerations, the template matching method performed very well in a broad range of microscopy conditions.     

Automatic CTF correction for single particles based upon multivariate statistical analysis of individual power spectra

Three-dimensional electron cryomicroscopy of randomly oriented single particles is a method that is suitable for the determination of three-dimensional structures of macromolecular complexes at molecular resolution. However, the electron-microscopical projection images are modulated by a contrast transfer function (CTF) that prevents the calculation of three-dimensional reconstructions of biological complexes at high resolution from uncorrected images. We describe here an automated method for the accurate determination and correction of the CTF parameters defocus, twofold astigmatism and amplitude-contrast proportion from single-particle images. At the same time, the method allows the frequency-dependent signal decrease (B factor) and the non-convoluted background signal to be estimated. The method involves the classification of the power spectra of single-particle images into groups with similar CTF parameters; this is done by multivariate statistical analysis (MSA) and hierarchically ascending classification (HAC). Averaging over several power spectra generates class averages with enhanced signal-to-noise ratios. The correct CTF parameters can be deduced from these class averages by applying an iterative correlation procedure with theoretical CTF functions; they are then used to correct the raw images. Furthermore, the method enables the tilt axis of the sample holder to be determined and allows the elimination of individual poor-quality images that show high drift or charging effects.     

EMAN: semiautomated software for high-resolution single-particle reconstructions

We present EMAN (Electron Micrograph ANalysis), a software package for performing semiautomated single-particle reconstructions from transmission electron micrographs. The goal of this project is to provide software capable of performing single-particle reconstructions beyond 10 A as such high-resolution data become available. A complete single-particle reconstruction algorithm is implemented. Options are available to generate an initial model for particles with no symmetry, a single axis of rotational symmetry, or icosahedral symmetry. Model refinement is an iterative process, which utilizes classification by model-based projection matching. CTF (contrast transfer function) parameters are determined using a new paradigm in which data from multiple micrographs are fit simultaneously. Amplitude and phase CTF correction is then performed automatically as part of the refinement loop. A graphical user interface is provided, so even those with little image processing experience will be able to begin performing reconstructions. Advanced users can directly use the lower level shell commands and even expand the package utilizing EMAN's extensive image-processing library. The package was written from scratch in C++ and is provided free of charge on our Web site. We present an overview of the package as well as several conformance tests with simulated data.     

PC12 Cells that Lack Synaptotagmin I Exhibit Loss of a Subpool of Small Dense Core Vesicles

Neurons communicate by releasing neurotransmitters that are stored in intracellular vesicular compartments. PC12 cells are frequently used as a model secretory cell line that is described to have two subpools of vesicles: small clear vesicles and dense core vesicles. We measured transmitter molecules released from vesicles in NGF-differentiated PC12 cells using carbon-fiber amperometry, and relative diameters of individual vesicles using electron microscopy. Both amperometry and electron micrograph data were analyzed by statistical and machine learning methods for Gaussian mixture models. An electron microscopy size correction algorithm was used to predict and correct for observation bias of vesicle size due to tangential slices through some vesicles. Expectation maximization algorithms were used to perform maximum likelihood estimation for the Gaussian parameters of different populations of vesicles, and were shown to be better than histogram and cumulative distribution function methods for analyzing mixed populations. The Bayesian information criterion was used to determine the most likely number of vesicle subpools observed in the amperometric and electron microscopy data. From this analysis, we show that there are three major subpools, not two, of vesicles stored and released from PC12 cells. The three subpools of vesicles include small clear vesicles and two subpools of dense core vesicles, a small and a large dense core vesicle subpool. Using PC12 cells stably transfected with short-hairpin RNA targeted to synaptotagmin I, an exocytotic Ca²⁺ sensor, we show that the presence and release of the small dense core vesicle subpool is dependent on synaptotagmin I. Furthermore, synaptotagmin I also plays a role in the formation and/or maintenance of the small dense core vesicle subpool in PC12 cells.     

FASTDEF: Fast defocus and astigmatism estimation for high-throughput transmission electron microscopy

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].     

Error analysis in the determination of the electron microscopical contrast transfer function parameters from experimental power Spectra

Background  

The transmission electron microscope is used to acquire structural information of macromolecular complexes. However, as any other imaging device, it introduces optical aberrations that must be corrected if high-resolution structural information is to be obtained. The set of all aberrations are usually modeled in Fourier space by the so-called Contrast Transfer Function (CTF). Before correcting for the CTF, we must first estimate it from the electron micrographs. This is usually done by estimating a number of parameters specifying a theoretical model of the CTF. This estimation is performed by minimizing some error measure between the theoretical Power Spectrum Density (PSD) and the experimentally observed PSD. The high noise present in the micrographs, the possible local minima of the error function for estimating the CTF parameters, and the cross-talking between CTF parameters may cause errors in the estimated CTF parameters.     

Accurate determination of local defocus and specimen tilt in electron microscopy

Accurate knowledge of defocus and tilt parameters is essential for the determination of three-dimensional protein structures at high resolution using electron microscopy. We present two computer programs, CTFFIND3 and CTFTILT, which determine defocus parameters from images of untilted specimens, as well as defocus and tilt parameters from images of tilted specimens, respectively. Both programs use a simple algorithm that fits the amplitude modulations visible in a power spectrum with a calculated contrast transfer function (CTF). The background present in the power spectrum is calculated using a low-pass filter. The background is then subtracted from the original power spectrum, allowing the fitting of only the oscillatory component of the CTF. CTFTILT determines specimen tilt parameters by measuring the defocus at a series of locations on the image while constraining them to a single plane. We tested the algorithm on images of two-dimensional crystals by comparing the results with those obtained using crystallographic methods. The images also contained contrast from carbon support film that added to the visibility of the CTF oscillations. The tests suggest that the fitting procedure is able to determine the image defocus with an error of about 10nm, whereas tilt axis and tilt angle are determined with an error of about 2 degrees and 1 degrees, respectively. Further tests were performed on images of single protein particles embedded in ice that were recorded from untilted or slightly tilted specimens. The visibility of the CTF oscillations from these images was reduced due to the lack of a carbon support film. Nevertheless, the test results suggest that the fitting procedure is able to determine image defocus and tilt angle with errors of about 100 nm and 6 degrees, respectively.     

Aβ(39-42) modulates Aβ oligomerization but not fibril formation

Recently, certain C-terminal fragments (CTFs) of Aβ42 have been shown to be effective inhibitors of Aβ42 toxicity. Here, we examine the interactions between the shortest CTF in the original series, Aβ(39-42), and full-length Aβ. Mass spectrometry results indicate that Aβ(39-42) binds directly to Aβ monomers and to the n = 2, 4, and 6 oligomers. The Aβ42:Aβ(39-42) complex is further probed using molecular dynamics simulations. Although the CTF was expected to bind to the hydrophobic C-terminus of Aβ42, the simulations show that Aβ(39-42) binds at several locations on Aβ42, including the C-terminus, other hydrophobic regions, and preferentially in the N-terminus. Ion mobility-mass spectrometry (IM-MS) and electron microscopy experiments indicate that Aβ(39-42) disrupts the early assembly of full-length Aβ. Specifically, the ion-mobility results show that Aβ(39-42) prevents the formation of large decamer/dodecamer Aβ42 species and, moreover, can remove these structures from solution. At the same time, thioflavin T fluorescence and electron microscopy results show that the CTF does not inhibit fibril formation, lending strong support to the hypothesis that oligomers and not amyloid fibrils are the Aβ form responsible for toxicity. The results emphasize the role of small, soluble assemblies in Aβ-induced toxicity and suggest that Aβ(39-42) inhibits Aβ-induced toxicity by a unique mechanism, modulating early assembly into nontoxic hetero-oligomers, without preventing fibril formation.     

Three-Dimensional Reconstruction with Contrast Transfer Function Correction from Energy-Filtered Cryoelectron Micrographs: Procedure and Application to the 70SEscherichia coliRibosome

Cryoelectron microscopy provides the means of studying macromolecules in their native state. However, the contrast transfer function (CTF) makes the images and the three-dimensional (3D) maps derived from them difficult to interpret. We developed methods to determine the CTF from experimental data and to obtain a CTF-corrected 3D reconstruction. The CTF correction and 3D reconstruction accomplished in one step make it easy to combine different defocus data sets and decrease the error accumulation in the computation. These methods were applied to energy-filtered images of the 70SEscherichia coliribosome, resulting in a distortion-free 3D map of the ribosome at 1/24.5 Å⁻¹resolution, as determined by the differential phase residual resolution criterion.     

Saccharomyces cerevisiae CTF18 and CTF4 are required for sister chromatid cohesion

CTF4 and CTF18 are required for high-fidelity chromosome segregation. Both exhibit genetic and physical ties to replication fork constituents. We find that absence of either CTF4 or CTF18 causes sister chromatid cohesion failure and leads to a preanaphase accumulation of cells that depends on the spindle assembly checkpoint. The physical and genetic interactions between CTF4, CTF18, and core components of replication fork complexes observed in this study and others suggest that both gene products act in association with the replication fork to facilitate sister chromatid cohesion. We find that Ctf18p, an RFC1-like protein, directly interacts with Rfc2p, Rfc3p, Rfc4p, and Rfc5p. However, Ctf18p is not a component of biochemically purified proliferating cell nuclear antigen loading RF-C, suggesting the presence of a discrete complex containing Ctf18p, Rfc2p, Rfc3p, Rfc4p, and Rfc5p. Recent identification and characterization of the budding yeast polymerase kappa, encoded by TRF4, strongly supports a hypothesis that the DNA replication machinery is required for proper sister chromatid cohesion. Analogous to the polymerase switching role of the bacterial and human RF-C complexes, we propose that budding yeast RF-C(CTF18) may be involved in a polymerase switch event that facilities sister chromatid cohesion. The requirement for CTF4 and CTF18 in robust cohesion identifies novel roles for replication accessory proteins in this process.     

TomoAlign: A novel approach to correcting sample motion and 3D CTF in CryoET

TomoAlign is a software package that integrates tools to mitigate two important resolution limiting factors in cryoET, namely the beam-induced sample motion and the contrast transfer function (CTF) of the microscope. The package is especially focused on cryoET of thick specimens where fiducial markers are required for accurate tilt-series alignment and sample motion estimation. TomoAlign models the beam-induced sample motion undergone during the tilt-series acquisition. The motion models are used to produce motion-corrected subtilt-series centered on the particles of interest. In addition, the defocus of each particle at each tilt image is determined and can be corrected, resulting in motion-corrected and CTF-corrected subtilt-series from which the subtomograms can be computed. Alternatively, the CTF information can be passed on so that CTF correction can be carried out entirely within external packages like Relion. TomoAlign serves as a versatile tool that can streamline the cryoET workflow from initial alignment of tilt-series to final subtomogram averaging during in situ structure determination.     

Auto-accumulation method using simulated annealing enables fully automatic particle pickup completely free from a matching template or learning data

Single-particle analysis is a 3-D structure determining method using electron microscopy (EM). In this method, a large number of projections is required to create 3-D reconstruction. In order to enable completely automatic pickup without a matching template or a training data set, we established a brand-new method in which the frames to pickup particles are randomly shifted and rotated over the electron micrograph and, using the total average image of the framed images as an index, each frame reaches a particle. In this process, shifts are selected to increase the contrast of the average. By iterated shifts and further selection of the shifts, the frames are induced to shift so as to surround particles. In this algorithm, hundreds of frames are initially distributed randomly over the electron micrograph in which multi-particle images are dispersed. Starting with these frames, one of them is selected and shifted randomly, and acceptance or non-acceptance of its new position is judged using the simulated annealing (SA) method in which the contrast score of the total average image is adopted as an index. After iteration of this process, the position of each frame converges so as to surround a particle and the framed images are picked up. This method is the first unsupervised fully automatic particle picking method which is applicable to EM of various kinds of proteins, especially to low-contrasted cryo-EM protein images.     

High-resolution cryo-EM performance comparison of two latest-generation cryo electron microscopes on the human ribosome

Recent technological advances in cryo electron microscopy (cryo-EM) have led to new opportunities in the structural biology field. Here we benchmark the performance of two 300 kV latest-generation cryo electron microscopes, Titan Krios G4 from Thermofisher Scientific and CRYO ARM 300 from Jeol, with regards to achieving high resolution single particle reconstructions on a real case sample. We compare potentially limiting factors such as drift rates, astigmatism & coma aberrations and performance during image processing and show that both microscopes, while comprising rather different technical setups & parameter settings and equipped with different types of energy filters & cameras, achieve a resolution of around 2 Å on the human ribosome, a non-symmetric object which constitutes a key drug target. Astigmatism correction, CTF refinement and correction of higher order aberrations through refinement in separate optics groups helped to account for astigmatism/coma caused by beam tilting during multi-spot and multi-hole acquisition in neighbouring holes without stage movement. The obtained maps resolve Mg²⁺ ions, water molecules, inhibitors and side-chains including chemical modifications. The fact that both instruments can resolve such detailed features will greatly facilitate understanding molecular mechanisms of various targets and helps in cryo-EM structure based drug design. The methods and analysis tools used here will be useful also to characterize existing instruments and optimize data acquisition settings and are applicable broadly to other drug targets in structural biology.     

Focal pair merging for contrast enhancement of single particles

A common technique in transmission electron microscopy is the collection of a focal pair, in which the first, close to focus image contains higher resolution information at lower contrast, and the second, far from focus image has high contrast but less reliable high-resolution information. Typically these second micrographs are used for purposes of particle selection, orientation estimate or micrograph evaluation. We introduce a technique for merging the information from both images, including signal to noise ratio weighting, contrast transfer function correction, and optional Weiner filtration. This produces a composite image with reduced contrast transfer function artifacts and optimized contrast. This technique is useful in numerous cases where low-contrast images are produced, such as small particles, proteins solubilized in detergent or projects with high-resolution goals when the first image is taken very close to focus.     

Accurate modeling of single-particle cryo-EM images quantitates the benefits expected from using Zernike phase contrast

The use of a Zernike-type phase plate in biologic cryo-electron microscopy allows the imaging, without using defocus, of what are predominantly phase objects. It is thought that such phase-plate implementations might result in higher quality images, free from the problems of CTF correction that occur when images must be recorded at extremely high values of defocus. In single-particle cryo-electron microscopy it is hoped that these improvements in image quality will facilitate work on structures that have proved difficult to study, either because of their relatively small size or because the structures are not completely homogeneous. There is still a need, however, to quantitate how much improvement can be gained by using a phase plate for single-particle cryo-electron microscopy. We present a method for quantitatively modeling the images recorded with 200keV electrons, for single particles embedded in vitreous ice. We then investigate what difference the use of a phase-plate device could have on the processing of single-particle data. We confirm that using a phase plate results in single-particle datasets in which smaller molecules can be detected, particles can be more accurately aligned and problems of heterogeneity can be more easily addressed.     

Precise beam-tilt alignment and collimation are required to minimize the phase error associated with coma in high-resolution cryo-EM

Electron microscopy at a resolution of 0.4nm or better requires more careful adjustment of the illumination than is the case at a resolution of 0.8nm. The use of current-axis alignment is not always sufficient, for example, to avoid the introduction of large phase errors, at higher resolution, due to axial coma. In addition, one must also ensure that off-axis coma does not corrupt the data quality at the higher resolution. We particularly emphasize that the standard CTF correction does not account for the phase error associated with coma. We explain the cause of both axial coma and the typically most troublesome component of off-axis coma in terms of the well-known shift of the electron diffraction pattern relative to the optical axis that occurs when the illumination is not parallel to the axis. We review the experimental conditions under which coma causes unacceptably large phase errors, and we discuss steps that can be taken when setting up the conditions of illumination, so as to ensure that neither axial nor off-axis coma is a problem.     

CTF determination and correction for low dose tomographic tilt series

The resolution of cryo-electron tomography can be limited by the first zero of the microscope’s contrast transfer function (CTF). To achieve higher resolution, it is critical to determine the CTF and correct its phase inversions. However, the extremely low signal-to-noise ratio (SNR) and the defocus gradient in the projections of tilted specimens make this process challenging. Two programs, CTFPLOTTER and CTFPHASEFLIP, have been developed to address these issues. CTFPLOTTER obtains a 1D power spectrum by periodogram averaging and rotational averaging and it estimates the noise background with a novel approach, which uses images taken with no specimen. The background-subtracted 1D power spectra from image regions at different defocus values are then shifted to align their first zeros and averaged together. This averaging improves the SNR sufficiently that it becomes possible to determine the defocus for subsets of the tilt series rather than just the entire series. CTFPHASEFLIP corrects images line-by-line by inverting phases appropriately in thin strips of the image at nearly constant defocus. CTF correction by these methods is shown to improve the resolution of aligned, averaged particles extracted from tomograms. However, some restoration of Fourier amplitudes at high frequencies is important for seeing the benefits from CTF correction.     

Optimal determination of particle orientation, absolute hand, and contrast loss in single-particle electron cryomicroscopy

A computational procedure is described for assigning the absolute hand of the structure of a protein or assembly determined by single-particle electron microscopy. The procedure requires a pair of micrographs of the same particle field recorded at two tilt angles of a single tilt-axis specimen holder together with the three-dimensional map whose hand is being determined. For orientations determined from particles on one micrograph using the map, the agreement (average phase residual) between particle images on the second micrograph and map projections is determined for all possible choices of tilt angle and axis. Whether the agreement is better at the known tilt angle and axis of the microscope or its inverse indicates whether the map is of correct or incorrect hand. An increased discrimination of correct from incorrect hand (free hand difference), as well as accurate identification of the known values for the tilt angle and axis, can be used as targets for rapidly optimizing the search or refinement procedures used to determine particle orientations. Optimized refinement reduces the tendency for the model to match noise in a single image, thus improving the accuracy of the orientation determination and therefore the quality of the resulting map. The hand determination and refinement optimization procedure is applied to image pairs of the dihydrolipoyl acetyltransferase (E2) catalytic core of the pyruvate dehydrogenase complex from Bacillus stearothermophilus taken by low-dose electron cryomicroscopy. Structure factor amplitudes of a three-dimensional map of the E2 catalytic core obtained by averaging untilted images of 3667 icosahedral particles are compared to a scattering reference using a Guinier plot. A noise-dependent structure factor weight is derived and used in conjunction with a temperature factor (B=-1000A(2)) to restore high-resolution contrast without amplifying noise and to visualize molecular features to 8.7A resolution, according to a new objective criterion for resolution assessment proposed here.     

Combining efficient conformational sampling with a deformable elastic network model facilitates structure refinement at low resolution

Structural studies of large proteins and protein assemblies are a difficult and pressing challenge in molecular biology. Experiments often yield only low-resolution or sparse data that are not sufficient to fully determine atomistic structures. We have developed a general geometry-based algorithm that efficiently samples conformational space under constraints imposed by low-resolution density maps obtained from electron microscopy or X-ray crystallography experiments. A deformable elastic network (DEN) is used to restrain the sampling to prior knowledge of an approximate structure. The DEN restraints dramatically reduce over-fitting, especially at low resolution. Cross-validation is used to optimally weight the structural information and experimental data. Our algorithm is robust even for noise-added density maps and has a large radius of convergence for our test case. The DEN restraints can also be used to enhance reciprocal space simulated annealing refinement.