High performance computing in structural determination by electron cryomicroscopy

Computational advances have significantly contributed to the current role of electron cryomicroscopy (cryoEM) in structural biology. The needs for computational power are constantly growing with the increasing complexity of algorithms and the amount of data needed to push the resolution limits. High performance computing (HPC) is becoming paramount in cryoEM to cope with those computational needs. Since the nineties, different HPC strategies have been proposed for some specific problems in cryoEM and, in fact, some of them are already available in common software packages. Nevertheless, the literature is scattered in the areas of computer science and structural biology. In this communication, the HPC approaches devised for the computation-intensive tasks in cryoEM (single particles and tomography) are retrospectively reviewed and the future trends are discussed. Moreover, the HPC capabilities available in the most common cryoEM packages are surveyed, as an evidence of the importance of HPC in addressing the future challenges.     

Structural characterization of components of protein assemblies by comparative modeling and electron cryo-microscopy

We explore structural characterization of protein assemblies by a combination of electron cryo-microscopy (cryoEM) and comparative protein structure modeling. Specifically, our method finds an optimal atomic model of a given assembly subunit and its position within an assembly by fitting alternative comparative models into a cryoEM map. The alternative models are calculated by MODELLER [J. Mol. Biol. 234 (1993) 313] from different sequence alignments between the modeled protein and its template structures. The fitting of these models into a cryoEM density map is performed either by FOLDHUNTER [J. Mol. Biol. 308 (2001) 1033] or by a new density fitting module of MODELLER (Mod-EM). Identification of the most accurate model is based on the correlation between the model accuracy and the quality of fit into the cryoEM density map. To quantify this correlation, we created a benchmark consisting of eight proteins of different structural folds with corresponding density maps simulated at five resolutions from 5 to 15 angstroms, with three noise levels each. Each of the proteins in the set was modeled based on 300 different alignments to their remotely related templates (12-32% sequence identity), spanning the range from entirely inaccurate to essentially accurate alignments. The benchmark revealed that one of the most accurate models can usually be identified by the quality of its fit into the cryoEM density map, even for noisy maps at 15 angstroms resolution. Therefore, a cryoEM density map can be helpful in improving the accuracy of a comparative model. Moreover, a pseudo-atomic model of a component in an assembly may be built better with comparative models of the native subunit sequences than with experimentally determined structures of their homologs.     

Web portal to an image database for high-resolution three-dimensional reconstruction

The exponential increase of image data in high-resolution reconstructions by electron cryomicroscopy (cryoEM) has posed a need for efficient data management solutions in addition to powerful data processing procedures. Although relational databases and web portals are commonly used to manage sequences and structures in biological research, their application in cryoEM has been limited due to the complexity in accomplishing the dual tasks of interacting with proprietary software and simultaneously providing data access to users without database knowledge. Here, we report our results in developing web portal to SQL image databases used by the Image Management and Icosahedral Reconstruction System (IMIRS) to manage cryoEM images for subnanometer-resolution reconstructions. Fundamental issues related to the design and deployment of web portals to image databases are described. A web browser-based user interface was designed to accomplish data reporting and other database-related services, including user authentication, data entry, graph-based data mining, and various query and reporting tasks with interactive image manipulation capabilities. With an integrated web portal, IMIRS represents the first cryoEM application that incorporates both web-based data reporting tools and a complete set of data processing modules. Our examples should thus provide general guidelines applicable to other cryoEM technology development efforts.     

On the structures of filamentous bacteriophage Ff (fd,f1, M13)

The filamentous bacteriophage (Inovirus) strain Ff (fd, f1, M13) is widely used in molecular biophysics as a simple model system. A low resolution molecular model of the fd protein coat has been reported, derived from iterative helical real space reconstruction of cryo-electron micrographs (cryoEM). This model is significantly different from the model previously derived from X-ray fibre diffraction and solid-state NMR. We show that the cryoEM model agrees neither with solid-state NMR data nor with X-ray fibre diffraction data of fd, and has some puzzling structural features, for instance nanometre holes through the protein coat. We refine the cryoEM model against the X-ray data, and find that the model after refinement closely approximates the model derived directly from X-ray fibre diffraction and solid-state NMR data. We suggest possible reasons for the differences between the models derived from cryoEM and X-ray diffraction.     

Automatic particle selection: results of a comparative study

Manual selection of single particles in images acquired using cryo-electron microscopy (cryoEM) will become a significant bottleneck when datasets of a hundred thousand or even a million particles are required for structure determination at near atomic resolution. Algorithm development of fully automated particle selection is thus an important research objective in the cryoEM field. A number of research groups are making promising new advances in this area. Evaluation of algorithms using a standard set of cryoEM images is an essential aspect of this algorithm development. With this goal in mind, a particle selection "bakeoff" was included in the program of the Multidisciplinary Workshop on Automatic Particle Selection for cryoEM. Twelve groups participated by submitting the results of testing their own algorithms on a common dataset. The dataset consisted of 82 defocus pairs of high-magnification micrographs, containing keyhole limpet hemocyanin particles, acquired using cryoEM. The results of the bakeoff are presented in this paper along with a summary of the discussion from the workshop. It was agreed that establishing benchmark particles and using bakeoffs to evaluate algorithms are useful in promoting algorithm development for fully automated particle selection, and that the infrastructure set up to support the bakeoff should be maintained and extended to include larger and more varied datasets, and more criteria for future evaluations.     

Virus-like particles of a fish nodavirus display a capsid subunit domain organization different from that of insect nodaviruses

The structure of recombinant virus-like particles of malabaricus grouper nervous necrosis virus (MGNNV), a fish nodavirus isolated from the grouper Epinephelus malabaricus, was determined by electron cryomicroscopy (cryoEM) and three-dimensional reconstruction at 23-A resolution. The cryoEM structure, sequence comparison, and protein fold recognition analysis indicate that the coat protein of MGNNV has two domains resembling those of tomato bushy stunt virus and Norwalk virus, rather than the expected single-domain coat protein of insect nodaviruses. The analysis implies that residues 83 to 216 fold as a beta-sandwich which forms the inner shell of the T=3 capsid and residues 217 to 308 form the trimeric surface protrusions observed in the cryoEM map. The structural similarities between fish nodaviruses and members of the tombusvirus and calicivirus groups provide significant new data for understanding the evolution of the nodavirus family.     

A machine learning approach for the identification of protein secondary structure elements from electron cryo-microscopy density maps

The accuracy of the secondary structure element (SSE) identification from volumetric protein density maps is critical for de-novo backbone structure derivation in electron cryo-microscopy (cryoEM). It is still challenging to detect the SSE automatically and accurately from the density maps at medium resolutions (～5-10 ?). We present a machine learning approach, SSELearner, to automatically identify helices and β-sheets by using the knowledge from existing volumetric maps in the Electron Microscopy Data Bank. We tested our approach using 10 simulated density maps. The averaged specificity and sensitivity for the helix detection are 94.9% and 95.8%, respectively, and those for the β-sheet detection are 86.7% and 96.4%, respectively. We have developed a secondary structure annotator, SSID, to predict the helices and β-strands from the backbone Cα trace. With the help of SSID, we tested our SSELearner using 13 experimentally derived cryo-EM density maps. The machine learning approach shows the specificity and sensitivity of 91.8% and 74.5%, respectively, for the helix detection and 85.2% and 86.5% respectively for the β-sheet detection in cryoEM maps of Electron Microscopy Data Bank. The reduced detection accuracy reveals the challenges in SSE detection when the cryoEM maps are used instead of the simulated maps. Our results suggest that it is effective to use one cryoEM map for learning to detect the SSE in another cryoEM map of similar quality.     

Fabs enable single particle cryoEM studies of small proteins

In spite of its recent achievements, the technique of single particle electron cryomicroscopy (cryoEM) has not been widely used to study proteins smaller than 100 kDa, although it is a highly desirable application of this technique. One fundamental limitation is that images of small proteins embedded in vitreous ice do not contain adequate features for accurate image alignment. We describe a general strategy to overcome this limitation by selecting a fragment antigen binding (Fab) to form a stable and rigid complex with a target protein, thus providing a defined feature for accurate image alignment. Using this approach, we determined a three-dimensional structure of an ～65 kDa protein by single particle cryoEM. Because Fabs can be readily generated against a wide range of proteins by phage display, this approach is generally applicable to study many small proteins by single particle cryoEM.     

Do's and Don'ts of Cryo-electron Microscopy: A Primer on Sample Preparation and High Quality Data Collection for Macromolecular 3D Reconstruction

Cryo-electron microscopy (cryoEM) entails flash-freezing a thin layer of sample on a support, and then visualizing the sample in its frozen hydrated state by transmission electron microscopy (TEM). This can be achieved with very low quantity of protein and in the buffer of choice, without the use of any stain, which is very useful to determine structure-function correlations of macromolecules. When combined with single-particle image processing, the technique has found widespread usefulness for 3D structural determination of purified macromolecules. The protocol presented here explains how to perform cryoEM and examines the causes of most commonly encountered problems for rational troubleshooting; following all these steps should lead to acquisition of high quality cryoEM images. The technique requires access to the electron microscope instrument and to a vitrification device. Knowledge of the 3D reconstruction concepts and software is also needed for computerized image processing. Importantly, high quality results depend on finding the right purification conditions leading to a uniform population of structurally intact macromolecules. The ability of cryoEM to visualize macromolecules combined with the versatility of single particle image processing has proven very successful for structural determination of large proteins and macromolecular machines in their near-native state, identification of their multiple components by 3D difference mapping, and creation of pseudo-atomic structures by docking of x-ray structures. The relentless development of cryoEM instrumentation and image processing techniques for the last 30 years has resulted in the possibility to generate de novo 3D reconstructions at atomic resolution level.     

Symmetry-adapted spherical harmonics method for high-resolution 3D single-particle reconstructions

Three-dimensional (3D) reconstruction is the last and an essential step toward high-resolution structural determination in single-particle cryo-electron microscopy (cryoEM). We have implemented a new algorithm for reconstructing 3D structures of macromolecular complexes with icosahedral symmetry from cryoEM images. Icosahedral symmetry-adapted functions (ISAFs) are used to interpolate structural factors in the reciprocal space to generate a 3D reconstruction in spherical coordinates. In our implementation, we introduced a recursive method for deriving higher order ISAFs from three lower order seed functions. We demonstrate improvements of our new method in both the noise suppression and the effective resolution in 3D reconstruction over the commonly used Fourier-Bessel synthesis method introduced by Crowther et al. three decades ago. Our 3D reconstruction method can be extended to macromolecular complexes with other symmetry types and is thus likely to impact future high-resolution cryoEM single-particle reconstruction efforts in general.     

Automated segmentation of molecular subunits in electron cryomicroscopy density maps

Electron cryomicroscopy (cryoEM) is capable of imaging large macromolecular machines composed of multiple components. However, it is currently only possible to achieve moderate resolution at which it may be possible to computationally extract the individual components in the machine. In this work, we present application details of an automated method for detecting and segmenting the components of a large machine in an experimentally determined density map. This method is applicable to object with and without symmetry and takes advantage of global and local symmetry axes if present. We have applied this segmentation algorithm to several cryoEM data sets already deposited in EMDB with various complexities, symmetries and resolutions and validated the results using manually segmented density and available structures of the components in the PDB. As such, automated segmentation could become a useful tool for the analysis of the ever-increasing number of structures of macromolecular machines derived from cryoEM.     

Interaction of coxsackievirus A21 with its cellular receptor, ICAM-1

Coxsackievirus A21 (CAV21), like human rhinoviruses (HRVs), is a causative agent of the common cold. It uses the same cellular receptor, intercellular adhesion molecule 1 (ICAM-1), as does the major group of HRVs; unlike HRVs, however, it is stable at acid pH. The cryoelectron microscopy (cryoEM) image reconstruction of CAV21 is consistent with the highly homologous crystal structure of poliovirus 1; like other enteroviruses and HRVs, CAV21 has a canyon-like depression around each of the 12 fivefold vertices. A cryoEM reconstruction of CAV21 complexed with ICAM-1 shows all five domains of the extracellular component of ICAM-1. The known atomic structure of the ICAM-1 amino-terminal domains D1 and D2 has been fitted into the cryoEM density of the complex. The site of ICAM-1 binding within the canyon of CAV21 overlaps the site of receptor recognition utilized by rhinoviruses and polioviruses. Interactions within this common region may be essential for triggering viral destabilization after attachment to susceptible cells.     

Ab initio modeling of the herpesvirus VP26 core domain assessed by CryoEM density

Efforts in structural biology have targeted the systematic determination of all protein structures through experimental determination or modeling. In recent years, 3-D electron cryomicroscopy (cryoEM) has assumed an increasingly important role in determining the structures of these large macromolecular assemblies to intermediate resolutions (6–10 Å). While these structures provide a snapshot of the assembly and its components in well-defined functional states, the resolution limits the ability to build accurate structural models. In contrast, sequence-based modeling techniques are capable of producing relatively robust structural models for isolated proteins or domains. In this work, we developed and applied a hybrid modeling approach, utilizing cryoEM density and ab initio modeling to produce a structural model for the core domain of a herpesvirus structural protein, VP26. Specifically, this method, first tested on simulated data, utilizes the cryoEM density map as a geometrical constraint in identifying the most native-like models from a gallery of models generated by ab initio modeling. The resulting model for the core domain of VP26, based on the 8.5-Å resolution herpes simplex virus type 1 (HSV-1) capsid cryoEM structure and mutational data, exhibited a novel fold. Additionally, the core domain of VP26 appeared to have a complementary interface to the known upper-domain structure of VP5, its cognate binding partner. While this new model provides for a better understanding of the assembly and interactions of VP26 in HSV-1, the approach itself may have broader applications in modeling the components of large macromolecular assemblies.     

Limiting factors in atomic resolution cryo electron microscopy: no simple tricks

To bring cryo electron microscopy (cryoEM) of large biological complexes to atomic resolution, several factors--in both cryoEM image acquisition and 3D reconstruction--that may be neglected at low resolution become significantly limiting. Here we present thorough analyses of four limiting factors: (a) electron-beam tilt, (b) inaccurate determination of defocus values, (c) focus gradient through particles, and (d) particularly for large particles, dynamic (multiple) scattering of electrons. We also propose strategies to cope with these factors: (a) the divergence and direction tilt components of electron-beam tilt could be reduced by maintaining parallel illumination and by using a coma-free alignment procedure, respectively. Moreover, the effect of all beam tilt components, including spiral tilt, could be eliminated by use of a spherical aberration corrector. (b) More accurate measurement of defocus value could be obtained by imaging areas adjacent to the target area at high electron dose and by measuring the image shift induced by tilting the electron beam. (c) Each known Fourier coefficient in the Fourier transform of a cryoEM image is the sum of two Fourier coefficients of the 3D structure, one on each of two curved 'characteristic surfaces' in 3D Fourier space. We describe a simple model-based iterative method that could recover these two Fourier coefficients on the two characteristic surfaces. (d) The effect of dynamic scattering could be corrected by deconvolution of a transfer function. These analyses and our proposed strategies offer useful guidance for future experimental designs targeting atomic resolution cryoEM reconstruction.     

Visualization of the externalized VP2 N termini of infectious human parvovirus B19

The structures of infectious human parvovirus B19 and empty wild-type particles were determined by cryoelectron microscopy (cryoEM) to 7.5-Å and 11.3-Å resolution, respectively, assuming icosahedral symmetry. Both of these, DNA filled and empty, wild-type particles contain a few copies of the minor capsid protein VP1. Comparison of wild-type B19 with the crystal structure and cryoEM reconstruction of recombinant B19 particles consisting of only the major capsid protein VP2 showed structural differences in the vicinity of the icosahedral fivefold axes. Although the unique N-terminal region of VP1 could not be visualized in the icosahedrally averaged maps, the N terminus of VP2 was shown to be exposed on the viral surface adjacent to the fivefold β-cylinder. The conserved glycine-rich region is positioned between two neighboring, fivefold-symmetrically related VP subunits and not in the fivefold channel as observed for other parvoviruses.     

Refinement of protein structures by iterative comparative modeling and CryoEM density fitting

We developed a method for structure characterization of assembly components by iterative comparative protein structure modeling and fitting into cryo-electron microscopy (cryoEM) density maps. Specifically, we calculate a comparative model of a given component by considering many alternative alignments between the target sequence and a related template structure while optimizing the fit of a model into the corresponding density map. The method relies on the previously developed Moulder protocol that iterates over alignment, model building, and model assessment. The protocol was benchmarked using 20 varied target-template pairs of known structures with less than 30% sequence identity and corresponding simulated density maps at resolutions from 5A to 25A. Relative to the models based on the best existing sequence profile alignment methods, the percentage of C(alpha) atoms that are within 5A of the corresponding C(alpha) atoms in the superposed native structure increases on average from 52% to 66%, which is half-way between the starting models and the models from the best possible alignments (82%). The test also reveals that despite the improvements in the accuracy of the fitness function, this function is still the bottleneck in reducing the remaining errors. To demonstrate the usefulness of the protocol, we applied it to the upper domain of the P8 capsid protein of rice dwarf virus that has been studied by cryoEM at 6.8A. The C(alpha) root-mean-square deviation of the model based on the remotely related template, bluetongue virus VP7, improved from 8.7A to 6.0A, while the best possible model has a C(alpha) RMSD value of 5.3A. Moreover, the resulting model fits better into the cryoEM density map than the initial template structure. The method is being implemented in our program MODELLER for protein structure modeling by satisfaction of spatial restraints and will be applicable to the rapidly increasing number of cryoEM density maps of macromolecular assemblies.     

Capturing pore-forming intermediates of MACPF and binary toxin assemblies by cryoEM

Deployed by both pathogenic bacteria and host immune systems, pore-forming proteins rupture target membranes and can serve as conduits for effector proteins. Understanding how these proteins work relies on capturing assembly intermediates. Advances in cryoEM allowing in silico purification of heterogeneous assemblies has led to new insights into two main classes of pore-forming proteins: membrane attack complex perforin (MACPF) proteins and binary toxins. The structure of an immune activation complex, sMAC, shows how pores form by sequential templating and insertion of β-hairpins. CryoEM structures of bacterial binary toxins present a series of transitions along the pore formation pathway and reveal a general mechanism of effector protein translocation. Future developments in time-resolved cryoEM could capture and place short-lived states along the trajectory of pore-formation.     

BCL::EM-Fit: rigid body fitting of atomic structures into density maps using geometric hashing and real space refinement

Cryo-electron microscopy (cryoEM) can visualize large macromolecular assemblies at resolutions often below 10? and recently as good as 3.8-4.5 ?. These density maps provide important insights into the biological functioning of molecular machineries such as viruses or the ribosome, in particular if atomic-resolution crystal structures or models of individual components of the assembly can be placed into the density map. The present work introduces a novel algorithm termed BCL::EM-Fit that accurately fits atomic-detail structural models into medium resolution density maps. In an initial step, a "geometric hashing" algorithm provides a short list of likely placements. In a follow up Monte Carlo/Metropolis refinement step, the initial placements are optimized by their cross correlation coefficient. The resolution of density maps for a reliable fit was determined to be 10 ? or better using tests with simulated density maps. The algorithm was applied to fitting of capsid proteins into an experimental cryoEM density map of human adenovirus at a resolution of 6.8 and 9.0 ?, and fitting of the GroEL protein at 5.4 ?. In the process, the handedness of the cryoEM density map was unambiguously identified. The BCL::EM-Fit algorithm offers an alternative to the established Fourier/Real space fitting programs. BCL::EM-Fit is free for academic use and available from a web server or as downloadable binary file at http://www.meilerlab.org.     

Resolution and Probabilistic Models of Components in CryoEM Maps of Mature P22 Bacteriophage

CryoEM continues to produce density maps of larger and more complex assemblies with multiple protein components of mixed symmetries. Resolution is not always uniform throughout a cryoEM map, and it can be useful to estimate the resolution in specific molecular components of a large assembly. In this study, we present procedures to 1) estimate the resolution in subcomponents by gold-standard Fourier shell correlation (FSC); 2) validate modeling procedures, particularly at medium resolutions, which can include loop modeling and flexible fitting; and 3) build probabilistic models that combine high-accuracy priors (such as crystallographic structures) with medium-resolution cryoEM densities. As an example, we apply these methods to new cryoEM maps of the mature bacteriophage P22, reconstructed without imposing icosahedral symmetry. Resolution estimates based on gold-standard FSC show the highest resolution in the coat region (7.6 Å), whereas other components are at slightly lower resolutions: portal (9.2 Å), hub (8.5 Å), tailspike (10.9 Å), and needle (10.5 Å). These differences are indicative of inherent structural heterogeneity and/or reconstruction accuracy in different subcomponents of the map. Probabilistic models for these subcomponents provide new insights, to our knowledge, and structural information when taking into account uncertainty given the limitations of the observed density.     

Visualization of alpha-helices in a 6-angstrom resolution cryoelectron microscopy structure of adenovirus allows refinement of capsid protein assignments

The structure of adenovirus was determined to a resolution of 6 A by cryoelectron microscopy (cryoEM) single-particle image reconstruction. Docking of the hexon and penton base crystal structures into the cryoEM density established that alpha-helices of 10 or more residues are resolved as rods. A difference map was calculated by subtracting a pseudoatomic capsid from the cryoEM reconstruction. The resulting density was analyzed in terms of observed alpha-helices and secondary structure predictions for the additional capsid proteins that currently lack atomic resolution structures (proteins IIIa, VI, VIII, and IX). Protein IIIa, which is predicted to be highly alpha-helical, is assigned to a cluster of helices observed below the penton base on the inner capsid surface. Protein VI is present in approximately 1.5 copies per hexon trimer and is predicted to have two long alpha-helices, one of which appears to lie inside the hexon cavity. Protein VIII is cleaved by the adenovirus protease into two fragments of 7.6 and 12.1 kDa, and the larger fragment is predicted to have one long alpha-helix, in agreement with the observed density for protein VIII on the inner capsid surface. Protein IX is predicted to have one long alpha-helix, which also has a strongly indicated propensity for coiled-coil formation. A region of density near the facet edge is now resolved as a four-helix bundle and is assigned to four copies of the C-terminal alpha-helix from protein IX.