Electron cryomicroscopy of biological machines at subnanometer resolution

Advances in electron cryomicroscopy (cryo-EM) have made possible the structural determination of large biological machines in the resolution range of 6-9 angstroms. Rice dwarf virus and the acrosomal bundle represent two distinct types of machines amenable to cryo-EM investigations at subnanometer resolutions. However, calculating the density map is only the first step, and much analysis remains to extract structural insights and the mechanism of action in these machines. This paper will review the computational and visualization methodologies necessary for analysis (structure mining) of the computed cryo-EM maps of these machines. These steps include component segmentation, averaging based on local symmetry among components, density connectivity trace, incorporation of bioinformatics analysis, and fitting of high-resolution component data, if available. The consequences of these analyses can not only identify accurately some of the secondary structure elements of the molecular components in machines but also suggest structural mechanisms related to their biological functions.     

Rapid increase of near atomic resolution virus capsid structures determined by cryo-electron microscopy

The recent technological advances in electron microscopes, detectors, as well as image processing and reconstruction software have brought single particle cryo-electron microscopy (cryo-EM) into prominence for determining structures of bio-molecules at near atomic resolution. This has been particularly true for virus capsids, ribosomes, and other large assemblies, which have been the ideal specimens for structural studies by cryo-EM approaches. An analysis of time series metadata of virus structures on the methods of structure determination, resolution of the structures, and size of the virus particles revealed a rapid increase in the virus structures determined by cryo-EM at near atomic resolution since 2010. In addition, the data highlight the median resolution (～3.0?Å) and size (～310.0?Å in diameter) of the virus particles determined by X-ray crystallography while no such limits exist for cryo-EM structures, which have a median diameter of 508?Å. Notably, cryo-EM virus structures in the last four years have a median resolution of 3.9?Å. Taken together with minimal sample requirements, not needing diffraction quality crystals, and being able to achieve similar resolutions of the crystal structures makes cryo-EM the method of choice for current and future virus capsid structure determinations.     

CryoRes: Local Resolution Estimation of Cryo-EM Density Maps by Deep Learning

Recent progress in cryo-EM research has ignited a revolution in biological macromolecule structure determination. Resolution is an essential parameter for quality assessment of a cryo-EM density map, and it is known that resolution varies in different regions of a map. Currently available methods for local resolution estimation require manual adjustment of parameters and in some cases necessitate acquisition or de novo generation of so-called “half maps”. Here, we developed CryoRes, a deep-learning algorithm to estimate local resolution directly from a single final cryo-EM density map, specifically by learning resolution-aware patterns of density map voxels through supervised training on a large dataset comprising 1,174 experimental cryo-EM density maps. CryoRes significantly outperforms all of the state-of-the-art competing resolution estimation methods, achieving an average RMSE of 2.26 Å for local resolution estimation relative to the currently most reliable FSC-based method blocres, yet requiring only the single final map as input. Further, CryoRes is able to generate a molecular mask for each map, with accuracy 12.12% higher than the masks generated by ResMap. CryoRes is ultra-fast, fully automatic, parameter-free, applicable to cryo-EM subtomogram data, and freely available at https://cryores.zhanglab.net.     

Cryo-EM of macromolecular assemblies at near-atomic resolution

With single-particle electron cryomicroscopy (cryo-EM), it is possible to visualize large, macromolecular assemblies in near-native states. Although subnanometer resolutions have been routinely achieved for many specimens, state of the art cryo-EM has pushed to near-atomic (3.3-4.6 ?) resolutions. At these resolutions, it is now possible to construct reliable atomic models directly from the cryo-EM density map. In this study, we describe our recently developed protocols for performing the three-dimensional reconstruction and modeling of Mm-cpn, a group II chaperonin, determined to 4.3 ? resolution. This protocol, utilizing the software tools EMAN, Gorgon and Coot, can be adapted for use with nearly all specimens imaged with cryo-EM that target beyond 5 ? resolution. Additionally, the feature recognition and computational modeling tools can be applied to any near-atomic resolution density maps, including those from X-ray crystallography.     

Symmetry-restrained flexible fitting for symmetric EM maps

Many large biological macromolecules have inherent structural symmetry, being composed of a few distinct subunits, repeated in a symmetric array. These complexes are often not amenable to traditional high-resolution structural determination methods, but can be imaged in functionally relevant states using cryo-electron microscopy (cryo-EM). A number of methods for fitting atomic-scale structures into cryo-EM maps have been developed, including the molecular dynamics flexible fitting (MDFF) method. However, quality and resolution of the cryo-EM map are the major determinants of a method's success. In order to incorporate knowledge of structural symmetry into the fitting procedure, we developed the symmetry-restrained MDFF method. The new method adds to the cryo-EM map-derived potential further restraints on the allowed conformations of a complex during fitting, thereby improving the quality of the resultant structure. The benefit of using symmetry-based restraints during fitting, particularly for medium to low-resolution data, is demonstrated for three different systems.     

Seeing GroEL at 6 A resolution by single particle electron cryomicroscopy

We present a reconstruction of native GroEL by electron cryomicroscopy (cryo-EM) and single particle analysis at 6 A resolution. alpha helices are clearly visible and beta sheet density is also visible at this resolution. While the overall conformation of this structure is quite consistent with the published X-ray data, a measurable shift in the positions of three alpha helices in the intermediate domain is observed, not consistent with any of the 7 monomeric structures in the Protein Data Bank model (1OEL). In addition, there is evidence for slight rearrangement or flexibility in parts of the apical domain. The 6 A resolution cryo-EM GroEL structure clearly demonstrates the veracity and expanding scope of cryo-EM and the single particle reconstruction technique for macromolecular machines.     

Applications of the molecular dynamics flexible fitting method

In recent years, cryo-electron microscopy (cryo-EM) has established itself as a key method in structural biology, permitting the structural characterization of large biomolecular complexes in various functional states. The data obtained through single-particle cryo-EM has recently seen a leap in resolution thanks to landmark advances in experimental and computational techniques, resulting in sub-nanometer resolution structures being obtained routinely. The remaining gap between these data and revealing the mechanisms of molecular function can be closed through hybrid modeling tools that incorporate known atomic structures into the cryo-EM data. One such tool, molecular dynamics flexible fitting (MDFF), uses molecular dynamics simulations to combine structures from X-ray crystallography with cryo-EM density maps to derive atomic models of large biomolecular complexes. The structures furnished by MDFF can be used subsequently in computational investigations aimed at revealing the dynamics of the complexes under study. In the present work, recent applications of MDFF are presented, including the interpretation of cryo-EM data of the ribosome at different stages of translation and the structure of a membrane-curvature-inducing photosynthetic complex.     

ATTRACT-EM: A New Method for the Computational Assembly of Large Molecular Machines Using Cryo-EM Maps

Many of the most important functions in the cell are carried out by proteins organized in large molecular machines. Cryo-electron microscopy (cryo-EM) is increasingly being used to obtain low resolution density maps of these large assemblies. A new method, ATTRACT-EM, for the computational assembly of molecular assemblies from their components has been developed. Based on concepts from the protein-protein docking field, it utilizes cryo-EM density maps to assemble molecular subunits at near atomic detail, starting from millions of initial subunit configurations. The search efficiency was further enhanced by recombining partial solutions, the inclusion of symmetry information, and refinement using a molecular force field. The approach was tested on the GroES-GroEL system, using an experimental cryo-EM map at 23.5 Å resolution, and on several smaller complexes. Inclusion of experimental information on the symmetry of the systems and the application of a new gradient vector matching algorithm allowed the efficient identification of docked assemblies in close agreement with experiment. Application to the GroES-GroEL complex resulted in a top ranked model with a deviation of 4.6 Å (and a 2.8 Å model within the top 10) from the GroES-GroEL crystal structure, a significant improvement over existing methods.     

Structure of HIV-1 Capsid Assemblies by Cryo-electron Microscopy and Iterative Helical Real-space Reconstruction

Cryo-electron microscopy (cryo-EM), combined with image processing, is an increasingly powerful tool for structure determination of macromolecular protein complexes and assemblies. In fact, single particle electron microscopy¹ and two-dimensional (2D) electron crystallography² have become relatively routine methodologies and a large number of structures have been solved using these methods. At the same time, image processing and three-dimensional (3D) reconstruction of helical objects has rapidly developed, especially, the iterative helical real-space reconstruction (IHRSR) method³, which uses single particle analysis tools in conjunction with helical symmetry. Many biological entities function in filamentous or helical forms, including actin filaments⁴, microtubules⁵, amyloid fibers⁶, tobacco mosaic viruses⁷, and bacteria flagella⁸, and, because a 3D density map of a helical entity can be attained from a single projection image, compared to the many images required for 3D reconstruction of a non-helical object, with the IHRSR method, structural analysis of such flexible and disordered helical assemblies is now attainable.In this video article, we provide detailed protocols for obtaining a 3D density map of a helical protein assembly (HIV-1 capsid⁹ is our example), including protocols for cryo-EM specimen preparation, low dose data collection by cryo-EM, indexing of helical diffraction patterns, and image processing and 3D reconstruction using IHRSR. Compared to other techniques, cryo-EM offers optimal specimen preservation under near native conditions. Samples are embedded in a thin layer of vitreous ice, by rapid freezing, and imaged in electron microscopes at liquid nitrogen temperature, under low dose conditions to minimize the radiation damage. Sample images are obtained under near native conditions at the expense of low signal and low contrast in the recorded micrographs. Fortunately, the process of helical reconstruction has largely been automated, with the exception of indexing the helical diffraction pattern. Here, we describe an approach to index helical structure and determine helical symmetries (helical parameters) from digitized micrographs, an essential step for 3D helical reconstruction. Briefly, we obtain an initial 3D density map by applying the IHRSR method. This initial map is then iteratively refined by introducing constraints for the alignment parameters of each segment, thus controlling their degrees of freedom. Further improvement is achieved by correcting for the contrast transfer function (CTF) of the electron microscope (amplitude and phase correction) and by optimizing the helical symmetry of the assembly.     

3.9 ? structure of the nucleosome core particle determined by phase-plate cryo-EM

The Volta phase plate is a recently developed electron cryo-microscopy (cryo-EM) device that enables contrast enhancement of biological samples. Here we have evaluated the potential of combining phase-plate imaging and single particle analysis to determine the structure of a small protein–DNA complex. To test the method, we made use of a 200 kDa Nucleosome Core Particle (NCP) reconstituted with 601 DNA for which a high-resolution X-ray crystal structure is known. We find that the phase plate provides a significant contrast enhancement that permits individual NCPs and DNA to be clearly identified in amorphous ice. The refined structure from 26,060 particles has an overall resolution of 3.9 Å and the density map exhibits structural features consistent with the estimated resolution, including clear density for amino acid side chains and DNA features such as the phosphate backbone. Our results demonstrate that phase-plate cryo-EM promises to become an important method to determine novel near-atomic resolution structures of small and challenging samples, such as nucleosomes in complex with nucleosome-binding factors.     

Near-atomic resolution reconstructions of icosahedral viruses from electron cryo-microscopy

Nine different near-atomic resolution structures of icosahedral viruses, determined by electron cryo-microscopy and published between early 2008 and late 2010, fulfil predictions made 15 years ago that single-particle cryo-EM techniques could visualize molecular detail at 3-4? resolution. This review summarizes technical developments, both in instrumentation and in computation, that have led to the new structures, which advance our understanding of virus assembly and cell entry.     

Structure of cyclodextrin glycosyltransferase refined at 2.0 A resolution.

The previously reported structural model of cyclodextrin glycosyltransferase (EC 2.4.1.19) from Bacillus circulans has been improved. For this purpose the known sequence was built into an electron density map established by multiple isomorphous replacement and subsequent solvent-flattening at 2.5 A resolution. The resulting model was refined at 2.0 A resolution using a simulated annealing refinement method. Based on 70,171 independent reflections in the range 7.0 to 2.0 A resolution, a final R-factor of 17.6% was obtained with a model obeying standard geometry within 0.013 A in bond lengths and 2.7 degrees in bond angles. The final model consists of all 684 amino acid residues, two calcium ions and 588 solvent molecules.     

The resolution dependence of optimal exposures in liquid nitrogen temperature electron cryomicroscopy of catalase crystals

Electron beam damage is the fundamental limit to resolution in electron cryomicroscopy (cryo-EM) of frozen, hydrated specimens. Radiation damage increases with the number of electrons used to obtain an image and affects information at higher spatial frequencies before low-resolution information. For the experimentalist, a balance exists between electron exposures sufficient to obtain a useful signal-to-noise ratio (SNR) in images and exposures that limit the damage to structural features. In single particle cryo-EM this balance is particularly delicate: low-resolution features must be imaged with a sufficient SNR to allow image alignment so that high-resolution features recorded below the noise level can be recovered by averaging independent images. By measuring the fading of Fourier components from images obtained at 200 kV of thin crystals of catalase embedded in ice, we have determined the electron exposures that will maximize the SNR at resolutions between 86 and 2.9 Å. These data allow for a rational choice of exposure for single particle cryo-EM. For example, for 20 Å resolution, the SNR is maximized at ～20 e^?/Ų, whereas for 3 Å resolution, it is maximized at ～10 e^?/Ų. We illustrate the effects of exposure in single particle cryo-EM with data collected at ～12–15 and ～24–30 e^?/Ų.     

An atomic model for the human septin hexamer by cryo-EM

In order to form functional filaments, human septins must assemble into hetero-oligomeric rod-like particles which polymerize end-to-end. The rules governing the assembly of these particles and the subsequent filaments are incompletely understood. Although crystallographic approaches have been successful in studying the separate components of the system, there has been difficulty in obtaining high resolution structures of the full particle. Here we report a first cryo-EM structure for a hexameric rod composed of human septins 2, 6 and 7 with a global resolution of ~3.6 Å and a local resolution of between ~3.0 Å and ~5.0 Å. By fitting the previously determined high-resolution crystal structures of the component subunits into the cryo-EM map, we are able to provide an essentially complete model for the particle. This exposes SEPT2 NC-interfaces at the termini of the hexamer and leaves internal cavities between the SEPT6-SEPT7 pairs. The floor of the cavity is formed by the two α₀ helices including their polybasic regions. These are locked into place between the two subunits by interactions made with the α₅ and α₆ helices of the neighbouring monomer together with its polyacidic region. The cavity may serve to provide space allowing the subunits to move with respect to one another. The elongated particle shows a tendency to bend at its centre where two copies of SEPT7 form a homodimeric G-interface. Such bending is almost certainly related to the ability of septin filaments to recognize and even induce membrane curvature.     

Cryo-EM density maps adjustment for subtraction,consensus and sharpening

Electron cryomicroscopy (cryo-EM) has emerged as a powerful structural biology instrument to solve near-atomic three-dimensional structures. Despite the fast growth in the number of density maps generated from cryo-EM data, comparison tools among these reconstructions are still lacking. Current proposals to compare cryo-EM data derived volumes perform map subtraction based on adjustment of each volume grey level to the same scale. We present here a more sophisticated way of adjusting the volumes before comparing, which implies adjustment of grey level scale and spectrum energy, but keeping phases intact inside a mask and imposing the results to be strictly positive. The adjustment that we propose leaves the volumes in the same numeric frame, allowing to perform operations among the adjusted volumes in a more reliable way. This adjustment can be a preliminary step for several applications such as comparison through subtraction, map sharpening, or combination of volumes through a consensus that selects the best resolved parts of each input map. Our development might also be used as a sharpening method using an atomic model as a reference. We illustrate the applicability of this algorithm with the reconstructions derived of several experimental examples. This algorithm is implemented in Xmipp software package and its applications are user-friendly accessible through the cryo-EM image processing framework Scipion.     

Flexible fitting of high-resolution x-ray structures into cryoelectron microscopy maps using biased molecular dynamics simulations

A methodology for flexible fitting of all-atom high-resolution structures into low-resolution cryoelectron microscopy (cryo-EM) maps is presented. Flexibility of the modeled structure is simulated by classical molecular dynamics and an additional effective potential is introduced to enhance the fitting process. The additional potential is proportional to the correlation coefficient between the experimental cryo-EM map and a synthetic map generated for an all-atom structure being fitted to the map. The additional forces are calculated as a gradient of the correlation coefficient. During the molecular dynamics simulations under the additional forces, the molecule undergoes a conformational transition that maximizes the correlation coefficient, which results in a high-accuracy fit of all-atom structure into a cryo-EM map. Using five test proteins that exhibit structural rearrangement during their biological activity, we demonstrate performance of our method. We also test our method on the experimental cryo-EM of elongation factor G and show that the model obtained is comparable to previous studies. In addition, we show that overfitting can be avoided by assessing the quality of the fitted model in terms of correlation coefficient and secondary structure preservation.     

Macromolecular electron microscopy in the era of structural genomics

Macromolecular machines carry out many cellular functions. Cryo-electron microscopy (cryo-EM) is emerging as a powerful method for studying the structure, assembly and dynamics of such macromolecules, and their interactions with substrates. With resolutions still improving, ‘single-particle’ analyses are already depicting secondary structure. Moreover, cryo-EM can be combined in several ways with X-ray diffraction to enhance the resolution of cryo-EM and the applicability of crystallography. Electron tomography holds promise for visualizing machines at work inside cells.     

Using a SMALP platform to determine a sub-nm single particle cryo-EM membrane protein structure

The field of membrane protein structural biology has been revolutionized over the last few years with a number of high profile structures being solved using cryo-EM including Piezo, Ryanodine receptor, TRPV1 and the Glutamate receptor. Further developments in the EM field hold the promise of even greater progress in terms of greater resolution, which for membrane proteins is still typically within the 4–7 Å range. One advantage of a cryo-EM approach is the ability to study membrane proteins in more “native” like environments for example proteoliposomes, amphipols and nanodiscs. Recently, styrene maleic acid co-polymers (SMA) have been used to extract membrane proteins surrounded by native lipids (SMALPs) maintaining a more natural environment. We report here the structure of the Escherichia coli multidrug efflux transporter AcrB in a SMALP scaffold to sub-nm resolution, with the resulting map being consistent with high resolution crystal structures and other EM derived maps. However, both the C-terminal helix (TM12) and TM7 are poorly defined in the map. These helices are at the exterior of the helical bundle and form the greater interaction with the native lipids and SMA polymer and may represent a more dynamic region of the protein. This work shows the promise of using an SMA approach for single particle cryo-EM studies to provide sub-nm structures.     

All-atom homology model of the Escherichia coli 30S ribosomal subunit

Understanding the structural basis of ribosomal function requires close comparison between biochemical and structural data. Although a large amount of biochemical data are available for the Escherichia coli ribosome, the structure has not been solved to atomic resolution. Using a new RNA homology procedure, we have modeled the all-atom structure of the E. coli 30S ribosomal subunit. We find that the tertiary structure of the ribosome core, including the A-, P- and E-sites, is highly conserved. The hypervariable regions in our structure, which differ from the structure of the 30S ribosomal subunit from Thermus thermophilus, are consistent with the cryo-EM map of the E. coli ribosome.     

Pf1 Inovirus. Electron density distribution calculated by a maximum entropy algorithm from native fibre diffraction data to 3 A resolution and single isomorphous replacement data to 5 A resolution

We have calculated the electron density distribution of the Pf1 strain of filamentous bacteriophage by a maximum entropy method. In the calculation we included native X-ray fibre diffraction data extending to 3 A resolution in the meridional direction on 60 layerlines that are resolved to 4 A in the equatorial direction, and lower resolution data from a single isomorphous derivative iodinated on the Tyr25 residue. The electron density map indicates that the 46-residue protein subunit is a single, gently curved stretch of alpha-helix with its axis at an angle of about 20 degrees to the axis of the virion. The alpha-helix subunit curves around the virion axis by about 1/6 turn, and decreases from about 27 A radius to about 13 A radius in the virion as the amino acid sequence of the subunit runs from the N terminus to the C terminus. Nearest-neighbour alpha-helical subunits are about 10 A apart along their length, and the axis of each subunit makes an unexpected negative angle with its nearest neighbours in the virion. To confirm the validity of the maximum entropy calculation, we have varied the constraints on the calculation. All variations result in either a map that is close to the original map or a map that cannot be interpreted in terms of secondary structure: we find only one map that makes structural sense.