Single-particle cryo-EM using alignment by classification (ABC): the structure of Lumbricus terrestris haemoglobin

Data-driven determination of number of discrete conformations in single-particle cryo-EM

BACKGROUND AND OBJECTIVE

One of the strengths of single-particle cryo-EM compared to other structural determination techniques is its ability to image heterogeneous samples containing multiple molecular species, different oligomeric states or distinct conformations. This is achieved using routines for in-silico 3D classification that are now well established in the field and have successfully been used to characterize the structural heterogeneity of important biomolecules. These techniques, however, rely on expert-user knowledge and trial-and-error experimentation to determine the correct number of conformations, making it a labor intensive, subjective, and difficult to reproduce procedure.

METHODS

We propose an approach to address the problem of automatically determining the number of discrete conformations present in heterogeneous single-particle cryo-EM datasets. We do this by systematically evaluating all possible partitions of the data and selecting the result that maximizes the average variance of similarities measured between particle images and the corresponding 3D reconstructions.

RESULTS

Using this strategy, we successfully analyzed datasets of heterogeneous protein complexes, including: 1) in-silico mixtures obtained by combining closely related antibody-bound HIV-1 Env trimers and other important membrane channels, and 2) naturally occurring mixtures from diverse and dynamic protein complexes representing varying degrees of structural heterogeneity and conformational plasticity.

CONCLUSIONS

The availability of unsupervised strategies for 3D classification combined with existing approaches for fully automatic pre-processing and 3D refinement, represents an important step towards converting single-particle cryo-EM into a high-throughput technique.

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

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

Cryo-EM structure of the mature and infective Mayaro virus at 4.4 Å resolution reveals features of arthritogenic alphaviruses

Mayaro virus (MAYV) is an emerging arbovirus of the Americas that may cause a debilitating arthritogenic disease. The biology of MAYV is not fully understood and largely inferred from related arthritogenic alphaviruses. Here, we present the structure of MAYV at 4.4 Å resolution, obtained from a preparation of mature, infective virions. MAYV presents typical alphavirus features and organization. Interactions between viral proteins that lead to particle formation are described together with a hydrophobic pocket formed between E1 and E2 spike proteins and conformational epitopes specific of MAYV. We also describe MAYV glycosylation residues in E1 and E2 that may affect MXRA8 host receptor binding, and a molecular “handshake” between MAYV spikes formed by N262 glycosylation in adjacent E2 proteins. The structure of MAYV is suggestive of structural and functional complexity among alphaviruses, which may be targeted for specificity or antiviral activity.

Mayaro virus (MAYV) is an emerging arbovirus in Central and South America that is transmitted by mosquitoes and causes arthritogenic disease. Here, the authors present the 4.4 Å resolution cryo-EM structure of MAYV and describe specific features of the virus, which could be exploited for the design of MAYV-specific diagnostics and therapeutics.

Improvement of cryo-EM maps by density modification

A density-modification procedure for improving maps from single-particle electron cryogenic microscopy (cryo-EM) is presented. The theoretical basis of the method is identical to that of maximum-likelihood density modification, previously used to improve maps from macromolecular X-ray crystallography. Key differences from applications in crystallography are that the errors in Fourier coefficients are largely in the phases in crystallography but in both phases and amplitudes in cryo-EM, and that half-maps with independent errors are available in cryo-EM. These differences lead to a distinct approach for combination of information from starting maps with information obtained in the density-modification process. The density-modification procedure was applied to a set of 104 datasets and improved map-model correlation and increased the visibility of details in many of the maps. The procedure requires two unmasked half-maps and a sequence file or other source of information on the volume of the macromolecule that has been imaged.

Validation tests for cryo-EM maps using an independent particle set

Cryo-electron microscopy (cryo-EM) has revolutionized structural biology by providing 3D density maps of biomolecules at near-atomic resolution. However, map validation is still an open issue. Despite several efforts from the community, it is possible to overfit 3D maps to noisy data. Here, we develop a novel methodology that uses a small independent particle set (not used during the 3D refinement) to validate the maps. The main idea is to monitor how the map probability evolves over the control set during the 3D refinement. The method is complementary to the gold-standard procedure, which generates two reconstructions at each iteration. We low-pass filter the two reconstructions for different frequency cutoffs, and we calculate the probability of each filtered map given the control set. For high-quality maps, the probability should increase as a function of the frequency cutoff and the refinement iteration. We also compute the similarity between the densities of probability distributions of the two reconstructions. As higher frequencies are included, the distributions become more dissimilar. We optimized the BioEM package to perform these calculations, and tested it over systems ranging from quality data to pure noise. Our results show that with our methodology, it possible to discriminate datasets that are constructed from noise particles. We conclude that validation against a control particle set provides a powerful tool to assess the quality of cryo-EM maps.

Cryo-EM structures from sub-nl volumes using pin-printing and jet vitrification

The increasing demand for cryo-electron microscopy (cryo-EM) reveals drawbacks in current sample preparation protocols, such as sample waste and lack of reproducibility. Here, we present several technical developments that provide efficient sample preparation for cryo-EM studies. Pin printing substantially reduces sample waste by depositing only a sub-nanoliter volume of sample on the carrier surface. Sample evaporation is mitigated by dewpoint control feedback loops. The deposited sample is vitrified by jets of cryogen followed by submersion into a cryogen bath. Because the cryogen jets cool the sample from the center, premounted autogrids can be used and loaded directly into automated cryo-EMs. We integrated these steps into a single device, named VitroJet. The device’s performance was validated by resolving four standard proteins (apoferritin, GroEL, worm hemoglobin, beta-galactosidase) to ~3 Å resolution using a 200-kV electron microscope. The VitroJet offers a promising solution for improved automated sample preparation in cryo-EM studies.

There is a need to further improve the automation of cryo-EM sample preparation to make it more easily accessible for non-specialists, reduce sample waste and increase reproducibility. Here, the authors present VitroJet, a single device, where sub-nl volumes of samples are deposited by pin printing thus eliminating the need for sample blotting, which is followed by jet vitrification, and they show that high-resolution structures can be obtained using four standard proteins.

Myriapod haemocyanin: the first three-dimensional reconstruction of Scolopendra subspinipes and preliminary structural analysis of S. viridicornis

Haemocyanins (Hcs) are copper-containing, respiratory proteins that occur in the haemolymph of many arthropod species. Here, we report the presence of Hcs in the chilopode Myriapoda, demonstrating that these proteins are more widespread among the Arthropoda than previously thought. The analysis of transcriptome of S. subspinipes subpinipes reveals the presence of two distinct subunits of Hc, where the signal peptide is present, and six of prophenoloxidase (PPO), where the signal peptide is absent, in the 75 kDa range. Size exclusion chromatography profiles indicate different quaternary organization for Hc of both species, which was corroborated by TEM analysis: S. viridicornis Hc is a 6 × 6-mer and S. subspinipes Hc is a 3 × 6-mer, which resembles the half-structure of the 6 × 6-mer but also includes the presence of phenoloxidases, since the 1 × 6-mer quaternary organization is commonly associated with hexamers of PPO. Studies with Chelicerata showed that PPO activity are exclusively associated with the Hcs. This study indicates that Scolopendra may have different proteins playing oxygen transport (Hc) and PO function, both following the hexameric oligomerization observed in Hcs.

Need for Cross-Validation of Single Particle Cryo-EM

Cross-validation is used to determine the validity of a model on unseen data by assessing if the model is overfitted to noise. It is widely used in many fields, from artificial intelligence to structural biology in X-ray crystallography and nuclear magnetic resonance. Although there are concerns of map overfitting in cryo-electron microscopy (cryo-EM), cross-validation is rarely used. The problem is that establishing a performance metric of the maps over unseen data (given by 2D-projection images) is difficult due to the low signal-to-noise ratios in the individual particles. Here, I present recent advances for cryo-EM map reconstruction. I highlight that the gold-standard procedure can fail to detect map overfitting in certain cases, showing the necessity of assessing the map quality on unbiased data. Finally, I describe the challenges and advantages of developing a robust cross-validation methodology for cryo-EM.

A revised order of subunits in mammalian septin complexes

Septins are GTP binding proteins considered to be novel components of the cytoskeleton. They polymerize into filaments based on hexameric or octameric core particles in which two copies of either three or four different septins, respectively, assemble into a specific sequence. Viable combinations of the 13 human septins are believed to obey substitution rules in which the different septins involved must come from distinct subgroups. The hexameric assembly, for example, has been reported to be SEPT7-SEPT6-SEPT2-SEPT2-SEPT6-SEPT7. Here, we have replaced SEPT2 by SEPT5 according to the substitution rules and used transmission electron microscopy to demonstrate that the resulting recombinant complex assembles into hexameric particles which are inverted with respect that predicted previously. MBP-SEPT5 constructs and immunostaining show that SEPT5 occupies the terminal positions of the hexamer. We further show that this is also true for the assembly including SEPT2, in direct contradiction with that reported previously. Consequently, both complexes expose an NC interface, as reported for yeast, which we show to be more susceptible to high salt concentrations. The correct assembly for the canonical combination of septins 2-6-7 is therefore established to be SEPT2-SEPT6-SEPT7-SEPT7-SEPT6-SEPT2, implying the need for revision of the mechanisms involved in filament assembly.

Cryo_fit: Democratization of flexible fitting for cryo-EM

Cryo-electron microscopy (cryo-EM) is becoming a method of choice for describing native conformations of biomolecular complexes at high resolution. The rapid growth of cryo-EM in recent years has created a high demand for automated solutions, both in hardware and software. Flexible fitting of atomic models to three-dimensional (3D) cryo-EM reconstructions by molecular dynamics (MD) simulation is a popular technique but often requires technical expertise in computer simulation. This work introduces cryo_fit, a package for the automatic flexible fitting of atomic models in cryo-EM maps using MD simulation. The package is integrated with the Phenix software suite. The module was designed to automate the multiple steps of MD simulation in a reproducible manner, as well as facilitate refinement and validation through Phenix. Through the use of cryo_fit, scientists with little experience in MD simulation can produce high quality atomic models automatically and better exploit the potential of cryo-EM.

Structure Determination by Single-Particle Cryo-Electron Microscopy: Only the Sky (and Intrinsic Disorder) is the Limit

Traditionally, X-ray crystallography and NMR spectroscopy represent major workhorses of structural biologists, with the lion share of protein structures reported in protein data bank (PDB) being generated by these powerful techniques. Despite their wide utilization in protein structure determination, these two techniques have logical limitations, with X-ray crystallography being unsuitable for the analysis of highly dynamic structures and with NMR spectroscopy being restricted to the analysis of relatively small proteins. In recent years, we have witnessed an explosive development of the techniques based on Cryo-electron microscopy (Cryo-EM) for structural characterization of biological molecules. In fact, single-particle Cryo-EM is a special niche as it is a technique of choice for the structural analysis of large, structurally heterogeneous, and dynamic complexes. Here, sub-nanometer atomic resolution can be achieved (i.e., resolution below 10 Å) via single-particle imaging of non-crystalline specimens, with accurate 3D reconstruction being generated based on the computational averaging of multiple 2D projection images of the same particle that was frozen rapidly in solution. We provide here a brief overview of single-particle Cryo-EM and show how Cryo-EM has revolutionized structural investigations of membrane proteins. We also show that the presence of intrinsically disordered or flexible regions in a target protein represents one of the major limitations of this promising technique.

Rational modulator design by exploitation of protein-protein complex structures

The horizon of drug discovery is currently expanding to target and modulate protein-protein interactions (PPIs) in globular proteins and intrinsically disordered proteins that are involved in various diseases. To either interrupt or stabilize PPIs, the 3D structure of target protein-protein (or protein-peptide) complexes can be exploited to rationally design PPI modulators (inhibitors or stabilizers) through structure-based molecular design. In this review, we present an overview of experimental and computational methods that can be used to determine 3D structures of protein-protein complexes. Several approaches including rational and in silico methods that can be applied to design peptides, peptidomimetics and small compounds by utilization of determined 3D protein-protein/peptide complexes are summarized and illustrated.

Accurate flexible refinement of atomic models against medium-resolution cryo-EM maps using damped dynamics

BACKGROUND

Dramatic progress has recently been made in cryo-electron microscopy technologies, which now make possible the reconstruction of a growing number of biomolecular structures to near-atomic resolution. However, the need persists for fitting and refinement approaches that address those cases that require modeling assistance.

METHODS

In this paper, we describe algorithms to optimize the performance of such medium-resolution refinement methods. These algorithms aim to automatically optimize the parameters that define the density shape of the flexibly fitted model, as well as the time-dependent damper cutoff distance. Atomic distance constraints can be prescribed for cases where extra containment of parts of the structure is helpful, such as in regions where the density map is poorly defined. Also, we propose a simple stopping criterion that estimates the probable onset of overfitting during the simulation.

RESULTS

The new set of algorithms produce more accurate fitting and refinement results, and yield a faster rate of convergence of the trajectory toward the fitted conformation. The latter is also more reliable due to the overfitting warning provided to the user.

CONCLUSIONS

The algorithms described here were implemented in the new Damped-Dynamics Flexible Fitting simulation tool "DDforge" in the Situs package.

New tools for the analysis and validation of cryo-EM maps and atomic models

Recent advances in the field of electron cryomicroscopy (cryo-EM) have resulted in a rapidly increasing number of atomic models of biomacromolecules that have been solved using this technique and deposited in the Protein Data Bank and the Electron Microscopy Data Bank. Similar to macromolecular crystallography, validation tools for these models and maps are required. While some of these validation tools may be borrowed from crystallography, new methods specifically designed for cryo-EM validation are required. Here, new computational methods and tools implemented in PHENIX are discussed, including d99 to estimate resolution, phenix.auto_sharpen to improve maps and phenix.mtriage to analyze cryo-EM maps. It is suggested that cryo-EM half-maps and masks should be deposited to facilitate the evaluation and validation of cryo-EM-derived atomic models and maps. The application of these tools to deposited cryo-EM atomic models and maps is also presented.

Mechanism for remodelling of the cell cycle checkpoint protein MAD2 by the ATPase TRIP13

The maintenance of genome stability during mitosis is coordinated by the spindle assembly checkpoint (SAC) through its effector the mitotic checkpoint complex (MCC), an inhibitor of the anaphase-promoting complex (APC/C, also known as the cyclosome)1,2. Unattached kinetochores control MCC assembly by catalysing a change in the topology of the β-sheet of MAD2 (an MCC subunit), thereby generating the active closed MAD2 (C-MAD2) conformer3-5. Disassembly of free MCC, which is required for SAC inactivation and chromosome segregation, is an ATP-dependent process driven by the AAA+ ATPase TRIP13. In combination with p31comet, an SAC antagonist6, TRIP13 remodels C-MAD2 into inactive open MAD2 (O-MAD2)7-10. Here, we present a mechanism that explains how TRIP13-p31comet disassembles the MCC. Cryo-electron microscopy structures of the TRIP13-p31comet-C-MAD2-CDC20 complex reveal that p31comet recruits C-MAD2 to a defined site on the TRIP13 hexameric ring, positioning the N terminus of C-MAD2 (MAD2NT) to insert into the axial pore of TRIP13 and distorting the TRIP13 ring to initiate remodelling. Molecular modelling suggests that by gripping MAD2NT within its axial pore, TRIP13 couples sequential ATP-driven translocation of its hexameric ring along MAD2NT to push upwards on, and simultaneously rotate, the globular domains of the p31comet-C-MAD2 complex. This unwinds a region of the αA helix of C-MAD2 that is required to stabilize the C-MAD2 β-sheet, thus destabilizing C-MAD2 in favour of O-MAD2 and dissociating MAD2 from p31comet. Our study provides insights into how specific substrates are recruited to AAA+ ATPases through adaptor proteins and suggests a model of how translocation through the axial pore of AAA+ ATPases is coupled to protein remodelling.