Subtomogram Averaging Research Articles

M-free: scoring the reference bias in sub-tomogram averaging and template matching.

Cryo-electron tomography provides a snapshot of the cellular proteome. With template matching, the spatial positions of various macromolecular complexes within their native cellular context can be detected. However, the growing awareness of the reference bias introduced by the cross-correlation based approaches, and more importantly the lack of a reliable confidence measurement in the selection of these macromolecular complexes, has restricted the use of these applications. Here we propose a heuristic, in which the reference bias is measured in real space in an analogous way to the R-free value in X-ray crystallography. We measure the reference bias within the mask used to outline the area of the template, and do not modify the template itself. The heuristic works by splitting the mask into a working and a testing area in a volume ratio of 9:1. While the working area is used during the calculation of the cross-correlation function, the information from both areas is explored to calculate the M-free score. We show using artificial data, that the M-free score gives a reliable measure for the reference bias. The heuristic can be applied in template matching and in sub-tomogram averaging. We further test the applicability of the heuristic in tomograms of purified macromolecules, and tomograms of whole Mycoplasma cells.

John Batten

<h3>Abstract</h3> Cryo-electron tomography (cryo-ET) is widely used to explore the 3D density of biomacromolecules. However, the heavy noise and missing wedge effect prevent directly visualizing and analyzing the 3D reconstructions....

Iterative reconstruction of cryo-electron tomograms using nonuniform fast Fourier transforms

Algorithms for three-dimensional (3D) reconstruction of objects based on their projections are essential in various biological and medical imaging modalities. In cryo-electron tomography (CET) a major challenge for reconstruction is the limited range of projection angles, which manifests itself as a "missing wedge" of data in Fourier space making the reconstruction problem ill-posed. Here, we apply an iterative reconstruction method that makes use of nonuniform fast Fourier transform (NUFFT) to the reconstruction of cryo-electron tomograms. According to several measures the reconstructions are superior to those obtained using conventional methods, most notably weighted backprojection. Most importantly, we show that it is possible to fill in partially the unsampled region in Fourier space with meaningful information without making assumptions about the data or applying prior knowledge. As a consequence, particles of known structure can be localized with higher confidence in cryotomograms and subtomogram averaging yields higher resolution densities.

Quantifying the local resolution of cryo-EM density maps

Using local sinusoidal features in a standard statistical testing framework, we propose a definition of local resolution for 3D density maps. The resulting algorithm has no free parameters and may be extended to other imaging modalities. Evaluating the local resolution of single particle reconstructions and sub-tomogram averages shows variable resolution across a 4 to 40Å range.

Determination of protein structure at 8.5 Å resolution using cryo-electron tomography and sub-tomogram averaging

Cryo-electron tomography combined with image processing by sub-tomogram averaging is unique in its power to resolve the structures of proteins and macromolecular complexes in situ. Limitations of the method, including the low signal to noise ratio within individual images from cryo-tomographic datasets and difficulties in determining the defocus at which the data was collected, mean that to date the very best structures obtained by sub-tomogram averaging are limited to a resolution of approximately 15 Å. Here, by optimizing data collection and defocus determination steps, we have determined the structure of assembled Mason-Pfizer monkey virus Gag protein using sub-tomogram averaging to a resolution of 8.5 Å. At this resolution alpha-helices can be directly and clearly visualized. These data demonstrate for the first time that high-resolution structural information can be obtained from cryo-electron tomograms using sub-tomogram averaging. Sub-tomogram averaging has the potential to allow detailed studies of unsolved and biologically relevant structures under biologically relevant conditions.

Up Close and Personal with an Internal-Membrane Virus

Up Close and Personal with an Internal-Membrane Virus

Optimal and fast rotational alignment of volumes with missing data in Fourier space

Electron tomography of intact cells has the potential to reveal the entire cellular content at a resolution corresponding to individual macromolecular complexes. Characterization of macromolecular complexes in tomograms is nevertheless an extremely challenging task due to the high level of noise, and due to the limited tilt angle that results in missing data in Fourier space. By identifying particles of the same type and averaging their 3D volumes, it is possible to obtain a structure at a more useful resolution for biological interpretation. Currently, classification and averaging of sub-tomograms is limited by the speed of computational methods that optimize alignment between two sub-tomographic volumes. The alignment optimization is hampered by the fact that the missing data in Fourier space has to be taken into account during the rotational search. A similar problem appears in single particle electron microscopy where the random conical tilt procedure may require averaging of volumes with a missing cone in Fourier space. We present a fast implementation of a method guaranteed to find an optimal rotational alignment that maximizes the constrained cross-correlation function (cCCF) computed over the actual overlap of data in Fourier space.

Novel structural labeling method using cryo-electron tomography and biotin–streptavidin system

There are a number of large macromolecular complexes that play important roles in the cell, and identifying the positions of their components is a key step to understanding their structure and function. Several structural labeling methods have been applied to electron microscopy in order to locate a specific component within a macromolecular complex, but each method is associated with problems in specificity, occupancy, signal intensity or precision. Here, we report a novel method for identifying the 3D locations of proteins using biotin-streptavidin labeling and cryo-electron tomography. We labeled a biotinylation-tagged intermediate chain of an axonemal dynein by streptavidin within the Chlamydomonas axoneme and visualized the 3D positions of the labels using subtomogram averaging. Increase of the density attributed to the bound streptavidin was validated by Student's t-test. In conclusion, the combination of the biotin-streptavidin system and cryo-electron tomography is a powerful method to investigate the structure of large macromolecular complexes.

Automated target segmentation and real space fast alignment methods for high-throughput classification and averaging of crowded cryo-electron subtomograms

Motivation: Cryo-electron tomography allows the imaging of macromolecular complexes in near living conditions. To enhance the nominal resolution of a structure it is necessary to align and average individual subtomograms each containing identical complexes. However, if the sample of complexes is heterogeneous, it is necessary to first classify subtomograms into groups of identical complexes. This task becomes challenging when tomograms contain mixtures of unknown complexes extracted from a crowded environment. Two main challenges must be overcomed: First, classification of subtomograms must be performed without knowledge of template structures. However, most alignment methods are too slow to perform reference-free classification of a large number of (e.g. tens of thousands) of subtomograms. Second, subtomograms extracted from crowded cellular environments, contain often fragments of other structures besides the target complex. However, alignment methods generally assume that each subtomogram only contains one complex. Automatic methods are needed to identify the target complexes in a subtomogram even when its shape is unknown.Results: In this article, we propose an automatic and systematic method for the isolation and masking of target complexes in subtomograms extracted from crowded environments. Moreover, we also propose a fast alignment method using fast rotational matching in real space. Our experiments show that, compared with our previously proposed fast alignment method in reciprocal space, our new method significantly improves the alignment accuracy for highly distorted and especially crowded subtomograms. Such improvements are important for achieving successful and unbiased high-throughput reference-free structural classification of complexes inside whole-cell tomograms.Contact: alber@usc.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Orthobunyavirus Ultrastructure and the Curious Tripodal Glycoprotein Spike

The genus Orthobunyavirus within the family Bunyaviridae constitutes an expanding group of emerging viruses, which threaten human and animal health. Despite the medical importance, little is known about orthobunyavirus structure, a prerequisite for understanding virus assembly and entry. Here, using electron cryo-tomography, we report the ultrastructure of Bunyamwera virus, the prototypic member of this genus. Whilst Bunyamwera virions are pleomorphic in shape, they display a locally ordered lattice of glycoprotein spikes. Each spike protrudes 18 nm from the viral membrane and becomes disordered upon introduction to an acidic environment. Using sub-tomogram averaging, we derived a three-dimensional model of the trimeric pre-fusion glycoprotein spike to 3-nm resolution. The glycoprotein spike consists mainly of the putative class-II fusion glycoprotein and exhibits a unique tripod-like arrangement. Protein–protein contacts between neighbouring spikes occur at membrane-proximal regions and intra-spike contacts at membrane-distal regions. This trimeric assembly deviates from previously observed fusion glycoprotein arrangements, suggesting a greater than anticipated repertoire of viral fusion glycoprotein oligomerization. Our study provides evidence of a pH-dependent conformational change that occurs during orthobunyaviral entry into host cells and a blueprint for the structure of this group of emerging pathogens.

Fast and accurate reference-free alignment of subtomograms

In cryoelectron tomography alignment and averaging of subtomograms, each dnepicting the same macromolecule, improves the resolution compared to the individual subtomogram. Major challenges of subtomogram alignment are noise enhancement due to overfitting, the bias of an initial reference in the iterative alignment process, and the computational cost of processing increasingly large amounts of data. Here, we propose an efficient and accurate alignment algorithm via a generalized convolution theorem, which allows computation of a constrained correlation function using spherical harmonics. This formulation increases computational speed of rotational matching dramatically compared to rotation search in Cartesian space without sacrificing accuracy in contrast to other spherical harmonic based approaches. Using this sampling method, a reference-free alignment procedure is proposed to tackle reference bias and overfitting, which also includes contrast transfer function correction by Wiener filtering. Application of the method to simulated data allowed us to obtain resolutions near the ground truth. For two experimental datasets, ribosomes from yeast lysate and purified 20S proteasomes, we achieved reconstructions of approximately 20Å and 16Å, respectively. The software is ready-to-use and made public to the community.

Structural biology in situ—the potential of subtomogram averaging

Cryo-electron tomography provides low-resolution 3D views of cells, organelles, or viruses. Macromolecular complexes present in multiple copies can be subsequently identified within the 3D reconstruction (the tomogram), computationally extracted, and averaged to obtain higher resolution 3D structures, as well as a map of their spatial distribution. This method, called subtomogram averaging or subvolume averaging, allows structures of macromolecular complexes to be resolved in situ. Recent applications have provided in situ structural data at resolutions of 2-4 nm on samples including polysomes, nuclear pores, vesicle coats, and viral surface proteins. Here I describe the method and discuss limitations, advances and recent applications. I speculate how the method will solve more structures at higher resolution, allowing in situ structural biology.

Computational Separation of Conformational Heterogeneity using Cryo-Electron Tomography and 3D Sub-Volume Averaging

We have previously used cryo-electron tomography combined with missing-wedge corrected, sub-volume averaging and classification to obtain 3D structures of macromolecular assemblies in cases where a single dominant species was present, and applied these methods to the analysis of a variety of trimeric HIV-1 and SIV envelope glycoproteins (Env). Here, we extend these studies by presenting a method for determining the distribution of conformational states found in a specimen and validate these procedures by successfully separating and reconstructing distinct 3D structures for unliganded and antibody-liganded as well as open and closed conformations of Env present simultaneously in mixtures. We show that identifying and removing spikes with the lowest signal-to-noise ratios improves the overall accuracy of alignment between individual Env sub-volumes, and that alignment accuracy, in turn, determines the success of image classification in assessing conformational heterogeneity in heterogeneous mixtures. This development extends the sub-tomogram averaging capabilities to heterogeneous samples. Furthermore, it turns cryo-electron tomography into a powerful analytical tool that can directly determine the relative amount of the different conformations found in the specimen.View Large Image | View Hi-Res Image | Download PowerPoint Slide

Conserved structural motifs in the central pair complex of eukaryotic flagella.

Cilia and flagella are conserved hair-like appendages of eukaryotic cells that function as sensing and motility generating organelles. Motility is driven by thousands of axonemal dyneins that require precise regulation. One essential motility regulator is the central pair complex (CPC) and many CPC defects cause paralysis of cilia/flagella. Several human diseases, such as immotile cilia syndrome, show CPC abnormalities, but little is known about the detailed three-dimensional (3D) structure and function of the CPC. The CPC is located in the center of typical [9+2] cilia/flagella and is composed of two singlet microtubules (MTs), each with a set of associated projections that extend toward the surrounding nine doublet MTs. Using cryo-electron tomography coupled with subtomogram averaging, we visualized and compared the 3D structures of the CPC in both the green alga Chlamydomonas and the sea urchin Strongylocentrotus at the highest resolution published to date. Despite the evolutionary distance between these species, their CPCs exhibit remarkable structural conservation. We identified several new projections, including those that form the elusive sheath, and show that the bridge has a more complex architecture than previously thought. Organism-specific differences include the presence of MT inner proteins in Chlamydomonas, but not Strongylocentrotus, and different overall outlines of the highly connected projection network, which forms a round-shaped cylinder in algae, but is more oval in sea urchin. These differences could be adaptations to the mechanical requirements of the rotating CPC in Chlamydomonas, compared to the Strongylocentrotus CPC which has a fixed orientation.

Single particle and molecular assembly analysis of polyribosomes by single- and double-tilt cryo electron tomography

Cryo electron tomography (cryo-ET) can provide cellular and molecular structural information on various biological samples. However, the detailed interpretation of tomograms reconstructed from single-tilt data tends to suffer from low signal-to-noise ratio and artefacts caused by some systematically missing angular views. While these can be overcome by sub-tomogram averaging, they remain limiting for the analysis of unique structures. Double-tilt ET can improve the tomogram quality by acquiring a second tilt series after an in-plane rotation, but its usage is not widespread yet because it is considered technically demanding and it is rarely used under cryo conditions. Here we show that double-tilt cryo-ET improves the quality of 3D reconstructions so significantly that even single particle analysis can be envisaged despite of the intrinsically low image contrast obtained from frozen-hydrated specimens. This is illustrated by the analysis of eukaryotic polyribosomes in which individual ribosomes were reconstructed using single-tilt, partial and full double-tilt geometries. The improved tomograms favour the faster convergence of iterative sub-tomogram averaging and allow a better 3D classification using multivariate statistical analysis. Our study of single particles and molecular assemblies within polysomes illustrates that the dual-axis approach is particularly useful for cryo applications of ET, both for unique objects and for structures that can be classified and averaged.

Unraveling the structure of membrane proteins in situ by transfer function corrected cryo-electron tomography

Cryo-electron tomography in combination with subtomogram averaging allows to investigate the structure of protein assemblies in their natural environment in a close to live state. To make full use of the structural information contained in tomograms it is necessary to analyze the contrast transfer function (CTF) of projections and to restore the phases of higher spatial frequencies. CTF correction is however hampered by the difficulty of determining the actual defocus values from tilt series data, which is due to the low signal-to-noise ratio of electron micrographs. In this study, an extended acquisition scheme is introduced that enables an independent CTF determination. Two high-dose images are recorded along the tilt axis on both sides of each projection, which allow an accurate determination of the defocus values of these images. These values are used to calculate the CTF for each image of the tilt series. We applied this scheme to the mycobacterial outer membrane protein MspA reconstituted in lipid vesicles and tested several variants of CTF estimation in combination with subtomogram averaging and correction of the modulation transfer function (MTF). The 3D electron density map of MspA was compared with a structure previously determined by X-ray crystallography. We were able to demonstrate that structural information up to a resolution of 16.8Å can be recovered using our CTF correction approach, whereas the uncorrected 3D map had a resolution of only 26.2Å.

Polarity and asymmetry in the arrangement of dynein and related structures in the Chlamydomonas axoneme

Understanding the molecular architecture of the flagellum is crucial to elucidate the bending mechanism produced by this complex organelle. The current known structure of the flagellum has not yet been fully correlated with the complex composition and localization of flagellar components. Using cryoelectron tomography and subtomogram averaging while distinguishing each one of the nine outer doublet microtubules, we systematically collected and reconstructed the three-dimensional structures in different regions of the Chlamydomonas flagellum. We visualized the radial and longitudinal differences in the flagellum. One doublet showed a distinct structure, whereas the other eight were similar but not identical to each other. In the proximal region, some dyneins were missing or replaced by minor dyneins, and outer-inner arm dynein links were variable among different microtubule doublets. These findings shed light on the intricate organization of Chlamydomonas flagella, provide clues to the mechanism that produces asymmetric flagellar beating, and pose a new challenge for the functional study of the flagella.

Visualization of the Two-Step Fusion Process of the Retrovirus Avian Sarcoma/Leukosis Virus by Cryo-Electron Tomography

Retrovirus infection starts with the binding of envelope glycoproteins to host cell receptors. Subsequently, conformational changes in the glycoproteins trigger fusion of the viral and cellular membranes. Some retroviruses, such as avian sarcoma/leukosis virus (ASLV), employ a two-step mechanism in which receptor binding precedes low-pH activation and fusion. We used cryo-electron tomography to study virion/receptor/liposome complexes that simulate the interactions of ASLV virions with cells. Binding the soluble receptor at neutral pH resulted in virions capable of binding liposomes tightly enough to alter their curvature. At virion-liposome interfaces, the glycoproteins are ∼3-fold more concentrated than elsewhere in the viral envelope, indicating specific recruitment to these sites. Subtomogram averaging showed that the oblate globular domain in the prehairpin intermediate (presumably the receptor-binding domain) is connected to both the target and the viral membrane by 2.5-nm-long stalks and is partially disordered, compared with its native conformation. Upon lowering the pH, fusion took place. Fusion is a stochastic process that, once initiated, must be rapid, as only final (postfusion) products were observed. These fusion products showed glycoprotein spikes on their surface, with their interiors occupied by patches of dense material but without capsids, implying their disassembly. In addition, some of the products presented a density layer underlying and resolved from the viral membrane, which may represent detachment of the matrix protein to facilitate the fusion process.

Structure of the yeast F 1 F o -ATP synthase dimer and its role in shaping the mitochondrial cristae

We used electron cryotomography of mitochondrial membranes from wild-type and mutant Saccharomyces cerevisiae to investigate the structure and organization of ATP synthase dimers in situ. Subtomogram averaging of the dimers to 3.7 nm resolution revealed a V-shaped structure of twofold symmetry, with an angle of 86° between monomers. The central and peripheral stalks are well resolved. The monomers interact within the membrane at the base of the peripheral stalks. In wild-type mitochondria ATP synthase dimers are found in rows along the highly curved cristae ridges, and appear to be crucial for membrane morphology. Strains deficient in the dimer-specific subunits e and g or the first transmembrane helix of subunit 4 lack both dimers and lamellar cristae. Instead, cristae are either absent or balloon-shaped, with ATP synthase monomers distributed randomly in the membrane. Computer simulations indicate that isolated dimers induce a plastic deformation in the lipid bilayer, which is partially relieved by their side-by-side association. We propose that the assembly of ATP synthase dimer rows is driven by the reduction in the membrane elastic energy, rather than by direct protein contacts, and that the dimer rows enable the formation of highly curved ridges in mitochondrial cristae.

GRecon: A Method for the Lipid Reconstitution of Membrane Proteins

Membrane proteins account for about 20–30 % of the protein-encoding genes in the genomes of all living organisms.2 Their fundamental importance in human health and disease is underlined by the fact that many drugs in current use act on them.3 Functional and structural studies of membrane proteins are therefore of increasing importance but more difficult and demanding than studies with soluble proteins, because membrane proteins function in the lipid bilayer of cell membranes. For in vitro studies, the proteins are first extracted from the membrane by detergent solubilization and then purified in detergent solution. Sensitive membrane proteins are often unstable in detergent solution, but quite stable once they are reconstituted into a lipid bilayer. Moreover, many membrane proteins require a lipid environment for activity. The reconstitution process, in which detergent is replaced by lipid, must be carefully controlled, as otherwise the proteins tend to denature and aggregate.4