High-resolution cryo-electron microscopy on macromolecular complexes and cell organelles

Cryo-electron microscopy techniques and computational 3-D reconstruction of macromolecular assemblies are tightly linked tools in modern structural biology. This symbiosis has produced vast amounts of detailed information on the structure and function of biological macromolecules. Typically, one of two fundamentally different strategies is used depending on the specimens and their environment. A: 3-D reconstruction based on repetitive and structurally identical unit cells that allow for averaging, and B: tomographic 3-D reconstructions where tilt-series between approximately ±60 and ±70° at small angular increments are collected from highly complex and flexible structures that are beyond averaging procedures, at least during the first round of 3-D reconstruction. Strategies of group A are averaging-based procedures and collect large number of 2-D projections at different angles that are computationally aligned, averaged together, and back-projected in 3-D space to reach a most complete 3-D dataset with high resolution, today often down to atomic detail. Evidently, success relies on structurally repetitive particles and an aligning procedure that unambiguously determines the angular relationship of all 2-D projections with respect to each other. The alignment procedure of small particles may rely on their packing into a regular array such as a 2-D crystal, an icosahedral (viral) particle, or a helical assembly. Critically important for cryo-methods, each particle will only be exposed once to the electron beam, making these procedures optimal for highest-resolution studies where beam-induced damage is a significant concern. In contrast, tomographic 3-D reconstruction procedures (group B) do not rely on averaging, but collect an entire dataset from the very same structure of interest. Data acquisition requires collecting a large series of tilted projections at angular increments of 1–2° or less and a tilt range of ±60° or more. Accordingly, tomographic data collection exposes its specimens to a large electron dose, which is particularly problematic for frozen-hydrated samples. Currently, cryo-electron tomography is a rapidly emerging technology, on one end driven by the newest developments of hardware such as super-stabile microscopy stages as well as the latest generation of direct electron detectors and cameras. On the other end, success also strongly depends on new software developments on all kinds of fronts such as tilt-series alignment and back-projection procedures that are all adapted to the very low-dose and therefore very noisy primary data. Here, we will review the status quo of cryo-electron microscopy and discuss the future of cellular cryo-electron tomography from data collection to data analysis, CTF-correction of tilt-series, post-tomographic sub-volume averaging, and 3-D particle classification. We will also discuss the pros and cons of plunge freezing of cellular specimens to vitrified sectioning procedures and their suitability for post-tomographic volume averaging despite multiple artifacts that may distort specimens to some degree.     

Corrim-based alignment for improved speed in single-particle image processing

The technique of single-particle electron cryomicroscopy is currently making possible the 3D structure determination of large macromolecular complexes at constantly increasing levels of resolution. Work at resolution now attainable requires many thousands of individual images to be processed computationally. The most time-consuming step of the image-processing procedure is usually the iterative alignment of individual particle images against a set of reference images derived from a preliminary 3-D structure. We have developed an improved multireference alignment procedure based on interpolated cross-correlation images (corrims) that results in an approximately 8-fold acceleration of the iterative alignment steps. These corrims can be used to restrict the number of image-alignment calculations by narrowing down the set of reference images. Another improvement in alignment speed has been achieved by optimising the software and its implementation on many parallel processors. This new corrim-based refinement has been found to work well with two different alignment algorithms, the commonly used "fast alignment by separate translational/rotational searches" and "exhaustive alignment by polar coordinates."     

Aligning multiple protein structures by deterministic annealing

Protein structure alignment plays a key role in protein structure prediction and fold family classification. An efficient method for multiple protein structure alignment in a mathematical manner is presented, based on deterministic annealing technique. The alignment problem is mapped onto a nonlinear continuous optimization problem (NCOP) with common consensus chain, matching assignment matrices and atomic coordinates as variables. At each step in the annealing procedure, the NCOP is decomposed into as many sub-problems as the number of protein chains, each of which is actually an independent pairwise structure alignment between a protein chain and the consensus chain and hence can be efficiently solved by the parallel computation technique. The proposed method is robust with respect to choice of iteration parameters for a wide range of proteins, and performs well in both multiple and pairwise structure alignment cases, compared with existing alignment methods.     

Electron tomography of fast frozen, stretched rigor fibers reveals elastic distortions in the myosin crossbridges

As a first step toward freeze-trapping and 3-D modeling of the very rapid load-induced structural responses of active myosin heads, we explored the conformational range of longer lasting force-dependent changes in rigor crossbridges of insect flight muscle (IFM). Rigor IFM fibers were slam-frozen after ramp stretch (1000 ms) of 1-2% and freeze-substituted. Tomograms were calculated from tilt series of 30 nm longitudinal sections of Araldite-embedded fibers. Modified procedures of alignment and correspondence analysis grouped self-similar crossbridge forms into 16 class averages with 4.5 nm resolution, revealing actin protomers and myosin S2 segments of some crossbridges for the first time in muscle thin sections. Acto-S1 atomic models manually fitted to crossbridge density required a range of lever arm adjustments to match variably distorted rigor crossbridges. Some lever arms were unchanged compared with low tension rigor, while others were bent and displaced M-ward by up to 4.5 nm. The average displacement was 1.6 +/- 1.0 nm. "Map back" images that replaced each unaveraged 39 nm crossbridge motif by its class average showed an ordered mix of distorted and unaltered crossbridges distributed along the 116 nm repeat that reflects differences in rigor myosin head loading even before stretch.     

Application of electron tomography to fungal ultrastructure studies

Access to structural information at the nanoscale enables fundamental insights into many complex biological systems. The development of the transmission electron microscope (TEM) has vastly increased our understanding of multiple biological systems. However, when attempting to visualize and understand the organizational and functional complexities that are typical of cells and tissues, the standard 2-D analyses that TEM affords often fall short. In recent years, high-resolution electron tomography methods, coupled with advances in specimen preparation and instrumentation and computational speed, have resulted in a revolution in the biological sciences. Electron tomography is analogous to medical computerized axial tomography (CAT-scan imaging) except at a far finer scale. It utilizes the TEM to assemble multiple projections of an object which are then combined for 3-D analyses. For biological specimens, tomography enables the highest 3-D resolution (5 nm spatial resolution) of internal structures in relatively thick slices of material (0.2-0.4 microm) without requiring the collection and alignment of large numbers of thin serial sections. Thus accurate and revealing 3-D reconstructions of complex cytoplasmic entities and architecture can be obtained. Electron tomography is now being applied to a variety of biological questions with great success. This review gives a brief introduction into cryopreservation and electron tomography relative to aspects of cytoplasmic organization in the hyphal tip of Aspergillus nidulans.     

The HSSP database of protein structure-sequence alignments.

HSSP is a derived database merging structural three dimensional (3-D) and sequence one dimensional(1-D) information. For each protein of known 3-D structure from the Protein Data Bank (PDB), the database has a multiple sequence alignment of all available homologues and a sequence profile characteristic of the family. The list of homologues is the result of a database search in Swissprot using a position-weighted dynamic programming method for sequence profile alignment (MaxHom). The database is updated frequently. The listed homologues are very likely to have the same 3-D structure as the PDB protein to which they have been aligned. As a result, the database is not only a database of aligned sequence families, but also a database of implied secondary and tertiary structures covering 27% of all Swissprot-stored sequences.     

The HSSP database of protein structure-sequence alignments.

HSSP is a derived database merging structural (3-D) and sequence (1-D) information. For each protein of known 3-D structure from the Protein Data Bank (PDB), the database has a multiple sequence alignment of all available homologues and a sequence profile characteristic of the family. The list of homologues is the result of a database search in SwissProt using a position-weighted dynamic programming method for sequence profile alignment (MaxHom). The database is updated frequently. The listed homologues are very likely to have the same 3-D structure as the PDB protein to which they have been aligned. As a result, the database is not only a database of aligned sequence families, but also a database of implied secondary and tertiary structures covering 29% of all SwissProt-stored sequences.     

The HSSP database of protein structure-sequence alignments and family profiles.

HSSP (http: //www.sander.embl-ebi.ac.uk/hssp/) is a derived database merging structure (3-D) and sequence (1-D) information. For each protein of known 3D structure from the Protein Data Bank (PDB), we provide a multiple sequence alignment of putative homologues and a sequence profile characteristic of the protein family, centered on the known structure. The list of homologues is the result of an iterative database search in SWISS-PROT using a position-weighted dynamic programming method for sequence profile alignment (MaxHom). The database is updated frequently. The listed putative homologues are very likely to have the same 3D structure as the PDB protein to which they have been aligned. As a result, the database not only provides aligned sequence families, but also implies secondary and tertiary structures covering 33% of all sequences in SWISS-PROT.     

Auto-accumulation method using simulated annealing enables fully automatic particle pickup completely free from a matching template or learning data

Single-particle analysis is a 3-D structure determining method using electron microscopy (EM). In this method, a large number of projections is required to create 3-D reconstruction. In order to enable completely automatic pickup without a matching template or a training data set, we established a brand-new method in which the frames to pickup particles are randomly shifted and rotated over the electron micrograph and, using the total average image of the framed images as an index, each frame reaches a particle. In this process, shifts are selected to increase the contrast of the average. By iterated shifts and further selection of the shifts, the frames are induced to shift so as to surround particles. In this algorithm, hundreds of frames are initially distributed randomly over the electron micrograph in which multi-particle images are dispersed. Starting with these frames, one of them is selected and shifted randomly, and acceptance or non-acceptance of its new position is judged using the simulated annealing (SA) method in which the contrast score of the total average image is adopted as an index. After iteration of this process, the position of each frame converges so as to surround a particle and the framed images are picked up. This method is the first unsupervised fully automatic particle picking method which is applicable to EM of various kinds of proteins, especially to low-contrasted cryo-EM protein images.     

Segmentation by classification: A novel and reliable approach for semi-automatic selection of HIV/SIV envelope spikes

We describe a semiautomated approach to segment Env spikes from the membrane envelope of Simian Immunodeficiency Virus visualized by cryoelectron tomography of frozen-hydrated specimens. Multivariate data analysis is applied to a large set of overlapping subvolumes extracted semiautomatically from the viral envelope and does not utilize a template of the target structure. The major manual step used in the method involves determination of six points that define an ellipsoid approximating the virion shape. The approach is robust to departures of the actual virion from this starting ellipsoid. A point cage of sufficient density is generated to ensure that any spike-like protein is identified multiple times. Subsequently translational alignment of class averages to a cylindrical reference on a curved surface separates subvolumes with spikes from those without. Spike containing subvolumes identified multiple times are removed by proximity analysis. Slightly different procedures segment spikes in the equatorial and the polar regions. Once all spikes are segmented, further alignment of class averages using separately the polar and spin angles produces recognizable spike images. Our approach localized 96% of the equatorial spikes and 85% of all spikes identified manually; it identifies a significant number of additional spikes missed by manual selection. Two types of spike shapes were segmented, one with near 3-fold symmetry resembling the conventional spike, the other had a T-shape resembling the spike structure obtained when antibodies such as PG9 bind to HIV Env. The approach should be applicable to segmentation of any protein spikes extending from a cellular or virion envelope.     

The Rhodopseudomonas viridis photosynthetic membrane: arrangement in situ

The organization of photosynthetic membranes in the cytoplasm of the photosynthetic bacterium Rh. viridis has been examined by several techniques for electron microscopy. Thin sections of membrane stacks show that the regular lattice of membrane subunits reported in other studies can be observed in thin section. Tilting of sections in the electron microscope shows that the regular lattices of several membranes overlap in a way that suggests they are in register with each other. This observation can be confirmed by freeze-fracture images in which a regular arrangement of membrane lattices can be observed, each perfectly aligned.Analysis of the spacings of membrane pairs shows that the photosynthetic membranes of Rh. viridis are very closely apposed. The mean diameter of two membranes is 160A, and the average space between two such membranes is only 42A. When a recently developed atomic level model of Rh. viridis reaction center is superimposed against these spacings, each reaction center extends from the surface of its respective membrane far enough to make contact with an apposing membrane. The limited free space between membranes and regular alignment of lattices has a number of implications for how this membrane is organized to carry out the process of energy transfer.     

Collective Motion of Spherical Bacteria

A large variety of motile bacterial species exhibit collective motions while inhabiting liquids or colonizing surfaces. These collective motions are often characterized by coherent dynamic clusters, where hundreds of cells move in correlated whirls and jets. Previously, all species that were known to form such motion had a rod-shaped structure, which enhances the order through steric and hydrodynamic interactions. Here we show that the spherical motile bacteria Serratia marcescens exhibit robust collective dynamics and correlated coherent motion while grown in suspensions. As cells migrate to the upper surface of a drop, they form a monolayer, and move collectively in whirls and jets. At all concentrations, the distribution of the bacterial speed was approximately Rayleigh with an average that depends on concentration in a non-monotonic way. Other dynamical parameters such as vorticity and correlation functions are also analyzed and compared to rod-shaped bacteria from the same strain. Our results demonstrate that self-propelled spherical objects do form complex ordered collective motion. This opens a door for a new perspective on the role of cell aspect ratio and alignment of cells with regards to collective motion in nature.     

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

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

YUP.SCX: coaxing atomic models into medium resolution electron density maps

The structures of large macromolecular complexes in different functional states can be determined by cryo-electron microscopy, which yields electron density maps of low to intermediate resolutions. The maps can be combined with high-resolution atomic structures of components of the complex, to produce a model for the complex that is more accurate than the formal resolution of the map. To this end, methods have been developed to dock atomic models into density maps rigidly or flexibly, and to refine a docked model so as to optimize the fit of the atomic model into the map. We have developed a new refinement method called YUP.SCX. The electron density map is converted into a component of the potential energy function to which terms for stereochemical restraints and volume exclusion are added. The potential energy function is then minimized (using simulated annealing) to yield a stereochemically-restrained atomic structure that fits into the electron density map optimally. We used this procedure to construct an atomic model of the 70S ribosome in the pre-accommodation state. Although some atoms are displaced by as much as 33 Å, they divide themselves into nearly rigid fragments along natural boundaries with smooth transitions between the fragments.     

Top-down approach in protein RDC data analysis: <Emphasis Type="Italic">de novo</Emphasis> estimation of the alignment tensor

In solution NMR spectroscopy the residual dipolar coupling (RDC) is invaluable in improving both the precision and accuracy of NMR structures during their structural refinement. The RDC also provides a potential to determine protein structure de novo. These procedures are only effective when an accurate estimate of the alignment tensor has already been made. Here we present a top–down approach, starting from the secondary structure elements and finishing at the residue level, for RDC data analysis in order to obtain a better estimate of the alignment tensor. Using only the RDCs from N–H bonds of residues in α-helices and CA–CO bonds in β-strands, we are able to determine the offset and the approximate amplitude of the RDC modulation-curve for each secondary structure element, which are subsequently used as targets for global minimization. The alignment order parameters and the orientation of the major principal axis of individual helix or strand, with respect to the alignment frame, can be determined in each of the eight quadrants of a sphere. The following minimization against RDC of all residues within the helix or strand segment can be carried out with fixed alignment order parameters to improve the accuracy of the orientation. For a helical protein Bax, the three components A _xx , A _yy and A _zz , of the alignment order can be determined with this method in average to within 2.3% deviation from the values calculated with the available atomic coordinates. Similarly for β-sheet protein Ubiquitin they agree in average to within 8.5%. The larger discrepancy in β-strand parameters comes from both the diversity of the β-sheet structure and the lower precision of CA–CO RDCs. This top-down approach is a robust method for alignment tensor estimation and also holds a promise for providing a protein topological fold using limited sets of RDCs.     

Fiducial-less alignment of cryo-sections

Cryo-electron tomography of vitreous sections is currently the most promising technique for visualizing arbitrary regions of eukaryotic cells or tissue at molecular resolution. Despite significant progress in the sample preparation techniques over the past few years, the three dimensional reconstruction using electron tomography is not as simple as in plunge frozen samples for various reasons, but mainly due to the effects of irradiation on the sections and the resulting poor alignment. Here, we present a new algorithm, which can provide a useful three-dimensional marker model after investigation of hundreds to thousands of observations calculated using local cross-correlation throughout the tilt series. The observations are chosen according to their coherence to a particular model and assigned to virtual markers. Through this type of measurement a merit figure can be calculated, precisely estimating the quality of the reconstruction. The merit figures of this alignment method are comparable to those obtained with plunge frozen samples using fiducial gold markers. An additional advantage of the algorithm is the implicit detection of areas in the sections that behave as rigid bodies and can thus be properly reconstructed.     

Reprint of "Fiducial-less alignment of cryo-sections" [J. Struct. Biol. 159 (2007) 413-423

Cryo-electron tomography of vitreous sections is currently the most promising technique for visualizing arbitrary regions of eukaryotic cells or tissue at molecular resolution. Despite significant progress in the sample preparation techniques over the past few years, the three dimensional reconstruction using electron tomography is not as simple as in plunge frozen samples for various reasons, but mainly due to the effects of irradiation on the sections and the resulting poor alignment. Here, we present a new algorithm, which can provide a useful three-dimensional marker model after investigation of hundreds to thousands of observations calculated using local cross-correlation throughout the tilt series. The observations are chosen according to their coherence to a particular model and assigned to virtual markers. Through this type of measurement a merit figure can be calculated, precisely estimating the quality of the reconstruction. The merit figures of this alignment method are comparable to those obtained with plunge frozen samples using fiducial gold markers. An additional advantage of the algorithm is the implicit detection of areas in the sections that behave as rigid bodies and can thus be properly reconstructed.     

Weak alignment offers new NMR opportunities to study protein structure and dynamics

Protein solution nuclear magnetic resonance (NMR) can be conducted in a slightly anisotropic environment, where the orientational distribution of the proteins is no longer random. In such an environment, the large one-bond internuclear dipolar interactions no longer average to zero and report on the average orientation of the corresponding vectors relative to the magnetic field. The desired very weak ordering, on the order of 10(-3), can be induced conveniently by the use of aqueous nematic liquid crystalline suspensions or by anisotropically compressed hydrogels. The resulting residual dipolar interactions are scaled down by three orders of magnitude relative to their static values, but nevertheless can be measured at high accuracy. They are very precise reporters on the average orientation of bonds relative to the molecular alignment frame, and they can be used in a variety of ways to enrich our understanding of protein structure and function. Applications to date have focused primarily on validation of structures, determined by NMR, X-ray crystallography, or homology modeling, and on refinement of structures determined by conventional NMR approaches. Although de novo structure determination on the basis of dipolar couplings suffers from a severe multiple minimum problem, related to the degeneracy of dipolar coupling relative to inversion of the internuclear vector, a number of approaches can address this problem and potentially can accelerate the NMR structure determination process considerably. In favorable cases, where large numbers of dipolar couplings can be measured, inconsistency between measured values can report on internal motions.     

Methods for identifying and averaging variable molecular conformations in tomograms of actively contracting insect flight muscle

During active muscle contraction, tension is generated through many simultaneous, independent interactions between the molecular motor myosin and the actin filaments. The ensemble of myosin motors displays heterogeneous conformations reflecting different mechanochemical steps of the ATPase pathway. We used electron tomography of actively contracting insect flight muscle fast-frozen, freeze substituted, Araldite embedded, thin-sectioned and stained, to obtain 3D snapshots of the multiplicity of actin-attached myosin structures. We describe procedures for alignment of the repeating lattice of sub-volumes (38.7 nm cross-bridge repeats bounded by troponin) and multivariate data analysis to identify self-similar repeats for computing class averages. Improvements in alignment and classification of repeat sub-volumes reveals (for the first time in active muscle images) the helix of actin subunits in the thin filament and the troponin density with sufficient clarity that a quasiatomic model of the thin filament can be built into the class averages independent of the myosin cross-bridges. We show how quasiatomic model building can identify both strong and weak myosin attachments to actin. We evaluate the accuracy of image classification to enumerate the different types of actin–myosin attachments.     

Comparative protein structure modeling by iterative alignment,model building and model assessment

Comparative or homology protein structure modeling is severely limited by errors in the alignment of a modeled sequence with related proteins of known three-dimensional structure. To ameliorate this problem, we have developed an automated method that optimizes both the alignment and the model implied by it. This task is achieved by a genetic algorithm protocol that starts with a set of initial alignments and then iterates through re-alignment, model building and model assessment to optimize a model assessment score. During this iterative process: (i) new alignments are constructed by application of a number of operators, such as alignment mutations and cross-overs; (ii) comparative models corresponding to these alignments are built by satisfaction of spatial restraints, as implemented in our program MODELLER; (iii) the models are assessed by a variety of criteria, partly depending on an atomic statistical potential. When testing the procedure on a very difficult set of 19 modeling targets sharing only 4–27% sequence identity with their template structures, the average final alignment accuracy increased from 37 to 45% relative to the initial alignment (the alignment accuracy was measured as the percentage of positions in the tested alignment that were identical to the reference structure-based alignment). Correspondingly, the average model accuracy increased from 43 to 54% (the model accuracy was measured as the percentage of the Cα atoms of the model that were within 5 Å of the corresponding Cα atoms in the superposed native structure). The present method also compares favorably with two of the most successful previously described methods, PSI-BLAST and SAM. The accuracy of the final models would be increased further if a better method for ranking of the models were available.