A Multi-model Approach to Assessing Local and Global Cryo-EM Map Quality

Elucidating structural variability in p53 conformers using combinatorial refinement strategies and molecular dynamics

Low molecular weight proteins and protein assemblies can now be investigated using cryo-electron microscopy (EM) as a complement to traditional structural biology techniques. It is important, however, to not lose sight of the dynamic information inherent in macromolecules that give rise to their exquisite functionality. As computational methods continue to advance the field of biomedical imaging, so must strategies to resolve the minute details of disease-related entities. Here, we employed combinatorial modeling approaches to assess flexible properties among low molecular weight proteins (~100 kDa or less). Through a blend of rigid body refinement and simulated annealing, we determined new hidden conformations for wild type p53 monomer and dimer forms. Structures for both states converged to yield new conformers, each revealing good stereochemistry and dynamic information about the protein. Based on these insights, we identified fluid parts of p53 that complement the stable central core of the protein responsible for engaging DNA. Molecular dynamics simulations corroborated the modeling results and helped pinpoint the more flexible residues in wild type p53. Overall, the new computational methods may be used to shed light on other small protein features in a vast ensemble of structural data that cannot be easily delineated by other algorithms.

Cryo-EM unveils kinesin KIF1A's processivity mechanism and the impact of its pathogenic variant P305L

Mutations in the microtubule-associated motor protein KIF1A lead to severe neurological conditions known as KIF1A-associated neurological disorders (KAND). Despite insights into its molecular mechanism, high-resolution structures of KIF1A-microtubule complexes remain undefined. Here, we present 2.7-3.5 Å resolution structures of dimeric microtubule-bound KIF1A, including the pathogenic P305L mutant, across various nucleotide states. Our structures reveal that KIF1A binds microtubules in one- and two-heads-bound configurations, with both heads exhibiting distinct conformations with tight inter-head connection. Notably, KIF1A's class-specific loop 12 (K-loop) forms electrostatic interactions with the C-terminal tails of both α- and β-tubulin. The P305L mutation does not disrupt these interactions but alters loop-12's conformation, impairing strong microtubule-binding. Structure-function analysis reveals the K-loop and head-head coordination as major determinants of KIF1A's superprocessive motility. Our findings advance the understanding of KIF1A's molecular mechanism and provide a basis for developing structure-guided therapeutics against KAND.

Structure and topography of the synaptic V-ATPase-synaptophysin complex

Synaptic vesicles are organelles with a precisely defined protein and lipid composition1,2, yet the molecular mechanisms for the biogenesis of synaptic vesicles are mainly unknown. Here we discovered a well-defined interface between the synaptic vesicle V-ATPase and synaptophysin by in situ cryo-electron tomography and single-particle cryo-electron microscopy of functional synaptic vesicles isolated from mouse brains3. The synaptic vesicle V-ATPase is an ATP-dependent proton pump that establishes the proton gradient across the synaptic vesicle, which in turn drives the uptake of neurotransmitters4,5. Synaptophysin6 and its paralogues synaptoporin7 and synaptogyrin8 belong to a family of abundant synaptic vesicle proteins whose function is still unclear. We performed structural and functional studies of synaptophysin-knockout mice, confirming the identity of synaptophysin as an interaction partner with the V-ATPase. Although there is little change in the conformation of the V-ATPase upon interaction with synaptophysin, the presence of synaptophysin in synaptic vesicles profoundly affects the copy number of V-ATPases. This effect on the topography of synaptic vesicles suggests that synaptophysin assists in their biogenesis. In support of this model, we observed that synaptophysin-knockout mice exhibit severe seizure susceptibility, suggesting an imbalance of neurotransmitter release as a physiological consequence of the absence of synaptophysin.

A structure-based mechanism for initiation of AP-3 coated vesicle formation

Adaptor protein complex 3 (AP-3) mediates cargo sorting from endosomes to lysosomes and lysosome-related organelles. Recently, it was shown that AP-3 is in a constitutively open, active conformation compared to the related AP-1 and AP-2 coat complexes, which are inactive until undergoing large conformational changes upon membrane recruitment. How AP-3 is regulated is therefore an open question. To understand the mechanism of AP-3 membrane recruitment and activation, we reconstituted the core of human AP-3 and determined multiple structures in the soluble and membrane-bound states using electron cryo-microscopy (cryo-EM). Similar to yeast AP-3, human AP-3 is in a constitutively open conformation, with the cargo-binding domain of the μ3 subunit conformationally free. To reconstitute AP-3 activation by the small GTPase Arf1, we used lipid nanodiscs to build Arf1-AP-3 complexes on membranes and determined three structures that show the stepwise conformational changes required for formation of AP-3 coated vesicles. First, membrane-recruitment is driven by one of two predicted Arf1 binding sites on AP-3. In this conformation, AP-3 is flexibly tethered to the membrane and its cargo binding domain remains conformationally dynamic. Second, cargo binding causes AP-3 to adopt a fixed position and rigidifies the complex, which stabilizes binding for a second Arf1 molecule. Finally, binding of the second Arf1 molecule provides the template for AP-3 dimerization, providing a glimpse into the first step of coat polymerization. We propose coat polymerization only occurs after cargo engagement, thereby linking cargo sorting with assembly of higher order coat structures. Additionally, we provide evidence for two amphipathic helices in AP-3, suggesting that AP-3 contributes to membrane deformation during coat assembly. In total, these data provide evidence for the first stages of AP-3 mediated vesicle coat assembly.

Cryo-EM structure and B-factor refinement with ensemble representation

Cryo-EM experiments produce images of macromolecular assemblies that are combined to produce three-dimensional density maps. Typically, atomic models of the constituent molecules are fitted into these maps, followed by a density-guided refinement. We introduce TEMPy-ReFF, a method for atomic structure refinement in cryo-EM density maps. Our method represents atomic positions as components of a Gaussian mixture model, utilising their variances as B-factors, which are used to derive an ensemble description. Extensively tested on a substantial dataset of 229 cryo-EM maps from EMDB ranging in resolution from 2.1-4.9 Å with corresponding PDB and CERES atomic models, our results demonstrate that TEMPy-ReFF ensembles provide a superior representation of cryo-EM maps. On a single-model basis, it performs similarly to the CERES re-refinement protocol, although there are cases where it provides a better fit to the map. Furthermore, our method enables the creation of composite maps free of boundary artefacts. TEMPy-ReFF is useful for better interpretation of flexible structures, such as those involving RNA, DNA or ligands.

Delineating Conformational Variability in Small Protein Structures Using Combinatorial Refinement Strategies

As small protein assemblies and even small proteins are becoming more amenable to cryo-Electron Microscopy (EM) structural studies, it is important to consider the complementary dynamic information present in the data. Current computational strategies are limited in their ability to resolve minute differences among low molecular weight entities. Here, we demonstrate a new combinatorial approach to delineate flexible conformations among small proteins using real-space refinement applications. We performed a meta-analysis of structural data for the SARS CoV-2 Nucleocapsid (N) protein using a combination of rigid-body refinement and simulated annealing methods. For the N protein monomer, we determined three new flexible conformers with good stereochemistry and quantitative comparisons provided new evidence of their dynamic properties. A similar analysis performed for the N protein dimer showed only minor structural differences among the flexible models. These results suggested a more stable view of the N protein dimer than the monomer structure. Taken together, the new computational strategies can delineate conformational changes in low molecular weight proteins that may go unnoticed by conventional assessments. The results also suggest that small proteins may be further stabilized for structural studies through the use of solution components that limit the movement of external flexible regions.

Adaptive Ensemble Refinement of Protein Structures in High Resolution Electron Microscopy Density Maps with Radical Augmented Molecular Dynamics Flexible Fitting

Recent advances in cryo-electron microscopy (cryo-EM) have enabled modeling macromolecular complexes that are essential components of the cellular machinery. The density maps derived from cryo-EM experiments are often integrated with manual, knowledge-driven or artificial intelligence-driven and physics-guided computational methods to build, fit, and refine molecular structures. Going beyond a single stationary-structure determination scheme, it is becoming more common to interpret the experimental data with an ensemble of models that contributes to an average observation. Hence, there is a need to decide on the quality of an ensemble of protein structures on-the-fly while refining them against the density maps. We introduce such an adaptive decision-making scheme during the molecular dynamics flexible fitting (MDFF) of biomolecules. Using RADICAL-Cybertools, the new RADICAL augmented MDFF implementation (R-MDFF) is examined in high-performance computing environments for refinement of two prototypical protein systems, adenylate kinase and carbon monoxide dehydrogenase. For these test cases, use of multiple replicas in flexible fitting with adaptive decision making in R-MDFF improves the overall correlation to the density by 40% relative to the refinements of the brute-force MDFF. The improvements are particularly significant at high, 2-3 Å map resolutions. More importantly, the ensemble model captures key features of biologically relevant molecular dynamics that are inaccessible to a single-model interpretation. Finally, the pipeline is applicable to systems of growing sizes, which is demonstrated using ensemble refinement of capsid proteins from the chimpanzee adenovirus. The overhead for decision making remains low and robust to computing environments. The software is publicly available on GitHub and includes a short user guide to install R-MDFF on different computing environments, from local Linux-based workstations to high-performance computing environments.

CryoFold 2.0: Cryo-EM Structure Determination with MELD

Cryo-electron microscopy data are becoming more prevalent and accessible at higher resolution levels, leading to the development of new computational tools to determine the atomic structure of macromolecules. However, while existing tools adapted from X-ray crystallography are suitable for the highest-resolution maps, new tools are needed for lower-resolution levels and to account for map heterogeneity. In this article, we introduce CryoFold 2.0, an integrative physics-based approach that combines Bayesian inference and the ability to handle multiple data sources with the molecular dynamics flexible fitting (MDFF) approach to determine the structures of macromolecules by using cryo-EM data. CryoFold 2.0 is incorporated into the MELD (modeling employing limited data) plugin, resulting in a pipeline that is more computationally efficient and accurate than running MELD or MDFF alone. The approach requires fewer computational resources and shorter simulation times than the original CryoFold, and it minimizes manual intervention. We demonstrate the effectiveness of the approach on eight different systems, highlighting its various benefits.

Computational pipeline provides mechanistic understanding of Omicron variant of concern neutralizing engineered ACE2 receptor traps

The SARS-CoV-2 Omicron variant, with 15 mutations in Spike receptor-binding domain (Spike-RBD), renders virtually all clinical monoclonal antibodies against WT SARS-CoV-2 ineffective. We recently engineered the SARS-CoV-2 host entry receptor, ACE2, to tightly bind WT-RBD and prevent viral entry into host cells ("receptor traps"). Here we determine cryo-EM structures of our receptor traps in complex with stabilized Spike ectodomain. We develop a multi-model pipeline combining Rosetta protein modeling software and cryo-EM to allow interface energy calculations even at limited resolution and identify interface side chains that allow for high-affinity interactions between our ACE2 receptor traps and Spike-RBD. Our structural analysis provides a mechanistic rationale for the high-affinity (0.53-4.2 nM) binding of our ACE2 receptor traps to Omicron-RBD confirmed with biolayer interferometry measurements. Finally, we show that ACE2 receptor traps potently neutralize Omicron and Delta pseudotyped viruses, providing alternative therapeutic routes to combat this evolving virus.

Image processing tools for the validation of CryoEM maps

The number of maps deposited in public databases (Electron Microscopy Data Bank, EMDB) determined by cryo-electron microscopy has quickly grown in recent years. With this rapid growth, it is critical to guarantee their quality. So far, map validation has primarily focused on the agreement between maps and models. From the image processing perspective, the validation has been mostly restricted to using two half-maps and the measurement of their internal consistency. In this article, we suggest that map validation can be taken much further from the point of view of image processing if 2D classes, particles, angles, coordinates, defoci, and micrographs are also provided. We present a progressive validation scheme that qualifies a result validation status from 0 to 5 and offers three optional qualifiers (A, W, and O) that can be added. The simplest validation state is 0, while the most complete would be 5AWO. This scheme has been implemented in a website https://biocomp.cnb.csic.es/EMValidationService/ to which reconstructed maps and their ESI can be uploaded.

Overview and applications of map and model validation tools in the CCP-EM software suite

Cryogenic electron microscopy (cryo-EM) has recently been established as a powerful technique for solving macromolecular structures. Although the best resolutions achievable are improving, a significant majority of data are still resolved at resolutions worse than 3 Å, where it is non-trivial to build or fit atomic models. The map reconstructions and atomic models derived from the maps are also prone to errors accumulated through the different stages of data processing. Here, we highlight the need to evaluate both model geometry and fit to data at different resolutions. Assessment of cryo-EM structures from SARS-CoV-2 highlights a bias towards optimising the model geometry to agree with the most common conformations, compared to the agreement with data. We present the CoVal web service which provides multiple validation metrics to reflect the quality of atomic models derived from cryo-EM data of structures from SARS-CoV-2. We demonstrate that further refinement can lead to improvement of the agreement with data without the loss of geometric quality. We also discuss the recent CCP-EM developments aimed at addressing some of the current shortcomings.

Molecular insights into antibody-mediated protection against the prototypic simian immunodeficiency virus

SIVmac239 infection of macaques is a favored model of human HIV infection. However, the SIVmac239 envelope (Env) trimer structure, glycan occupancy, and the targets and ability of neutralizing antibodies (nAbs) to protect against SIVmac239 remain unknown. Here, we report the isolation of SIVmac239 nAbs that recognize a glycan hole and the V1/V4 loop. A high-resolution structure of a SIVmac239 Env trimer-nAb complex shows many similarities to HIV and SIVcpz Envs, but with distinct V4 features and an extended V1 loop. Moreover, SIVmac239 Env has a higher glycan shield density than HIV Env that may contribute to poor or delayed nAb responses in SIVmac239-infected macaques. Passive transfer of a nAb protects macaques from repeated intravenous SIVmac239 challenge at serum titers comparable to those described for protection of humans against HIV infection. Our results provide structural insights for vaccine design and shed light on antibody-mediated protection in the SIV model.

Computational pipeline provides mechanistic understanding of Omicron variant of concern neutralizing engineered ACE2 receptor traps

The SARS-CoV-2 Omicron variant, with 15 mutations in Spike receptor binding domain (Spike-RBD), renders virtually all clinical monoclonal antibodies against WT SARS-CoV-2 ineffective. We recently engineered the SARS-CoV-2 host entry receptor, ACE2, to tightly bind WT-Spike-RBD and prevent viral entry into host cells ("receptor traps"). Here we determine cryo-EM structures of our receptor traps in complex with full length Spike. We develop a multi-model pipeline combining Rosetta protein modeling software and cryo-EM to allow interface energy calculations even at limited resolution and identify interface side chains that allow for high affinity interactions between our ACE2 receptor traps and Spike-RBD. Our structural analysis provides a mechanistic rationale for the high affinity (0.53 - 4.2nM) binding of our ACE2 receptor traps to Omicron-RBD confirmed with biolayer interferometry measurements. Finally, we show that ACE2 receptor traps potently neutralize Omicron- and Delta-pseudotyped viruses, providing alternative therapeutic routes to combat this evolving virus.

Exploring cryo-electron microscopy with molecular dynamics

Single particle analysis cryo-electron microscopy (EM) and molecular dynamics (MD) have been complimentary methods since cryo-EM was first applied to the field of structural biology. The relationship started by biasing structural models to fit low-resolution cryo-EM maps of large macromolecular complexes not amenable to crystallization. The connection between cryo-EM and MD evolved as cryo-EM maps improved in resolution, allowing advanced sampling algorithms to simultaneously refine backbone and sidechains. Moving beyond a single static snapshot, modern inferencing approaches integrate cryo-EM and MD to generate structural ensembles from cryo-EM map data or directly from the particle images themselves. We summarize the recent history of MD innovations in the area of cryo-EM modeling. The merits for the myriad of MD based cryo-EM modeling methods are discussed, as well as, the discoveries that were made possible by the integration of molecular modeling with cryo-EM. Lastly, current challenges and potential opportunities are reviewed.

CryoFold: determining protein structures and data-guided ensembles from cryo-EM density maps

Cryo-electron microscopy (EM) requires molecular modeling to refine structural details from data. Ensemble models arrive at low free-energy molecular structures, but are computationally expensive and limited to resolving only small proteins that cannot be resolved by cryo-EM. Here, we introduce CryoFold - a pipeline of molecular dynamics simulations that determines ensembles of protein structures directly from sequence by integrating density data of varying sparsity at 3-5 Å resolution with coarse-grained topological knowledge of the protein folds. We present six examples showing its broad applicability for folding proteins between 72 to 2000 residues, including large membrane and multi-domain systems, and results from two EMDB competitions. Driven by data from a single state, CryoFold discovers ensembles of common low-energy models together with rare low-probability structures that capture the equilibrium distribution of proteins constrained by the density maps. Many of these conformations, unseen by traditional methods, are experimentally validated and functionally relevant. We arrive at a set of best practices for data-guided protein folding that are controlled using a Python GUI.

Validation, analysis and annotation of cryo-EM structures

The process of turning 2D micrographs into 3D atomic models of the imaged macromolecules has been under rapid development and scrutiny in the field of cryo-EM. Here, some important methods for validation at several stages in this process are described. Firstly, how Fourier shell correlation of two independent maps and phase randomization beyond a certain frequency address the assessment of map resolution is reviewed. Techniques for local resolution estimation and map sharpening are also touched upon. The topic of validating models which are either built de novo or based on a known atomic structure fitted into a cryo-EM map is then approached. Map-model comparison using Q-scores and Fourier shell correlation plots is used to assure the agreement of the model with the observed map density. The importance of annotating the model with B factors to account for the resolvability of individual atoms in the map is illustrated. Finally, the timely topic of detecting and validating water molecules and metal ions in maps that have surpassed ∼2 Å resolution is described.

Structural basis of mechano-chemical coupling by the mitotic kinesin KIF14

KIF14 is a mitotic kinesin whose malfunction is associated with cerebral and renal developmental defects and several cancers. Like other kinesins, KIF14 couples ATP hydrolysis and microtubule binding to the generation of mechanical work, but the coupling mechanism between these processes is still not fully clear. Here we report 20 high-resolution (2.7-3.9 Å) cryo-electron microscopy KIF14-microtubule structures with complementary functional assays. Analysis procedures were implemented to separate coexisting conformations of microtubule-bound monomeric and dimeric KIF14 constructs. The data provide a comprehensive view of the microtubule and nucleotide induced KIF14 conformational changes. It shows that: 1) microtubule binding, the nucleotide species, and the neck-linker domain govern the transition between three major conformations of the motor domain; 2) an undocked neck-linker prevents the nucleotide-binding pocket to fully close and dampens ATP hydrolysis; 3) 13 neck-linker residues are required to assume a stable docked conformation; 4) the neck-linker position controls the hydrolysis rather than the nucleotide binding step; 5) the two motor domains of KIF14 dimers adopt distinct conformations when bound to the microtubule; and 6) the formation of the two-heads-bound-state introduces structural changes in both motor domains of KIF14 dimers. These observations provide the structural basis for a coordinated chemo-mechanical kinesin translocation model.

Molecular Insights Into the Gating Kinetics of the Cardiac hERG Channel, Illuminated by Structure and Molecular Dynamics

The rapidly activating delayed rectifier K⁺ current generated by the cardiac hERG potassium channel encoded by KCNH2 is the most important reserve current for cardiac repolarization. The unique inward rectification characteristics of the hERG channel depend on the gating regulation, which involves crucial structural domains and key single amino acid residues in the full-length hERG channel. Identifying critical molecules involved in the regulation of gating kinetics for the hERG channel requires high-resolution structures and molecular dynamics simulation models. Based on the latest progress in hERG structure and molecular dynamics simulation research, summarizing the molecules involved in the changes in the channel state helps to elucidate the unique gating characteristics of the channel and the reason for its high affinity to cardiotoxic drugs. In this review, we aim to summarize the significant advances in understanding the voltage gating regulation of the hERG channel based on its structure obtained from cryo-electron microscopy and computer simulations, which reveal the critical roles of several specific structural domains and amino acid residues.

Conformational Plasticity of Hepatitis B Core Protein Spikes Promotes Peptide Binding Independent of the Secretion Phenotype

Hepatitis B virus is a major human pathogen, which forms enveloped virus particles. During viral maturation, membrane-bound hepatitis B surface proteins package hepatitis B core protein capsids. This process is intercepted by certain peptides with an “LLGRMKG” motif that binds to the capsids at the tips of dimeric spikes. With microcalorimetry, electron cryo microscopy and peptide microarray-based screens, we have characterized the structural and thermodynamic properties of peptide binding to hepatitis B core protein capsids with different secretion phenotypes. The peptide “GSLLGRMKGA” binds weakly to hepatitis B core protein capsids and mutant capsids with a premature (F97L) or low-secretion phenotype (L60V and P5T). With electron cryo microscopy, we provide novel structures for L60V and P5T and demonstrate that binding occurs at the tips of the spikes at the dimer interface, splaying the helices apart independent of the secretion phenotype. Peptide array screening identifies “SLLGRM” as the core binding motif. This shortened motif binds only to one of the two spikes in the asymmetric unit of the capsid and induces a much smaller conformational change. Altogether, these comprehensive studies suggest that the tips of the spikes act as an autonomous binding platform that is unaffected by mutations that affect secretion phenotypes.

Below 3 Å structure of apoferritin using a multipurpose TEM with a side entry cryoholder

Recently, the structural analysis of protein complexes by cryo-electron microscopy (cryo-EM) single particle analysis (SPA) has had great impact as a biophysical method. Many results of cryo-EM SPA are based on data acquired on state-of-the-art cryo-electron microscopes customized for SPA. These are currently only available in limited locations around the world, where securing machine time is highly competitive. One potential solution for this time-competitive situation is to reuse existing multi-purpose equipment, although this comes with performance limitations. Here, a multi-purpose TEM with a side entry cryo-holder was used to evaluate the potential of high-resolution SPA, resulting in a 3 Å resolution map of apoferritin with local resolution extending to 2.6 Å. This map clearly showed two positions of an aromatic side chain. Further, examination of optimal imaging conditions depending on two different multi-purpose electron microscope and camera combinations was carried out, demonstrating that higher magnifications are not always necessary or desirable. Since automation is effectively a requirement for large-scale data collection, and augmenting the multi-purpose equipment is possible, we expanded testing by acquiring data with SerialEM using a β-galactosidase test sample. This study demonstrates the possibilities of more widely available and established electron microscopes, and their applications for cryo-EM SPA.

State-dependent sequential allostery exhibited by chaperonin TRiC/CCT revealed by network analysis of Cryo-EM maps

The eukaryotic chaperonin TRiC/CCT plays a major role in assisting the folding of many proteins through an ATP-driven allosteric cycle. Recent structures elucidated by cryo-electron microscopy provide a broad view of the conformations visited at various stages of the chaperonin cycle, including a sequential activation of its subunits in response to nucleotide binding. But we lack a thorough mechanistic understanding of the structure-based dynamics and communication properties that underlie the TRiC/CCT machinery. In this study, we present a computational methodology based on elastic network models adapted to cryo-EM density maps to gain a deeper understanding of the structure-encoded allosteric dynamics of this hexadecameric machine. We have analysed several structures of the chaperonin resolved in different states toward mapping its conformational landscape. Our study indicates that the overall architecture intrinsically favours cooperative movements that comply with the structural variabilities observed in experiments. Furthermore, the individual subunits CCT1-CCT8 exhibit state-dependent sequential events at different states of the allosteric cycle. For example, in the ATP-bound state, subunits CCT5 and CCT4 selectively initiate the lid closure motions favoured by the overall architecture; whereas in the apo form of the heteromer, the subunit CCT7 exhibits the highest predisposition to structural change. The changes then propagate through parallel fluxes of allosteric signals to neighbours on both rings. The predicted state-dependent mechanisms of sequential activation provide new insights into TRiC/CCT intra- and inter-ring signal transduction events.

The accuracy of protein models automatically built into cryo-EM maps with ARP/wARP

A new module of the ARP/wARP suite for automated model building into cryo-EM maps is presented.

Recent developments in cryogenic electron microscopy (cryo-EM) have enabled structural studies of large macromolecular complexes at resolutions previously only attainable using macromolecular crystallography. Although a number of methods can already assist in de novo building of models into high-resolution cryo-EM maps, automated and reliable map interpretation remains a challenge. Presented here is a systematic study of the accuracy of models built into cryo-EM maps using ARP/wARP. It is demonstrated that the local resolution is a good indicator of map interpretability, and for the majority of the test cases ARP/wARP correctly builds 90% of main-chain fragments in regions where the local resolution is 4.0 Å or better. It is also demonstrated that the coordinate accuracy for models built into cryo-EM maps is comparable to that of X-ray crystallographic models at similar local cryo-EM and crystallographic resolutions. The model accuracy also correlates with the refined atomic displacement parameters.

Present and Emerging Methodologies in Cryo-EM Single-Particle Analysis

Over the past decade, technical and methodological improvements in cryogenic electron microscopy (cryo-EM) single-particle analysis have enabled routine high-resolution structural analyses of biological macromolecules, resulting in a flood of new molecular insights into protracted biological questions. However, despite the tremendous progress and success of the field in recent years, opportunities for improvement remain in various aspects of the cryo-EM single-particle analysis workflow (e.g., sample preparation, image acquisition and processing, and structure validation). Here, we review recent advances that have contributed to the principal methods in cryo-EM and identify persisting challenges and bottlenecks that will require further methodological and hardware development.

Density modification of cryo-EM maps

Density modification uses expectations about features of a map such as a flat solvent and expected distributions of density in the region of the macromolecule to improve individual Fourier terms representing the map. This process transfers information from one part of a map to another and can improve the accuracy of a map. Here, the assumptions behind density modification for maps from electron cryomicroscopy are examined and a procedure is presented that allows the incorporation of model-based information. Density modification works best in cases where unfiltered, unmasked maps with clear boundaries between the macromolecule and solvent are visible, and where there is substantial noise in the map, both in the region of the macromolecule and the solvent. It also is most effective if the characteristics of the map are relatively constant within regions of the macromolecule and the solvent. Model-based information can be used to improve density modification, but model bias can in principle occur. Here, model bias is reduced by using ensemble models that allow an estimation of model uncertainty. A test of model bias is presented that suggests that even if the expected density in a region of a map is specified incorrectly by using an incorrect model, the incorrect expectations do not strongly affect the final map.

Advances in methods for atomic resolution macromolecular structure determination

Recent technical advances have dramatically increased the power and scope of structural biology. New developments in high-resolution cryo-electron microscopy, serial X-ray crystallography, and electron diffraction have been especially transformative. Here we highlight some of the latest advances and current challenges at the frontiers of atomic resolution methods for elucidating the structures and dynamical properties of macromolecules and their complexes.

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.

Structures Illuminate Cardiac Ion Channel Functions in Health and in Long QT Syndrome

The cardiac action potential is critical to the production of a synchronized heartbeat. This electrical impulse is governed by the intricate activity of cardiac ion channels, among them the cardiac voltage-gated potassium (Kv) channels KCNQ1 and hERG as well as the voltage-gated sodium (Nav) channel encoded by SCN5A. Each channel performs a highly distinct function, despite sharing a common topology and structural components. These three channels are also the primary proteins mutated in congenital long QT syndrome (LQTS), a genetic condition that predisposes to cardiac arrhythmia and sudden cardiac death due to impaired repolarization of the action potential and has a particular proclivity for reentrant ventricular arrhythmias. Recent cryo-electron microscopy structures of human KCNQ1 and hERG, along with the rat homolog of SCN5A and other mammalian sodium channels, provide atomic-level insight into the structure and function of these proteins that advance our understanding of their distinct functions in the cardiac action potential, as well as the molecular basis of LQTS. In this review, the gating, regulation, LQTS mechanisms, and pharmacological properties of KCNQ1, hERG, and SCN5A are discussed in light of these recent structural findings.

Structural basis for distinct operational modes and protease activation in AAA+ protease Lon

Substrate-bound structures of AAA+ protein translocases reveal a conserved asymmetric spiral staircase architecture wherein a sequential ATP hydrolysis cycle drives hand-over-hand substrate translocation. However, this configuration is unlikely to represent the full conformational landscape of these enzymes, as biochemical studies suggest distinct conformational states depending on the presence or absence of substrate. Here, we used cryo-electron microscopy to determine structures of the Yersinia pestis Lon AAA+ protease in the absence and presence of substrate, uncovering the mechanistic basis for two distinct operational modes. In the absence of substrate, Lon adopts a left-handed, "open" spiral organization with autoinhibited proteolytic active sites. Upon the addition of substrate, Lon undergoes a reorganization to assemble an enzymatically active, right-handed "closed" conformer with active protease sites. These findings define the mechanistic principles underlying the operational plasticity required for processing diverse protein substrates.

What Could Go Wrong? A Practical Guide to Single-Particle Cryo-EM: From Biochemistry to Atomic Models

Cryo-electron microscopy (cryo-EM) has enjoyed explosive recent growth due to revolutionary advances in hardware and software, resulting in a steady stream of long-awaited, high-resolution structures with unprecedented atomic detail. With this comes an increased number of microscopes, cryo-EM facilities, and scientists eager to leverage the ability to determine protein structures without crystallization. However, numerous pitfalls and considerations beset the path toward high-resolution structures and are not necessarily obvious from literature surveys. Here, we detail the most common misconceptions when initiating a cryo-EM project and common technical hurdles, as well as their solutions, and we conclude with a vision for the future of this exciting field.

Cryo-EM map interpretation and protein model-building using iterative map segmentation

A procedure for building protein chains into maps produced by single-particle electron cryo-microscopy (cryo-EM) is described. The procedure is similar to the way an experienced structural biologist might analyze a map, focusing first on secondary structure elements such as helices and sheets, then varying the contour level to identify connections between these elements. Since the high density in a map typically follows the main-chain of the protein, the main-chain connection between secondary structure elements can often be identified as the unbranched path between them with the highest minimum value along the path. This chain-tracing procedure is then combined with finding side-chain positions based on the presence of density extending away from the main path of the chain, allowing generation of a Cα model. The Cα model is converted to an all-atom model and is refined against the map. We show that this procedure is as effective as other existing methods for interpretation of cryo-EM maps and that it is considerably faster and produces models with fewer chain breaks than our previous methods that were based on approaches developed for crystallographic maps.

Evolving data standards for cryo-EM structures

Electron cryo-microscopy (cryo-EM) is increasingly being used to determine 3D structures of a broad spectrum of biological specimens from molecules to cells. Anticipating this progress in the early 2000s, an international collaboration of scientists with expertise in both cryo-EM and structure data archiving was established (EMDataResource, previously known as EMDataBank). The major goals of the collaboration have been twofold: to develop the necessary infrastructure for archiving cryo-EM-derived density maps and models, and to promote development of cryo-EM structure validation standards. We describe how cryo-EM data archiving and validation have been developed and jointly coordinated for the Electron Microscopy Data Bank and Protein Data Bank archives over the past two decades, as well as the impact of evolving technology on data standards. Just as for X-ray crystallography and nuclear magnetic resonance, engaging the scientific community via workshops and challenging activities has played a central role in developing recommendations and requirements for the cryo-EM structure data archives.

Modelling structures in cryo-EM maps

Recent advances in structure determination of sub-cellular structures using cryo-electron microscopy and tomography have enabled us to understand their architecture in a more detailed manner and gain insight into their function. The choice of approach to use for atomic model building, fitting, refinement and validation in the 3D map resulting from these experiments depends primarily on the resolution of the map and the prior information on the corresponding model. Here, we survey some of such methods and approaches and highlight their uses in specific recent examples.

Unique Structural Features of the Mitochondrial AAA+ Protease AFG3L2 Reveal the Molecular Basis for Activity in Health and Disease

Mitochondrial AAA+ quality-control proteases regulate diverse aspects of mitochondrial biology through specialized protein degradation, but the underlying mechanisms of these enzymes remain poorly defined. The mitochondrial AAA+ protease AFG3L2 is of particular interest, as genetic mutations localized throughout AFG3L2 are linked to diverse neurodegenerative disorders. However, a lack of structural data has limited our understanding of how mutations impact enzymatic function. Here, we used cryoelectron microscopy (cryo-EM) to determine a substrate-bound structure of the catalytic core of human AFG3L2. This structure identifies multiple specialized structural features that integrate with conserved motifs required for ATP-dependent translocation to unfold and degrade targeted proteins. Many disease-relevant mutations localize to these unique structural features of AFG3L2 and distinctly influence its activity and stability. Our results provide a molecular basis for neurological phenotypes associated with different AFG3L2 mutations and establish a structural framework to understand how different members of the AAA+ superfamily achieve specialized biological functions.

An allosteric network in spastin couples multiple activities required for microtubule severing

The AAA+ ATPase spastin remodels microtubule arrays through severing and its mutation is the most common cause of hereditary spastic paraplegias (HSP). Polyglutamylation of the tubulin C-terminal tail recruits spastin to microtubules and modulates severing activity. Here, we present a ~3.2 Å resolution cryo-EM structure of the Drosophila melanogaster spastin hexamer with a polyglutamate peptide bound in its central pore. Two electropositive loops arranged in a double-helical staircase coordinate the substrate sidechains. The structure reveals how concurrent nucleotide and substrate binding organizes the conserved spastin pore loops into an ordered network that is allosterically coupled to oligomerization, and suggests how tubulin tail engagement activates spastin for microtubule disassembly. This allosteric coupling may apply generally in organizing AAA+ protein translocases into their active conformations. We show that this allosteric network is essential for severing and is a hotspot for HSP mutations.

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.

Cryo-Electron Microscopy Methodology: Current Aspects and Future Directions

Cryo-electron microscopy (cryo-EM) has emerged as a powerful structure determination technique. Its most prolific branch is single particle analysis (SPA), a method being used in a growing number of laboratories worldwide to determine high-resolution protein structures. Cryo-electron tomography (cryo-ET) is another powerful approach that enables visualization of protein complexes in their native cellular environment. Despite the wide-ranging success of cryo-EM, there are many methodological aspects that could be improved. Those include sample preparation, sample screening, data acquisition, image processing, and structure validation. Future developments will increase the reliability and throughput of the technique and reduce the cost and skill level barrier for its adoption.

Structure Reveals a Mechanism of CRISPR-RNA-Guided Nuclease Recruitment and Anti-CRISPR Viral Mimicry

Bacteria and archaea have evolved sophisticated adaptive immune systems that rely on CRISPR RNA (crRNA)-guided detection and nuclease-mediated elimination of invading nucleic acids. Here, we present the cryo-electron microscopy (cryo-EM) structure of the type I-F crRNA-guided surveillance complex (Csy complex) from Pseudomonas aeruginosa bound to a double-stranded DNA target. Comparison of this structure to previously determined structures of this complex reveals a ∼180-degree rotation of the C-terminal helical bundle on the "large" Cas8f subunit. We show that the double-stranded DNA (dsDNA)-induced conformational change in Cas8f exposes a Cas2/3 "nuclease recruitment helix" that is structurally homologous to a virally encoded anti-CRISPR protein (AcrIF3). Structural homology between Cas8f and AcrIF3 suggests that AcrIF3 is a mimic of the Cas8f nuclease recruitment helix.

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.

Automated cryo-EM structure refinement using correlation-driven molecular dynamics

We present a correlation-driven molecular dynamics (CDMD) method for automated refinement of atomistic models into cryo-electron microscopy (cryo-EM) maps at resolutions ranging from near-atomic to subnanometer. It utilizes a chemically accurate force field and thermodynamic sampling to improve the real-space correlation between the modeled structure and the cryo-EM map. Our framework employs a gradual increase in resolution and map-model agreement as well as simulated annealing, and allows fully automated refinement without manual intervention or any additional rotamer- and backbone-specific restraints. Using multiple challenging systems covering a wide range of map resolutions, system sizes, starting model geometries and distances from the target state, we assess the quality of generated models in terms of both model accuracy and potential of overfitting. To provide an objective comparison, we apply several well-established methods across all examples and demonstrate that CDMD performs best in most cases.

Exploring applications of crowdsourcing to cryo-EM

Extraction of particles from cryo-electron microscopy (cryo-EM) micrographs is a crucial step in processing single-particle datasets. Although algorithms have been developed for automatic particle picking, these algorithms generally rely on two-dimensional templates for particle identification, which may exhibit biases that can propagate artifacts through the reconstruction pipeline. Manual picking is viewed as a gold-standard solution for particle selection, but it is too time-consuming to perform on data sets of thousands of images. In recent years, crowdsourcing has proven effective at leveraging the open web to manually curate datasets. In particular, citizen science projects such as Galaxy Zoo have shown the power of appealing to users’ scientific interests to process enormous amounts of data. To this end, we explored the possible applications of crowdsourcing in cryo-EM particle picking, presenting a variety of novel experiments including the production of a fully annotated particle set from untrained citizen scientists. We show the possibilities and limitations of crowdsourcing particle selection tasks, and explore further options for crowdsourcing cryo-EM data processing.

Achieving better-than-3-Å resolution by single-particle cryo-EM at 200 keV

Nearly 98% of single particle cryo-EM structures resolved to better than 4 Å resolution have been determined using 300 keV transmission electron microscopes. We demonstrate that it is possible to obtain reconstructions of macromolecular complexes at a range of sizes to better than 3 Å resolution using a 200 keV transmission electron microscope. These structures are of sufficient quality to unambiguously assign amino acid rotameric conformations and identify ordered water molecules.