Towards atomic resolution structural determination by single-particle cryo-electron microscopy

Recent advances in cryo-electron microscopy and single-particle reconstruction (collectively referred to as 'cryoEM') have made it possible to determine the three-dimensional (3D) structures of several macromolecular complexes at near-atomic resolution ( approximately 3.8-4.5A). These achievements were accomplished by overcoming the challenges in sample handling, instrumentation, image processing, and model building. At near-atomic resolution, many detailed structural features can be resolved, such as the turns and deep grooves of helices, strand separation in beta sheets, and densities for loops and bulky amino acid side chains. Such structural data of the cytoplasmic polyhedrosis virus (CPV), the Epsilon 15 bacteriophage and the GroEL complex have provided valuable constraints for atomic model building using integrative tools, thus significantly enhancing the value of the cryoEM structures. The CPV structure revealed a drastic conformational change from a helix to a beta hairpin associated with RNA packaging and replication, coupling of RNA processing and release, and the long sought-after polyhedrin-binding domain. These latest advances in single-particle cryoEM provide exciting opportunities for the 3D structural determination of viruses and macromolecular complexes that are either too large or too heterogeneous to be investigated by conventional X-ray crystallography or nuclear magnetic resonance (NMR) methods.     

Do's and Don'ts of Cryo-electron Microscopy: A Primer on Sample Preparation and High Quality Data Collection for Macromolecular 3D Reconstruction

Cryo-electron microscopy (cryoEM) entails flash-freezing a thin layer of sample on a support, and then visualizing the sample in its frozen hydrated state by transmission electron microscopy (TEM). This can be achieved with very low quantity of protein and in the buffer of choice, without the use of any stain, which is very useful to determine structure-function correlations of macromolecules. When combined with single-particle image processing, the technique has found widespread usefulness for 3D structural determination of purified macromolecules. The protocol presented here explains how to perform cryoEM and examines the causes of most commonly encountered problems for rational troubleshooting; following all these steps should lead to acquisition of high quality cryoEM images. The technique requires access to the electron microscope instrument and to a vitrification device. Knowledge of the 3D reconstruction concepts and software is also needed for computerized image processing. Importantly, high quality results depend on finding the right purification conditions leading to a uniform population of structurally intact macromolecules. The ability of cryoEM to visualize macromolecules combined with the versatility of single particle image processing has proven very successful for structural determination of large proteins and macromolecular machines in their near-native state, identification of their multiple components by 3D difference mapping, and creation of pseudo-atomic structures by docking of x-ray structures. The relentless development of cryoEM instrumentation and image processing techniques for the last 30 years has resulted in the possibility to generate de novo 3D reconstructions at atomic resolution level.     

SPRING – An image processing package for single-particle based helical reconstruction from electron cryomicrographs

Helical reconstruction from electron cryomicrographs has become a routine technique for macromolecular structure determination of helical assemblies since the first days of Fourier-based three-dimensional image reconstruction. In the past decade, the single-particle technique has had an important impact on the advancement of helical reconstruction. Here, we present the software package SPRING that combines Fourier based symmetry analysis and real-space helical processing into a single workflow. One of the most time-consuming steps in helical reconstruction is the determination of the initial symmetry parameters. First, we propose a class-based helical reconstruction approach that enables the simultaneous exploration and evaluation of many symmetry combinations at low resolution. Second, multiple symmetry solutions can be further assessed and refined by single-particle based helical reconstruction using the correlation of simulated and experimental power spectra. Finally, the 3D structure can be determined to high resolution. In order to validate the procedure, we use the reference specimen Tobacco Mosaic Virus (TMV). After refinement of the helical symmetry, a total of 50,000 asymmetric units from two micrographs are sufficient to reconstruct a subnanometer 3D structure of TMV at 6.4 Å resolution. Furthermore, we introduce the individual programs of the software and discuss enhancements of the helical reconstruction workflow. Thanks to its user-friendly interface and documentation, SPRING can be utilized by the novice as well as the expert user. In addition to the study of well-ordered helical structures, the development of a streamlined workflow for single-particle based helical reconstruction opens new possibilities to analyze specimens that are heterogeneous in symmetries.     

Single particle macromolecular structure determination via electron microscopy

Three-dimensional structure determination of macromolecules and macromolecular complexes is an integral part of understanding biological functions. For large protein and macromolecular complexes structure determination is often performed using electron cryomicroscopy where projection images of individual macromolecular complexes are combined to produce a three-dimensional reconstruction. Single particle methods have been devised to perform this structure determination for macromolecular complexes with little or no underlying symmetry. These computational methods generally involve an iterative process of aligning unique views of the macromolecular images followed by determination of the angular components that define those views. In this review, this structure determination process is described with the aim of clarifying a seemingly complex structural method.     

Chemically functionalized carbon films for single molecule imaging

Many biological complexes are naturally low in abundance and pose a significant challenge to their structural and functional studies. Here we describe a new method that utilizes strong oxidation and chemical linkage to introduce a high density of bioactive ligands onto nanometer-thick carbon films and enable selective enrichment of individual macromolecular complexes at subnanogram levels. The introduced ligands are physically separated. Ni-NTA, Protein G and DNA/RNA oligonucleotides were covalently linked to the carbon surface. They embody negligible mass and their stability makes the functionalized films able to survive long-term storage and tolerate variations in pH, temperature, salts, detergents, and solvents. We demonstrated the application of the new method to the electron microscopic imaging of the substrate-bound C3PO, an RNA-processing enzyme important for the RNA interference pathway. On the ssRNA-linked carbon surface, the formation of C3PO oligomers at subnanomolar concentrations likely mimics their assembly onto ssRNA substrates presented by their native partners. Interestingly, the 3D reconstructions by negative stain EM reveal a side port in the C3PO/ssRNA complex, and the 15 Å cryoEM map showed extra density right above the side port, which probably represents the ssRNA. These results suggest a new way for ssRNAs to interact with the active sites of the complex. Together our data demonstrate that the surface-engineered carbon films are suitable for selectively enriching low-abundance biological complexes at nanomolar level and for developing novel applications on a large number of surface-presented molecules.     

Cryoelectron microscopy reveals new features in the three-dimensional structure of phosphorylase kinase

Phosphorylase kinase (PhK), a regulatory enzyme in the cascade activation of glycogenolysis, is a 1.3-MDa hexadecameric complex, (alphabetagammadelta)(4). PhK comprises two arched octameric (alphabetagammadelta)(2) lobes that are oriented back-to-back with overall D(2) symmetry and connected by small bridges. These interlobal bridges, arguably the most questionable structural component of PhK, are one of several structural features that potentially are artifactually generated or altered by conventional sample preparation techniques for electron microscopy (EM). To minimize such artifacts, we have solved by cryoEM the first three-dimensional (3D) structure of nonactivated PhK from images of frozen hydrated molecules of the kinase. Minimal dose electron micrographs of PhK in vitreous ice revealed particles in a multitude of orientations. A simple model was used to orient the individual images for 3D reconstruction, followed by multiple rounds of refinement. Three-dimensional reconstruction of nonactivated PhK from approximately 5000 particles revealed a bridged, bilobal molecule with a resolution estimated by Fourier shell correlation analysis at 25 A. This new structure suggests that several prominent features observed in the structure of PhK derived from negatively stained particles arise as artifacts of specimen preparation. In comparison to the structure from negative staining, the cryoEM structure shows three important differences: (1) a dihedral angle between the two lobes of approximately 90 degrees instead of 68 degrees, (2) a compact rather than extended structure for the lobes, and (3) the presence of four, rather than two, connecting bridges, which provides the first direct evidence for these components as authentic elements of the kinase solution structure.     

IMIRS: a high-resolution 3D reconstruction package integrated with a relational image database

Recent advances in electron cryomicroscopy instrumentation and single particle reconstruction have created opportunities for high-throughput and high-resolution three-dimensional (3D) structure determination of macromolecular complexes. However, it has become impractical and inefficient to rely on conventional text file data management and command-line programs to organize and process the increasing numbers of image data required in high-resolution studies. Here, we present a distributed relational database for managing complex datasets and its integration into our high-resolution software package IMIRS (Image Management and Icosahedral Reconstruction System). IMIRS consists of a complete set of modular programs for icosahedral reconstruction organized under a graphical user interface and provides options for user-friendly, step-by-step data processing as well as automatic reconstruction. We show that the integration of data management with processing in IMIRS automates the tedious tasks of data management, enables data coherence, and facilitates information sharing in a distributed computer and user environment without significantly increasing the time of program execution. We demonstrate the applicability of IMIRS in icosahedral reconstruction toward high resolution by using it to obtain an 8-A 3D structure of an intermediate-sized dsRNA virus.     

Cryo-electron tomography: gaining insight into cellular processes by structural approaches

Visualization of cellular processes at a resolution of the individual protein should involve integrative and complementary approaches that can eventually draw realistic functional and cellular landscapes. Electron tomography of vitrified but otherwise unaltered cells emerges as a central method for three-dimensional reconstruction of cellular architecture at a resolution of 2-6 nm. While a combination of correlative light-based microscopy with cryo-electron tomography (cryo-ET) provides medium-resolution insight into pivotal cellular processes, fitting high-resolution structural approaches, for example, X-ray crystallography, into reconstructed macromolecular assemblies provides unprecedented information on native protein assemblies. Thus, cryo-ET bridges the resolution gap between cellular and structural biology. In this article, we focus on the study of eukaryotic cells and macromolecular complexes in a close-to-life-state. We discuss recent developments and structural findings enabling major strides to be made in understanding complex physiological functions.     

Automated segmentation of molecular subunits in electron cryomicroscopy density maps

Electron cryomicroscopy (cryoEM) is capable of imaging large macromolecular machines composed of multiple components. However, it is currently only possible to achieve moderate resolution at which it may be possible to computationally extract the individual components in the machine. In this work, we present application details of an automated method for detecting and segmenting the components of a large machine in an experimentally determined density map. This method is applicable to object with and without symmetry and takes advantage of global and local symmetry axes if present. We have applied this segmentation algorithm to several cryoEM data sets already deposited in EMDB with various complexities, symmetries and resolutions and validated the results using manually segmented density and available structures of the components in the PDB. As such, automated segmentation could become a useful tool for the analysis of the ever-increasing number of structures of macromolecular machines derived from cryoEM.     

Efficient implementation of a filtered back-projection algorithm using a voxel-by-voxel approach

The large amount of image data necessary for high-resolution 3D reconstruction of macromolecular assemblies leads to significant increases in the computational time. One of the most time consuming operations is 3D density map reconstruction, and software optimization can greatly reduce the time required for any given structural study. The majority of algorithms proposed for improving the computational effectiveness of a 3D reconstruction are based on a ray-by-ray projection of each image into the reconstructed volume. In this paper, we propose a novel fast implementation of the "filtered back-projection" algorithm based on a voxel-by-voxel principle. Our version of this implementation has been exhaustively tested using both model and real data. We compared 3D reconstructions obtained by the new approach with results obtained by the filtered Back-Projections algorithm and the Fourier-Bessel algorithm commonly used for reconstructing icosahedral viruses. These computational experiments demonstrate the robustness, reliability, and efficiency of this approach.     

Limiting factors in atomic resolution cryo electron microscopy: no simple tricks

To bring cryo electron microscopy (cryoEM) of large biological complexes to atomic resolution, several factors--in both cryoEM image acquisition and 3D reconstruction--that may be neglected at low resolution become significantly limiting. Here we present thorough analyses of four limiting factors: (a) electron-beam tilt, (b) inaccurate determination of defocus values, (c) focus gradient through particles, and (d) particularly for large particles, dynamic (multiple) scattering of electrons. We also propose strategies to cope with these factors: (a) the divergence and direction tilt components of electron-beam tilt could be reduced by maintaining parallel illumination and by using a coma-free alignment procedure, respectively. Moreover, the effect of all beam tilt components, including spiral tilt, could be eliminated by use of a spherical aberration corrector. (b) More accurate measurement of defocus value could be obtained by imaging areas adjacent to the target area at high electron dose and by measuring the image shift induced by tilting the electron beam. (c) Each known Fourier coefficient in the Fourier transform of a cryoEM image is the sum of two Fourier coefficients of the 3D structure, one on each of two curved 'characteristic surfaces' in 3D Fourier space. We describe a simple model-based iterative method that could recover these two Fourier coefficients on the two characteristic surfaces. (d) The effect of dynamic scattering could be corrected by deconvolution of a transfer function. These analyses and our proposed strategies offer useful guidance for future experimental designs targeting atomic resolution cryoEM reconstruction.     

A multiresolution approach to orientation assignment in 3D electron microscopy of single particles

Three-dimensional (3D) electron microscopy (3DEM) aims at the determination of the spatial distribution of the Coulomb potential of macromolecular complexes. The 3D reconstruction of a macromolecule using single-particle techniques involves thousands of 2D projections. One of the key parameters required to perform such a 3D reconstruction is the orientation of each projection image as well as its in-plane orientation. This information is unknown experimentally and must be determined using image-processing techniques. We propose the use of wavelets to match the experimental projections with those obtained from a reference 3D model. The wavelet decomposition of the projection images provides a framework for a multiscale matching algorithm in which speed and robustness to noise are gained. Furthermore, this multiresolution approach is combined with a novel orientation selection strategy. Results obtained from computer simulations as well as experimental data encourage the use of this approach.     

Scoring functions for cryoEM density fitting

In fitting atomic structures into cryoEM density maps of macromolecular assemblies, the cross-correlation function (CCF) is the most prevalent method of scoring the goodness-of-fit. However, there are still many possible, less studied ways of scoring fits. In this paper, we introduce four scores new to cryoEM fitting and compare their performance to three known scores. Our benchmark consists of (a) 4 protein assemblies with simulated maps at 5-20 ? resolution, including the heptameric ring of GroEL; and (b) 4 experimental maps of GroEL at ～6-23 ? resolution with corresponding fitted atomic models. We perturb each fit 1000 times and assess each new fit with each score. The correlation between a score and the Cα RMSD of each fit from the "correctly" fitted structure shows that the CCF is one of the best scores, but in certain situations could be augmented or even replaced by other scores. For instance, our implementation of a score based on mutual information outperforms or is comparable to the CCF in almost all test cases, and our new "envelope score" works as well as the CCF at sub-nanometer resolution but is an order of magnitude faster to calculate. The results also suggest that the width of the Gaussian function used to blur the atomic structure into a density map can significantly affect the fitting process. Finally, we show that our score-testing method, when combined with the Laplacian CCF or the mutual information scores, can be used as a statistical tool for improving cryoEM density fitting.     

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.     

SwarmPS: rapid, semi-automated single particle selection software

Single particle analysis (SPA) coupled with high-resolution electron cryo-microscopy is emerging as a powerful technique for the structure determination of membrane protein complexes and soluble macromolecular assemblies. Current estimates suggest that approximately 10(4)-10(5) particle projections are required to attain a 3A resolution 3D reconstruction (symmetry dependent). Selecting this number of molecular projections differing in size, shape and symmetry is a rate-limiting step for the automation of 3D image reconstruction. Here, we present Swarm(PS), a feature rich GUI based software package to manage large scale, semi-automated particle picking projects. The software provides cross-correlation and edge-detection algorithms. Algorithm-specific parameters are transparently and automatically determined through user interaction with the image, rather than by trial and error. Other features include multiple image handling (approximately 10(2)), local and global particle selection options, interactive image freezing, automatic particle centering, and full manual override to correct false positives and negatives. Swarm(PS) is user friendly, flexible, extensible, fast, and capable of exporting boxed out projection images, or particle coordinates, compatible with downstream image processing suites.     

Visualization of macromolecular complexes using cryo-electron microscopy with FEI Tecnai transmission electron microscopes

This protocol details the steps used for visualizing the frozen-hydrated grids as prepared following the accompanying protocol entitled 'Preparation of macromolecular complexes for visualization using cryo-electron microscopy.' This protocol describes how to transfer the grid to the microscope using a standard cryo-transfer holder or, alternatively, using a cryo-cartridge loading system, and how to collect low-dose data using an FEI Tecnai transmission electron microscope. This protocol also summarizes and compares the various options that are available in data collection for three-dimensional (3D) single-particle reconstruction. These options include microscope settings, choice of detectors and data collection strategies both in situations where a 3D reference is available and in the absence of such a reference (random-conical and common lines).     

A 11.5 A single particle reconstruction of GroEL using EMAN.

Single-particle analysis has become an increasingly important method for structural determination of large macromolecular assemblies. GroEL is an 800 kDa molecular chaperone, which, along with its co-chaperonin GroES, promotes protein folding both in vitro and in the bacterial cell. EMAN is a single-particle analysis software package, which was first publicly distributed in 2000. We present a three-dimensional reconstruction of native naked GroEL to approximately 11.5 A performed entirely with EMAN. We demonstrate that the single-particle reconstruction, X-ray scattering data and X-ray crystal structure all agree well at this resolution. These results validate the specific methods of image restoration, reconstruction and evaluation techniques implemented in EMAN. It also demonstrates that the single-particle reconstruction technique and X-ray crystallography will yield consistent structure factors, even at low resolution, when image restoration is performed correctly. A detailed comparison of the single-particle and X-ray structures exhibits some small variations in the equatorial domain of the molecule, likely due to the absence of crystal packing forces in the single-particle reconstruction.     

A measure for the angle between projections based on the extent of correlation between corresponding central sections

A pre-condition for the ab initio assignment of Euler angles to a set of projections from an asymmetric object is that at least three of the available projections correspond to rotations about different axes. For symmetric objects this condition may be relaxed. There are some applications of single-particle electron microscopy, such as the reconstruction of filamentous macromolecular assemblies, where all available projections more-or-less correspond to rotations about a common rotation axis making it difficult to satisfy this condition. Here, a method has been developed to overcome this problem, based on the fact that the correlation between two central sections of the Fourier transform of a compact object will not be limited to an infinitesimal central line but will have a finite extent, which is related to the angle between the corresponding projections. Projections from model filaments, with different degrees of rotational symmetry about the long axis, have been used to test the methodology. The results show that angle determination is robust down to signal-to-noise ratios as low as 2 and that, in general, the error decreases as the degree of symmetry increases. The method has been used to assign angles to a set of negatively stained muscle thick filament projections to obtain an initial 3D reconstruction. The main features of the projections are seen to be faithfully reproduced in the reprojections from the reconstruction. A real-space adaptation of this method is also discussed.     

Automated three-dimensional reconstruction of keyhole limpet hemocyanin type 1

We have reconstructed a three-dimensional map of keyhole limpet hemocyanin isoform 1 (KLH1), using our automated data collection software, Leginon, integrated with particle selection algorithms, and the SPIDER reconstruction package. KLH1, a 7.9 MDa macromolecule, is an extracellular respiratory pigment composed of two asymmetric decamers, and presents an overall D(5) point-group symmetry. The reconstruction is in agreement with previous data published on molluscan hemocyanins. The reconstructed map (11.3A resolution, 3sigma criterion) was used to fit an available X-ray crystallography structure of Octopus dofleini Odg, solved at 2.3A [J. Mol. Biol. 278 (4) (1998) 855], with satisfactory results. The results validate the approach of automating the cryoEM process and demonstrate that the quality of the images acquired and the particles selected is comparable to those obtained using manual methods. Several problems remain to be solved however before these results can be generalized.