Validation of the orthogonal tilt reconstruction method with a biological test sample

Electron microscopy of frozen-hydrated samples (cryo-EM) can yield high resolution structures of macromolecular complexes by accurately determining the orientation of large numbers of experimental views of the sample relative to an existing 3D model. The "initial model problem", the challenge of obtaining these orientations ab initio, remains a major bottleneck in determining the structure of novel macromolecules, chiefly those lacking internal symmetry. We previously proposed a method for the generation of initial models--orthogonal tilt reconstruction (OTR)--that bypasses limitations inherent to the other two existing methods, random conical tilt (RCT) and angular reconstitution (AR). Here we present a validation of OTR with a biological test sample whose structure was previously solved by RCT: the complex between the yeast exosome and the subunit Rrp44. We show that, as originally demonstrated with synthetic data, OTR generates initial models that do not exhibit the "missing cone" artifacts associated with RCT and show an isotropic distribution of information when compared with the known structure. This eliminates the need for further user intervention to solve these artifacts and makes OTR ideal for automation and the analysis of heterogeneous samples. With the former in mind, we propose a set of simple quantitative criteria that can be used, in combination, to select from a large set of initial reconstructions a subset that can be used as reliable references for refinement to higher resolution.     

The orthogonal tilt reconstruction method: an approach to generating single-class volumes with no missing cone for ab initio reconstruction of asymmetric particles

Generating reliable initial models is a critical step in the reconstruction of asymmetric single-particles by 3D electron microscopy. This is particularly difficult to do if heterogeneity is present in the sample. The Random Conical Tilt (RCT) method, arguably the most robust presently to accomplish this task, requires significant user intervention to solve the "missing cone" problem. We present here a novel approach, termed the orthogonal tilt reconstruction method, that eliminates the missing cone altogether, making it possible for single-class volumes to be used directly as initial references in refinement without further processing. The method involves collecting data at +45 degrees and -45 degrees tilts and only requires that particles adopt a relatively large number of orientations on the grid. One tilted data set is used for alignment and classification and the other set--which provides views orthogonal to those in the first--is used for reconstruction, resulting in the absence of a missing cone. We have tested this method with synthetic data and compared its performance to that of the RCT method. We also propose a way of increasing the level of homogeneity in individual 2D classes (and volumes) in a heterogeneous data set and identifying the most homogeneous volumes.     

Automated acquisition of electron microscopic random conical tilt sets

Single particle reconstruction using the random conical tilt data collection geometry is a robust method for the initial determination of macromolecular structures by electron microscopy. Unfortunately, the broad adoption of this powerful approach has been limited by the practical challenges inherent in manual data collection of the required pairs of matching high and low tilt images (typically 60 degrees and 0 degrees). The microscopist is obliged to keep the imaging area centered during tilting as well as to maintain accurate focus in the tilted image while minimizing the overall electron dose, a challenging and time consuming process. To help solve these problems, we have developed an automated system for the rapid acquisition of accurately aligned and focused tilt pairs. The system has been designed to minimize the dose incurred during alignment and focusing, making it useful in both negative stain and cryo-electron microscopy. The system includes a feature for montaging untilted images to ensure that all of the particles in the tilted image may be used in the reconstruction.     

Initial evaluation of a direct detection device detector for single particle cryo-electron microscopy

We report on initial results of using a new direct detection device (DDD) for single particle reconstruction of vitreous ice embedded specimens. Images were acquired on a Tecnai F20 at 200 keV and a nominal magnification of 29,000×. This camera has a significantly improved signal to noise ratio and modulation transfer function (MTF) at 200 keV compared to a standard CCD camera installed on the same microscope. Control of the DDD has been integrated into Leginon, an automated data collection system. Using GroEL as a test specimen, we obtained images of ∼30 K particles with the CCD and the DDD from the same specimen sample using essentially identical imaging conditions. Comparison of the maps reconstructed from the CCD images and the DDD images demonstrates the improved performance of the DDD. We also obtained a 3D reconstruction from ∼70 K GroEL particles acquired using the DDD; the quality of the density map demonstrates the potential of this new recording device for cryoEM data acquisition.     

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.     

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.     

Electron tomography of vitreous sections from cultured mammalian cells

Cryo-electron tomography of appropriately thin, frozen-hydrated biological specimens has excellent potential for investigating the 3D macromolecular architecture of eukaryotic cells and tissues. Since cardiomyocytes are too thick to be visualised in an intact state, we grew immortalised cell line HL-1 to sub-confluency and harvested the cells by enzymatic detachment prior to hyperbaric freezing, ultramicrotomy, and tomography. We improved the efficiency of tomographic acquisition from vitreous cryosections by implementing two new features: (1) fluorescence microscopy at cryogenic temperatures to search for features of interest without expending any of the tolerable electron dose on secondary (non-imaging) tasks, and (2) the use of colloidal gold as fiducial markers. Vital fluorescent staining and subsequent cryo-fluorescence microscopy of vitreous sections were used to localise mitochondria lying in positions suitable for acquiring tilt series, taking into account section flatness, presence of contamination and proximity to grid bars. To provide a simple and robust means of aligning tomograms, we developed a universally applicable protocol for depositing colloidal gold onto vitreous sections, analogous to the method for applying quantum dots described by Masich et al. [Masich, S., Östberg, T., Norlén, L., Shupliakov, O., Daneholt, B., 2006. A procedure to deposit fiducial markers on vitreous cryo-sections for cellular tomography. J. Struct. Biol. 156, 461–468]. Tomograms of thin sections (nominal thickness 65–85 nm) of cardiac mitochondria revealed the interconnectivity of cristae and junctions with the inner mitochondrial membrane. In some cases, ATP synthases could be identified without ambiguity. These findings confirm the feasibility of investigating the structural biology of mammalian cells in three dimensions and at a resolution of 6–8 nm.     

Reprint of "On the feasibility of visualizing ultrasmall gold labels in biological specimens by STEM tomography" [J. Struct. Biol. 159 (2007) 507-522

Labeling with heavy atom clusters attached to antibody fragments is an attractive technique for determining the 3D distribution of specific proteins in cells using electron tomography. However, the small size of the labels makes them very difficult to detect by conventional bright-field electron tomography. Here, we evaluate quantitative scanning transmission electron microscopy (STEM) at a beam voltage of 300 kV for detecting 11-gold atom clusters (Undecagold) and 1.4 nm-diameter nanoparticles (Nanogold) for a variety of specimens and imaging conditions. STEM images as well as tomographic tilt series are simulated by means of the NIST Elastic-Scattering Cross-Section Database for gold clusters embedded in carbon. The simulations indicate that the visibility in 2D of Undecagold clusters in a homogeneous matrix is maximized for low inner collection semi-angles of the STEM annular dark-field detector (15–20 mrad). Furthermore, our calculations show that the visibility of Undecagold in 3D reconstructions is significantly higher than in 2D images for an inhomogeneous matrix corresponding to fluctuations in local density. The measurements demonstrate that it is possible to detect Nanogold particles in plastic sections of tissue freeze-substituted in the presence of osmium. STEM tomography has the potential to localize specific proteins in permeabilized cells using antibody fragments tagged with small heavy atom clusters. Our quantitative analysis provides a framework for determining the detection limits and optimal experimental conditions for localizing these small clusters.     

On the feasibility of visualizing ultrasmall gold labels in biological specimens by STEM tomography

Labeling with heavy atom clusters attached to antibody fragments is an attractive technique for determining the 3D distribution of specific proteins in cells using electron tomography. However, the small size of the labels makes them very difficult to detect by conventional bright-field electron tomography. Here, we evaluate quantitative scanning transmission electron microscopy (STEM) at a beam voltage of 300 kV for detecting 11-gold atom clusters (Undecagold) and 1.4 nm-diameter nanoparticles (Nanogold) for a variety of specimens and imaging conditions. STEM images as well as tomographic tilt series are simulated by means of the NIST Elastic-Scattering Cross-Section Database for gold clusters embedded in carbon. The simulations indicate that the visibility in 2D of Undecagold clusters in a homogeneous matrix is maximized for low inner collection semi-angles of the STEM annular dark-field detector (15–20 mrad). Furthermore, our calculations show that the visibility of Undecagold in 3D reconstructions is significantly higher than in 2D images for an inhomogeneous matrix corresponding to fluctuations in local density. The measurements demonstrate that it is possible to detect Nanogold particles in plastic sections of tissue freeze-substituted in the presence of osmium. STEM tomography has the potential to localize specific proteins in permeabilized cells using antibody fragments tagged with small heavy atom clusters. Our quantitative analysis provides a framework for determining the detection limits and optimal experimental conditions for localizing these small clusters.     

Parakeet: a digital twin software pipeline to assess the impact of experimental parameters on tomographic reconstructions for cryo-electron tomography

In cryo-electron tomography (cryo-ET) of biological samples, the quality of tomographic reconstructions can vary depending on the transmission electron microscope (TEM) instrument and data acquisition parameters. In this paper, we present Parakeet, a ‘digital twin’ software pipeline for the assessment of the impact of various TEM experiment parameters on the quality of three-dimensional tomographic reconstructions. The Parakeet digital twin is a digital model that can be used to optimize the performance and utilization of a physical instrument to enable in silico optimization of sample geometries, data acquisition schemes and instrument parameters. The digital twin performs virtual sample generation, TEM image simulation, and tilt series reconstruction and analysis within a convenient software framework. As well as being able to produce physically realistic simulated cryo-ET datasets to aid the development of tomographic reconstruction and subtomogram averaging programs, Parakeet aims to enable convenient assessment of the effects of different microscope parameters and data acquisition parameters on reconstruction quality. To illustrate the use of the software, we present the example of a quantitative analysis of missing wedge artefacts on simulated planar and cylindrical biological samples and discuss how data collection parameters can be modified for cylindrical samples where a full 180° tilt range might be measured.     

Automated correlation of single particle tilt pairs for Random Conical Tilt and Orthogonal Tilt Reconstructions

One of the major methodological challenges in single particle electron microscopy is obtaining initial reconstructions which represent the structural heterogeneity of the dataset. Random Conical Tilt and Orthogonal Tilt Reconstruction techniques in combination with 3D alignment and classification can be used to obtain initial low-resolution reconstructions which represent the full range of structural heterogeneity of the dataset. In order to achieve statistical significance, however, a large number of 3D reconstructions, and, in turn, a large number of tilted image pairs are required. The extraction of single particle tilted image pairs from micrographs can be tedious and time-consuming, as it requires intensive user input even for semi-automated approaches. To overcome the bottleneck of manual selection of a large number of tilt pairs, we developed an algorithm for the correlation of single particle images from tilted image pairs in a fully automated and user-independent manner. The algorithm reliably correlates correct pairs even from noisy micrographs. We further demonstrate the applicability of the algorithm by using it to obtain initial references both from negative stain and unstained cryo datasets.     

Retrospective: radiation damage and its associated "information limitations"

The fact that radiation damage would limit the usefulness of electron microscopy with biological specimens was a concern in the earliest days of the field. Good estimates of what that limitation must be can be made by using Rose's empirical relationship between the inherent image contrast, the exposure used to record an image, and the smallest feature size that is detectable. Such estimates show that it is necessary to average many images in order to obtain statistically well-defined data at high resolution. Structures are now routinely obtained by averaging large numbers of shot-noise limited images, and some of these extend to atomic resolution. The signal level in current images is nevertheless far below what physics would allow it to be. A possible explanation is that beam-induced movement limits the quality of images recorded by electron microscopy. For specimens embedded in vitreous ice, beam-induced movement can even be severe enough to limit the resolution achieved during tomographic reconstruction. The fact that very high-quality images can nevertheless be obtained, although only unpredictably, suggests that it may be possible to devise new techniques of specimen preparation and/or data collection that at least partially overcome beam-induced movement. If so, the need for image averaging would be correspondingly reduced.     

Automated cryoEM data acquisition and analysis of 284742 particles of GroEL

One of the goals in developing our automated electron microscopy data acquisition system, Leginon, was to improve both the ease of use and the throughput of the process of acquiring low dose images of macromolecular specimens embedded in vitreous ice. In this article, we demonstrate the potential of the Leginon system for high-throughput data acquisition by describing an experiment in which we acquired images of more than 280,000 particles of GroEL in a single 25 h session at the microscope. We also demonstrate the potential for an automated pipeline for molecular microscopy by showing that these particles can be subjected to completely automated procedures to reconstruct a three-dimensional (3D) density map to a resolution better than 8 A. In generating the 3D maps, we used a variety of metadata associated with the data acquisition and processing steps to sort and select the particles. These metadata provide a number of insights into factors that affect the quality of the acquired images and the resulting reconstructions. In particular, we show that the resolution of the reconstructed 3D density maps improves with decreasing ice thickness. These data provide a basis for assessing the capabilities of high-throughput macromolecular microscopy.     

A procedure to deposit fiducial markers on vitreous cryo-sections for cellular tomography

We describe a novel approach for the accurate alignment of images in electron tomography of vitreous cryo-sections. Quantum dots, suspended in organic solvents at cryo-temperatures, are applied directly onto the sections and are subsequently used as fiducial markers to align the tilt series. Data collection can be performed from different regions of the vitreous sections, even when the sections touch the grid only at a few places. We present high-resolution tomograms of some organelles in cryo-sections of human skin cells using this method. The average error in image alignment was about 1nm and the resolution was estimated to be 5-7nm. Thus, the use of section-attached quantum dots as fiducial markers in electron tomography of vitreous cryo-sections facilitates high-resolution in situ 3D imaging of organelles and macromolecular complexes in their native hydrated state.     

Automated Identification of Filaments in Cryoelectron Microscopy Images

Since the foundation for the three-dimensional image reconstruction of helical objects from electron micrographs was laid more than 30 years ago, there have been sustained developments in specimen preparation, data acquisition, image analysis, and interpretation of results. However, the boxing of filaments in large numbers of images--one of the critical steps toward the reconstruction at high resolution--is still constrained by manual processing even though interactive interfaces have been built to aid the tedious and sometimes inaccurate boxing process. This article describes an accurate approach for automated detection of filamentous structures in low-contrast images acquired in defocus pairs using cryoelectron microscopy. The performance of the approach has been evaluated across various magnifications and at a series of defocus values using tobacco mosaic virus (TMV) preserved in vitreous ice as a test specimen. By integrating the proposed approach into our automated data acquisition and reconstruction system, we are now able to generate a three-dimensional map of TMV to approximately 10-A resolution within 24 h of inserting the specimen grid into the microscope.     

AFM characterization of tilt and intrinsic flexibility of Rhodobacter sphaeroides light harvesting complex 2 (LH2)

Atomic force microscopy (AFM) has developed into a powerful tool to investigate membrane protein surfaces in a close-to-native environment. Here we report on the surface topography of Rhodobacter sphaeroides light harvesting complex 2 (LH2) reconstituted into two-dimensional crystals. These photosynthetic trans-membrane proteins formed cylindrical oligomeric complexes, which inserted tilted into the lipid membrane. This peculiar packing of an integral membrane protein allowed us to determine oligomerization and tilt of the LH2 complexes, but also protrusion height and intrinsic flexibility of their individual subunits. Furthermore the surface contouring reliability and limits of the atomic force microscopy could be studied. The two-dimensional crystals examined had sizes of up to 5 microm and, as revealed by a 10 A cryo electron microscopy projection map, p22(1)2(1) crystal symmetry. The unit cell had dimensions of a = b = 150 A and gamma = 90 degrees, and housed four nonameric complexes, two pointing up and two pointing down. AFM topographs of these 2D crystals had a lateral resolution of 10 A. Further, the high vertical resolution of approximately 1 A, allowed the protrusion height of the cylindrical LH2 complexes over the membrane to be determined. This was maximally 13.1 A on one side and 3.8 A on the other. Interestingly, the protrusion height varied across the LH2 complexes, showing the complexes to be inserted with a 6.2 degree tilt with respect to the membrane plane. A detailed analysis of the individual subunits showed the intrinsic flexibility of the membrane protruding peptide stretches to be equal and independent of their protrusion height. Furthermore, our analysis of membrane proteins within this peculiar packing confirmed the high vertical resolution of the atomic force microscopy on biological samples, and led us to conclude that the image acquisition function was equally accurate for contouring protrusions with heights up to approximately 15 A.     

Automatic alignment of transmission electron microscope tilt series without fiducial markers

Accurate image alignment is needed for computing three-dimensional reconstructions from transmission electron microscope tilt series. So far, the best results have been obtained by using colloidal gold beads as fiducial markers. If their use has not been possible for some reason, the only option has been the automatic cross-correlation-based registration methods. However, the latter methods are inaccurate and, as we will show, inappropriate for the whole problem. Conversely, we propose a novel method that uses the actual 3D motion model but works without any fiducial markers in the images. The method is based on matching and tracking some interest points of the intensity surface by first solving the underlying geometrical constraint of consecutive images in the tilt series. The results show that our method is near the gold marker alignment in the level of accuracy and hence opens the way for new opportunities in the analysis of electron tomography reconstructions, especially when markers cannot be used.     

Structure of HIV-1 Capsid Assemblies by Cryo-electron Microscopy and Iterative Helical Real-space Reconstruction

Cryo-electron microscopy (cryo-EM), combined with image processing, is an increasingly powerful tool for structure determination of macromolecular protein complexes and assemblies. In fact, single particle electron microscopy¹ and two-dimensional (2D) electron crystallography² have become relatively routine methodologies and a large number of structures have been solved using these methods. At the same time, image processing and three-dimensional (3D) reconstruction of helical objects has rapidly developed, especially, the iterative helical real-space reconstruction (IHRSR) method³, which uses single particle analysis tools in conjunction with helical symmetry. Many biological entities function in filamentous or helical forms, including actin filaments⁴, microtubules⁵, amyloid fibers⁶, tobacco mosaic viruses⁷, and bacteria flagella⁸, and, because a 3D density map of a helical entity can be attained from a single projection image, compared to the many images required for 3D reconstruction of a non-helical object, with the IHRSR method, structural analysis of such flexible and disordered helical assemblies is now attainable.In this video article, we provide detailed protocols for obtaining a 3D density map of a helical protein assembly (HIV-1 capsid⁹ is our example), including protocols for cryo-EM specimen preparation, low dose data collection by cryo-EM, indexing of helical diffraction patterns, and image processing and 3D reconstruction using IHRSR. Compared to other techniques, cryo-EM offers optimal specimen preservation under near native conditions. Samples are embedded in a thin layer of vitreous ice, by rapid freezing, and imaged in electron microscopes at liquid nitrogen temperature, under low dose conditions to minimize the radiation damage. Sample images are obtained under near native conditions at the expense of low signal and low contrast in the recorded micrographs. Fortunately, the process of helical reconstruction has largely been automated, with the exception of indexing the helical diffraction pattern. Here, we describe an approach to index helical structure and determine helical symmetries (helical parameters) from digitized micrographs, an essential step for 3D helical reconstruction. Briefly, we obtain an initial 3D density map by applying the IHRSR method. This initial map is then iteratively refined by introducing constraints for the alignment parameters of each segment, thus controlling their degrees of freedom. Further improvement is achieved by correcting for the contrast transfer function (CTF) of the electron microscope (amplitude and phase correction) and by optimizing the helical symmetry of the assembly.