Comparing Cryo-EM Reconstructions and Validating Atomic Model Fit Using Difference Maps

Alignment of density maps in Wasserstein distance

In this article, we propose an algorithm for aligning three-dimensional objects when represented as density maps, motivated by applications in cryogenic electron microscopy. The algorithm is based on minimizing the 1-Wasserstein distance between the density maps after a rigid transformation. The induced loss function enjoys a more benign landscape than its Euclidean counterpart and Bayesian optimization is employed for computation. Numerical experiments show improved accuracy and efficiency over existing algorithms on the alignment of real protein molecules. In the context of aligning heterogeneous pairs, we illustrate a potential need for new distance functions.

Accurate global and local 3D alignment of cryo-EM density maps using local spatial structural features

Advances in cryo-electron microscopy (cryo-EM) imaging technologies have led to a rapidly increasing number of cryo-EM density maps. Alignment and comparison of density maps play a crucial role in interpreting structural information, such as conformational heterogeneity analysis using global alignment and atomic model assembly through local alignment. Here, we present a fast and accurate global and local cryo-EM density map alignment method called CryoAlign, that leverages local density feature descriptors to capture spatial structure similarities. CryoAlign is a feature-based cryo-EM map alignment tool, in which the employment of feature-based architecture enables the rapid establishment of point pair correspondences and robust estimation of alignment parameters. Extensive experimental evaluations demonstrate the superiority of CryoAlign over the existing methods in terms of both alignment accuracy and speed.

ChemEM: Flexible Docking of Small Molecules in Cryo-EM Structures

Cryo-electron microscopy (cryo-EM), through resolution advancements, has become pivotal in structure-based drug discovery. However, most cryo-EM structures are solved at 3-4 Å resolution, posing challenges for small-molecule docking and structure-based virtual screening due to issues in the precise positioning of ligands and the surrounding side chains. We present ChemEM, a software package that employs cryo-EM data for the accurate docking of one or multiple ligands in a protein-binding site. Validated against a highly curated benchmark of high- and medium-resolution cryo-EM structures and the corresponding high-resolution controls, ChemEM displayed impressive performance, accurately placing ligands in all but one case, often surpassing cryo-EM PDB-deposited solutions. Even without including the cryo-EM density, the ChemEM scoring function outperformed the well-established AutoDock Vina score. Using ChemEM, we illustrate that valuable information can be extracted from maps at medium resolution and underline the utility of cryo-EM structures for drug discovery.

In-cell structural insight into the stability of sperm microtubule doublet

The propulsion for mammalian sperm swimming is generated by flagella beating. Microtubule doublets (DMTs) along with microtubule inner proteins (MIPs) are essential structural blocks of flagella. However, the intricate molecular architecture of intact sperm DMT remains elusive. Here, by in situ cryo-electron tomography, we solved the in-cell structure of mouse sperm DMT at 4.5-7.5 Å resolutions, and built its model with 36 kinds of MIPs in 48 nm periodicity. We identified multiple copies of Tektin5 that reinforce Tektin bundle, and multiple MIPs with different periodicities that anchor the Tektin bundle to tubulin wall. This architecture contributes to a superior stability of A-tubule than B-tubule of DMT, which was revealed by structural comparison of DMTs from the intact and deformed axonemes. Our work provides an overall molecular picture of intact sperm DMT in 48 nm periodicity that is essential to understand the molecular mechanism of sperm motility as well as the related ciliopathies.

ArpC5 isoforms regulate Arp2/3 complex-dependent protrusion through differential Ena/VASP positioning

Regulation of the Arp2/3 complex is required for productive nucleation of branched actin networks. An emerging aspect of regulation is the incorporation of subunit isoforms into the Arp2/3 complex. Specifically, both ArpC5 subunit isoforms, ArpC5 and ArpC5L, have been reported to fine-tune nucleation activity and branch junction stability. We have combined reverse genetics and cellular structural biology to describe how ArpC5 and ArpC5L differentially affect cell migration. Both define the structural stability of ArpC1 in branch junctions and, in turn, by determining protrusion characteristics, affect protein dynamics and actin network ultrastructure. ArpC5 isoforms also affect the positioning of members of the Ena/Vasodilator-stimulated phosphoprotein (VASP) family of actin filament elongators, which mediate ArpC5 isoform-specific effects on the actin assembly level. Our results suggest that ArpC5 and Ena/VASP proteins are part of a signaling pathway enhancing cell migration.

How to correct relative voxel scale factors for calculations of vector-difference Fourier maps in cryo-EM

The atomic coordinates derived from cryo-electron microscopy (cryo-EM) maps can be inaccurate when the voxel scaling factors are not properly calibrated. Here, we describe a method for correcting relative voxel scaling factors between pairs of cryo-EM maps for the same or similar structures that are expanded or contracted relative to each other. We find that the correction of scaling factors reduces the amplitude differences of Fourier-inverted structure factors from voxel-rescaled maps by up to 20-30%, as shown by two cryo-EM maps of the SARS-CoV-2 spike protein measured at pH 4.0 and pH 8.0. This allows for the calculation of the difference map after properly scaling, revealing differences between the two structures for individual amino acid residues. Unexpectedly, the analysis uncovers two previously overlooked differences of amino acid residues in structures and their local structural changes. Furthermore, we demonstrate the method as applied to two cryo-EM maps of monomeric apo-photosystem II from the cyanobacteria Synechocystis sp. PCC 6803 and Thermosynechococcus elongatus. The resulting difference maps reveal many changes in the peripheral transmembrane PsbX subunit between the two species.

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.

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.

Cryo-EM reveals mechanisms of angiotensin I-converting enzyme allostery and dimerization

Hypertension (high blood pressure) is a major risk factor for cardiovascular disease, which is the leading cause of death worldwide. The somatic isoform of angiotensin I‐converting enzyme (sACE) plays a critical role in blood pressure regulation, and ACE inhibitors are thus widely used to treat hypertension and cardiovascular disease. Our current understanding of sACE structure, dynamics, function, and inhibition has been limited because truncated, minimally glycosylated forms of sACE are typically used for X‐ray crystallography and molecular dynamics simulations. Here, we report the first cryo‐EM structures of full‐length, glycosylated, soluble sACE (sACE^S1211). Both monomeric and dimeric forms of the highly flexible apo enzyme were reconstructed from a single dataset. The N‐ and C‐terminal domains of monomeric sACE^S1211 were resolved at 3.7 and 4.1 Å, respectively, while the interacting N‐terminal domains responsible for dimer formation were resolved at 3.8 Å. Mechanisms are proposed for intradomain hinging, cooperativity, and homodimerization. Furthermore, the observation that both domains were in the open conformation has implications for the design of sACE modulators.

Adeno-associated virus (AAV) cell entry: structural insights

Adeno-associated virus (AAV) is the leading vector in emerging treatments of inherited diseases. Higher transduction efficiencies and cellular specificity are required for broader clinical application, motivating investigations of virus-host molecular interactions during cell entry. High-throughput methods are identifying host proteins more comprehensively, with subsequent molecular studies revealing unanticipated complexity and serotype specificity. Cryogenic electron microscopy (cryo-EM) provides a path towards structural details of these sometimes heterogeneous virus-host complexes, and is poised to illuminate more fully the steps in entry. Here presented, is progress in understanding the distinct steps of glycan attachment, and receptor-mediated entry/trafficking. Comparison with structures of antibody complexes provides new insights on immune neutralization with implications for the design of improved gene therapy vectors.

Structure of the decoy module of human glycoprotein 2 and uromodulin and its interaction with bacterial adhesin FimH

Glycoprotein 2 (GP2) and uromodulin (UMOD) filaments protect against gastrointestinal and urinary tract infections by acting as decoys for bacterial fimbrial lectin FimH. By combining AlphaFold2 predictions with X-ray crystallography and cryo-EM, we show that these proteins contain a bipartite decoy module whose new fold presents the high-mannose glycan recognized by FimH. The structure rationalizes UMOD mutations associated with kidney diseases and visualizes a key epitope implicated in cast nephropathy.

AlphaFold2 predictions, X-ray crystallography and cryo-EM analyses reveal how related human glycoproteins GP2 and uromodulin catch pathogenic bacteria by presenting a high-mannose glycan that acts as a decoy for fimbrial adhesin FimH.

Cryo-EM single-particle structure refinement and map calculation using Servalcat

In 2020, cryo-EM single-particle analysis achieved true atomic resolution thanks to technological developments in hardware and software. The number of high-resolution reconstructions continues to grow, increasing the importance of the accurate determination of atomic coordinates. Here, a new Python package and program called Servalcat is presented that is designed to facilitate atomic model refinement. Servalcat implements a refinement pipeline using the program REFMAC5 from the CCP4 package. After the refinement, Servalcat calculates a weighted Fo - Fc difference map, which is derived from Bayesian statistics. This map helps manual and automatic model building in real space, as is common practice in crystallography. The Fo - Fc map helps in the visualization of weak features including hydrogen densities. Although hydrogen densities are weak, they are stronger than in the electron-density maps produced by X-ray crystallography, and some H atoms are even visible at ∼1.8 Å resolution. Servalcat also facilitates atomic model refinement under symmetry constraints. If point-group symmetry has been applied to the map during reconstruction, the asymmetric unit model is refined with the appropriate symmetry constraints.

Cryo-EM structure and kinetics reveal electron transfer by 2D diffusion of cytochrome c in the yeast III-IV respiratory supercomplex

Energy conversion in aerobic organisms involves an electron current from low-potential donors, such as NADH and succinate, to dioxygen through the membrane-bound respiratory chain. Electron transfer is coupled to transmembrane proton transport, which maintains the electrochemical proton gradient used to produce ATP and drive other cellular processes. Electrons are transferred from respiratory complexes III to IV (CIII and CIV) by water-soluble cytochrome (cyt.) c In Saccharomyces cerevisiae and some other organisms, these complexes assemble into larger CIII₂CIV_1/2 supercomplexes, the functional significance of which has remained enigmatic. In this work, we measured the kinetics of the S. cerevisiae supercomplex cyt. c-mediated QH₂:O₂ oxidoreductase activity under various conditions. The data indicate that the electronic link between CIII and CIV is confined to the surface of the supercomplex. Single-particle electron cryomicroscopy (cryo-EM) structures of the supercomplex with cyt. c show the positively charged cyt. c bound to either CIII or CIV or along a continuum of intermediate positions. Collectively, the structural and kinetic data indicate that cyt. c travels along a negatively charged patch on the supercomplex surface. Thus, rather than enhancing electron transfer rates by decreasing the distance that cyt. c must diffuse in three dimensions, formation of the CIII₂CIV_1/2 supercomplex facilitates electron transfer by two-dimensional (2D) diffusion of cyt. c This mechanism enables the CIII₂CIV_1/2 supercomplex to increase QH₂:O₂ oxidoreductase activity and suggests a possible regulatory role for supercomplex formation in the respiratory chain.

Cryo-EM density maps adjustment for subtraction, consensus and sharpening

Electron cryomicroscopy (cryo-EM) has emerged as a powerful structural biology instrument to solve near-atomic three-dimensional structures. Despite the fast growth in the number of density maps generated from cryo-EM data, comparison tools among these reconstructions are still lacking. Current proposals to compare cryo-EM data derived volumes perform map subtraction based on adjustment of each volume grey level to the same scale. We present here a more sophisticated way of adjusting the volumes before comparing, which implies adjustment of grey level scale and spectrum energy, but keeping phases intact inside a mask and imposing the results to be strictly positive. The adjustment that we propose leaves the volumes in the same numeric frame, allowing to perform operations among the adjusted volumes in a more reliable way. This adjustment can be a preliminary step for several applications such as comparison through subtraction, map sharpening, or combination of volumes through a consensus that selects the best resolved parts of each input map. Our development might also be used as a sharpening method using an atomic model as a reference. We illustrate the applicability of this algorithm with the reconstructions derived of several experimental examples. This algorithm is implemented in Xmipp software package and its applications are user-friendly accessible through the cryo-EM image processing framework Scipion.

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.

Full-length de novo protein structure determination from cryo-EM maps using deep learning

MOTIVATION

Advances in microscopy instruments and image processing algorithms have led to an increasing number of cryo-EM maps. However, building accurate models for the EM maps at 3-5 Å resolution remains a challenging and time-consuming process. With the rapid growth of deposited EM maps, there is an increasing gap between the maps and reconstructed/modeled 3-dimensional (3D) structures. Therefore, automatic reconstruction of atomic-accuracy full-atomstructures fromEMmaps is pressingly needed.

RESULTS

We present a semi-automatic de novo structure determination method using a deep learningbased framework, named as DeepMM, which builds atomic-accuracy all-atom models from cryo-EM maps at near-atomic resolution. In our method, the main-chain and Cα positions as well as their amino acid and secondary structure types are predicted in the EM map using Densely Connected Convolutional Networks. DeepMM was extensively validated on 40 simulated maps at 5 Å resolution and 30 experimental maps at 2.6-4.8 Å resolution as well as an EMDB-wide data set of 2931 experimental maps at 2.6-4.9 Å resolution, and compared with state-of-the-art algorithms including RosettaES, MAINMAST, and Phenix. Overall, our DeepMM algorithm obtained a significant improvement over existing methods in terms of both accuracy and coverage in building full-length protein structures on all test sets, demonstrating the efficacy and general applicability of DeepMM.

AVAILABILITY

http://huanglab.phys.hust.edu.cn/DeepMM.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

VESPER: global and local cryo-EM map alignment using local density vectors

An increasing number of density maps of biological macromolecules have been determined by cryo-electron microscopy (cryo-EM) and stored in the public database, EMDB. To interpret the structural information contained in EM density maps, alignment of maps is an essential step for structure modeling, comparison of maps, and for database search. Here, we developed VESPER, which captures the similarity of underlying molecular structures embedded in density maps by taking local gradient directions into consideration. Compared to existing methods, VESPER achieved substantially more accurate global and local alignment of maps as well as database retrieval.

Macromolecular refinement of X-ray and cryoelectron microscopy structures with Phenix/OPLS3e for improved structure and ligand quality

With the advent of the resolution revolution in cryoelectron microscopy (cryo-EM), low-resolution refinement is common, and likewise increases the need for a reliable force field. Here, we report on the incorporation of the OPLS3e force field with the VSGB2.1 solvation model in the structure determination package Phenix. Our results show significantly improved structure quality and reduced ligand strain at lower resolution for X-ray refinement. For refinement of cryo-EM-based structures, we find comparable quality structures, goodness-of-fit, and reduced ligand strain. We also show how structure quality and ligand strain are related to the map-model cross-correlation as a function of data weight, and how that can detect overfitting. Signs of overfitting are found in over half of our cryo-EM dataset, which can be remedied by a re-refinement at a lower data weight. Finally, a start-to-end script for refining structures with Phenix/OPLS3e is available in the Schrödinger 2020-3 distribution.

TEMPy2: a Python library with improved 3D electron microscopy density-fitting and validation workflows

Structural determination of molecular complexes by cryo-EM requires large, often complex processing of the image data that are initially obtained. Here, TEMPy2, an update of the TEMPy package to process, optimize and assess cryo-EM maps and the structures fitted to them, is described. New optimization routines, comprehensive automated checks and workflows to perform these tasks are described.