The ribosome at improved resolution: new techniques for merging and orientation refinement in 3D cryo-electron microscopy of biological particles

Estimating alignment errors in sets of 2-D images

We describe a robust and accurate method for the estimation of alignment errors for a set of two-dimensional images, in the case where the true pattern is unknown. The intended application of the proposed method is cryo-electron microscopy, where two-dimensional views of individual proteins in random orientations are observed in the electron microscope at low signal-to-noise ratio. By representing images in the basis of Fourier-harmonic coordinates and constructing averages and average intensities, we demonstrate that the variances of translation and rotational errors as well as of the Gaussian noise can be recovered. This machinery therefore allows one to isolate the various categories of errors that impede the quality of results in single particle reconstructions into constituent parts: translational errors, rotational errors, and additive noise.

XMIPP: a new generation of an open-source image processing package for electron microscopy

X-windows based microscopy image processing package (Xmipp) is a specialized suit of image processing programs, primarily aimed at obtaining the 3D reconstruction of biological specimens from large sets of projection images acquired by transmission electron microscopy. This public-domain software package was introduced to the electron microscopy field eight years ago, and since then it has changed drastically. New methodologies for the analysis of single-particle projection images have been added to classification, contrast transfer function correction, angular assignment, 3D reconstruction, reconstruction of crystals, etc. In addition, the package has been extended with functionalities for 2D crystal and electron tomography data. Furthermore, its current implementation in C++, with a highly modular design of well-documented data structures and functions, offers a convenient environment for the development of novel algorithms. In this paper, we present a general overview of a new generation of Xmipp that has been re-engineered to maximize flexibility and modularity, potentially facilitating its integration in future standardization efforts in the field. Moreover, by focusing on those developments that distinguish Xmipp from other packages available, we illustrate its added value to the electron microscopy community.

3D reconstruction of the Mu transposase and the Type 1 transpososome: a structural framework for Mu DNA transposition

Mu DNA transposition proceeds through a series of higher-order nucleoprotein complexes called transpososomes. The structural core of the transpososome is a tetramer of the transposase, Mu A, bound to the two transposon ends. High-resolution structural analysis of the intact transposase and the transpososome has not been successful to date. Here we report the structure of Mu A at 16-angstroms and the Type 1 transpososome at 34-angstroms resolution, by 3D reconstruction of images obtained by scanning transmission electron microscopy (STEM) at cryo-temperatures. Electron spectroscopic imaging (ESI) of the DNA-phosphorus was performed in conjunction with the structural investigation to derive the path of the DNA through the transpososome and to define the DNA-binding surface in the transposase. Our model of the transpososome fits well with the accumulated biochemical literature for this intricate transposition system, and lays a structural foundation for biochemical function, including catalysis in trans and the complex circuit of macromolecular interactions underlying Mu DNA transposition.

Dynamics of EF-G interaction with the ribosome explored by classification of a heterogeneous cryo-EM dataset

A method of supervised classification using two available structure templates was applied to investigate the possible heterogeneity existing in a large cryo-EM dataset of an Escherichia coli 70S ribosome-EF-G complex. Two subpopulations showing the ribosome in distinct conformational states, related by a ratchet-like rotation of the 30S subunit with respect to the 50S subunit, were extracted from the original dataset. The possible presence of additional intermediate states is discussed.

Structure and molecular organization of mammalian fatty acid synthase

De novo synthesis of fatty acids in the cytosol of animal cells is carried out by the multifunctional, homodimeric fatty acid synthase (FAS). Cryo-EM analysis of single FAS particles imaged under conditions that limit conformational variability, combined with gold labeling of the N termini and structural analysis of the FAS monomers, reveals two coiled monomers in an overlapping arrangement. Comparison of dimeric FAS structures related to different steps in the fatty acid synthesis process indicates that only limited local rearrangements are required for catalytic interaction among different functional domains. Monomer coiling probably contributes to FAS efficiency and provides a structural explanation for the reported activity of a FAS monomer dimerized to a catalytically inactive partner. The new FAS structure provides a new paradigm for understanding the architecture of FAS and the related modular polyketide synthases.

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

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

Noise bias in the refinement of structures derived from single particles

One of the main goals in the determination of three-dimensional macromolecular structures from electron microscope images of individual molecules and complexes (single particles) is a sufficiently high spatial resolution, about 4 A, at which the interpretation with an atomic model becomes possible. To reach high resolution, an iterative refinement procedure using an expectation maximization algorithm is often used that leads to a more accurate alignment of the positional and orientational parameters for each particle. We show here the results of refinement algorithms that use a phase residual, a linear correlation coefficient, or a weighted correlation coefficient to align individual particles. The algorithms were applied to computer-generated data sets that contained projections from model structures, as well as noise. The algorithms show different degrees of over-fitting, especially at high resolution where the signal is weak. We demonstrate that the degree of over-fitting is reduced with a weighting scheme that depends on the signal-to-noise ratio in the data. The weighting also improves the accuracy of resolution measurement by the commonly used Fourier shell correlation. The performance of the refinement algorithms is compared to that using a maximum likelihood approach. The weighted correlation coefficient was implemented in the computer program FREALIGN.

Unified 3-D structure and projection orientation refinement using quasi-Newton algorithm

We describe an algorithm for simultaneous refinement of a three-dimensional (3-D) density map and of the orientation parameters of two-dimensional (2-D) projections that are used to reconstruct this map. The application is in electron microscopy, where the 3-D structure of a protein has to be determined from a set of 2-D projections collected at random but initially unknown angles. The design of the algorithm is based on the assumption that initial low resolution approximation of the density map and reasonable guesses for orientation parameters are available. Thus, the algorithm is applicable in final stages of the structure refinement, when the quality of the results is of main concern. We define the objective function to be minimized in real space and solve the resulting nonlinear optimization problem using a Quasi-Newton algorithm. We calculate analytical derivatives with respect to density distribution and the finite difference approximations of derivatives with respect to orientation parameters. We demonstrate that calculation of derivatives is robust with respect to noise in the data. This is due to the fact that noise is annihilated by the back-projection operations. Our algorithm is distinguished from other orientation refinement methods (i) by the simultaneous update of the density map and orientation parameters resulting in a highly efficient computational scheme and (ii) by the high quality of the results produced by a direct minimization of the discrepancy between the 2-D data and the projected views of the reconstructed 3-D structure. We demonstrate the speed and accuracy of our method by using simulated data.

Rearrangement of the 16S precursor subunits is essential for the formation of the active 20S proteasome

Proteasome-dependent proteolysis is essential for a number of key cellular processes and requires a sophisticated biogenesis pathway to function. Here, we have arrested the assembly process in its dynamic progression at the short-lived 16S state. Structural analysis of the 16S proteasome precursor intermediates by electron microscopy, and single particle analysis reveals major conformational changes in the structure of the beta-ring in comparison with one-half of the 20S proteasome. The individual beta-subunits in the 16S precursor complex rotate with respect to their positions in the x-ray crystallographic structure of the fully assembled 20S. This rearrangement results in a movement of the catalytic residue threonine-1 from the protected location in 16S precursor complexes to a more exposed position in the 20S structure. Thereby, our findings provide a molecular explanation for the structural rearrangements necessary for the dimerization of two 16S precursor complexes and the subsequent final maturation to active 20S proteasomes.

Structural Studies of Arthrin: Monoubiquitinated Actin

Here, we report on the structure and in situ location of arthrin (monoubiquitinated actin). Labelling of insect muscle thin filaments with a ubiquitin antibody reveals that every seventh subunit along the filament long-pitch helices is ubiquitinated. A three-dimensional reconstruction of frozen-hydrated arthrin filaments was produced. This was based on a novel algorithm that divides filament images into short segments that are used for single-particle image processing. Difference maps with an actin filament reconstruction locate ubiquitin at the side of actin sub-domain 1 opposite where myosin binds. Consistent with the reconstructions, peptide mapping places the ubiquitin linkage on lysine 118 in actin. Molecular modelling was used to generate arthrin monomers from ubiquitin and actin crystal structures. Filament models constructed from these monomers were compared with the arthrin reconstruction. The reconstruction suggests ubiquitin attached to Lys118 adopts one or a few conformers, stabilized by a small interface with actin. The function of actin ubiquitination is not known, but may involve regulation of muscle contractile activity.

Structure and function of CRSP/Med2; a promoter-selective transcriptional coactivator complex

The multi-subunit, human CRSP coactivator-also known as Mediator (Med)-regulates transcription by mediating signals between enhancer-bound factors (activators) and the core transcriptional machinery. Interestingly, different activators are known to bind distinct subunits within the CRSP/Med complex. We have isolated a stable, endogenous CRSP/Med complex (CRSP/Med2) that specifically lacks both the Med220 and the Med70 subunits. The three-dimensional structure of CRSP/Med2 was determined to 31 A resolution using electron microscopy and single-particle reconstruction techniques. Despite lacking both Med220 and Med70, CRSP/Med2 displays potent, activator-dependent transcriptional coactivator function in response to VP16, Sp1, and Sp1/SREBP-1a in vitro using chromatin templates. However, CRSP/Med2 is unable to potentiate activated transcription from a vitamin D receptor-responsive promoter, which requires interaction with Med220 for coactivator recruitment, whereas VDR-directed activation by CRSP/Med occurs normally. Thus, it appears that CRSP/Med may be regulated by a combinatorial assembly mechanism that allows promoter-selective function upon exchange of specific coactivator targets.

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

Distinct conformational states of nuclear receptor-bound CRSP-Med complexes

The human CRSP-Med coactivator complex is targeted by a diverse array of sequence-specific regulatory proteins. Using EM and single-particle reconstruction techniques, we recently completed a structural analysis of CRSP-Med bound to VP16 and SREBP-1a. Notably, these activators induced distinct conformational states upon binding the coactivator. Ostensibly, these different conformational states result from VP16 and SREBP-1a targeting distinct subunits in the CRSP-Med complex. To test this, we conducted a structural analysis of CRSP-Med bound to either thyroid hormone receptor (TR) or vitamin D receptor (VDR), both of which interact with the same subunit (Med220) of CRSP-Med. Structural comparison of TR- and VDR-bound complexes (at a resolution of 29 A) indeed reveals a shared conformational feature that is distinct from other known CRSP- Med structures. Importantly, this nuclear receptor-induced structural shift seems largely dependent on the movement of Med220 within the complex.

Location of divergent region 2 on the three-dimensional structure of cardiac muscle ryanodine receptor/calcium release channel

Ryanodine receptors (RyRs) are a family of calcium release channels found on intracellular calcium-handing organelles. Molecular cloning studies have identified three different RyR isoforms, which are 66-70% identical in amino acid sequence. In mammals, the three isoforms are encoded by three separate genes located on different chromosomes. The major variations among the isoforms occur in three regions, known as divergent regions 1, 2, and 3 (DR1, DR2, and DR3). In the present study, a modified RyR2 (cardiac isoform) cDNA was constructed, into which was inserted a green fluorescent protein (GFP)-encoding cDNA within DR2, specifically after amino acid residue Thr1366 (RyR2(T1366-GFP)). HEK293 cells expressing RyR2(T1366-GFP) cDNAs showed caffeine-sensitive and ryanodine-sensitive calcium release, demonstrating that RyR2(T1366-GFP) forms functional calcium release channels. Cells expressing RyR2(T1366-GFP) were identified readily by the characteristic fluorescence of GFP, indicating that the overall structure of the inserted GFP was retained. Cryo-electron microscopy (cryo-EM) of purified RyR2(T1366-GFP) showed structurally intact receptors, and a three-dimensional reconstruction was obtained by single-particle image processing. The location of the inserted GFP was obtained by comparing this three-dimensional reconstruction to one obtained for wild-type RyR2. The inserted GFP and, consequently Thr1366 within DR2, was mapped on the three-dimensional structure of RyR2 to domain 6, one of the characteristic cytoplasmic domains that form part of the multi-domain "clamp" regions of RyR2. The three-dimensional location of DR2 suggests that it plays roles in the RyR conformational changes that occur during channel gating, and possibly in RyR's interaction with the dihydropyridine receptor in excitation-contraction coupling. This study further demonstrates the feasibility and reliability of the GFP insertion/cryo-EM approach for correlating RyR's amino acid sequence with its three-dimensional structure, thereby enhancing our understanding of the structural basis of RyR function.

Gridding-based direct Fourier inversion of the three-dimensional ray transform

We describe a fast and accurate direct Fourier method for reconstructing a function f of three variables from a number of its parallel beam projections. The main application of our method is in single particle analysis, where the goal is to reconstruct the mass density of a biological macromolecule. Typically, the number of projections is extremely large, and each projection is extremely noisy. The projection directions are random and initially unknown. However, it is possible to determine both the directions and f by an iterative procedure; during each stage of the iteration, one has to solve a reconstruction problem of the type considered here. Our reconstruction algorithm is distinguished from other direct Fourier methods by the use of gridding techniques that provide an efficient means to compute a uniformly sampled version of a function g from a nonuniformly sampled version of Fg, the Fourier transform of g, or vice versa. We apply the two-dimensional reverse gridding method to each available projection of f, the function to be reconstructed, in order to obtain Ff on a special spherical grid. Then we use the three-dimensional gridding method to reconstruct f from this sampled version of Ff. This stage requires a proper weighting of the samples of Ff to compensate for their nonuniform distribution. We use a fast method for computing appropriate weights that exploits the special properties of the spherical sampling grid for Ff and involves the computation of a Voronoi diagram on the unit sphere. We demonstrate the excellent speed and accuracy of our method by using simulated data.

Recognition and separation of single particles with size variation by statistical analysis of their images

Macromolecules may occupy conformations with structural differences that cannot be resolved biochemically. The separation of mixed molecular populations is a pressing problem in single-particle analysis. Until recently, the task of distinguishing small structural variations was intractable, but developments in cryo-electron microscopy hardware and software now make it possible to address this problem. We have developed a general strategy for recognizing and separating structures of variable size from cryo-electron micrographs of single particles. The method uses a combination of statistical analysis and projection matching to multiple models. Identification of size variations by multivariate statistical analysis was used to do an initial separation of the data and generate starting models by angular reconstitution. Refinement was performed using alternate projection matching to models and angular reconstitution of the separated subsets. The approach has been successful at intermediate resolution, taking it within range of resolving secondary structure elements of proteins. Analysis of simulated and real data sets is used to illustrate the problems encountered and possible solutions. The strategy developed was used to resolve the structures of two forms of a small heat shock protein (Hsp26) that vary slightly in diameter and subunit packing.

Automated three-dimensional reconstruction of keyhole limpet hemocyanin type 1

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

Structure of the mammalian mitochondrial ribosome reveals an expanded functional role for its component proteins

The mitochondrial ribosome is responsible for the biosynthesis of protein components crucial to the generation of ATP in the eukaryotic cell. Because the protein:RNA ratio in the mitochondrial ribosome (approximately 69:approximately 31) is the inverse of that of its prokaryotic counterpart (approximately 33:approximately 67), it was thought that the additional and/or larger proteins of the mitochondrial ribosome must compensate for the shortened rRNAs. Here, we present a three-dimensional cryo-electron microscopic map of the mammalian mitochondrial 55S ribosome carrying a tRNA at its P site, and we find that instead, many of the proteins occupy new positions in the ribosome. Furthermore, unlike cytoplasmic ribosomes, the mitochondrial ribosome possesses intersubunit bridges composed largely of proteins; it has a gatelike structure at its mRNA entrance, perhaps involved in recruiting unique mitochondrial mRNAs; and it has a polypeptide exit tunnel that allows access to the solvent before the exit site, suggesting a unique nascent-polypeptide exit mechanism.

RNA Polymerase II/TFIIF Structure and Conserved Organization of the Initiation Complex

The structure of an RNA polymerase II/general transcription factor TFIIF complex was determined by cryo-electron microscopy and single particle analysis. Density due to TFIIF was not concentrated in one area but rather was widely distributed across the surface of the polymerase. The largest subunit of TFIIF interacted with the dissociable Rpb4/Rpb7 polymerase subunit complex and with the mobile "clamp." The distribution of the second largest subunit of TFIIF was very similar to that previously reported for the sigma subunit in the bacterial RNA polymerase holoenzyme, consisting of a series of globular domains extending along the polymerase active site cleft. This result indicates that the second TFIIF subunit is a true structural homolog of the bacterial sigma factor and reveals an important similarity of the transcription initiation mechanism between bacteria and eukaryotes. The structure of the RNAPII/TFIIF complex suggests a model for the organization of a minimal transcription initiation complex.

Automatic particle detection through efficient Hough transforms

Manual selection of single particles in images acquired using cryo-electron microscopy (cryoEM) will become a significant bottleneck when a very large number of images are required to achieve three-dimensional reconstructions at near atomic resolution. Investigation of fast, accurate approaches for automatic particle detection has become one of the current challenges in the cryoEM community. At the same time, the investigation is hampered by the fact that few benchmark particles or image datasets exist in the community. The unavailability of such data makes it difficult to evaluate newly developed algorithms and to leverage expertise from other disciplines. The paper presents our recent contribution to this effort. It also describes our newly developed computational framework for particle detection, through the application of edge detection and a sequence of ordered Hough transforms. Experimental results using keyhole limpet hemocyanin (KLH) as a model particle are very promising. In addition, it introduces a newly established web site, designed to support the investigation of automatic particle detection by providing an annotated image dataset of KLH available to the general scientific community.

Revisiting the structure of Alvinella pompejana hemoglobin at 20A resolution by cryoelectron microscopy

The hemoglobin of the polychaete worm Alvinella pompejana was reconstructed at 20A resolution from frozen-hydrated samples observed by electron microscopy according to the random conical tilt series method. This three-dimensional reconstruction was mirror-inverted with respect to a previous volume published by de Haas et al. in 1996. In order to explain this handedness discrepancy, various 3D reconstructions using different reference volumes were carried out showing that the choice of the first volume was the keystone during the refinement process. The 3D reconstruction volume of A. pompejana Hb presented structural features characteristic of annelid Hbs with two hexagonal layers each comprising six hollow globular subassemblies and a complex of non-heme linker chains. Moreover, the eclipsed conformation of the two hexagonal layers and a HGS architecture similar to that described for Arenicola marina Hb led to the conclusion that A. pompejana Hb belonged to the architectural type II according to the definition of Jouan et al. (2001). A comparison between this cryo-EM volume and X-ray crystallography density maps of Lumbricus terrestris type-I Hb (Royer et al., 2000) showed that the triple stranded coiled coil structures of linker chains were different. Based on this observation, a model was proposed to explain the eclipsed conformation of the two hexagonal layers of type-II Hbs.

An approach to examining model dependence in EM reconstructions using cross-validation

Reference bias refers to a common problem in fitting experimental data to an initial model. Given enough free parameters, a good fit of any experimental data to the model can be obtained, even if the experimental data contain only noise. Reference-based alignment methods used in electron microscopy (EM) are subject to this type of bias, in that images containing pure noise can regenerate the reference. Cross-validation is based on the idea that the experimental data used to assess the validity of the fitting should not be the same data as were used to do the fitting. Here we present the application of cross-validation to one form of reference-based alignment: 3D-projection matching in single-particle reconstructions. Our results show that reference bias is indeed present in reconstructions, but that the effect is small for real data compared to that for random noise, and that this difference in behavior is magnified, rather than diminished, during iterative refinement.

Three-dimensional localization of divergent region 3 of the ryanodine receptor to the clamp-shaped structures adjacent to the FKBP binding sites

Of the three divergent regions of ryanodine receptors (RyRs), divergent region 3 (DR3) is the best studied and is believed to be involved in excitation-contraction coupling as well as in channel regulation by Ca(2+) and Mg(2+). To gain insight into the structural basis of DR3 function, we have determined the location of DR3 in the three-dimensional structure of RyR2. We inserted green fluorescent protein (GFP) into the middle of the DR3 region after Thr-1874 in the sequence. HEK293 cells expressing this GFP-RyR2 fusion protein, RyR2(T1874-GFP,) were readily detected by their green fluorescence, indicating proper folding of the inserted GFP. RyR2(T1874-GFP) was further characterized functionally by assays of Ca(2+) release and [(3)H]ryanodine binding. These analyses revealed that RyR2(T1874-GFP) functions as a caffeine- and ryanodine-sensitive Ca(2+) release channel and displays Ca(2+) dependence and [(3)H]ryanodine binding properties similar to those of the wild type RyR2. RyR2(T1874-GFP) was purified from cell lysates in a single step by affinity chromatography using GST-FKBP12.6 as the affinity ligand. The three-dimensional structure of the purified RyR2(T1874-GFP) was then reconstructed using cryoelectron microscopy and single particle image analysis. Comparison of the three-dimensional reconstructions of wild type RyR2 and RyR2(T1874-GFP) revealed the location of the inserted GFP, and hence the DR3 region, in one of the characteristic domains of RyR, domain 9, in the clamp-shaped structure adjacent to the FKBP12 and FKBP12.6 binding sites. COOH-terminal truncation analysis demonstrated that a region between 1815 and 1855 near DR3 is essential for GST-FKBP12.6 binding. These results provide a structural basis for the role of the DR3 region in excitation-contraction coupling and in channel regulation.

Three-dimensional reconstruction of the Saccharomyces cerevisiae multidrug resistance protein Pdr5p

Pdr5p, the major multidrug exporter in Saccharomyces cerevisiae, is a member of the ATP-binding cassette (ABC) superfamily. Pdr5p shares similar mechanisms of substrate recognition and transport with the human MDR1-Pgp, despite an inverted topology of transmembrane and ATP-binding domains. The hexahistidine-tagged Pdr5p multidrug transporter was highly overexpressed in yeast strains where other ABC genes have been deleted. After solubilization and purification, the 160-kDa recombinant Pdr5p has been reconstituted into a lipid bilayer. Controlled detergent removal from Pdr5p-lipid-detergent micelles allowed the production of peculiar square-shaped particles coexisting with liposomes and proteoliposomes. These particles having 11 nm in side were well suited for single particle analysis by electron microscopy. From such analysis, a computed volume has been determined at 25-A resolution, giving insight into the structural organization of Pdr5p. Comparison with the reported structures of different bacterial ABC transporters was consistent with a dimeric organization of Pdr5p in the square particles. Each monomer was composed of three subregions corresponding to a membrane region of about 50 A in height that joins two well separated protruding stalks of about 40 A in height, ending each one with a cytoplasmic nucleotide-binding domain (NBD) lobe of about 50-60 A in diameter. The three-dimensional reconstruction of Pdr5p revealed a close arrangement and a structural asymmetric organization of the two NBDs that appeared oriented perpendicularly within a monomer. The existence of different angular positions of the NBDs, with respect to the stalks, suggest rotational movements during the catalytic cycle.

The variance of icosahedral virus models is a key indicator in the structure determination: a model-free reconstruction of viruses, suitable for refractory particles

A model-free method to determine the three-dimensional structure of icosahedral viruses is described. The novel strategy is based upon the approximate principle that correct virus structures have high variance as do all other well-detailed structures, even wrong ones. The original projections of individual particles are reduced to a radius of 25 pixels and are used to compute single particle reconstruction models by assigning them 1800 different Euler triads. The variance of the models obtained from all projections is stored in maps and a decimation process is carried out. In a first stage, thresholds are adopted for the variance values, and in a second stage, carried out by correspondence analysis and classification, 30 clusters of models are sorted out. The clusters are refined to yield models contained in boxes of 64(3) voxels. The refined models with highest variance and closest similarity represent the correct solution. Once enlarged, these models can be used to align all available projections in their original scale in a customary projection-matching process. The method has proved successful in determining the structures of poliovirus, of the empty and filled capsids of L-A virus, and of a modified capsid of hepatitis B virus.

The prepower stroke conformation of myosin V

We have used electron microscopy and single-particle image processing to study head conformation in myosin V molecules. We find that in the presence of ATP, many heads have a sharply angled conformation that is rare in its absence. The sharply angled conformation is similar to a myosin II atomic structure proposed to mimic the prepower stroke state. The leading head in molecules attached to actin by both heads has a similar conformation, but is also sharply angled in a second plane by tethering through the trail head. The lead head lever joins the motor domain approximately 5 nm axially from where it joins the trail motor. These positions locate the converter subdomain and show the lead motor is in the prepower stroke conformation. Tethering by the trail head places the lead head motor domain at the correct axial position along the actin for binding, but at the wrong orientation. Attachment is achieved either by bending the lead head lever throughout its length or at the pliant point. The microscopy shows that most of the walking stride is produced by changes in lever angle brought about by converter movement, but is augmented by distortion produced by thermal energy.

Three-dimensional reconstruction of the recombinant type 2 ryanodine receptor and localization of its divergent region 1

Isoform 2 of the ryanodine receptor (RyR2) is the major calcium release channel in cardiac muscle. In the present study, two kinds of RyR2 cDNA were constructed, one encoding the wild type mouse RyR2 (RyR2(wt)) and the other encoding modified RyR2, into which was inserted a cDNA encoding green fluorescent protein (GFP). GFP was inserted into the divergent region 1 (DR1) of RyR2, after the Asp-4365 (RyR2(D4365-GFP)). HEK293 cells expressing both RyR2(wt) and RyR2(D4365-GFP) cDNAs showed caffeine- and ryanodine-sensitive calcium release, demonstrating that both wild type and modified RyR2s form functional calcium release channels. Cells expressing the fusion protein, RyR2(D4365-GFP), were readily identified by their fluorescence due to the presence of GFP, indicating that the inserted GFP folded properly. Both expressed RyR2s were purified from cell lysates in a single step by affinity chromatography using a GST-FKBP12.6 as the affinity ligand. Cryoelectron microscopy of purified RyR2s showed structurally intact receptors, and three-dimensional reconstructions were obtained by single particle image processing. The three-dimensional reconstruction of RyR2(wt) appeared very similar to that of the native RyR2 purified from dog heart. The location of the inserted GFP, and consequently of DR1, was mapped on the three-dimensional structure of RyR2 to one of the subunit's characteristic domains, domain 3, also known as the "handle" domain. This study describes the first internal fusion of a protein into a ryanodine receptor, and it demonstrates the potential of this technology for localizing functional and structural domains on the three-dimensional structure of RyR.

Structural analysis of the RSC chromatin-remodeling complex

Electron microscopy of the RSC chromatin-remodeling complex reveals a ring of protein densities around a central cavity. The size and shape of the cavity correspond closely to those of a nucleosome. Results of nuclease protection analysis are consistent with nucleosome binding in the cavity. Such binding could explain the ability of RSC to expose nucleosomal DNA in the presence of ATP without loss of associated histones.

Microtubule structure at 8 A resolution

We have obtained a 3D reconstruction of intact microtubules, using cryoelectron microscopy and image processing, at a resolution of about 8 A, sufficient to resolve much of the secondary structure. The level of detail in the map allows docking of the tubulin structure previously determined by electron crystallography, with very strong constraints, providing several important insights not previously available through docking tubulin into lower-resolution maps. This work provides an improved picture of the interactions between adjacent protofilaments, which are responsible for microtubule stability, and also suggests that some structural features are different in microtubules from those in the zinc sheets with which the tubulin structure was determined.

Single-particle imaging of macromolecules by cryo-electron microscopy

Cryo-electron microscopy (cryo-EM) of biological molecules in single-particle (i.e., unordered, nonaggregated) form is a new approach to the study of molecular assemblies, which are often too large and flexible to be amenable to X-ray crystallography. New insights into biological function on the molecular level are expected from cryo-EM applied to the study of such complexes "trapped" at different stages of their conformational changes and dynamical interactions. Important molecular machines involved in the fundamental processes of transcription, mRNA splicing, and translation are examples for successful applications of the new technique, combined with structural knowledge gained by conventional techniques of structure determination, such as X-ray crystallography and NMR.

Cryo-EM reveals an active role for aminoacyl-tRNA in the accommodation process

During the elongation cycle of protein biosynthesis, the specific amino acid coded for by the mRNA is delivered by a complex that is comprised of the cognate aminoacyl-tRNA, elongation factor Tu and GTP. As this ternary complex binds to the ribosome, the anticodon end of the tRNA reaches the decoding center in the 30S subunit. Here we present the cryo- electron microscopy (EM) study of an Escherichia coli 70S ribosome-bound ternary complex stalled with an antibiotic, kirromycin. In the cryo-EM map the anticodon arm of the tRNA presents a new conformation that appears to facilitate the initial codon-anticodon interaction. Furthermore, the elbow region of the tRNA is seen to contact the GTPase-associated center on the 50S subunit of the ribosome, suggesting an active role of the tRNA in the transmission of the signal prompting the GTP hydrolysis upon codon recognition.

Efficiency of 2D alignment methods

In single particle analysis, the alignment of two-dimensional images is a fundamental step aimed at bringing into register various particle views of biological macromolecules observed with the electron microscope. The computational efficiency of this step is a deciding factor in design of alignment strategies for large sets of noisy data and in development of three-dimensional structure refinement methods. In addition, the accuracy of the alignment method varies depending on the numerical solutions adopted to efficiently perform exhaustive searches for three orientation parameters. The selected alignment methods are analyzed in terms of their computational complexity and the estimates of numbers of arithmetic operations for each method are given. The tests of alignment accuracy are performed using images simulated in accordance with the linear theory of image formation in the electron microscope. It is demonstrated that the efficiency of the alignment methods can be improved if approximate centers of gravity of particle views are known. The accuracy of the methods considered is largely affected, particularly for high noise levels, by the order in which interpolation steps are applied.

Human CRSP interacts with RNA polymerase II CTD and adopts a specific CTD-bound conformation

Activation of gene transcription in mammalian cells requires several classes of coactivators that participate in different steps of the activation cascade. Using conventional and affinity chromatography, we have isolated a human coactivator complex that interacts directly with the C-terminal domain (CTD) of RNA polymerase II (Pol II). The CTD-binding complex is structurally and functionally indistinguishable from our previously isolated CRSP coactivator complex. The closely related, but transcriptionally inactive, ARC-L complex failed to interact with the CTD, indicating a significant biochemical difference between CRSP and ARC-L that may, in part, explain their functional divergence. Electron microscopy and three-dimensional single-particle reconstruction reveals a conformation for CTD-CRSP that is structurally distinct from unliganded CRSP or CRSP bound to SREBP-1a, but highly similar to CRSP bound to the VP16 activator. Together, our findings suggest that the human CRSP coactivator functions, at least in part, by mediating activator-dependent recruitment of RNA Pol II via the CTD.

Cryo-negative staining reduces electron-beam sensitivity of vitrified biological particles

Beam damage is the main resolution-limiting factor when biological particles are observed by cryoelectron microscopy in a thin vitrified solution film. Furthermore, the low contrast of the specimen frequently makes observation difficult and limits the possibility of image processing. Cryo-negative staining, in which the particles are vitrified in a thin layer of concentrated ammonium molybdate solution, makes it possible to visualize the particles with a much better signal-to-noise ratio (SNR) while keeping the specimen in a good state of preservation. We have observed the Escherichia coli GroEL chaperonin, prepared in a native vitrified solution and by cryo-negative staining after electron exposure from 1000 to 3000e(-)/nm(2). We have compared the resulting three-dimensional models obtained from these different conditions and have tested their fit with the atomic model of the protein subunit obtained from X-ray crystallography. It is found that, down to 1.5-nm resolution, the particles appear to be faithfully represented in the cryo-negatively stained preparation, but there is an approximately 10-fold increase of SNR compared with the native vitrified preparation. Furthermore, for the same range of irradiation and down to the same resolution, the particles seem unaffected by beam damage, whereas the damage is severe in the native vitrified particles.

Three-dimensional spectral signal-to-noise ratio for a class of reconstruction algorithms

A three-dimensional (3D) version of the spectral signal-to-noise ratio (SSNR)-based resolution measure is introduced. The measure is defined for a class of 3D reconstruction algorithms that use interpolation in Fourier space. The statistical properties of the SSNR are discussed and related to the properties of another resolution measure, the Fourier shell correlation (FSC). The new measure was tested on 3D structures calculated from a simulated set of quasi-evenly spaced 2D projections using a nearest-neighbor interpolation and a gridding algorithm. In the latter case, the results agree very well with the FSC-based estimate, with the exception of very high SSNR values. The main applicability of the 3D SSNR is tomography, where due to the small number of projections collected, FSC cannot be used. The new measure was applied to three sets of tomographic data. It was demonstrated that the measure is sufficiently sensitive to yield theoretically expected results. Therefore, the 3D SSNR opens up the possibility of evaluating the quality of tomographic reconstructions in an objective manner. The 3D distribution of SSNR is of major interest in single-particle analysis. It is shown that the new measure can be used to evaluate the anisotropy of 3D reconstructions. The distribution of SSNR is characterized by three anisotropy indices derived from principal axes of the 3D inertia covariance matrix of the SSNR. These indices are used to construct a 3D Fourier filter which, when applied to a 3D reconstruction of a macromolecule, maximizes the SNR in real space and minimizes real-space artifacts caused by uneven distribution of 2D projections.

Striking conformational change suspected within the phosphoribulokinase dimer induced by interaction with GAPDH

A multitechnique approach was used to study the [glyceraldehyde-3-phosphate dehydrogenase](2 x 4)-[phosphoribulokinase](2 x 2) multienzymatic complex of the alga Chlamydomonas reinhardtii. On the one hand, each component of the complex was compared with known atomic structures of related enzymes or of similar enzymes originating from different organisms. On the other hand, the overall low resolution architecture of the whole complex was studied using cryoelectron microscopy and image processing techniques. The dimers of phosphoribulokinase are suspected to undergo a dramatic change in activity during a cycle of binding and detaching from tetramers of glyceraldehyde-3-phosphate dehydrogenase. This is likely supported by strong structural differences between the modeled phosphoribulokinase dimers and the counterpart in the three-dimensional reconstruction volume of the whole complex obtained from cryoelectron microscope images.

Three-dimensional structure by cryo-electron microscopy of YvcC, an homodimeric ATP-binding cassette transporter from Bacillus subtilis

YvcC, a multidrug transporter from Bacillus subtilis, is a member of the ATP-binding cassette superfamily, highly homologous to each half of human multidrug-resistance P-glycoprotein and to several other bacterial half-ABC transporters. Here, the purified recombinant histidine-tagged YvcC has been reconstituted into a lipid bilayer. Controlled and partial detergent removal from YvcC-lipid micelles allowed the production of particularly interesting lipid-detergent-YvcC ring-shaped particles, about 40 nm in diameter, well suited for single particle analysis by cryo-electron microscopy. Furthermore, binding of these histidine-tagged ring-shaped particles to lipid layers functionalized with a Ni(2+)-chelating head group generated a preferential perpendicular orientation, eliminating the missing cone in the final three-dimensional reconstruction. From such analysis, a computed volume has been determined to 2.5 nm resolution giving a detailed insight into the structural organization of this half-ABC transporter within a membrane. The repetitive unit in the ring-shaped particles is consistent with a homodimeric organization of YvcC. Each subunit was composed of three domains: a 5 nm height transmembrane region, a stalk of about 4 nm in height and 2 nm in diameter, and a cytoplasmic lobe of about 5-6 nm in diameter. The latest domain, which fitted with the reported X-ray structure of HisP, was identified as the nucleotide-binding domain (NBD). The 3D reconstruction of the YvcC homodimer well compared with the very recent X-ray crystallographic data on the MsbA homodimer from Escherichia coli, supporting the existence of a central open chamber between the two subunits constituting the homodimer. In addition, the 3D reconstruction of YvcC embedded in a membrane revealed an asymmetric organization of the two NBDs sites within the homodimer, as well as a dimeric interaction between two homodimers.

Apocalmodulin and Ca2+-calmodulin bind to neighboring locations on the ryanodine receptor

Calmodulin (CaM) binds to the ryanodine receptor/calcium release channel of skeletal muscle (RyR1), both in the absence and presence of Ca(2+), and regulates the activity of the channel activity by activating and inhibiting it, respectively. Using cryo-electron microscopy and three-dimensional reconstruction, we found that one apoCaM binds per RyR1 subunit along the sides of the cytoplasmic assembly of the receptor. This location is distinct from but close to the location found for Ca(2+)-CaM, providing a structural basis for efficient switching of CaM between these two positions with the oscillating intracellular Ca(2+) concentration that generates muscle relaxation/contraction cycles. The locations of apoCaM and Ca(2+)-CaM at a critical region for RYR1-dihydropyridine receptor interaction are suggestive of a direct role for CaM in the mechanism of excitation-contraction coupling.

Three-dimensional structure of phosphorylase kinase at 22 A resolution and its complex with glycogen phosphorylase b

Phosphorylase kinase (PhK) integrates hormonal and neuronal signals and is a key enzyme in the control of glycogen metabolism. PhK is one of the largest of the protein kinases and is composed of four types of subunit, with stoichiometry (alphabetagammadelta)(4) and a total MW of 1.3 x 10(6). PhK catalyzes the phosphorylation of inactive glycogen phosphorylase b (GPb), resulting in the formation of active glycogen phosphorylase a (GPa) and the stimulation of glycogenolysis. We have determined the three-dimensional structure of PhK at 22 A resolution by electron microscopy with the random conical tilt method. We have also determined the structure of PhK decorated with GPb at 28 A resolution. GPb is bound toward the ends of each of the lobes with an apparent stoichiometry of four GPb dimers per (alphabetagammadelta)(4) PhK. The PhK/GPb model provides an explanation for the formation of hybrid GPab intermediates in the PhK-catalyzed phosphorylation of GPb.

Single particle macromolecular structure determination via electron microscopy

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

Localization of L11 protein on the ribosome and elucidation of its involvement in EF-G-dependent translocation

L11 protein is located at the base of the L7/L12 stalk of the 50 S subunit of the Escherichia coli ribosome. Because of the flexible nature of the region, recent X-ray crystallographic studies of the 50 S subunit failed to locate the N-terminal domain of the protein. We have determined the position of the complete L11 protein by comparing a three-dimensional cryo-EM reconstruction of the 70 S ribosome, isolated from a mutant lacking ribosomal protein L11, with the three-dimensional map of the wild-type ribosome. Fitting of the X-ray coordinates of L11-23 S RNA complex and EF-G into the cryo-EM maps combined with molecular modeling, reveals that, following EF-G-dependent GTP hydrolysis, domain V of EF-G intrudes into the cleft between the 23 S ribosomal RNA and the N-terminal domain of L11 (where the antibiotic thiostrepton binds), causing the N-terminal domain to move and thereby inducing the formation of the arc-like connection with the G' domain of EF-G. The results provide a new insight into the mechanism of EF-G-dependent translocation.

Cryo-electron microscopy as an investigative tool: the ribosome as an example

Cryo-electron microscopy allows the visualization of macromolecules in their native state. Combined with techniques of three-dimensional reconstruction, cryo-EM images of single molecules can be used to study macromolecular interactions. The ribosome, a large RNA-protein complex with multiple binding interactions, is an excellent test case illustrating the power of these new techniques. Conformational changes during the binding of tRNA and protein factors to the ribosome can now be studied without the interference of crystal packing. Now that the first X-ray structures of ribosomal subunits have become available, conformational changes observed by cryo-EM in different functional states can be traced back to internal rearrangements of the underlying structural framework. Electron microscopy, X-ray crystallography, and modeling should be used together in the endeavor to understand the functioning of the translational machinery.

Three-dimensional reconstruction of the recombinant type 3 ryanodine receptor and localization of its amino terminus

Recombinant type 3 ryanodine receptor (RyR3) has been purified in quantities sufficient for structural characterization by cryoelectron microscopy and three-dimensional (3D) reconstruction. Two cDNAs were prepared and expressed in HEK293 cells, one encoding the wild-type RyR3 and the other encoding RyR3 containing glutathione S-transferase (GST) fused to its amino terminus (GST-RyR3). RyR3 was purified from detergent-solubilized transfected cells by affinity chromatography using 12.6-kDa FK506-binding protein in the form of a GST fusion as the affinity ligand. Purification of GST-RyR3 was achieved by affinity chromatography by using glutathione-Sepharose. Purified recombinant RyR3 and GST-RyR3 proteins exhibited high-affinity [(3)H]ryanodine binding that was sensitive to activation by Ca(2+) and caffeine and to inhibition by Mg(2+). 3D reconstructions of both recombinant RyR3 and GST-RyR3 appeared very similar to that of the native RyR3 purified from bovine diaphragm. Comparison of the 3D reconstructions of RyR3 and GST-RyR3 revealed that the GST domains and, hence, the amino termini of the RyR3 subunits are located in the "clamp" structures that form the corners of the square-shaped cytoplasmic region of homotetrameric RyR3. This study describes the 3D reconstruction of a recombinant ryanodine receptor and it demonstrates the potential of this technology for characterizing functional and structural perturbations introduced by site-directed mutagenesis.

The structure of the multidrug resistance protein 1 (MRP1/ABCC1). crystallization and single-particle analysis

Multidrug resistance protein 1 (MRP1/ABCC1) is an ATP-binding cassette (ABC) polytopic membrane transporter of considerable clinical importance that confers multidrug resistance on tumor cells by reducing drug accumulation by active efflux. MRP1 is also an efficient transporter of conjugated organic anions. Like other ABC proteins, including the drug resistance conferring 170-kDa P-glycoprotein (ABCB1), the 190-kDa MRP1 has a core structure consisting of two membrane-spanning domains (MSDs), each followed by a nucleotide binding domain (NBD). However, unlike P-glycoprotein and most other ABC superfamily members, MRP1 contains a third MSD with five predicted transmembrane segments with an extracytosolic NH(2) terminus. Moreover, the two nucleotide-binding domains of MRP1 are considerably more divergent than those of P-glycoprotein. In the present study, the first structural details of MRP1 purified from drug-resistant lung cancer cells have been obtained by electron microscopy of negatively stained single particles and two-dimensional crystals formed after reconstitution of purified protein with lipids. The crystals display p2 symmetry with a single dimer of MRP1 in the unit cell. The overall dimensions of the MRP1 monomer are approximately 80 x 100 A. The MRP1 monomer shows some pseudo-2-fold symmetry in projection, and in some orientations of the detergent-solubilized particles, displays a stain filled depression (putative pore) appearing toward the center of the molecule, presumably to enable transport of substrates. These data represent the first structural information of this transporter to approximately 22-A resolution and provide direct structural evidence for a dimeric association of the transporter in a reconstituted lipid bilayer.

Strul--a method for 3D alignment of single-particle projections based on common line correlation in Fourier space

A central problem of 3D reconstruction in single-particle electron microscopy is the determination of relative orientations of the individual projections contributing to the reconstruction. This article describes an implementation of the method of common lines correlation in Fourier space that allows generation of common lines between an arbitrary number of projections which might posses an arbitrary point group symmetry. Based on this method, it is possible to optimize rotational and translational alignment parameters for individual single-particle projections. The underlying philosophy and details of implementation are discussed, and as an illustration a 3D reconstruction in ice of peroxisomal alcohol oxidase from Pichia pastoris, an octameric assembly with 422-symmetry and a molecular weight of 592 kDa is presented.

Classification and reconstruction of a heterogeneous set of electron microscopic images: a case study of GroEL-substrate complexes

Image analysis methods were used to separate images of a large macromolecular complex, the chaperonin GroEL, in a preparation in which it is partially liganded to a nonnative protein substrate, glutamine synthetase. The relatively small difference ( approximately 6%) in size between the chaperonin in its free and complexed forms, and the absence of gross changes in overall conformation, made separation of the two types of particles challenging. Different approaches were evaluated and used for alignment and classification of images, both in two common projections and in three dimensions, yielding 2D averages and a 3D reconstruction. The results of 3D analysis describe the conformational changes effected by binding of this particular protein substrate and demonstrate the utility of 2D analysis as an indicator of structural change in this system.

Lumbricus terrestris hemoglobin--the architecture of linker chains and structural variation of the central toroid

The extracellular giant hemoglobin from the earthworm Lumbricus terrestris was reconstructed at 14.9-A resolution from cryo-electron microscope images, using a new procedure for estimating parameters of the contrast transfer (CTF) function. In this approach, two important CTF parameters, defocus and amplitude contrast ratio, can be refined iteratively within the framework of 3D projection alignment procedure, using minimization of sign disagreement between theoretical CTF and cross-resolution curves. The 3D cryo-EM map is in overall good agreement with the recent X-ray crystallography map of Royer et al. (2000, Proc. Natl. Acad. Sci. USA 97, 7107-7111), and it reveals the local threefold arrangement of the three linker chains present within each 1/12 of the complex. The 144 globin chains and 36 linker chains within the complex are clearly visible, and the interdigitation of the 12 coiled-coil helical spokes forming the central toroidal piece is confirmed. Based on these findings, two mechanisms of the dodecameric unit assembly are proposed and termed "zigzag" and "pairwise" polymerizations. However, the detection by cryo-EM of 12 additional rod-like bodies within the toroid raises the possibility that the architecture of the toroid is more complex than previously thought or that yet unknown ligands or allosteric effectors for this oxygen carrier are present.

Occurrence of two architectural types of hexagonal bilayer hemoglobin in annelids: comparison of 3D reconstruction volumes of Arenicola marina and Lumbricus terrestris hemoglobins

A 3D reconstruction at 25 A resolution of native hemoglobin of the polychaete worm Arenicola marina was carried out from frozen-hydrated specimens examined in the electron microscope. The reconstruction volume of this large extracellular multimeric respiratory pigment appears as a hexagonal bilayer structure with eclipsed vertices in its upper and lower hexagonal layers. Conversely, in hemoglobins of oligochaetes, achaetes, and vestimentiferans and in chlorocruorins of the Sabellidae (polychaete) family, the vertices of the upper layer are 16 degrees clockwise rotated with respect to those of the lower layer. The fact that two other polychaete hemoglobins (Alvinella pompejana and Tylorrhynchus heterochaetus) have the same architecture as Arenicola led us to define two types of hexagonal bilayer hemoglobins/chlorocruorins: (i) type-I present in oligochaete, achaete, and vestimentiferan hemoglobins and in Sabellidae chlorocruorins; and (ii) type-II present in polychaete hemoglobins. A comparative study of the hemoglobins of Lumbricus terrestris (type-I) and Arenicola marina (type-II) showed that only two small differences located in the c4 and c5 linking units are responsible of the important architectural difference present in oligomers. A likely scheme proposed to explain the phylogenic distribution of the two types suggests that Clitellata, Sabellida (polychaete), and vestimentiferan hemoglobins and chlorocruorins derive from a type-I ancestral molecule, while Terebellida (Alvinella), Phyllodocida (Tylorrhynchus), and Scolecida (Arenicola) and possibly other polychaetes derive from an ancestor molecule with type-II hemoglobin. The architectures of the hollow globular substructures are highly similar in Arenicola and Lumbricus hemoglobins, with 12 globin chains and three linking units (c3a, c3b, and c4). The central piece of Arenicola hemoglobin is an ellipsoid while that of Lumbricus is a toroid. No phylogenic correlation could be found between the structure of the central pieces and the architecture type.