Atomic Resolution Cryo-EM Structure of β-Galactosidase

Could Egg White Lysozyme be Solved by Single Particle Cryo-EM?

The combination of high-end cryogenic transmission electron microscopes (cryo-EM), direct electron detectors, and advanced image algorithms allows researchers to obtain the 3D structures of much smaller macromolecules than years ago. However, there are still major challenges for the single-particle cryo-EM method to achieve routine structure determinations for macromolecules much smaller than 100 kDa, which are the majority of all plant and animal proteins. These challenges include sample characteristics such as sample heterogeneity, beam damage, ice layer thickness, stability, and quality, as well as hardware limitations such as detector performance, beam, and phase plate quality. Here, single particle data sets were simulated for samples that were ideal in terms of homogeneity, distribution, and stability, but with realistic parameters for ice layer, dose, detector performance, and beam characteristics. Reference data were calculated for human apo-ferritin using identical parameters reported for an experimental data set downloaded from EMPIAR. Processing of the simulated data set resulted in a value of 1.86 Å from 20 214 particles, similar to a 2 Å density map obtained from 29 224 particles selected from real micrographs. Simulated data sets were then generated for a 14 kDa protein, hen egg white lysozyme (HEWL), with and without an ideal phase plate (PP). Whereas we could not obtain a high-resolution 3D reconstruction of HEWL for the data set without PP, the one with PP resulted in a 2.78 Å resolution density map from 225 751 particles. Our simulator and simulations could help in pushing the size limits of cryo-EM.

GemSpot: A Pipeline for Robust Modeling of Ligands into Cryo-EM Maps

Producing an accurate atomic model of biomolecule-ligand interactions from maps generated by cryoelectron microscopy (cryo-EM) often presents challenges inherent to the methodology and the dynamic nature of ligand binding. Here, we present GemSpot, an automated pipeline of computational chemistry methods that take into account EM map potentials, quantum mechanics energy calculations, and water molecule site prediction to generate candidate poses and provide a measure of the degree of confidence. The pipeline is validated through several published cryo-EM structures of complexes in different resolution ranges and various types of ligands. In all cases, at least one identified pose produced both excellent interactions with the target and agreement with the map. GemSpot will be valuable for the robust identification of ligand poses and drug discovery efforts through cryo-EM.

Confidence maps: statistical inference of cryo-EM maps

Confidence maps provide complementary information for interpreting cryo-EM densities as they indicate statistical significance with respect to background noise. They can be thresholded by specifying the expected false-discovery rate (FDR), and the displayed volume shows the parts of the map that have the corresponding level of significance. Here, the basic statistical concepts of confidence maps are reviewed and practical guidance is provided for their interpretation and usage inside the CCP-EM suite. Limitations of the approach are discussed and extensions towards other error criteria such as the family-wise error rate are presented. The observed map features can be rendered at a common isosurface threshold, which is particularly beneficial for the interpretation of weak and noisy densities. In the current article, a practical guide is provided to the recommended usage of confidence maps.

Outlier Profiles of Atomic Structures Derived from X-ray Crystallography and from Cryo-Electron Microscopy

Background: As more protein atomic structures are determined from cryo-electron microscopy (cryo-EM) density maps, validation of such structures is an important task. Methods: We applied a histogram-based outlier score (HBOS) to six sets of cryo-EM atomic structures and five sets of X-ray atomic structures, including one derived from X-ray data with better than 1.5 Å resolution. Cryo-EM data sets contain structures released by December 2016 and those released between 2017 and 2019, derived from resolution ranges 0-4 Å and 4-6 Å respectively. Results: The distribution of HBOS values in five sets of X-ray structures show that HBOS is sensitive distinguishing sets of X-ray structures derived from different resolution ranges-higher than 1.5 Å, 1.5-2.0 Å, 2.0-2.5 Å, 2.5-3.0 Å, and 3.0-3.5 Å. The overall quality of cryo-EM structures is likely improved, as shown in a comparison of cryo-EM structures released before the end of 2016, those between 2017 and 2018, and those between 2018 and 2019. Our investigation shows that leucine (LEU) has a significantly higher rate of HBOS outliers than that of the reference data set (X-ray-1.5) and of other residue types in the cryo-EM data sets. HBOS was able to detect outliers for those residues that are currently marked as green in PDB validation reports. Conclusions: The HBOS profile of a dataset is a potential method to characterize the overall structural quality of the set. Residue LEU deserves special attention since it has a significantly higher HBOS outlier rate in sets of cryo-EM structures and those X-ray structures derived from X-ray data of lower than 2.5 Å resolutions. Most HBOS outlier residues from the EM-0-4-2019 set are located on loops for most types of residues.

Modeling of a 14 kDa RUVBL2-Binding Domain with Medium Resolution Cryo-EM Density

The number of high-resolution structures of protein complexes obtained using cryo-electron microscopy (cryo-EM) is increasing rapidly. Cryo-EM maps of large macromolecular complexes frequently contain regions resolved at different resolution levels, and modeling atomic structures de novo can be difficult for domains determined at worse than 5 Å in the absence of atomic information from other structures. Here we describe the details and step-by-step decisions in the strategy we followed to model the RUVBL2-binding domain (RBD), a 14 kDa domain at the C-terminus of RNA Polymerase II associated protein 3 (RPAP3) for which atomic information was not available. Modeling was performed on a cryo-EM map at 4.0-5.5 Å resolution, integrating information from secondary structure predictions, homology modeling, restraints from cross-linked mass spectrometry, and molecular dynamics (MD) in AMBER. Here, we compare our model with the structure of RBD determined by NMR to evaluate our strategy. We also perform new MD simulations to describe important residues mediating the interaction of RBD with RUVBL2 and analyze their conservation in RBD homologous domains. Our approach and its evaluation can serve as an example to address the analysis of medium resolution regions in cryo-EM maps.

Deep Learning to Predict Protein Backbone Structure from High-Resolution Cryo-EM Density Maps

Cryo-electron microscopy (cryo-EM) has become a leading technology for determining protein structures. Recent advances in this field have allowed for atomic resolution. However, predicting the backbone trace of a protein has remained a challenge on all but the most pristine density maps (<2.5 Å resolution). Here we introduce a deep learning model that uses a set of cascaded convolutional neural networks (CNNs) to predict Cα atoms along a protein's backbone structure. The cascaded-CNN (C-CNN) is a novel deep learning architecture comprised of multiple CNNs, each predicting a specific aspect of a protein's structure. This model predicts secondary structure elements (SSEs), backbone structure, and Cα atoms, combining the results of each to produce a complete prediction map. The cascaded-CNN is a semantic segmentation image classifier and was trained using thousands of simulated density maps. This method is largely automatic and only requires a recommended threshold value for each protein density map. A specialized tabu-search path walking algorithm was used to produce an initial backbone trace with Cα placements. A helix-refinement algorithm made further improvements to the α-helix SSEs of the backbone trace. Finally, a novel quality assessment-based combinatorial algorithm was used to effectively map protein sequences onto Cα traces to obtain full-atom protein structures. This method was tested on 50 experimental maps between 2.6 Å and 4.4 Å resolution. It outperformed several state-of-the-art prediction methods including Rosetta de-novo, MAINMAST, and a Phenix based method by producing the most complete predicted protein structures, as measured by percentage of found Cα atoms. This method accurately predicted 88.9% (mean) of the Cα atoms within 3 Å of a protein's backbone structure surpassing the 66.8% mark achieved by the leading alternate method (Phenix based fully automatic method) on the same set of density maps. The C-CNN also achieved an average root-mean-square deviation (RMSD) of 1.24 Å on a set of 50 experimental density maps which was tested by the Phenix based fully automatic method. The source code and demo of this research has been published at https://github.com/DrDongSi/Ca-Backbone-Prediction.

Single-particle cryo-electron microscopy: Mathematical theory, computational challenges, and opportunities

In recent years, an abundance of new molecular structures have been elucidated using cryo-electron microscopy (cryo-EM), largely due to advances in hardware technology and data processing techniques. Owing to these new exciting developments, cryo-EM was selected by Nature Methods as Method of the Year 2015, and the Nobel Prize in Chemistry 2017 was awarded to three pioneers in the field. The main goal of this article is to introduce the challenging and exciting computational tasks involved in reconstructing 3-D molecular structures by cryo-EM. Determining molecular structures requires a wide range of computational tools in a variety of fields, including signal processing, estimation and detection theory, high-dimensional statistics, convex and non-convex optimization, spectral algorithms, dimensionality reduction, and machine learning. The tools from these fields must be adapted to work under exceptionally challenging conditions, including extreme noise levels, the presence of missing data, and massively large datasets as large as several Terabytes. In addition, we present two statistical models: multi-reference alignment and multi-target detection, that abstract away much of the intricacies of cryo-EM, while retaining some of its essential features. Based on these abstractions, we discuss some recent intriguing results in the mathematical theory of cryo-EM, and delineate relations with group theory, invariant theory, and information theory.

Fragment-based drug discovery using cryo-EM

Recent advances in electron cryo-microscopy (cryo-EM) structure determination have pushed the resolutions obtainable by the method into the range widely considered to be of utility for drug discovery. Here, we review the use of cryo-EM in fragment-based drug discovery (FBDD) based on in-house method development. We demonstrate not only that cryo-EM can reveal details of the molecular interactions between fragments and a protein, but also that the current reproducibility, quality, and throughput are compatible with FBDD. We exemplify this using the test system β-galactosidase (Bgal) and the oncology target pyruvate kinase 2 (PKM2).

Sub-2 Angstrom resolution structure determination using single-particle cryo-EM at 200 keV

Although the advent of direct electron detectors (DEDs) and software developments have enabled the routine use of single-particle cryogenic electron microscopy (cryo-EM) for structure determination of well-behaved specimens to high-resolution, there nonetheless remains a discrepancy between the resolutions attained for biological specimens and the information limits of modern transmission electron microscopes (TEMs). Instruments operating at 300 kV equipped with DEDs are the current paradigm for high-resolution single-particle cryo-EM, while 200 kV TEMs remain comparatively underutilized for purposes beyond sample screening. Here, we expand upon our prior work and demonstrate that one such 200 kV microscope, the Talos Arctica, equipped with a K2 DED is capable of determining structures of macromolecules to as high as ∼1.7 Å resolution. At this resolution, ordered water molecules are readily assigned and holes in aromatic residues can be clearly distinguished in the reconstructions. This work emphasizes the utility of 200 kV electrons for high-resolution single-particle cryo-EM and applications such as structure-based drug design.

Computational Protocol for Assessing the Optimal Pixel Size to Improve the Accuracy of Single-particle Cryo-electron Microscopy Maps

Cryo-electron microscopy (cryo-EM) single-particle analysis has come a long way in achieving atomic-level resolution when imaging biomolecules. To obtain the best possible three-dimensional (3D) structure in cryo-EM, many parameters have to be carefully considered. Here we address the often-overlooked parameter of the pixel size, which describes the magnification of the image produced by the experiment. While efforts are made to refine and validate this parameter in the analysis of cryo-EM experimental data, there is no systematic protocol in place. Since the pixel size parameter can have an impact on the resolution and accuracy of a cryo-EM map, and the atomic resolution 3D structure models derived from it, we propose a computational protocol to estimate the appropriate pixel size parameter. In our protocol, we fit and refine atomic structures against cryo-EM maps at multiple pixel sizes. The resulting fitted and refined structures are evaluated using the GOAP (generalized orientation-dependent, all-atom statistical potential) score, which we found to perform better than other commonly used functions, such as Molprobity and the correlation coefficient from refinement. Finally, we describe the efficacy of this protocol in retrieving appropriate pixel sizes for several examples; simulated data based on yeast elongation factor 2 and experimental data from Gro-EL chaperone, beta-galactosidase, and the TRPV1 ion channel.

In situ structure determination at nanometer resolution using TYGRESS

The resolution of subtomogram averages calculated from cryo-electron tomograms (cryo-ET) of crowded cellular environments is often limited owing to signal loss in, and misalignment of, the subtomograms. By contrast, single-particle cryo-electron microscopy (SP-cryo-EM) routinely reaches near-atomic resolution of isolated complexes. We report a method called 'tomography-guided 3D reconstruction of subcellular structures' (TYGRESS) that is a hybrid of cryo-ET and SP-cryo-EM, and is able to achieve close-to-nanometer resolution of complexes inside crowded cellular environments. TYGRESS combines the advantages of SP-cryo-EM (images with good signal-to-noise ratio and contrast, as well as minimal radiation damage) and subtomogram averaging (three-dimensional alignment of macromolecules in a complex sample). Using TYGRESS, we determined the structure of the intact ciliary axoneme with up to resolution of 12 Å. These results reveal many structural details that were not visible by cryo-ET alone. TYGRESS is generally applicable to cellular complexes that are amenable to subtomogram averaging.

Cryo-EM reconstruction of continuous heterogeneity by Laplacian spectral volumes

Single-particle electron cryomicroscopy is an essential tool for high-resolution 3D reconstruction of proteins and other biological macromolecules. An important challenge in cryo-EM is the reconstruction of non-rigid molecules with parts that move and deform. Traditional reconstruction methods fail in these cases, resulting in smeared reconstructions of the moving parts. This poses a major obstacle for structural biologists, who need high-resolution reconstructions of entire macromolecules, moving parts included. To address this challenge, we present a new method for the reconstruction of macromolecules exhibiting continuous heterogeneity. The proposed method uses projection images from multiple viewing directions to construct a graph Laplacian through which the manifold of three-dimensional conformations is analyzed. The 3D molecular structures are then expanded in a basis of Laplacian eigenvectors, using a novel generalized tomographic reconstruction algorithm to compute the expansion coefficients. These coefficients, which we name spectral volumes, provide a high-resolution visualization of the molecular dynamics. We provide a theoretical analysis and evaluate the method empirically on several simulated data sets.

Disruption of the HIV-1 Envelope allosteric network blocks CD4-induced rearrangements

The trimeric HIV-1 Envelope protein (Env) mediates viral-host cell fusion via a network of conformational transitions, with allosteric elements in each protomer orchestrating host receptor-induced exposure of the co-receptor binding site and fusion elements. To understand the molecular details of this allostery, here, we introduce Env mutations aimed to prevent CD4-induced rearrangements in the HIV-1 BG505 Env trimer. Binding analysis and single-molecule Förster Resonance Energy Transfer confirm that these mutations prevent CD4-induced transitions of the HIV-1 Env. Structural analysis by single-particle cryo-electron microscopy performed on the BG505 SOSIP mutant Env proteins shows rearrangements in the gp120 topological layer contacts with gp41. Displacement of a conserved tryptophan (W571) from its typical pocket in these Env mutants renders the Env insensitive to CD4 binding. These results reveal the critical function of W571 as a conformational switch in Env allostery and receptor-mediated viral entry and provide insights on Env conformation that are relevant for vaccine design.

Why recombinant antibodies - benefits and applications

Antibodies (Abs) are ubiquitous reagents for biological and biochemical research and are rapidly expanding into new therapeutic areas. They are one of the most important probes for determining how proteins function under normal and pathophysiological conditions. Abs are required for the quantification of targets, detection of temporal and spatial patterns of protein expression in cells and tissues, and identification of interacting partners and their biological activities. Their remarkable specificity and unique binding properties can facilitate three-dimensional structure determination using X-ray crystallography and electron cryomicroscopy. While hybridoma technology that involves animal immunization is often productive, many antigen targets do not generate useful Abs. This is particularly true if unique states of the target or critical non-immunogenic target sequences need to be recognized by the Abs. By using the methods of recombinant antibody generation, identification, and engineering, these 'hybridoma-refractory' antigens can be readily targeted. Specific, reproducible, and renewable recombinant Abs are proving to be invaluable reagents in applications ranging from biological discovery to structure determination of challenging macromolecules.

Virus Structures by X-Ray Free-Electron Lasers

Until recently X-ray crystallography has been the standard technique for virus structure determinations. Available X-ray sources have continuously improved over the decades, leading to the realization of X-ray free-electron lasers (XFELs). They provide high-intensity femtosecond X-ray pulses, which allow for new kinds of experiments by making use of the diffraction-before-destruction principle. By overcoming classical dose constraints, they at least in principle allow researchers to perform X-ray virus structure determination for single particles at room temperature. Simultaneously, the availability of XFELs led to the development of the method of serial femtosecond crystallography, where a crystal structure is determined from the measurement of hundreds to thousands of microcrystals. In the case of virus crystallography this method does not require freezing of the crystals and allows researchers to perform experiments under non-equilibrium conditions (e.g., by laser-induced temperature jumps or rapid chemical mixing), which is currently not possible with electron microscopy.

Non-uniformity of projection distributions attenuates resolution in Cryo-EM

Virtually all single-particle cryo-EM experiments currently suffer from specimen adherence to the air-water interface, leading to a non-uniform distribution in the set of projection views. Whereas it is well accepted that uniform projection distributions can lead to high-resolution reconstructions, non-uniform (anisotropic) distributions can negatively affect map quality, elongate structural features, and in some cases, prohibit interpretation altogether. Although some consequences of non-uniform sampling have been described qualitatively, we know little about how sampling quantitatively affects resolution in cryo-EM. Here, we show how inhomogeneity in any projection distribution scheme attenuates the global Fourier Shell Correlation (FSC) in relation to the number of particles and a single geometrical parameter, which we term the sampling compensation factor (SCF). The reciprocal of the SCF is defined as the average over Fourier shells of the reciprocal of the per-particle sampling and normalized to unity for uniform distributions. The SCF therefore ranges from one to zero, with values close to the latter implying large regions of poorly sampled or completely missing data in Fourier space. Using two synthetic test cases, influenza hemagglutinin and human apoferritin, we demonstrate how any amount of sampling inhomogeneity always attenuates the FSC compared to a uniform distribution. We advocate quantitative evaluation of the SCF criterion to approximate the effect of non-uniform sampling on resolution within experimental single-particle cryo-EM reconstructions.

CryoEM maps are full of potential

Electron microscopy is based on elastic scattering due to Coulomb forces between the incident electrons and the sample; thus, electron scattering is dependent on the charge distribution in the sample. Unlike atomic scattering factors for X-rays, electron scattering factors for some atoms are strongly dependent on scattering angle, and the scattering factor for ionic oxygen is negative at low scattering angle. This phenomenon can result in a significant negative contribution to Coulomb potential maps by oxygen and can result in deviations in the positions of positive map features from atomic centers. An important factor that can also complicate the interpretation of cryoEM maps is the exquisite sensitivity of macromolecules to damage from electron irradiation, especially the carboxylates of acidic amino acids. Ideally, when compared with electron density maps derived by X-ray crystallography, Coulomb potential maps can provide additional details about the electrostatic environment and charge state of atoms. Enhancements in model building, refinement and computational simulation will be required to realize the full potential of EM-derived maps to reveal deeper insight into the electronic structure and functional properties of macromolecular complexes and their interactions with binding partners, ligands, cofactors, and drugs.

Analysis of discrete local variability and structural covariance in macromolecular assemblies using Cryo-EM and focused classification

Single-particle electron cryo-microscopy and computational image classification can be used to analyze structural variability in macromolecules and their assemblies. In some cases, a particle may contain different regions that each display a range of distinct conformations. We have developed strategies, implemented within the Frealign and cisTEM image processing packages, to focus-classify on specific regions of a particle and detect potential covariance. The strategies are based on masking the region of interest using either a 2-D mask applied to reference projections and particle images, or a 3-D mask applied to the 3-D volume. We show that focused classification approaches can be used to study structural covariance, a concept that is likely to gain more importance as datasets grow in size, allowing the distinction of more structural states and smaller differences between states. Finally, we apply the approaches to an experimental dataset containing the HIV-1 Transactivation Response (TAR) element RNA fused into the large bacterial ribosomal subunit to deconvolve structural mobility within localized regions of interest, and to a dataset containing assembly intermediates of the large subunit to measure structural covariance.

Interpretation of medium resolution cryoEM maps of multi-protein complexes

•

CryoEM maps at medium (3.5–6 Å) resolution can be challenging to interpret.

•

Integration of multiple methods can inform cryoEM studies.

•

Mass spectrometry and biochemistry facilitate map interpretation and model building.

Electron cryo-microscopy (cryoEM) is used to determine structures of biological molecules, including multi-protein complexes. Maps at better than 3.0 Å resolution are relatively straightforward to interpret since atomic models of proteins and nucleic acids can be built directly. Still, these resolutions are often difficult to achieve, and map quality frequently varies within a structure. This results in data that are challenging to interpret, especially when crystal structures or suitable homology models are not available. Recent advances in mass spectrometry techniques, computational methods and model building tools facilitate subunit/domain fitting into maps, elucidation of protein contacts, and de novo generation of atomic models. Here, we review techniques for map interpretation and provide examples from recent studies of multi-protein complexes.

X-ray radiation damage to biological samples: recent progress

With the continuing development of beamlines for macromolecular crystallography (MX) over the last few years providing ever higher X-ray flux densities, it has become even more important to be aware of the effects of radiation damage on the resulting structures. Nine papers in this issue cover a range of aspects related to the physics and chemistry of the manifestations of this damage, as observed in both MX and small-angle X-ray scattering (SAXS) on crystals, solutions and tissue samples. The reports include measurements of the heating caused by X-ray irradiation in ruby microcrystals, low-dose experiments examining damage rates as a function of incident X-ray energy up to 30 keV on a metallo-enzyme using a CdTe detector of high quantum efficiency as well as a theoretical analysis of the gains predicted in diffraction efficiency using these detectors, a SAXS examination of low-dose radiation exposure effects on the dissociation of a protein complex related to human health, theoretical calculations describing radiation chemistry pathways which aim to explain the specific structural damage widely observed in proteins, investigation of radiation-induced damage effects in a DNA crystal, a case study on a metallo-enzyme where structural movements thought to be mechanism related might actually be radiation-damage-induced changes, and finally a review describing what X-ray radiation-induced cysteine modifications can teach us about protein dynamics and catalysis. These papers, along with some other relevant literature published since the last Journal of Synchrotron Radiation Radiation Damage special issue in 2017, are briefly summarized below.

Single particle cryo-EM reconstruction of 52 kDa streptavidin at 3.2 Angstrom resolution

The fast development of single-particle cryogenic electron microscopy (cryo-EM) has made it more feasible to obtain the 3D structure of well-behaved macromolecules with a molecular weight higher than 300 kDa at ~3 Å resolution. However, it remains a challenge to obtain the high-resolution structures of molecules smaller than 200 kDa using single-particle cryo-EM. In this work, we apply the Cs-corrector-VPP-coupled cryo-EM to study the 52 kDa streptavidin (SA) protein supported on a thin layer of graphene and embedded in vitreous ice. We are able to solve both the apo-SA and biotin-bound SA structures at near-atomic resolution using single-particle cryo-EM. We demonstrate that the method has the potential to determine the structures of molecules as small as 39 kDa.

It remains a challenge to obtain high-resolution structures of molecules smaller than 200 kDa using single particle cryo-EM. Here, the authors apply the Cs-corrector-VPP coupled cryo-EM and solve structures of the 52 kDa streptavidin (SA) protein at near-atomic resolution.

Cryo-EM structures of GPCRs coupled to Gs, Gi and Go

Advances in electron cryo-microscopy (cryo-EM) now permit the structure determination of G protein-coupled receptors (GPCRs) coupled to heterotrimeric G proteins by single-particle imaging. A combination of G protein engineering and the development of antibodies that stabilise the heterotrimeric G protein facilitate the formation of stable GPCR-G protein complexes suitable for structural biology. Structures have been determined of GPCRs coupled to either heterotrimeric G proteins (G_s, G_i or G_o) or mini-G proteins (mini-G_s or mini-G_o) by single-particle cryo-EM and X-ray crystallography, respectively. This review describes the technical breakthroughs allowing their structure determination and compares the different techniques. In addition, we compare the structures of GPCRs coupled either to G_s, G_i or G_o and analyse the contributions of amino acid residues to the GPCR-G protein interface. There is no unique set of interactions that specifies coupling either to G_s, G_i or G_o. Instead, there is a common core of interactions between the C-terminal α-helix of the G protein α-subunit and helices H3, H5 and H6 of the receptor. In addition, there are varying degrees of interaction between all the other GPCR helices and intracellular loops to five regions of the α-subunit and four regions of the β-subunit. These data support the contention that there is not a simple linear barcode that defines the specificity of G protein coupling and thus how a G protein couples to a GPCR cannot currently be determined from simply analysing amino acid sequences. Although the overall architecture of GPCR-G protein complexes is conserved, there are significant differences in the molecular details. The number and type of molecular interactions between amino acid residues at the interfaces varies, resulting in subtly different orientation and position of the G protein with respect to the GPCR. This in turn affects the interface surface area that varies between 845 Ų and 1490 Ų, which could impact upon the lifetime of signalling complexes in the cell.

Challenges and opportunities in cryo-EM single-particle analysis

Cryogenic electron microscopy (cryo-EM) enables structure determination of macromolecular objects and their assemblies. Although the techniques have been developing for nearly four decades, they have gained widespread attention in recent years due to technical advances on numerous fronts, enabling traditional microscopists to break into the world of molecular structural biology. Many samples can now be routinely analyzed at near-atomic resolution using standard imaging and image analysis techniques. However, numerous challenges to conventional workflows remain, and continued technical advances open entirely novel opportunities for discovery and exploration. Here, I will review some of the main methods surrounding cryo-EM with an emphasis specifically on single-particle analysis, and I will highlight challenges, open questions, and opportunities for methodology development.

Frontiers in Cryo Electron Microscopy of Complex Macromolecular Assemblies

In recent years, cryo electron microscopy (cryo-EM) technology has been transformed with the development of better instrumentation, direct electron detectors, improved methods for specimen preparation, and improved software for data analysis. Analyses using single-particle cryo-EM methods have enabled determination of structures of proteins with sizes smaller than 100 kDa and resolutions of ∼2 Å in some cases. The use of electron tomography combined with subvolume averaging is beginning to allow the visualization of macromolecular complexes in their native environment in unprecedented detail. As a result of these advances, solutions to many intractable challenges in structural and cell biology, such as analysis of highly dynamic soluble and membrane-embedded protein complexes or partially ordered protein aggregates, are now within reach. Recent reports of structural studies of G protein-coupled receptors, spliceosomes, and fibrillar specimens illustrate the progress that has been made using cryo-EM methods, and are the main focus of this review.

High-resolution structure determination of sub-100 kDa complexes using conventional cryo-EM

Determining high-resolution structures of biological macromolecules amassing less than 100 kilodaltons (kDa) has been a longstanding goal of the cryo-electron microscopy (cryo-EM) community. While the Volta phase plate has enabled visualization of specimens in this size range, this instrumentation is not yet fully automated and can present technical challenges. Here, we show that conventional defocus-based cryo-EM methodologies can be used to determine high-resolution structures of specimens amassing less than 100 kDa using a transmission electron microscope operating at 200 keV coupled with a direct electron detector. Our ~2.7 Å structure of alcohol dehydrogenase (82 kDa) proves that bound ligands can be resolved with high fidelity to enable investigation of drug-target interactions. Our ~2.8 Å and ~3.2 Å structures of methemoglobin demonstrate that distinct conformational states can be identified within a dataset for proteins as small as 64 kDa. Furthermore, we provide the sub-nanometer cryo-EM structure of a sub-50 kDa protein.

Extension of the transferable aspherical pseudoatom data bank for the comparison of molecular electrostatic potentials in structure-activity studies

The transferable aspherical pseudoatom data bank, UBDB2018, is extended with over 130 new atom types present in small and biological molecules of great importance in biology and chemistry. UBDB2018 can be applied either as a source of aspherical atomic scattering factors in a standard X-ray experiment (d_min ≃ 0.8 Å) instead of the independent atom model (IAM), and can therefore enhance the final crystal structure geometry and refinement parameters; or as a tool to reconstruct the molecular charge-density distribution and derive the electrostatic properties of chemical systems for which 3D structural data are available. The extended data bank has been extensively tested, with the focus being on the accuracy of the molecular electrostatic potential computed for important drug-like molecules, namely the HIV-1 protease inhibitors. The UBDB allows the reconstruction of the reference B3LYP/6-31G** potentials, with a root-mean-squared error of 0.015 e bohr^-1 computed for entire potential grids which span values from ca 200 e bohr^-1 to ca -0.1 e bohr^-1 and encompass both the inside and outside regions of a molecule. UBDB2018 is shown to be applicable to enhancing the physical meaning of the molecular electrostatic potential descriptors used to construct predictive quantitative structure-activity relationship/quantitative structure-property relationship (QSAR/QSPR) models for drug discovery studies. In addition, it is suggested that electron structure factors computed from UBDB2018 may significantly improve the interpretation of electrostatic potential maps measured experimentally by means of electron diffraction or single-particle cryo-EM methods.

Capabilities of the Falcon III detector for single-particle structure determination

Direct electron detectors are an essential asset for the resolution revolution in electron cryo microscopy of biological objects. The direct detectors provide two modes of data acquisition; the counting mode in which single electrons are counted, and the integrating mode in which the signal that arises from the incident electrons is integrated. While counting mode leads to far higher detective quantum efficiency at all spatial frequencies, the integrating mode enables faster data acquisition at higher exposure rates. For optimal throughput at best possible resolution it is important to understand when the better performance in counting mode becomes essential for solving a structure and when the lower detective quantum efficiency in integrating mode can be compensated by increasing the number of particles in the data set. Here, we provide a case study of the Falcon III camera, which has counting mode capability at exposure rates of <0.9 e^-/Px² and integrating mode capability at exposure rates above 10 e^-/Px². We found that counting mode gives better resolution for medium sized complexes such as the β-galactosidase (465 kDa) (2.2 Å, 97% of Nyquist vs. 2.4 Å, 89% of Nyquist) with data sets of similar size. However, for larger particles such as Hepatitis B virus capsid like particles (4.8 MDa) we did not find any resolution gain in counting mode.

Cryo-EM in drug discovery

The impact of structural biology on drug discovery is well documented, and the workhorse technique for the past 30 years or so has been X-ray crystallography. With the advent of several technological improvements, including direct electron detectors, automation, better microscope vacuums and lenses, phase plates and improvements in computing power enabled by GPUs, it is now possible to record and analyse images of protein structures containing high-resolution information. This review, from a pharmaceutical perspective, highlights some of the most relevant and interesting protein structures for the pharmaceutical industry and shows examples of how ligand-binding sites, membrane proteins, both big and small, pseudo symmetry and complexes are being addressed by this technique.

Thresholding of cryo-EM density maps by false discovery rate control

Cryo-EM now commonly generates close-to-atomic resolution as well as intermediate resolution maps from macromolecules observed in isolation and in situ. Interpreting these maps remains a challenging task owing to poor signal in the highest resolution shells and the necessity to select a threshold for density analysis. In order to facilitate this process, a statistical framework for the generation of confidence maps by multiple hypothesis testing and false discovery rate (FDR) control has been developed. In this way, three-dimensional confidence maps contain signal separated from background noise in the form of local detection rates of EM density values. It is demonstrated that confidence maps and FDR-based thresholding can be used for the interpretation of near-atomic resolution single-particle structures as well as lower resolution maps determined by subtomogram averaging. Confidence maps represent a conservative way of interpreting molecular structures owing to minimized noise. At the same time they provide a detection error with respect to background noise, which is associated with the density and is particularly beneficial for the interpretation of weaker cryo-EM densities in cases of conformational flexibility and lower occupancy of bound molecules and ions in the structure.

Active site flexibility revealed in crystal structures of Parabacteroides merdae β-glucuronidase from the human gut microbiome

β-Glucuronidase (GUS) enzymes in the gastrointestinal tract are involved in maintaining mammalian-microbial symbiosis and can play key roles in drug efficacy and toxicity. Parabacteroides merdae GUS was identified as an abundant mini-Loop 2 (mL2) type GUS enzyme in the Human Microbiome Project gut metagenomic database. Here, we report the crystal structure of P. merdae GUS and highlight the differences between this enzyme and extant structures of gut microbial GUS proteins. We find that P. merdae GUS exhibits a distinct tetrameric quaternary structure and that the mL2 motif traces a unique path within the active site, which also includes two arginines distinctive to this GUS. We observe two states of the P. merdae GUS active site; a loop repositions itself by more than 50 Å to place a functionally-relevant residue into the enzyme's catalytic site. Finally, we find that P. merdae GUS is able to bind to homo and heteropolymers of the polysaccharide alginic acid. Together, these data broaden our understanding of the structural and functional diversity in the GUS family of enzymes present in the human gut microbiome and point to specialization as an important feature of microbial GUS orthologs.

New tools for automated high-resolution cryo-EM structure determination in RELION-3

Here, we describe the third major release of RELION. CPU-based vector acceleration has been added in addition to GPU support, which provides flexibility in use of resources and avoids memory limitations. Reference-free autopicking with Laplacian-of-Gaussian filtering and execution of jobs from python allows non-interactive processing during acquisition, including 2D-classification, de novo model generation and 3D-classification. Per-particle refinement of CTF parameters and correction of estimated beam tilt provides higher resolution reconstructions when particles are at different heights in the ice, and/or coma-free alignment has not been optimal. Ewald sphere curvature correction improves resolution for large particles. We illustrate these developments with publicly available data sets: together with a Bayesian approach to beam-induced motion correction it leads to resolution improvements of 0.2–0.7 Å compared to previous RELION versions.

Sub-2 Å Ewald curvature corrected structure of an AAV2 capsid variant

Single-particle cryogenic electron microscopy (cryo-EM) provides a powerful methodology for structural biologists, but the resolutions typically attained with experimentally determined structures have lagged behind microscope capabilities. Here, we exploit several technical advances to improve resolution, including per-particle contrast transfer function (CTF) refinement and correction for Ewald sphere curvature. The latter is demonstrated with several experimental samples and should become more standard as resolutions increase or at lower microscope accelerating voltages. The combined application of the described methods to micrographs recorded on a Titan Krios enables structure determination at ~1.86-Å resolution of an adeno-associated virus serotype 2 variant (AAV2), an important gene-delivery vehicle. The resulting structural details provide an improved model for understanding the biology of AAV that will guide future vector development for gene therapy.

Single-particle cryo-EM is a powerful method for macromolecular structure determination. Here the authors demonstrate that Ewald sphere curvature correction, sub-Angstrom pixilation and per-particle CTF refinement can improve map quality and resolution and present the 1.86 Å cryo-EM structure of an adeno-associated virus serotype 2 variant.

Pushing the Limits of Computational Structure-Based Drug Design with a Cryo-EM Structure: The Ca2+ Channel α2δ-1 Subunit as a Test Case

Cryo-electron microscopy (cryo-EM) is emerging as a real alternative for structural elucidation. In spite of this, very few cryo-EM structures have been described so far as successful platforms for in silico drug design. Gabapentin and pregabalin are some of the most successful drugs in the treatment of epilepsy and neuropathic pain. Although both are in clinical use and are known to exert their effects by binding to the regulatory α2δ subunit of voltage gated calcium channels, their binding modes have never been characterized. We describe here the successful use of an exhaustive protein-ligand sampling algorithm on the α2δ-1 subunit of the recently published cryo-EM structure, with the goal of characterizing the ligand entry path and binding mode for gabapentin, pregabalin, and several other amino acidic α2δ-1 ligands. Our studies indicate that (i) all simulated drugs explore the same path for accessing the occluded binding site on the interior of the α2δ-1 subunit; (ii) they all roughly occupy the same pocket; (iii) the plasticity of the binding site allows the accommodation of a variety of amino acidic modulators, driven by the flexible "capping loop" delineated by residues Tyr426-Val435 and the floppy nature of Arg217; (iv) the predicted binding modes are in line with previously available mutagenesis data, confirming Arg217 as key for binding, with Asp428 and Asp467 highlighted as additional anchoring points for all amino acidic drugs. The study is one of the first proofs that latest-generation cryo-EM structures combined with exhaustive computational methods can be exploited in early drug discovery.

The first single particle analysis Map Challenge: A summary of the assessments

The recent successes of cryo-electron microscopy fostered great expectation of solving many new and previously recalcitrant biomolecular structures. However, it also brings with it the danger of compromising the validity of the outcomes if not done properly. The Map Challenge is a first step in assessing the state of the art and to shape future developments in data processing. The organizers presented seven cases for single particle reconstruction, and 27 members of the community responded with 66 submissions. Seven groups analyzed these submissions, resulting in several assessment reports, summarized here. We devised a range of analyses to evaluate the submitted maps, including visual impressions, Fourier shell correlation, pairwise similarity and interpretation through modeling. Unfortunately, we did not find strong trends. We ascribe this to the complexity of the challenge, dealing with multiple cases, software packages and processing approaches. This puts the user in the spotlight, where his/her choices becomes the determinant of map quality. The future focus should therefore be on promulgating best practices and encapsulating these in the software. Such practices include adherence to validation principles, most notably the processing of independent sets, proper resolution-limited alignment, appropriate masking and map sharpening. We consider the Map Challenge to be a highly valuable exercise that should be repeated frequently or on an ongoing basis.