cryo-EM Map Research Articles

Solving the α-helix correspondence problem at medium-resolution Cryo-EM maps through modeling and 3D matching

Cryo-electron microscopy (cryo-EM) has recently emerged as a prominent biophysical method for macromolecular structure determination. Many research efforts have been devoted to produce cryo-EM images, density maps, at near-atomic resolution. Despite many advances in technology, the resolution of the generated density maps may not be sufficiently adequate and informative to directly construct the atomic structure of proteins. At medium-resolution (∼4-10 Å), secondary structure elements (α-helices and β-sheets) are discernible, whereas finding the correspondence of secondary structure elements detected in the density map with those on the sequence remains a challenging problem. In this paper, an automatic framework is proposed to solve α-helix correspondence problem in three-dimensional space. Through modeling of the sequence with the aid of a novel strategy, the α-helix correspondence problem is initially transformed into a complete weighted bipartite graph matching problem. An innovative correlation-based scoring function based on a well-known and robust statistical method is proposed for weighting the graph. Moreover, two local optimization algorithms, which are Greedy and Improved Greedy algorithms, have been presented to find α-helix correspondence. A widely used data set including 16 reconstructed and 4 experimental cryo-EM maps were chosen to verify the accuracy and reliability of the proposed automatic method. The experimental results demonstrate that the automatic method is highly efficient (86.25% accuracy), robust (11.3% error rate), fast (∼1.4 s), and works independently from cryo-EM skeleton.

Real-space quantum-based refinement for cryo-EM: Q|R#3.

Electron cryo-microscopy (cryo-EM) is rapidly becoming a major competitor toX-ray crystallography, especially for large structures that are difficult or impossible to crystallize. While recent spectacular technological improvements have led to significantly higher resolution three-dimensional reconstructions, the average quality of cryo-EM maps is still at the low-resolution end of the range compared with crystallography. A long-standing challenge for atomic model refinement has been the production of stereochemically meaningful models for this resolution regime. Here, it is demonstrated that including accurate model geometry restraints derived from ab initio quantum-chemical calculations (HF-D3/6-31G) can improve the refinement of an example structure (chain A of PDB entry 3j63). The robustness of the procedure is tested for additional structures with up to 7000 atoms (PDB entry 3a5x and chain C of PDB entry 5fn5) using the less expensive semi-empirical (GFN1-xTB) model. The necessary algorithms enabling real-space quantum refinement have been implemented in the latest version of qr.refine and are described here.

Multiplexed Functional Assessment of Genetic Variants in CARD11

SummaryGenetic testing has increased the number of variants identified in disease genes, but the diagnostic utility is limited by lack of understanding variant function. CARD11 encodes an adaptor protein that expresses dominant-negative and gain-of-function variants associated with distinct immunodeficiencies. Here, we used a “cloning-free” saturation genome editing approach in a diploid cell line to simultaneously score 2,542 variants for decreased or increased function in the region of CARD11 associated with immunodeficiency. We also described an exon-skipping mechanism for CARD11 dominant-negative activity. The classification of reported clinical variants was sensitive (94.6%) and specific (88.9%), which rendered the data immediately useful for interpretation of seven coding and splicing variants implicated in immunodeficiency found in our clinic. This approach is generalizable for variant interpretation in many other clinically actionable genes, in any relevant cell type.

Structural basis for polarized elongation of actin filaments

Actin filaments elongate and shorten much faster at their barbed end than their pointed end, but the molecular basis of this difference has not been understood. We use all-atom molecular dynamics simulations to investigate the properties of subunits at both ends of the filament. The terminal subunits tend toward conformations that resemble actin monomers in solution, while contacts with neighboring subunits progressively flatten the conformation of internal subunits. At the barbed end the terminal subunit is loosely tethered by its DNase-1 loop to the third subunit, because its monomer-like conformation precludes stabilizing contacts with the penultimate subunit. The motions of the terminal subunit make the partially flattened penultimate subunit accessible for binding monomers. At the pointed end, unique contacts between the penultimate and terminal subunits are consistent with existing cryogenic electron microscopic (cryo-EM) maps, limit binding to incoming monomers, and flatten the terminal subunit, which likely promotes ATP hydrolysis and rapid phosphate release. These structures explain the distinct polymerization kinetics of the two ends.

Structural basis for SARM1 inhibition and activation under energetic stress.

SARM1, an executor of axonal degeneration, displays NADase activity that depletes the key cellular metabolite, NAD+, in response to nerve injury. The basis of SARM1 inhibition and its activation under stress conditions are still unknown. Here, we present cryo-EM maps of SARM1 at 2.9 and 2.7 Å resolutions. These indicate that SARM1 homo-octamer avoids premature activation by assuming a packed conformation, with ordered inner and peripheral rings, that prevents dimerization and activation of the catalytic domains. This inactive conformation is stabilized by binding of SARM1's own substrate NAD+ in an allosteric location, away from the catalytic sites. This model was validated by mutagenesis of the allosteric site, which led to constitutively active SARM1. We propose that the reduction of cellular NAD+ concentration contributes to the disassembly of SARM1's peripheral ring, which allows formation of active NADase domain dimers, thereby further depleting NAD+ to cause an energetic catastrophe and cell death.

Structure of Nora virus at 2.7\xa0\xc5 resolution and implications for receptor binding, capsid stability and taxonomy

Nora virus, a virus of Drosophila, encapsidates one of the largest single-stranded RNA virus genomes known. Its taxonomic affinity is uncertain as it has a picornavirus-like cassette of enzymes for virus replication, but the capsid structure was at the time for genome publication unknown. By solving the structure of the virus, and through sequence comparison, we clear up this taxonomic ambiguity in the invertebrate RNA virosphere. Despite the lack of detectable similarity in the amino acid sequences, the 2.7 Å resolution cryoEM map showed Nora virus to have T = 1 symmetry with the characteristic capsid protein β-barrels found in all the viruses in the Picornavirales order. Strikingly, α-helical bundles formed from the extended C-termini of capsid protein VP4B and VP4C protrude from the capsid surface. They are similar to signalling molecule folds and implicated in virus entry. Unlike other viruses of Picornavirales, no intra-pentamer stabilizing annulus was seen, instead the intra-pentamer stability comes from the interaction of VP4C and VP4B N-termini. Finally, intertwining of the N-termini of two-fold symmetry-related VP4A capsid proteins and RNA, provides inter-pentamer stability. Based on its distinct structural elements and the genetic distance to other picorna-like viruses we propose that Nora virus, and a small group of related viruses, should have its own family within the order Picornavirales.

Crystal structure of the Cenp-HIKHead-TW sub-module of the inner kinetochore CCAN complex.

Kinetochores are large multi-subunit complexes that attach centromeric chromatin to microtubules of the mitotic spindle, enabling sister chromatid segregation in mitosis. The inner kinetochore constitutive centromere associated network (CCAN) complex assembles onto the centromere-specific Cenp-A nucleosome (Cenp-ANuc), thereby coupling the centromere to the microtubule-binding outer kinetochore. CCAN is a conserved 14–16 subunit complex composed of discrete modules. Here, we determined the crystal structure of the Saccharomyces cerevisiae Cenp-HIKHead-TW sub-module, revealing how Cenp-HIK and Cenp-TW interact at the conserved Cenp-HIKHead–Cenp-TW interface. A major interface is formed by the C-terminal anti-parallel α-helices of the histone fold extension (HFE) of the Cenp-T histone fold domain (HFD) combining with α-helix H3 of Cenp-K to create a compact three α-helical bundle. We fitted the Cenp-HIKHead-TW sub-module to the previously determined cryo-EM map of the S. cerevisiae CCAN–Cenp-ANuc complex. This showed that the HEAT repeat domain of Cenp-IHead and C-terminal HFD of Cenp-T of the Cenp-HIKHead-TW sub-module interact with the nucleosome DNA gyre at a site close to the Cenp-ANuc dyad axis. Our structure provides a framework for understanding how Cenp-T links centromeric Cenp-ANuc to the outer kinetochore through its HFD and N-terminal Ndc80-binding motif, respectively.

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 ofa 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.

Combine Cryo-EM Density Map and Residue Contact for Protein Structure Prediction - A Case Study.

Cryo-electron microscopy is a major structure determination technique for large molecular machines and membrane-associated complexes. Although atomic structures have been determined directly from cryo-EM density maps with high resolutions, current structure determination methods for medium resolution (5 to 10 Å) cryo-EM maps are limited by the availability of structure templates. Secondary structure traces are lines detected from a cryo-EM density map for α-helices and β-strands of a protein. When combined with secondary structure sequence segments predicted from a protein sequence, it is possible to generate a set of likely topologies of α-traces and β-sheet traces. A topology describes the overall folding relationship among secondary structures; it is a critical piece of information for deriving the corresponding atomic structure. We propose a method for protein structure prediction that combines three sources of information: the secondary structure traces detected from the cryo-EM density map, predicted secondary structure sequence segments, and amino acid contact pairs predicted using MULTICOM. A case study shows that using amino acid contact prediction from MULTICOM improves the ranking of the true topology. Our observations convey that using a small set of highly voted secondary structure contact pairs enhances the ranking in all experiments conducted for this case.

Cryo-EM structure of the entire mammalian F-type ATP synthase.

The majority of adenosine triphosphate (ATP) powering cellular processes in eukaryotes is produced by the mitochondrial F1Fo ATP synthase. Here, we present the atomic models of the membrane Fo domain and the entire mammalian (ovine) F1Fo, determined by cryo-electron microscopy. Subunits in the membrane domain are arranged in the 'proton translocation cluster' attached to the c-ring and a more distant 'hook apparatus' holding subunit e. Unexpectedly, this subunit is anchored to a lipid 'plug' capping the c-ring. We present a detailed proton translocation pathway in mammalian Fo and key inter-monomer contacts in F1Fo multimers. Cryo-EM maps of F1Fo exposed to calcium reveal a retracted subunit e and a disassembled c-ring, suggesting permeability transition pore opening. We propose a model for the permeability transition pore opening, whereby subunit e pulls the lipid plug out of the c-ring. Our structure will allow the design of drugs for many emerging applications in medicine.

Apple latent spherical virus structure with stable capsid frame supports quasi-stable protrusions expediting genome release

Picorna-like plant viruses are non-enveloped RNA spherical viruses of ~30 nm. Part of the survival of these viruses depends on their capsid being stable enough to harbour the viral genome and yet malleable enough to allow its release. However, molecular mechanisms remain obscure. Here, we report a structure of a picorna-like plant virus, apple latent spherical virus, at 2.87 Å resolution by single-particle cryo-electron microscopy (cryo-EM) with a cold-field emission beam. The cryo-EM map reveals a unique structure composed of three capsid proteins Vp25, Vp20, and Vp24. Strikingly Vp25 has a long N-terminal extension, which substantially stabilises the capsid frame of Vp25 and Vp20 subunits. Cryo-EM images also resolve RNA genome leaking from a pentameric protrusion of Vp24 subunits. The structures and observations suggest that genome release occurs through occasional opening of the Vp24 subunits, possibly suppressed to a low frequency by the rigid frame of the other subunits.

Reformulation of the self-guided molecular simulation method.

Self-guided molecular/Langevin dynamics (SGMD/SGLD) simulation methods were developed to enhance conformational sampling through promoting low frequency motion of molecular systems and have been successfully applied in many simulation studies. Quantitative understanding of conformational distribution in SGLD has been achieved by separating microscopic properties according to frequency. However, a missing link between the guiding factors and conformational distributions makes it highly empirical and system dependent when choosing the values of the guiding parameters. Based on the understanding that molecular interactions are the source of energy barriers and diffusion friction, this work reformulates the equation of the low frequency motion to resemble Langevin dynamics. This reformulation leads to new forms of guiding forces and establishes a relation between the guiding factors and conformational distributions. We call simulations with these new guiding forces the generalized self-guided molecular/Langevin dynamics (SGMDg/SGLDg). In addition, we present a new way to calculate low frequency properties and an efficient algorithm to implement SGMDg/SGLDg that minimizes memory usage and inter-processor communication. Through example simulations with a skewed double well system, an argon fluid, and a cryo-EM map flexible fitting case, we demonstrate the guiding effects on conformational distributions and conformational searching.

Connexin-46/50 in a dynamic lipid environment resolved by CryoEM at 1.9\u2009\xc5

Gap junctions establish direct pathways for cells to transfer metabolic and electrical messages. The local lipid environment is known to affect the structure, stability and intercellular channel activity of gap junctions; however, the molecular basis for these effects remains unknown. Here, we incorporate native connexin-46/50 (Cx46/50) intercellular channels into a dual lipid nanodisc system, mimicking a native cell-to-cell junction. Structural characterization by CryoEM reveals a lipid-induced stabilization to the channel, resulting in a 3D reconstruction at 1.9 Å resolution. Together with all-atom molecular dynamics simulations, it is shown that Cx46/50 in turn imparts long-range stabilization to the dynamic local lipid environment that is specific to the extracellular lipid leaflet. In addition, ~400 water molecules are resolved in the CryoEM map, localized throughout the intercellular permeation pathway and contributing to the channel architecture. These results illustrate how the aqueous-lipid environment is integrated with the architectural stability, structure and function of gap junction communication channels.

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.

Assessment of Force Field Accuracy Using Cryogenic Electron Microscopy Data of Hyper-thermostable Glutamate Dehydrogenase

Molecular dynamics (MD) simulations in biophysically relevant time scales of microseconds is a powerful tool for studying biomolecular processes, but results often display force field dependency. Therefore, assessment of force field accuracy using experimental data of biomolecules in solution is essential for simulation studies. Here, we propose the use of structural models obtained via cryo-electron microscopy (cryoEM), which provides biomolecular structures in vitreous ice mimicking the environment in solution. The accuracy of the AMBER (ff99SB-ILDN-NMR, ff14SB, ff15ipq, and ff15FB) and CHARMM (CHARMM22 and CHARMM36m) force fields was assessed by comparing their MD trajectories with the cryoEM data of thermostable hexameric glutamate dehydrogenase (GDH), which included a cryoEM map at a resolution of approximately 3 Å and structure models of subunits reflecting metastable conformations in domain motion occurring in GDH. In the assessment, we validated the force fields with respect to the reproducibility and stability of secondary structures and intersubunit interactions in the cryoEM data. Furthermore, we evaluated the force fields regarding the reproducibility of the energy landscape in the domain motion expected from the cryoEM data. As a result, among the six force fields, ff15FB and ff99SB-ILDN-NMR displayed good agreement with the experiment. The present study demonstrated the advantages of the high-resolution cryoEM map and suggested the optimal force field to reproduce experimentally observed protein structures.

Hierarchical natural move Monte Carlo refines flexible RNA structures into cryo-EM densities.

Ribonucleic acids (RNAs) play essential roles in living cells. Many of them fold into defined three-dimensional (3D) structures to perform functions. Recent advances in single-particle cryo-electron microscopy (cryo-EM) have enabled structure determinations of RNA to atomic resolutions. However, most RNA molecules are structurally flexible, limiting the resolution of their structures solved by cryo-EM. In modeling these molecules, several computational methods are limited by the requirement of massive computational resources and/or the low efficiency in exploring large-scale structural variations. Here we use hierarchical natural move Monte Carlo (HNMMC), which takes advantage of collective motions for groups of nucleic acid residues, to refine RNA structures into their cryo-EM maps, preserving atomic details in the models. After validating the method on a simulated density map of tRNA, we applied it to objectively obtain the model of the folding intermediate for the specificity domain of ribonuclease P from Bacillus subtilis and refine a flexible ribosomal RNA (rRNA) expansion segment from the Mycobacterium tuberculosis (Mtb) ribosome in different conformational states. Finally, we used HNMMC to model atomic details and flexibility for two distinct conformations of the complete genomic RNA (gRNA) inside MS2, a single-stranded RNA virus, revealing multiple pathways for its capsid assembly.

Deep learning enables the atomic structure determination of the Fanconi Anemia core complex from cryoEM.

Cryo-electron microscopy of protein complexes often leads to moderate resolution maps (4-8 Å), with visible secondary-structure elements but poorly resolved loops, making model building challenging. In the absence of high-resolution structures of homologues, only coarse-grained structural features are typically inferred from these maps, and it is often impossible to assign specific regions of density to individual protein subunits. This paper describes a new method for overcoming these difficulties that integrates predicted residue distance distributions from a deep-learned convolutional neural network, computational protein folding using Rosetta, and automated EM-map-guided complex assembly. We apply this method to a 4.6 Å resolution cryoEM map of Fanconi Anemia core complex (FAcc), an E3 ubiquitin ligase required for DNA interstrand crosslink repair, which was previously challenging to interpret as it comprises 6557 residues, only 1897 of which are covered by homology models. In the published model built from this map, only 387 residues could be assigned to the specific subunits with confidence. By building and placing into density 42 deep-learning-guided models containing 4795 residues not included in the previously published structure, we are able to determine an almost-complete atomic model of FAcc, in which 5182 of the 6557 residues were placed. The resulting model is consistent with previously published biochemical data, and facilitates interpretation of disease-related mutational data. We anticipate that our approach will be broadly useful for cryoEM structure determination of large complexes containing many subunits for which there are no homologues of known structure.

Improvement of cryo-EM maps by density modification.

A density modification procedure for improving maps from single-particle electron cryo-microscopy 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 electron cryo-microscopy, and that half-maps with independent errors are available in electron cryo-microscopy. 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.

The de Rham-Hodge Analysis and Modeling of Biomolecules.

Biological macromolecules have intricate structures that underpin their biological functions. Understanding their structure-function relationships remains a challenge due to their structural complexity and functional variability. Although de Rham-Hodge theory, a landmark of twentieth-century mathematics, has had a tremendous impact on mathematics and physics, it has not been devised for macromolecular modeling and analysis. In this work, we introduce de Rham-Hodge theory as a unified paradigm for analyzing the geometry, topology, flexibility, and Hodge mode analysis of biological macromolecules. Geometric characteristics and topological invariants are obtained either from the Helmholtz-Hodge decomposition of the scalar, vector, and/or tensor fields of a macromolecule or from the spectral analysis of various Laplace-de Rham operators defined on the molecular manifolds. We propose Laplace-de Rham spectral-based models for predicting macromolecular flexibility. We further construct a Laplace-de Rham-Helfrich operator for revealing cryo-EM natural frequencies. Extensive experiments are carried out to demonstrate that the proposed de Rham-Hodge paradigm is one of the most versatile tools for the multiscale modeling and analysis of biological macromolecules and subcellular organelles. Accurate, reliable, and topological structure-preserving algorithms for implementing discrete exterior calculus (DEC) have been developed to facilitate the aforementioned modeling and analysis of biological macromolecules. The proposed de Rham-Hodge paradigm has potential applications to subcellular organelles and the structure construction from medium- or low-resolution cryo-EM maps, and functional predictions from massive biomolecular datasets.

Measurement of Atom Resolvability in CryoEM Maps with Q-scores

Measurement of Atom Resolvability in CryoEM Maps with Q-scores