Cryo-EM Model Research Articles

Yeast ribosome biogenesis factors Puf6 and Nog2 and ribosomal proteins uL2 and eL43 act in concert to facilitate the release of nascent large ribosomal subunits from the nucleolus.

Large ribosomal subunit precursors (pre-LSUs) are primarily synthesized in the nucleolus. At an undetermined step in their assembly, they are released into the nucleoplasm. Structural models of yeast pre-LSUs at various stages of assembly have been collected using cryo-EM. However, which cryo-EM model is closest to the final nucleolar intermediate of the LSU has yet to be determined. To elucidate the mechanisms of the release of pre-LSUs from the nucleolus, we assayed effects of depleting or knocking out two yeast ribosome biogenesis factors (RiBi factors), Puf6 and Nog2, and two ribosomal proteins, uL2 and eL43. These proteins function during or stabilize onto pre-LSUs between the late nucleolar stages to early nucleoplasmic stages of ribosome biogenesis. By characterizing the phenotype of these four mutants, we determined that a particle that is intermediate between the cryo-EM model State NE1 and State NE2 likely represents the final nucleolar assembly intermediate of the LSU. We conclude that the release of the RiBi factors Nip7, Nop2 and Spb1 and the subsequent stabilization of rRNA domains IV and V may be key triggers for the release of pre-LSUs from the nucleolus.

RIBFIND2: Identifying rigid bodies in protein and nucleic acid structures.

Molecular structures are often fitted into cryo-EM maps by flexible fitting. When this requires large conformational changes, identifying rigid bodies can help optimize the model-map fit. Tools for identifying rigid bodies in protein structures exist, however an equivalent for nucleic acid structures is lacking. With the increase in cryo-EM maps containing RNA and progress in RNA structure prediction, there is a need for such tools. We previously developed RIBFIND, a program for clustering protein secondary structures into rigid bodies. In RIBFIND2, this approach is extended to nucleic acid structures. RIBFIND2 can identify biologically relevant rigid bodies in important groups of complex RNA structures, capturing a wide range of dynamics, including large rigid-body movements. The usefulness of RIBFIND2-assigned rigid bodies in cryo-EM model refinement was demonstrated on three examples, with two conformations each: Group II Intron complexed IEP, Internal Ribosome Entry Site and the Processome, using cryo-EM maps at 2.7-5 Å resolution. A hierarchical refinement approach, performed on progressively smaller sets of RIBFIND2 rigid bodies, was clearly shown to have an advantage over classical all-atom refinement. RIBFIND2 is available via a web server with structure visualization and as a standalone tool.

Particle fusion of super-resolution data reveals the unit structure of Nup96 in Nuclear Pore Complex

Single molecule localization microscopy offers resolution nearly down to the molecular level with specific molecular labelling, and is thereby a promising tool for structural biology. In practice, however, the actual value to this field is limited primarily by incomplete fluorescent labelling of the structure. This missing information can be completed by merging information from many structurally identical particles in a particle fusion approach similar to cryo-EM single-particle analysis. In this paper, we present a data analysis of particle fusion results of fluorescently labelled Nup96 nucleoporins in the Nuclear Pore Complex to show that Nup96 occurs in a spatial arrangement of two rings of 8 units with two Nup96 copies per unit giving a total of 32 Nup96 copies per pore. We use Artificial Intelligence assisted modeling in Alphafold to extend the existing cryo-EM model of Nup96 to accurately pinpoint the positions of the fluorescent labels and show the accuracy of the match between fluorescent and cryo-EM data to be better than 3 nm in-plane and 5 nm out-of-plane.

Capsid structure of a fungal dsRNA megabirnavirus reveals its previously unidentified surface architecture.

Rosellinia necatrix megabirnavirus 1-W779 (RnMBV1) is a non-enveloped icosahedral double-stranded (ds)RNA virus that infects the ascomycete fungus Rosellinia necatrix, a causative agent that induces a lethal plant disease white root rot. Herein, we have first resolved the atomic structure of the RnMBV1 capsid at 3.2 Å resolution using cryo-electron microscopy (cryo-EM) single-particle analysis. Compared with other non-enveloped icosahedral dsRNA viruses, the RnMBV1 capsid protein structure exhibits an extra-long C-terminal arm and a surface protrusion domain. In addition, the previously unrecognized crown proteins are identified in a symmetry-expanded cryo-EM model and are present over the 3-fold axes. These exclusive structural features of the RnMBV1 capsid could have been acquired for playing essential roles in transmission and/or particle assembly of the megabirnaviruses. Our findings, therefore, will reinforce the understanding of how the structural and molecular machineries of the megabirnaviruses influence the virulence of the disease-related ascomycete fungus.

Ring-stacked capsids of white spot syndrome virus and structural transitions with genome ejection.

White spot syndrome virus (WSSV) is one of the largest DNA viruses and the major pathogen responsible for white spot syndrome in crustaceans. The WSSV capsid is critical for genome encapsulation and ejection and exhibits the rod-shaped and oval-shaped structures during the viral life cycle. However, the detailed architecture of the capsid and the structural transition mechanism remain unclear. Here, using cryo-electron microscopy (cryo-EM), we obtained a cryo-EM model of the rod-shaped WSSV capsid and were able to characterize its ring-stacked assembly mechanism. Furthermore, we identified an oval-shaped WSSV capsid from intact WSSV virions and analyzed the structural transition mechanism from the oval-shaped to rod-shaped capsids induced by high salinity. These transitions, which decrease internal capsid pressure, always accompany DNA release and mostly eliminate the infection of the host cells. Our results demonstrate an unusual assembly mechanism of the WSSV capsid and offer structural insights into the pressure-driven genome release.

Structural studies of human serum albumin using cryo-EM up to 0.38 nm resolution.

Single particle cryo-electron microscopy (cryoEM) offers for large proteins the resolution close to the crystallographic resolution. Unfortunately, the current limitation for cryoEM is the molecular weight, and the smallest proteins, tested individually by the cryoEM method, rarely are less than 100 kDa in size. The aim of our research was to study the structure of human serum albumin and its oligomeric forms using a combination of cryoEM, small angle X-ray scattering (SAXS) and dynamic light scattering (DLS). HSA is a monomeric (66.4 kDa) protein, which for years has been a model system for the transport of drugs and other substances, because in the human body it binds and carries lipids or metal cations (including zinc and copper ions) as well as for basic biophysical research. Currently, over 130 HSA structures in complexes with various ligands and alone at the resolution ranging from 0.189 nm to 0.39 nm have been deposited in Protein Data Bank. Using cryoEM, we were able to obtain the HSA structure with a 0.38 nm resolution. At present, it is probably one of the smallest protein structures determined using cryoEM at near-crystallographic resolution. Using SAXS and DLS we analyzed the HSA monomer structure and its oligomers in solution and compared with the cryoEM model. This research was funded by National Science Centre (Poland), grant number 2017/27/B/ST4/00485.

Fast and automated protein-DNA/RNA macromolecular complex modeling from cryo-EM maps.

Cryo-electron microscopy (cryo-EM) allows a macromolecular structure such as protein-DNA/RNA complexes to be reconstructed in a three-dimensional coulomb potential map. The structural information of these macromolecular complexes forms the foundation for understanding the molecular mechanism including many human diseases. However, the model building of large macromolecular complexes is often difficult and time-consuming. We recently developed DeepTracer-2.0, an artificial-intelligence-based pipeline that can build amino acid and nucleic acid backbones from a single cryo-EM map, and even predict the best-fitting residues according to the density of side chains. The experiments showed improved accuracy and efficiency when benchmarking the performance on independent experimental maps of protein-DNA/RNA complexes and demonstrated the promising future of macromolecular modeling from cryo-EM maps. Our method and pipeline could benefit researchers worldwide who work in molecular biomedicine and drug discovery, and substantially increase the throughput of the cryo-EM model building. The pipeline has been integrated into the web portal https://deeptracer.uw.edu/.

Direct calculation of cryo-EM and crystallographic model maps for real-space refinement.

This work addresses the problem of the calculation of limited-resolution maps from an atomic model in cryo-electron microscopy and in X-ray and neutron crystallography, including cases where the resolution varies from one molecular region to another. Such maps are necessary in real-space refinement for comparison with the experimental maps. For an appropriate numeric comparison, the calculated maps should reproduce not only the structural features contained in the experimental maps but also the principal map distortions. These model maps can be obtained with no use of Fourier transforms but, similar to density distributions, as a sum of individual atomic contributions. Such contributions, referred to as atomic density images, are atomic densities morphed to reflect distortions of the experimental map, in particular the loss of resolution. They are described by functions composed of a central peak surrounded by Fourier ripples. For practical calculations, atomic images should be cut at some distance. It is shown that to reach a reasonable accuracy such a distance should be significantly larger than the distance customarily applied when calculating density distributions. This is a consequence of the slow rate with which the amplitude of the Fourier ripples decreases. Such a large distance means that at least a few ripples should be included in calculations in order to obtain a map that is sufficiently accurate. Oscillating functions describing these atomic contributions depend, for a given atomic type, on the resolution and on the atomic displacement parameter values. To express both the central peak and the Fourier ripples of the atomic images, these functions are represented by the sums of especially designed terms, each concentrated in a spherical shell and depending analytically on the atomic parameters. In this work, the strength of the dependence of the accuracy of resulting map on the accuracy of the atomic displacement parameters and on the truncation distance, i.e. the number of ripples included in atomic density images, is analyzed. This analysis is completed by practical aspects of the calculation of maps of inhomogeneous resolution. Tests show that the calculation of limited-resolution maps from an atomic model as a sum of atomic contributions requires a large truncation radius extending beyond the central peak of an atomic image and the first Fourier ripples. The article discusses the practical details of such calculations expressing atomic contributions as analytic functions of the atomic coordinates, the atomic displacement parameters and the local resolution.

Structural basis for recognition of N-formyl peptides as pathogen-associated molecular patterns

The formyl peptide receptor 1 (FPR1) is primarily responsible for detection of short peptides bearing N-formylated methionine (fMet) that are characteristic of protein synthesis in bacteria and mitochondria. As a result, FPR1 is critical to phagocyte migration and activation in bacterial infection, tissue injury and inflammation. How FPR1 distinguishes between formyl peptides and non-formyl peptides remains elusive. Here we report cryo-EM structures of human FPR1-Gi protein complex bound to S. aureus-derived peptide fMet-Ile-Phe-Leu (fMIFL) and E. coli-derived peptide fMet-Leu-Phe (fMLF). Both structures of FPR1 adopt an active conformation and exhibit a binding pocket containing the R2015.38XXXR2055.42 (RGIIR) motif for formyl group interaction and receptor activation. This motif works together with D1063.33 for hydrogen bond formation with the N-formyl group and with fMet, a model supported by MD simulation and functional assays of mutant receptors with key residues for recognition substituted by alanine. The cryo-EM model of agonist-bound FPR1 provides a structural basis for recognition of bacteria-derived chemotactic peptides with potential applications in developing FPR1-targeting agents.

Comparative analysis of human and bovine thyroglobulin structures

In biology, evolutionary conserved protein sequences show homologous physiological phenotypes in their structures and functions. If the protein has a vital function, its sequence is usually conserved across the species. However, in highly conserved protein there still remains small differences across the species. Upon protein–protein interaction (PPI), it is observed that the conserved proteins can have different binding partners that are considered to be caused by the small sequence variations in a specific domain. Thyroglobulin (TG) is the most commonly found protein in the thyroid gland of vertebrates and serves as the precursor of the thyroid hormones, tetraiodothyronine and triiodothyronine that are critical for growth, development and metabolism in vertebrates. In this study, we comparatively analyzed the sequences and structures of the highly conserved regions of TG from two different species in relation to their PPIs. In order to do so, we employed SIM for sequence alignment, STRING for PPI analysis and cryo-electron microscopy for 3D structural analysis. Our Cryo-EM model for TG of Bos taurus determined at 7.1 Å resolution fitted well with the previously published Cryo-EM model for Homo sapiens TG. By demonstrating overall structural homology between TGs from different species, we address that local amino acid sequence variation is sufficient to alter PPIs specific for the organism. We predict that our result will contribute to a deeper understanding in the evolutionary pattern applicable to many other proteins.

Energy landscape of the SARS-CoV-2 reveals extensive conformational heterogeneity.

Cryo-electron microscopy (cryo-EM) has produced a number of structural models of the SARS-CoV-2 spike, already prompting biomedical outcomes. However, these reported models and their associated electrostatic potential maps represent an unknown admixture of conformations stemming from the underlying energy landscape of the spike protein. As with any protein, some of the spike's conformational motions are expected to be biophysically relevant, but cannot be interpreted only by static models. Using experimental cryo-EM images, we present the energy landscape of the glycosylated spike protein, and identify the diversity of low-energy conformations in the vicinity of its open (so called 1RBD-up) state. The resulting atomic refinement reveal global and local molecular rearrangements that cannot be inferred from an average 1RBD-up cryo-EM model. Here we report varied degrees of “openness” in global conformations of the 1RBD-up state, not revealed in the single-model interpretations of the density maps, together with conformations that overlap with the reported models. We discover how the glycan shield contributes to the stability of these low-energy conformations. Five out of six binding sites we analyzed, including those for engaging ACE2, therapeutic mini-proteins, linoleic acid, two different kinds of antibodies, switch conformations between their known apo- and holo-conformations, even when the global spike conformation is 1RBD-up. This apo-to-holo switching is reminiscent of a conformational preequilibrium. We found only one binding site, namely that of AB-C135 remains in apo state within all the sampled free energy-minimizing models, suggesting an induced fit mechanism for the docking of this antibody to the spike.

EMDA: A Python package for Electron Microscopy Data Analysis

An open-source Python library EMDA for cryo-EM map and model manipulation is presented with a specific focus on validation. The use of several functionalities in the library is presented through several examples. The utility of local correlation as a metric for identifying map-model differences and unmodeled regions in maps, and how it is used as a metric of map-model validation is demonstrated. The mapping of local correlation to individual atoms, and its use to draw insights on local signal variations are discussed. EMDA’s likelihood-based map overlay is demonstrated by carrying out a superposition of two domains in two related structures. The overlay is carried out first to bring both maps into the same coordinate frame and then to estimate the relative movement of domains. Finally, the map magnification refinement in EMDA is presented with an example to highlight the importance of adjusting the map magnification in structural comparison studies.

Cryo-EM structure provides insights into the dimer arrangement of the O-linked \u03b2-N-acetylglucosamine transferase OGT

The O-linked β-N-acetylglucosamine modification is a core signalling mechanism, with erroneous patterns leading to cancer and neurodegeneration. Although thousands of proteins are subject to this modification, only a single essential glycosyltransferase catalyses its installation, the O-GlcNAc transferase, OGT. Previous studies have provided truncated structures of OGT through X-ray crystallography, but the full-length protein has never been observed. Here, we report a 5.3 Å cryo-EM model of OGT. We show OGT is a dimer, providing a structural basis for how some X-linked intellectual disability mutations at the interface may contribute to disease. We observe that the catalytic section of OGT abuts a 13.5 tetratricopeptide repeat unit region and find the relative positioning of these sections deviate from the previously proposed, X-ray crystallography-based model. We also note that OGT exhibits considerable heterogeneity in tetratricopeptide repeat units N-terminal to the dimer interface with repercussions for how OGT binds protein ligands and partners.

What is a single cryo-EM image worth? A survey of high-resolution cryo-EM model system datasets

What is a single cryo-EM image worth? A survey of high-resolution cryo-EM model system datasets

Cryo-EM model validation recommendations based on outcomes of the 2019 EMDataResource challenge

This paper describes outcomes of the 2019 Cryo-EM Model Challenge. The goals were to (1) assess the quality of models that can be produced from cryogenic electron microscopy (cryo-EM) maps using current modeling software, (2) evaluate reproducibility of modeling results from different software developers and users and (3) compare performance of current metrics used for model evaluation, particularly Fit-to-Map metrics, with focus on near-atomic resolution. Our findings demonstrate the relatively high accuracy and reproducibility of cryo-EM models derived by 13 participating teams from four benchmark maps, including three forming a resolution series (1.8 to 3.1 Å). The results permit specific recommendations to be made about validating near-atomic cryo-EM structures both in the context of individual experiments and structure data archives such as the Protein Data Bank. We recommend the adoption of multiple scoring parameters to provide full and objective annotation and assessment of the model, reflective of the observed cryo-EM map density.

Structural basis of fast- and slow-severing actin–cofilactin boundaries

Members of the ADF/cofilin family of regulatory proteins bind actin filaments cooperatively, locally change actin subunit conformation and orientation, and sever filaments at “boundaries” between bare and cofilin-occupied segments. A cluster of bound cofilin introduces two distinct classes of boundaries due to the intrinsic polarity of actin filaments, one at the “pointed” end side and the other at the “barbed” end-side of the cluster; severing occurs more readily at the pointed end side of the cluster (“fast-severing” boundary) than the barbed end side (“slow-severing” boundary). A recent electron-cryomicroscopy (cryo-EM) model of the slow-severing boundary revealed structural “defects” at the interface that potentially contribute to severing. However, the structure of the fast-severing boundary remains uncertain. Here, we use extensive molecular dynamics simulations to produce atomic resolution models of both severing boundaries. Our equilibrated simulation model of the slow-severing boundary is consistent with the cryo-EM structural model. Simulations indicate that actin subunits at both boundaries adopt structures intermediate between those of bare and cofilin-bound actin subunits. These “intermediate” states have compromised intersubunit contacts, but those at the slow-severing boundary are stabilized by cofilin bridging interactions, accounting for its lower fragmentation probability. Simulations where cofilin proteins are removed from cofilactin filaments favor a mechanism in which a cluster of two contiguously bound cofilins is needed to fully stabilize the cofilactin conformation, promote cooperative binding interactions, and accelerate filament severing. Together, these studies provide a molecular-scale foundation for developing coarse-grained and theoretical descriptions of cofilin-mediated actin filament severing.

Cryo-EM Structure Determination and Model Fitting of the Proton-Gated Ligand-Gated Ion Channel glic at Multiple pH States

Cryo-EM Structure Determination and Model Fitting of the Proton-Gated Ligand-Gated Ion Channel glic at Multiple pH States

Outcomes of the Cryo-EM Map and Model Challenges

The next step in enabling cryo-EM to continue its growth as a leading method for structure detemination is to create best practices in reconstruction and model building and robust community-approved structure validation methods. We have sponsored two challenges to highlight the need for map and model validation standards, and to expedite development of quantitative tools for cryo-EM structure comparison and assessment with community participation (challenges.emdataresource.org). The scientific community responded enthusiastically with more than 90 scientists worldwide participating. Following a participants' workshop held in October 2017 at Stanford University/SLAC National Accelerator Laboratory, key results and recommendations are now collated in a 2018 special issue of the Journal of Structural Biology. In this presentation we will describe how this round of Cryo-EM Challenges has stimulated important discussions that will help propel this field forward.

Domain Requirements and Genetic Interactions of the Mud1 Subunit of the Saccharomyces cerevisiae U1 snRNP.

Mud1 is an inessential 298-amino acid protein subunit of the Saccharomyces cerevisiae U1 snRNP. Mud1 consists of N-terminal and C-terminal RRM domains (RRM1 and RRM2) separated by a linker domain. Synthetic lethal interactions of mud1∆ with deletions of inessential spliceosome components Nam8, Mud2, and Msl1, or missense mutations in the branchpoint-binding protein Msl5 enabled us to dissect genetically the domain requirements for Mud1 function. We find that the biological activities of Mud1 can be complemented by co-expressing separately the RRM1 (aa 1-127) and linker-RRM2 (aa 128-298) modules. Whereas RRM1 and RRM2 (aa 197-298) per se are inactive in all tests of functional complementation, the linker-RRM2 by itself partially complements a subset of synthetic lethal mud1∆ interactions. Linker segment aa 155 to 196 contains a nuclear localization signal rich in basic amino acids that is necessary for RRM2 activity in mud1∆ complementation. Alanine scanning mutagenesis indicates that none of the individual RRM1 amino acid contacts to U1 snRNA in the cryo-EM model of the yeast U1 snRNP is necessary for mud1∆ complementation activity.

Assessment of detailed conformations suggests strategies for improving cryoEM models: Helix at lower resolution, ensembles, pre-refinement fixups, and validation at multi-residue length scale.

We find that the overall quite good methods used in the CryoEM Model Challenge could still benefit greatly from several strategies for improving local conformations. Our assessments primarily use validation criteria from the MolProbity web service. Those criteria include MolProbity's all-atom contact analysis, updated versions of standard conformational validations for protein and RNA, plus two recent additions: first, flags for cis-nonPro and twisted peptides, and second, the CaBLAM system for diagnosing secondary structure, validating Cα backbone, and validating adjacent peptide CO orientations in the context of the Cα trace. In general, automated ab initio building of starting models is quite good at backbone connectivity but often fails at local conformation or sequence register, especially at poorer than 3.5 Å resolution. However, we show that even if criteria (such as Ramachandran or rotamer) are explicitly restrained to improve refinement behavior and overall validation scores, automated optimization of a deposited structure seldom corrects specific misfittings that start in the wrong local minimum, but just hides them. Therefore, local problems should be identified, and as many as possible corrected, before starting refinement. Secondary structures are confusing at 3-4 Å but can be better recognized at 6-8 Å. In future model challenges, specific steps being tested (such as segmentation) and the required documentation (such as PDB code of starting model) should each be explicitly defined, so competing methods on a given task can be meaningfully compared. Individual local examples are presented here, to understand what local mistakes and corrections look like in 3D, how they probably arise, and what possible improvements to methodology might help avoid them. At these resolutions, both structural biologists and end-users need meaningful estimates of local uncertainty, perhaps through explicit ensembles. Fitting problems can best be diagnosed by validation that spans multiple residues; CaBLAM is such a multi-residue tool, and its effectiveness is demonstrated.