Three-dimensional structure of bovine NADH:ubiquinone oxidoreductase of the mitochondrial respiratory chain

We have studied the structure of bovine heart mitochondrial NADH:ubiquinone (Q) oxidoreductase (EC 1.6.99.3) by image analysis of electron micrographs. A three-dimensional reconstruction was calculated from a tilt-series of a two-dimensional crystal of the molecule. Our interpretation of the position of the molecule in the unit cell of the crystal is supported by additional (low-resolution) analysis of images of single molecules. The three-dimensional reconstruction was calculated with the aid of an iterative real-space reconstruction algorithm. The various projections used as input to the algorithm were obtained by averaging the images of the tilted crystal through a Fourier-space peak-filtering procedure. The reconstructed unit cell measures 15.2 X 15.2 nm in the plane of the two-dimensional crystal and has a height of 10-11 nm. The unit cell contains one molecule consisting of four large subunits. At the present resolution of about 1.3 nm in the untilted projection, these four monomers are seen as two dimers related by a two-fold axis. Two views of the single particles have been recognized; they are the top and side view of the building block of the crystal. After computer image alignment and correspondence analysis, clusters of similar particles have been averaged. In the averages an uneven stain distribution is seen around the molecules, which may result from preferential staining of hydrophilic parts of the molecule. The molecular mass of the whole molecule was determined from scanning transmission electron microscopy measurements as (1.6 +/- 0.2) X 10(6) daltons.     

Model for the structure of bacteriorhodopsin based on high-resolution electron cryo-microscopy

The light-driven proton pump bacteriorhodopsin occurs naturally as two-dimensional crystals. A three-dimensional density map of the structure, at near-atomic resolution, has been obtained by studying the crystals using electron cryo-microscopy to obtain electron diffraction patterns and high-resolution micrographs. New methods were developed for analysing micrographs from tilted specimens, incorporating methods previously developed for untilted specimens that enable large areas to be analysed and corrected for distortions. Data from 72 images, from both tilted and untilted specimens, were analysed to produce the phases of 2700 independent Fourier components of the structure. The amplitudes of these components were accurately measured from 150 diffraction patterns. Together, these data represent about half of the full three-dimensional transform to 3.5 A. The map of the structure has a resolution of 3.5 A in a direction parallel to the membrane plane but lower than this in the perpendicular direction. It shows many features in the density that are resolved from the main density of the seven alpha-helices. We interpret these features as the bulky aromatic side-chains of phenylalanine, tyrosine and tryptophan residues. There is also a very dense feature, which is the beta-ionone ring of the retinal chromophore. Using these bulky side-chains as guide points and taking account of bulges in the helices that indicate smaller side-chains such as leucine, a complete atomic model for bacteriorhodopsin between amino acid residues 8 and 225 has been built. There are 21 amino acid residues, contributed by all seven helices, surrounding the retinal and 26 residues, contributed by five helices, forming the proton pathway or channel. Ten of the amino acid residues in the middle of the proton channel are also part of the retinal binding site. The model also provides a useful basis for consideration of the mechanism of proton pumping and allows a consistent interpretation of a great deal of other experimental data. In particular, the structure suggests that pK changes in the Schiff base must act as the means by which light energy is converted into proton pumping pressure in the channel. Asp96 is on the pathway from the cytoplasm to the Schiff base and Asp85 is on the pathway from the Schiff base to the extracellular surface.     

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.     

High-resolution spot-scan electron microscopy of microcrystals of an alpha-helical coiled-coil protein

We describe the electron microscopy of a crystalline assembly of an alpha-helical coiled-coil protein extracted from the ootheca of the praying mantis. Electron diffraction patterns of unstained crystals show crystal lattice sampling of the coiled-coil molecular transform to a resolution beyond 1.5 A. Using a "spot-scan" method of electron imaging, micrographs of unstained crystals have been obtained that visibly diffract laser light from crystal spacings as small as 4.3 A. A projection map was calculated to 4 A using electron diffraction amplitudes and phases from computer-processed images. The projection map clearly shows modulations in density arising from the 5.1 A alpha-helical repeat, the first time this type of modulation has been revealed by electron microscopy. The crystals have p2 plane group symmetry with a = 92.4 A, b = 150.7 A, y = 92.4 degrees. Examination of tilted specimens shows that c is approximately 18 A, indicating that the unit cell is only one molecule thick. A preliminary interpretation shows tightly packed molecules some 400 A long lying with their long axes in the plane of the projection. The molecules have a coiled-coil configuration for most of their length. The possible modes of packing of the molecules in three dimensions are discussed.     

A comparison of negatively stained electron micrographs and projections obtained from single crystal X-ray studies: The lipovitellin complex from lamprey

The analysis of single crystals of the lipovitellin complex from lamprey has made it possible to compare electron density projections derived from X-ray diffraction and density modification methods with previously published electron micrographs. The close correlation between the images obtained by the two methods demonstrates that the fidelity of images obtained by electron microscopy is excellent despite the loss of specimen order. These correlations also attest to the ability of density modification methods, currently used in macromolecular crystallography, to estimate the phases of small angle X-ray reflections.     

Spatial disorders and computational cures

Image averaging provides a powerful method for enhancing the yield of interpretable information from electron micrographs of biological macromolecules. However, as originally conceived, the full benefit of averaging is achieved only with perfectly ordered two-dimensional crystals. More recent developments, reviewed here, allow one to rectify disordered lattices, straighten randomly bent filaments, and combine multiple images of free-standing particles, thus extending the advantages of image averaging to virtually every class of macromolecular specimen.     

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

Focus gradient correction applied to tilt series image data used in electron tomography

The resolution in 3D reconstructions from tilt series is limited to the information below the first zero of the contrast transfer function unless the signal is corrected computationally. The restoration is usually based on the assumption of a linear space-invariant system and a linear relationship between object mass density and observed image contrast. The space-invariant model is no longer valid when applied to tilted micrographs because the defocus varies in a direction perpendicular to the tilt axis and with it the shape of the associated point spread function. In this paper, a method is presented for determining the defocus gradient in thin specimens such as sections and 2D crystals, and for restoration of the images subsequently used for 3D reconstruction. The alignment procedure for 3D reconstruction includes area matching and tilt geometry refinement. A map with limited resolution computed from uncorrected micrographs is compared to a volume computed from corrected micrographs with extended resolution.     

Three-dimensional structure of T3 connector purified from overexpressing bacteria.

The bacteriophage T3 connector has been purified from overexpressed protein in Escherichia coli, harboring a plasmid containing the gene encoding p8 protein. The connector, which is composed of 12 copies of p8, has been crystallized in two-dimensional sheets and studied by electron microscopy from negatively stained specimens. A two-dimensional Fourier filtering and averaging procedure was performed with crystalline specimens. In addition, single particle averaging techniques were used with other preparations. The average images obtained from these two approaches gave similar results. A three-dimensional reconstruction from two-dimensional crystals of T3 connectors was obtained by collecting several sets of tilted views and using standard Fourier procedures. The resolution of the three-dimensional map was 1.65 nm. The reconstructed connector shows two main domains: a wider one with 12 small units in the periphery and with an external diameter of 14.9 nm, and a smaller one with 8.5 nm diameter. The height of the reconstructed connector has been determined to be around 8.5 nm. The reconstruction clearly shows an internal open channel running along the longitudinal axis of the particle and having an average diameter of 3.7 nm.     

Structure of light-harvesting chlorophyll a/b protein complex from plant photosynthetic membranes at 7 A resolution in projection

The structure of thin three-dimensional crystals of the light-harvesting chlorophyll a/b protein complex, an integral membrane protein from the photosynthetic membrane of chloroplasts, has been determined at 7 A (1 A = 0.1 nm) resolution in projection. The structure analysis was carried out by image processing of low-dose electron micrographs, and electron diffraction of thin three-dimensional crystals preserved in tannin. The three-dimensional crystals appeared to be stacks of two-dimensional crystals having p321 symmetry. Results of the image analysis indicated that the crystals were disordered, due to random translational displacement of stacked layers. This was established by a translation search routine that used the low-resolution projection of a single layer as a reference. The reference map was derived from the symmetrized average of two images that showed features consistent with the projected structure of negatively stained two-dimensional crystals. The phase shift resulting from the displacement of each layer was corrected. Phase shifts were then refined by minimizing the phase residual, bringing all layers to the same phase origin. Refined phases from different images were in agreement and reliable to 7 A resolution. A projection map was generated from the averaged phases and electron diffraction amplitudes. The map showed that the complex was a trimer composed of three protein monomers related by 3-fold symmetry. The projected density within the protein monomer suggested membrane-spanning alpha-helices roughly perpendicular to the crystal plane. The density in the centre and on the periphery of the trimeric complex was lower than that of the protein, indicating that this region contained low-density matter, such as lipids and antenna chlorophylls.     

Electron crystallography of macromolecular periodic arrays on phospholipid monolayers

Electron crystallography has the potential of yielding structural information equivalent to x-ray diffraction. The major difficulty has been preparing specimens with the required structural order and size for diffraction and imaging in the electron microscope. 2D crystallization on phospholipid monolayers is capable of fulfilling both of these requirements. Crystals can form as a result of specific interactions with a protein's ligand or an analog, suitably linked to a lipid tail; or on a surface of complementary head-group charge. With such choices, the availability of a suitable lipid is limited only by synthetic chemistry. Ultimately, it is the quality and regularity of the protein-protein interactions that determine the crystalline order, as it is with any protein crystal. In the case of streptavidin, the monolayer crystal diffracts beyond 2.5 Å. A 3 Å projection map reconstructed from electron diffraction amplitudes and phases from images shows density which can be interpreted as β-sheets and clusters of side chains. It remains to be shown that the monolayer crystals are flat and diffract as well at high tilt angle as untilted. Technological issues such as charging must be resolved. With parallel advances in data collection and processing, electron crystallography of monolayer macromolecular crystals will eventually take its place beside x-ray crystallography and NMR as a routine and efficient structural technique.     

Crystals of Lumbricus erythrocruorin

Lumbricus terrestris erythrocruorin, a 3.9 X 10(6) Mr respiratory protein, has been crystallized in four different forms. Despite the high molecular symmetry apparent from images in electron micrographs, only one crystal form expresses any molecular symmetry as crystallographic symmetry. The lattice parameters provide upper limits on the molecular dimensions of 267 A X 308 A X 172 A (1 A = 0.1 nm), which agree well with dimensions obtained from electron micrographs of negatively stained molecules. We have collected diffraction data to 5.5 A from type III crystals and have begun a structural analysis.     

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.     

Three-dimensional reconstruction of tubulin in zinc-induced sheets.

The three-dimensional structure of porcine brain tubulin in planar sheets formed in the presence of zinc has been determined to a resolution of approximately 20 Å by electron microscopy and image reconstruction on negatively stained samples. The samples were prepared with a mica floatation technique, which yields tubulin sheets with 36 reciprocal space maxima on lattice lines at 21, 28, 42 and 84 Å^?1 in Fourier transforms of digitized images. In order to obtain three-dimensional data, sheets were tilted with the goniometer stage of the electron microscope to provide images at various angles between 0 ° and ± 60 °. Transforms of 33 tilted images plus the transform of untilted sheets based on an average of nine untilted images were combined to give the third dimension of reciprocal space ($\text{z}{}^1$). These data, were expressed in terms of the phases and amplitudes along the $\text{z}{}^1$ lattice line for each of the 36 maxima observed in untilted samples, as well as five additional lattice lines which have zero-amplitudes in the non-tilted central section of the three-dimensional transform. Home of these zero-amplitudes arise from systematic absences which are due to a 2-fold screw axis relating adjacent protofilaments of tubulin in the zinc-induced sheets. Thus in the three-dimensional reconstructions of the sheets a polarity of the protofilaments is apparent, with adjacent protofilaments aligned in opposite directions to give an antiparallel pattern, in contrast to normal microtubules composed of protofilaments in parallel alignment. Two classes of morphological units, each with a mass corresponding to a molecular weight of about 55,000, are found to alternate along the protofilaments. These distinct morphological units are identified as the α and β subunits of tubulin, confirming the representation of tubulin as an αβ heterodimer. Furthermore, the extensive internal contact between subunits within a dimer can readily be distinguished from the less extensive contact between dimer units. Such differences in contacts were not apparent in the earlier two-dimensional reconstructions. In addition, areas of excluded stain joining one class of subunits to the subunits of the other class in adjacent protofilaments have been resolved for tubulin polymerized in zinc-induced sheets. Of the two classes of subunits one is distinguished by a prominent cleft. Identification of which class of subunits is α and which is β is not yet possible.     

Preparation and analysis of large, flat crystals of Ca(2+)-ATPase for electron crystallography.

Obtaining large, flat, well ordered crystals represents the key to structure determination by electron crystallography. Multilamellar crystals of Ca(2+)-ATPase are a good candidate for this methodology, and we have optimized methods of crystallization and of preparation for cryoelectron microscopy. In particular, high concentrations of glycerol were found to prevent nucleation and to reduce stacking; thus, by seeding solutions containing 40% glycerol, we obtained thin crystals that were 5-30 microns in diameter and 2-10 unit cells thick. We found that removing vesicles and minimizing concentrations of divalent cations were critical to preparing flat crystals in the frozen-hydrated state. Finally, we developed two methods for determining the number of lamellae composing individual crystals, information that is required for structure determination of this crystal form. The first method, using low magnification images of freeze-dried crystals, is more practical in our case. Nevertheless, the alternative method, involving analysis of Laue zones from electron diffraction patterns of slightly tilted crystals, may be of general use in structure determination from thin, three-dimensional crystals.     

Analysis of symmetry and three-dimensional reconstruction of thin gp32*I crystals.

Thin, multilayered crystals of gp32*I were analyzed by negative stain electron microscopy and image processing. Images of untilted crystals exhibited different projection symmetries and structural motifs. Systematic analysis of these images categorized the projections into four types. Areas producing the type 1 projection were reconstructed in three-dimensions from four tilt series containing 111 images. The three-dimensional data has excellent p121 plane group symmetry and reveals that the gp32*I molecule contains two large domains linked together by a small domain. Computer simulations utilizing projection data suggested that the type 2 and 3 projections arise from superposition of type 1 projections related by a 21 screw axis along the projection axis. The three-dimensional reconstruction was utilized in a final simulation that explained the occurrence of the fourth type of projection. This work provides a firm foundation for future high-resolution analysis of the crystal by electron cryomicroscopy.     

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.     

Improved specimen preparation for cryo-electron microscopy using a symmetric carbon sandwich technique

Image shift due to beam-induced specimen charging has become the most severe problem in electron microscopy for imaging two-dimensional (2D) crystals of biological macromolecules, especially in the case of highly tilted specimens. Image shift causes diffraction spots perpendicular to the tilt axis to disappear even at medium or low resolution. The yield of good images from tilted specimens prepared on a single layer of continuous carbon support film is therefore very low. In this paper, we have used 2D crystals of aquaporin-4 to investigate the effect of a carbon sandwich preparation method on specimen charging. We find that a larger number of images show sharp diffraction spots perpendicular to the tilt axis if crystals are placed in between two sheets of carbon film as compared to images taken from specimens prepared by the conventional single carbon support film technique. Our results demonstrate that the reproducible carbon sandwich preparation technique overcomes the severe specimen charging problem and thus has the potential to significantly speed up structure analysis by electron crystallography.     

Correlation analysis of gap junction lattice images.

Fourier averages of connexon images computed from low-irradiation electron micrographs of isolated negatively stained gap junction domains exhibited differences in stain distribution and connexon orientation. To analyze these polymorphic structures, correlation averaging methods were applied to images from negatively stained and frozen-hydrated specimens. For the negatively stained specimens, separate averages over two subsets of connexons with differing degrees of stain accumulation in the axial channel were obtained. Two populations of connexons with opposite skew orientations were distinguishable within a single junctional domain of a frozen-hydrated specimen. Correlation maps calculated using the left- and right-skewed references showed that the selected connexons tend to locally cluster. Using correlation methods to analyze packing disorder in a typical connexon lattice, we estimated the root-mean-square variation in the nearest neighbor pair separation to be approximately 11% of the lattice constant. Displacements of the connexons relative to each other increased with increasing pair separation in the lattice, rather like a liquid, although long-range orientation order was conserved as in a crystal. These results support the hypothesis that the hexagonal ordering of the connexons results from short-range repulsive forces.