EMDA: A Python package for Electron Microscopy Data Analysis

X-rays, electrons, and neutrons as probes of atomic matter

X-rays, electrons, and neutrons probe different properties of matter. X-rays feel electron density (ED). Electrons sense the electrostatic potential (ESP) of electrons and nuclei. Neutrons are sensitive to nuclear coherent scattering length (NCSL). While NCSL maps are widely understood to be different, ED and ESP maps are tacitly assumed to be similar. Here, I show that the belief in ED and ESP map equivalence is mistaken, but contains a grain of truth. Using density functional theory (DFT), the Bethe-Mott (BM) relation, and the Thomas-Fermi (TF) and Cromer-Mann (CM) atomic models, I show that ED and ESP maps are indeed more similar to each other than to NCSL maps. Nevertheless, peak and integrated map values depend differently on the atomic order number and on the contributions from electrons in the inner and outer CM shells. ED and ESP maps also differ in the sign and relative magnitude of excess charge effects.

The open gate of the AMPA receptor forms a Ca2+ binding site critical in regulating ion transport

Alpha-amino-3-hydroxyl-5-methyl-4-isoxazole-propionic acid receptors (AMPARs) are cation-selective ion channels that mediate most fast excitatory neurotransmission in the brain. Although their gating mechanism has been studied extensively, understanding how cations traverse the pore has remained elusive. Here we investigated putative ion and water densities in the open pore of Ca2+-permeable AMPARs (rat GRIA2 flip-Q isoform) at 2.3-2.6 Å resolution. We show that the ion permeation pathway attains an extracellular Ca2+ binding site (site-G) when the channel gate moves into the open configuration. Site-G is highly selective for Ca2+ over Na+, favoring the movement of Ca2+ into the selectivity filter of the pore. Seizure-related N619K mutation, adjacent to site-G, promotes channel opening but attenuates Ca2+ binding and thus diminishes Ca2+ permeability. Our work identifies the importance of site-G, which coordinates with the Q/R site of the selectivity filter to ensure the preferential transport of Ca2+ through the channel pore.

Community recommendations on cryoEM data archiving and validation

In January 2020, a workshop was held at EMBL-EBI (Hinxton, UK) to discuss data requirements for the deposition and validation of cryoEM structures, with a focus on single-particle analysis. The meeting was attended by 47 experts in data processing, model building and refinement, validation, and archiving of such structures. This report describes the workshop's motivation and history, the topics discussed, and the resulting consensus recommendations. Some challenges for future methods-development efforts in this area are also highlighted, as is the implementation to date of some of the recommendations.

Community recommendations on cryoEM data archiving and validation: Outcomes of a wwPDB/EMDB workshop on cryoEM data management, deposition and validation

In January 2020, a workshop was held at EMBL-EBI (Hinxton, UK) to discuss data requirements for deposition and validation of cryoEM structures, with a focus on single-particle analysis. The meeting was attended by 47 experts in data processing, model building and refinement, validation, and archiving of such structures. This report describes the workshop's motivation and history, the topics discussed, and consensus recommendations resulting from the workshop. Some challenges for future methods-development efforts in this area are also highlighted, as is the implementation to date of some of the recommendations.

Accurate magnification determination for cryoEM using gold

Determining the correct magnified pixel size of single-particle cryoEM micrographs is necessary to maximize resolution and enable accurate model building. Here we describe a simple and rapid procedure for determining the absolute magnification in an electron cryomicroscope to a precision of <0.5%. We show how to use the atomic lattice spacings of crystals of thin and readily available test specimens, such as gold, as an absolute reference to determine magnification for both room temperature and cryogenic imaging. We compare this method to other commonly used methods, and show that it provides comparable accuracy in spite of its simplicity. This magnification calibration method provides a definitive reference quantity for data analysis and processing, simplifies the combination of multiple datasets from different microscopes and detectors, and improves the accuracy with which the contrast transfer function of the microscope can be determined. We also provide an open source program, magCalEM, which can be used to accurately estimate the magnified pixel size of a cryoEM dataset ex post facto.

Outcomes of the EMDataResource Cryo-EM Ligand Modeling Challenge

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein/nucleic-acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: E. coli beta-galactosidase with inhibitor, SARS-CoV-2 RNA-dependent RNA polymerase with covalently bound nucleotide analog, and SARS-CoV-2 ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. We found that (1) the quality of submitted ligand models and surrounding atoms varied, as judged by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics, and contact scores, and (2) a composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Architecture of the biofilm-associated archaic Chaperone-Usher pilus CupE from Pseudomonas aeruginosa

Chaperone-Usher Pathway (CUP) pili are major adhesins in Gram-negative bacteria, mediating bacterial adherence to biotic and abiotic surfaces. While classical CUP pili have been extensively characterized, little is known about so-called archaic CUP pili, which are phylogenetically widespread and promote biofilm formation by several human pathogens. In this study, we present the electron cryomicroscopy structure of the archaic CupE pilus from the opportunistic human pathogen Pseudomonas aeruginosa. We show that CupE1 subunits within the pilus are arranged in a zigzag architecture, containing an N-terminal donor β-strand extending from each subunit into the next, where it is anchored by hydrophobic interactions, with comparatively weaker interactions at the rest of the inter-subunit interface. Imaging CupE pili on the surface of P. aeruginosa cells using electron cryotomography shows that CupE pili adopt variable curvatures in response to their environment, which might facilitate their role in promoting cellular attachment. Finally, bioinformatic analysis shows the widespread abundance of cupE genes in isolates of P. aeruginosa and the co-occurrence of cupE with other cup clusters, suggesting interdependence of cup pili in regulating bacterial adherence within biofilms. Taken together, our study provides insights into the architecture of archaic CUP pili, providing a structural basis for understanding their role in promoting cellular adhesion and biofilm formation in P. aeruginosa.

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 asymmetry of a photosynthetic supercomplex from green sulfur bacteria

The photochemical reaction center (RC) features a dimeric architecture for charge separation across the membrane. In green sulfur bacteria (GSB), the trimeric Fenna-Matthews-Olson (FMO) complex mediates the transfer of light energy from the chlorosome antenna complex to the RC. Here we determine the structure of the photosynthetic supercomplex from the GSB Chlorobaculum tepidum using single-particle cryogenic electron microscopy (cryo-EM) and identify the cytochrome c subunit (PscC), two accessory protein subunits (PscE and PscF), a second FMO trimeric complex, and a linker pigment between FMO and the RC core. The protein subunits that are assembled with the symmetric RC core generate an asymmetric photosynthetic supercomplex. One linker bacteriochlorophyll (BChl) is located in one of the two FMO-PscA interfaces, leading to differential efficiencies of the two energy transfer branches. The two FMO trimeric complexes establish two different binding interfaces with the RC cytoplasmic surface, driven by the associated accessory subunits. This structure of the GSB photosynthetic supercomplex provides mechanistic insight into the light excitation energy transfer routes and a possible evolutionary transition intermediate of the bacterial photosynthetic supercomplex from the primitive homodimeric RC.

Structure of dynein-dynactin on microtubules shows tandem adaptor binding

Cytoplasmic dynein is a microtubule motor that is activated by its cofactor dynactin and a coiled-coil cargo adaptor1-3. Up to two dynein dimers can be recruited per dynactin, and interactions between them affect their combined motile behaviour4-6. Different coiled-coil adaptors are linked to different cargos7,8, and some share motifs known to contact sites on dynein and dynactin4,9-13. There is limited structural information on how the resulting complex interacts with microtubules and how adaptors are recruited. Here we develop a cryo-electron microscopy processing pipeline to solve the high-resolution structure of dynein-dynactin and the adaptor BICDR1 bound to microtubules. This reveals the asymmetric interactions between neighbouring dynein motor domains and how they relate to motile behaviour. We found that two adaptors occupy the complex. Both adaptors make similar interactions with the dyneins but diverge in their contacts with each other and dynactin. Our structure has implications for the stability and stoichiometry of motor recruitment by cargos.

Pathogen-sugar interactions revealed by universal saturation transfer analysis

Many pathogens exploit host cell-surface glycans. However, precise analyses of glycan ligands binding with heavily modified pathogen proteins can be confounded by overlapping sugar signals and/or compounded with known experimental constraints. Universal saturation transfer analysis (uSTA) builds on existing nuclear magnetic resonance spectroscopy to provide an automated workflow for quantitating protein-ligand interactions. uSTA reveals that early-pandemic, B-origin-lineage severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike trimer binds sialoside sugars in an "end-on" manner. uSTA-guided modeling and a high-resolution cryo-electron microscopy structure implicate the spike N-terminal domain (NTD) and confirm end-on binding. This finding rationalizes the effect of NTD mutations that abolish sugar binding in SARS-CoV-2 variants of concern. Together with genetic variance analyses in early pandemic patient cohorts, this binding implicates a sialylated polylactosamine motif found on tetraantennary N-linked glycoproteins deep in the human lung as potentially relevant to virulence and/or zoonosis.

A molecular prior distribution for Bayesian inference based on Wilson statistics

BACKGROUND AND OBJECTIVE

Wilson statistics describe well the power spectrum of proteins at high frequencies. Therefore, it has found several applications in structural biology, e.g., it is the basis for sharpening steps used in cryogenic electron microscopy (cryo-EM). A recent paper gave the first rigorous proof of Wilson statistics based on a formalism of Wilson's original argument. This new analysis also leads to statistical estimates of the scattering potential of proteins that reveal a correlation between neighboring Fourier coefficients. Here we exploit these estimates to craft a novel prior that can be used for Bayesian inference of molecular structures.

METHODS

We describe the properties of the prior and the computation of its hyperparameters. We then evaluate the prior on two synthetic linear inverse problems, and compare against a popular prior in cryo-EM reconstruction at a range of SNRs.

RESULTS

We show that the new prior effectively suppresses noise and fills-in low SNR regions in the spectral domain. Furthermore, it improves the resolution of estimates on the problems considered for a wide range of SNR and produces Fourier Shell Correlation curves that are insensitive to masking effects.

CONCLUSIONS

We analyze the assumptions in the model, discuss relations to other regularization strategies, and postulate on potential implications for structure determination in cryo-EM.

Cryo-electron Microscopy of Adeno-associated Virus

Adeno-associated virus (AAV) has a single-stranded DNA genome encapsidated in a small icosahedrally symmetric protein shell with 60 subunits. AAV is the leading delivery vector in emerging gene therapy treatments for inherited disorders, so its structure and molecular interactions with human hosts are of intense interest. A wide array of electron microscopic approaches have been used to visualize the virus and its complexes, depending on the scientific question, technology available, and amenability of the sample. Approaches range from subvolume tomographic analyses of complexes with large and flexible host proteins to detailed analysis of atomic interactions within the virus and with small ligands at resolutions as high as 1.6 Å. Analyses have led to the reclassification of glycan receptors as attachment factors, to structures with a new-found receptor protein, to identification of the epitopes of antibodies, and a new understanding of possible neutralization mechanisms. AAV is now well-enough characterized that it has also become a model system for EM methods development. Heralding a new era, cryo-EM is now also being deployed as an analytic tool in the process development and production quality control of high value pharmaceutical biologics, namely AAV vectors.

Validation analysis of EMDB entries

The Electron Microscopy Data Bank (EMDB) is the central archive of the electron cryo-microscopy (cryo-EM) community for storing and disseminating volume maps and tomograms. With input from the community, EMDB has developed new resources for the validation of cryo-EM structures, focusing on the quality of the volume data alone and that of the fit of any models, themselves archived in the Protein Data Bank (PDB), to the volume data. Based on recommendations from community experts, the validation resources are developed in a three-tiered system. Tier 1 covers an extensive and evolving set of validation metrics, including tried and tested metrics as well as more experimental ones, which are calculated for all EMDB entries and presented in the Validation Analysis (VA) web resource. This system is particularly useful for cryo-EM experts, both to validate individual structures and to assess the utility of new validation metrics. Tier 2 comprises a subset of the validation metrics covered by the VA resource that have been subjected to extensive testing and are considered to be useful for specialists as well as nonspecialists. These metrics are presented on the entry-specific web pages for the entire archive on the EMDB website. As more experience is gained with the metrics included in the VA resource, it is expected that consensus will emerge in the community regarding a subset that is suitable for inclusion in the tier 2 system. Tier 3, finally, consists of the validation reports and servers that are produced by the Worldwide Protein Data Bank (wwPDB) Consortium. Successful metrics from tier 2 will be proposed for inclusion in the wwPDB validation pipeline and reports. The details of the new resource are described, with an emphasis on the tier 1 system. The output of all three tiers is publicly available, either through the EMDB website (tiers 1 and 2) or through the wwPDB ftp sites (tier 3), although the content of all three will evolve over time (fastest for tier 1 and slowest for tier 3). It is our hope that these validation resources will help the cryo-EM community to obtain a better understanding of the quality and of the best ways to assess the quality of cryo-EM structures in EMDB and PDB.

Assessing the Mobility of Severe Acute Respiratory Syndrome Coronavirus-2 Spike Protein Glycans by Structural and Computational Methods

Two years after its emergence, the coronavirus disease-2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) remains difficult to control despite the availability of several vaccines. The extensively glycosylated SARS-CoV-2 spike (S) protein, which mediates host cell entry by binding to the angiotensin converting enzyme 2 (ACE2) through its receptor binding domain (RBD), is the major target of neutralizing antibodies. Like to many other viral fusion proteins, the SARS-CoV-2 spike protein utilizes a glycan shield to thwart the host immune response. To grasp the influence of chemical signatures on carbohydrate mobility and reconcile the cryo-EM density of specific glycans we combined our cryo-EM map of the S ectodomain to 4.1 Å resolution, reconstructed from a limited number of particles, and all-atom molecular dynamics simulations. Chemical modifications modeled on representative glycans (defucosylation, sialylation and addition of terminal LacNAc units) show no significant influence on either protein shielding or glycan flexibility. By estimating at selected sites the local correlation between the full density map and atomic model-based maps derived from molecular dynamics simulations, we provide insight into the geometries of the α-Man-(1→3)-[α-Man-(1→6)-]-β-Man-(1→4)-β-GlcNAc(1→4)-β-GlcNAc core common to all N-glycosylation sites.