A binary segmentation approach for boxing ribosome particles in cryo EM micrographs

Three-dimensional reconstruction of ribosome particles from electron micrographs requires selection of many single-particle images. Roughly 100,000 particles are required to achieve approximately 10 A resolution. Manual selection of particles, by visual observation of the micrographs on a computer screen, is recognized as a bottleneck in automated single-particle reconstruction. This paper describes an efficient approach for automated boxing of ribosome particles in micrographs. Use of a fast, anisotropic non-linear reaction-diffusion method to pre-process micrographs and rank-leveling to enhance the contrast between particles and the background, followed by binary and morphological segmentation constitute the core of this technique. Modifying the shape of the particles to facilitate segmentation of individual particles within clusters and boxing the isolated particles is successfully attempted. Tests on a limited number of micrographs have shown that over 80% success is achieved in automatic particle picking.     

Nucleosome mass distribution using image averaging

Scanning transmission electron microscopy has been used to obtain mass distribution data for unstained, randomly oriented dinucleosomes. These data were then subjected to computer analysis to increase the signal/noise ratio, and define the major sources of image variation within the population. Average images of groups of similar particles have a mass distribution consistent with current models of nucleosome structure, and a resolution approaching 3 nm. The potential value of this method for analyzing structural variants of nucleosomes is discussed.     

Conical tomography of freeze-fracture replicas: a method for the study of integral membrane proteins inserted in phospholipid bilayers

We have used conical tomography to study the structure of integral proteins in their phospholipid bilayer environments. Complete conical series were collected from replicas of the water channel aquaporin-0 (AQP0), a 6.6 nm side tetramer with a molecular weight of approximately 120 kDa that was purified and reconstituted in liposomes. The replicas were tilted at 38 degrees , 50 degrees or 55 degrees and rotated by 2.5 degrees , 4 degrees , or 5 degrees increments until completing 360 degrees turns. The elliptical paths of between 6 and 12 freeze-fracture particles aligned the images to a common coordinate system. Using the weighted back projection algorithm, small volumes of the replicas were independently reconstructed to reconstitute the field. Using the Fourier Shell Correlation computed from reconstructions of even and odd projections of the series, we estimated a resolution of 2-3 nm, a value that was close to the thickness of the replica (approximately 1.5 nm). The 3D reconstructions exhibited isotropic resolution along the x-y plane, which simplified the analysis of particles oriented randomly in the membrane plane. In contrast to reconstructions from single particles imaged using random conical tilt [J. Mol. Biol. 325 (2003) 210], the reconstructions using conical tomography allowed the size and shape of individual particles representing the AQP0 channel to be identified without averaging or imposing symmetry. In conclusion, the reconstruction of freeze-fracture replicas with electron tomography has provided a novel experimental approach for the study of integral proteins inserted in phospholipid bilayers.     

Studies of the structure of the T4 bacteriophage tail sheath: I. The recovery of three-dimensional structural information from the extended sheath

The extended tail sheath of bacteriophage T4 has been used to study the transfer of information from an electron micrograph to the three-dimensional reconstruction obtained from it. Two methods have been developed to assess micrograph images of helical particles and their reconstructions. First, a filter has been designed which eliminates all structure in the image inconsistent with the symmetry and assumed radius of the helical particle. Individual micrographs can therefore be assessed with respect to their consistency with the assumed symmetry and radius, before reconstruction. Second, a map of the root-meansquare deviation of individual reconstructions from their average provides a quantitative measure of the consistency of the individual sets of tail data and allows the regions in the average reconstruction which are most sensitive to differences between the particles to be identified.The averaged reconstruction is used to examine the problems related to resolution and reproducibility of the structural information and to define the extent of the different components of the extended sheath.     

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.     

Fourier amplitude decay of electron cryomicroscopic images of single particles and effects on structure determination

Several factors, including spatial and temporal coherence of the electron microscope, specimen movement, recording medium, and scanner optics, contribute to the decay of the measured Fourier amplitude in electron image intensities. We approximate the combination of these factors as a single Gaussian envelope function, the width of which is described by a single experimental B-factor. We present an improved method for estimating this B-factor from individual micrographs by combining the use of X-ray solution scattering and numerical fitting to the average power spectrum of particle images. A statistical estimation from over 200 micrographs of herpes simplex virus type-1 capsids was used to estimate the spread in the experimental B-factor of the data set. The B-factor is experimentally shown to be dependent on the objective lens defocus setting of the microscope. The average B-factor, the X-ray scattering intensity of the specimen, and the number of particles required to determine the structure at a lower resolution can be used to estimate the minimum fold increase in the number of particles that would be required to extend a single particle reconstruction to a specified higher resolution. We conclude that microscope and imaging improvements to reduce the experimental B-factor will be critical for obtaining an atomic resolution structure.     

The structure of the F-actin filament and the actin molecule.

A consensus view on the three-dimensional structure of the F-actin filament and the relative strength of the intersubunit contacts in the filament has been established from an atomic filament model and recent three-dimensional reconstructions from electron micrographs of F-actin filaments. Functional implications of recent structural and biochemical data indicating a rather dynamic filament structure are discussed.     

A two step approach for semi-automated particle selection from low contrast cryo-electron micrographs

Over recent years advances in cryo-electron microscopy for the study of macromolecular structure have resulted in resolutions in the range 10-15 A becoming routine. With this drive for increased resolution comes the need to collect larger datasets, commonly >10,000 particle images. Manual selection of particles from micrographs is often difficult and with such large numbers of particles now involved it is also laborious and a common bottleneck. Automated methods do exist but are normally restricted to specific samples or data, i.e., spherical particles, no aggregation, high contrast, and low noise. A two step approach has been developed that remains general and can be applied to low contrast, high noise micrographs of small molecules. Specifically, application of the approach is presented using micrographs of Escherichia coli RNA polymerase, which due to low contrast and the relatively small size of the molecule prove difficult to pick manually. To test the automated approach, independent reconstructions of RNA polymerase were carried out using manual and automatically picked data. The two reconstructions are shown to be comparable and the reconstruction from the automatically picked dataset is at a higher resolution, due to an increase in the number of particles picked.     

Computer-aided three-dimensional reconstruction of nematode embryos from EM serial sections

Embryos of the nematode Caenorhabditis elegans were serially sectioned and photographed in the electron microscope (EM). The micrographs were used to produce three-dimensional (3D) reconstructions. Size and position of each nucleus were entered into a computer, displayed as spheres, and were color-coded to indicate lineage membership. Location in space and position in the cell cycle are generally adequate criteria to identify cells. The reconstructions allow visualization of lineage-related topographic patterns and ultrastructural analysis of differently determined cells.     

Electron microscopic visualisation of the 5S rRNA-YL3 complex from Saccharomyces cerevisiae

The complex comprising 5S ribosomal RNA and the ribosomal protein YL3 (5S rRNP) was isolated from yeast (Saccharomyces cerevisiae), and positively contrasted preparations were imaged by transmission electron microscopy. The overall dimensions of the 5S rRNP complex in the micrographs were 10 nm by 6 min. Three predominant projections were selected from several hundred putative particles for digitisation and computer averaging to yield two-dimensional constructions with reproducible spatial resolutions exceeding 2 run. The enhanced projection images were compatible with structural models of this complex based on biochemical studies.     

Toward understanding the structure and interactions of microtubules and motor proteins

To obtain an overall three-dimensional picture of the interaction between microtubules and the motor proteins of the kinesin family it will be necessary to take account of both atomic resolution structures obtained by X-ray crystallography and medium resolution reconstructions obtained by electron cryomicroscopy. We examine the problems associated with obtaining the required structural information from electron micrographs of vitreous ice-embedded microtubules decorated with motor domains. We find that the minus-end directed motor, ncd, decorates microtubules with an 80 Å periodicity as for kinesin. Our theoretical analysis and experiments with ncd illustrate the difficulty in determining unambiguously the surface lattice organization by diffraction analysis of micrographs. 3D reconstructions of decorated microtubules are required to accurately locate the motor domains. Helical diffraction theory is not usually applicable because microtubules are cylindrical structures that rarely have complete helical symmetry. We propose using a back-projection method based on the long pitch helices formed by individual protofilaments. Model reconstructions show that this approach is feasible. © 1995 Wiley-Liss, Inc.     

Confocal scanning light microscopy of the Escherichia coli nucleoid: comparison with phase-contrast and electron microscope images.

The nucleoid of living and OsO4- or glutaraldehyde-fixed cells of Escherichia coli strains was studied with a phase-contrast microscope, a confocal scanning light microscope, and an electron microscope. The trustworthiness of the images obtained with the confocal scanning light microscope was investigated by comparison with phase-contrast micrographs and reconstructions based on serially sectioned material of DNA-containing and DNA-less cells. This comparison showed higher resolution of the confocal scanning light microscope as compared with the phase-contrast microscope, and agreement with results obtained with the electron microscope. The effects of fixation on the structure of the nucleoid were studied in E. coli B/r H266. Confocal scanning light micrographs and electron microscopic reconstructions showed that the shape of the nucleoid remained similar after OsO4 or glutaraldehyde fixation; however, the OsO4 nucleoid appeared to be somewhat smaller and more centralized within the cell.     

Deducing fast electron density changes in randomly orientated uncrystallized biomolecules in a pump–probe experiment

We propose a method for deducing time-resolved structural changes in uncrystallized biomolecules in solution. The method relies on measuring the angular correlations of the intensities, when averaged over a large number of diffraction patterns from randomly oriented biomolecules in solution in a liquid solvent. The experiment is somewhat like a pump–probe version of an experiment on small angle X-ray scattering, except that the data expected by the algorithm are not just the radial variation of the averaged intensities. The differences of these correlation functions as measured from a photoexcited and dark structure enable the direct calculation of the difference electron density with a knowledge of only the dark structure. We exploit a linear relation we derive between the difference in these correlation functions and the difference electron density, applicable for small structural changes.     

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.     

A two-exposure technique for ice-embedded samples successfully reconstructs the chlorocruorin pigment of Sabella spallanzanii at 2. 1 Nm resolution.

A technique for reconstructing ice-embedded macromolecules from electron micrographs taken at two specimen tilts (+/-23 degrees ) has been used to determine the structure of chlorocruorin isolated from the Polychaete annelid Sabella spallanzanii. Images of individual molecules were extracted in couples from two micrographs of the same field of view so each couple consists of two projections of the same molecule. One couple was used as a fixed reference for alignment. Different references yielded reconstructions with different orientations. These were merged to give a model against which the orientation of 1624 first-exposure images was refined to give a final reconstruction at 2.1 nm resolution. The structure of this hematic pigment, essentially the same as that for Lumbricus terrestris, is a bilayer structure with overall symmetry D6, containing six hollow groups per layer. A hollow group is formed by six globular masses and has approximate threefold symmetry. Other structural elements connect the two layers and the hollow groups in a layer. This non-globin material occupies about 15% of the total molecular volume. The results show that the double-exposure strategy, previously described by some of the authors and tested in computer simulations, performs well in real experiments and could be used to obtain preliminary reconstructions in a semiautomatic way.     

Structural analysis of tropomyosin tactoids.

A refined computer graphics approach to correlation of molecular sequence with electron micrographs of tropomyosin tactoids is presented. It is shown that antiparallel structures with molecular chains in phase, 21 or 14 residue overlap and C or N terminal overlap are consistent with the morphology. The C terminal overlap structure previously postulated gives the best direct correlation but chemical evidence appears to support the N terminal overlap structure. The parallel tactoid form appears more complicated and no adequate structure has yet been elucidated.     

Reconstructive three-dimensional electron microscopy. A routine biologic tool

Twenty years ago a laboratory could devote an entire year or more to the collection and analysis of a single set of serial electron micrographs. In contrast, simple technical improvements have now made it possible to take embedded material and have in hand complete computer reconstructions of cells' organelles, microtubules, etc., in less than a week. With a few additional minor improvements, this time could be reduced to only two or three days. Experience in our laboratory suggests that almost without exception these reconstructions provide new insights into both the structure and function of cells. We illustrate this point by presenting a new, unpublished anatomic feature of mammalian nuclei, the "nuclear tube." This example is typical of many other unpublished incidental findings we have made over the last five years using serial electron microscopy as a routine tool, and we believe it represents only the tip of a largely unexplored world of three-dimensional cytoarchitecture.     

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.     

ResLog plots as an empirical metric of the quality of cryo-EM reconstructions

Compared to the field of X-ray crystallography, the field of single particle three-dimensional electron microscopy has few reliable metrics for assessing the quality of 3D reconstructions. New metrics are needed that can determine whether a given 3D reconstruction accurately reflects the structure of the particles from which it was derived or instead depicts a plausible though incorrect structure due to coarse misalignment of particles. Here an empirical procedure is presented for differentiating between a reconstruction with well-aligned particles and a reconstruction with grossly misclassified particles. For a given dataset, 3D reconstructions are computed from subsets of particles with decreasing numbers of particles contributing to the reconstruction. A plot of inverse resolution vs. the logarithm of the number of particles (a “ResLog” plot) provides metrics for the reliability of the reconstruction and the overall quality of the dataset and processing. Specifically, the y-intercept of a regression line provides a measure of the relative accuracy of the particle alignment and classification, and the slope is an indicator of the overall data quality including the imaging conditions and processing steps. ResLog plots can also be used to optimize conditions for data collection and reconstruction parameters. Although resolution estimates can vary by method of calculation, ResLog-derived parameters are consistent whether calculated by Fourier shell correlation or Fourier neighbor correlation, or a new coordinate-based metric that serves as a yardstick for structures where atomic coordinates are available. ResLog plots could become part of a standard set of parameters to be included in 3D reconstruction reports.