EM Density Map Research Articles

Building Protein Atomic Models from Cryo-EM Density Maps and Residue Co-Evolution.

Electron cryo-microscopy (cryo-EM) has emerged as a powerful method by which to obtain three-dimensional (3D) structures of macromolecular complexes at atomic or near-atomic resolution. However, de novo building of atomic models from near-atomic resolution (3–5 Å) cryo-EM density maps is a challenging task, in particular because poorly resolved side-chain densities hamper sequence assignment by automatic procedures at a lower resolution. Furthermore, segmentation of EM density maps into individual subunits remains a difficult problem when the structure of the subunits is not known, or when significant conformational rearrangement occurs between the isolated and associated form of the subunits. To tackle these issues, we have developed a graph-based method to thread most of the C- trace of the protein backbone into the EM density map. The EM density is described as a weighted graph such that the resulting minimum spanning tree encompasses the high-density regions of the map. A pruning algorithm cleans the tree and finds the most probable positions of the C- atoms, by using side-chain density when available, as a collection of C- trace fragments. By complementing experimental EM maps with contact predictions from sequence co-evolutionary information, we demonstrate that this approach can correctly segment EM maps into individual subunits and assign amino acid sequences to backbone traces to generate atomic models.

Mass-selective and ice-free electron cryomicroscopy protein sample preparation via native electrospray ion-beam deposition.

Despite tremendous advances in sample preparation and classification algorithms for electron cryomicroscopy (cryo-EM) and single-particle analysis (SPA), sample heterogeneity remains a major challenge and can prevent access to high-resolution structures. In addition, optimization of preparation conditions for a given sample can be time-consuming. In the current work, it is demonstrated that native electrospray ion-beam deposition (native ES-IBD) is an alternative, reliable approach for the preparation of extremely high-purity samples, based on mass selection in vacuum. Folded protein ions are generated by native electrospray ionization, separated from other proteins, contaminants, aggregates, and fragments, gently deposited on cryo-EM grids, frozen in liquid nitrogen, and subsequently imaged by cryo-EM. We demonstrate homogeneous coverage of ice-free cryo-EM grids with mass-selected protein complexes. SPA reveals that the complexes remain folded and assembled, but variations in secondary and tertiary structures are currently limiting information in 2D classes and 3D EM density maps. We identify and discuss challenges that need to be addressed to obtain a resolution comparable to that of the established cryo-EM workflow. Our results show the potential of native ES-IBD to increase the scope and throughput of cryo-EM for protein structure determination and provide an essential link between gas-phase and solution-phase protein structures.

Single Particle Cryo-Electron Microscopy: From Sample to Structure.

Cryo-electron microscopy (cryoEM) is a powerful technique for structure determination of macromolecular complexes, via single particle analysis (SPA). The overall process involves i) vitrifying the specimen in a thin film supported on a cryoEM grid; ii) screening the specimen to assess particle distribution and ice quality; iii) if the grid is suitable, collecting a single particle dataset for analysis; and iv) image processing to yield an EM density map. In this protocol, an overview for each of these steps is provided, with a focus on the variables which a user can modify during the workflow and the troubleshooting of common issues. With remote microscope operation becoming standard in many facilities, variations on imaging protocols to assist users in efficient operation and imaging when physical access to the microscope is limited will be described.

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.

Identification of a conserved virion-stabilizing network inside the interprotomer pocket of enteroviruses

Enteroviruses pose a persistent and widespread threat to human physical health, with no specific treatments available. Small molecule capsid binders have the potential to be developed as antivirals that prevent virus attachment and entry into host cells. To aid with broad-range drug development, we report here structures of coxsackieviruses B3 and B4 bound to different interprotomer-targeting capsid binders using single-particle cryo-EM. The EM density maps are beyond 3 Å resolution, providing detailed information about interactions in the ligand-binding pocket. Comparative analysis revealed the residues that form a conserved virion-stabilizing network at the interprotomer site, and showed the small molecule properties that allow anchoring in the pocket to inhibit virus disassembly.

A New Protocol for Atomic-Level Protein Structure Modeling and Refinement Using Low-to-Medium Resolution Cryo-EM Density Maps

The rapid progress of cryo-electron microscopy (cryo-EM) in structural biology has raised an urgent need for robust methods to create and refine atomic-level structural models using low-resolution EM density maps. We propose a new protocol to create initial models using I-TASSER protein structure prediction, followed by EM density map-based rigid-body structure fitting, flexible fragment adjustment and atomic-level structure refinement simulations. The protocol was tested on a large set of 285 non-homologous proteins and generated structural models with correct folds for 260 proteins, where 28% had RMSDs below 2 Å. Compared to other state-of-the-art methods, the major advantage of the proposed pipeline lies in the uniform structure prediction and refinement protocol, as well as the extensive structural re-assembly simulations, which allow for low-to-medium resolution EM density map-guided structure modeling starting from amino acid sequences. Interestingly, the quality of both the image fitting and subsequent structure refinement was found to be strongly correlated with the correctness of the initial I-TASSER models; this is mainly due to the different correlation patterns observed between force field and structural quality for the models with template modeling score (or TM-score, a metric quantifying the similarity of models to the native) above and below a threshold of 0.5. Overall, the results demonstrate a new avenue that is ready to use for large-scale cryo-EM-based structure modeling and atomic-level density map-guided structure refinement.

MAINMASTseg: Automated Map Segmentation Method for Cryo-EM Density Maps with Symmetry.

For structural interpretation of cryo-electron microscopy (cryo-EM) density maps that contain multiple chains, map segmentation is an important step. If a map is segmented accurately into regions of individual protein components, the structure of each protein can be separately modeled using an existing modeling tool. Here, we developed new software, MAINMASTseg, for segmenting maps with symmetry. MAINMASTseg is an extension of the MAINMAST de novo cryo-EM protein structure modeling tool, which builds protein structures from a graph structure that captures the distribution of salient density points in the map. MAINMASTseg uses this graph and segments the map by considering symmetry corresponding density points in the graph. We tested MAINMASTseg on a data set of 38 experimentally determined EM density maps. MAINMASTseg successfully identified an individual protein unit for the majority of the maps, which was significantly better than two other popular existing methods, Segger and Phenix. The software is made freely available for academic users at http://kiharalab.org/mainmast_seg.

Molecular Basis for Hormone Recognition and Activation of Corticotropin-Releasing Factor Receptors.

Corticotropin-releasing factor (CRF) and the three related peptides urocortins 1-3 (UCN1-UCN3) are endocrine hormones that control the stress responses by activating CRF1R and CRF2R, two members of class B G-protein-coupled receptors (GPCRs). Here, we present two cryoelectron microscopy (cryo-EM) structures of UCN1-bound CRF1R and CRF2R with the stimulatory G protein. In both structures, UCN1 adopts a single straight helix with its N terminus dipped into the receptor transmembrane bundle. Although the peptide-binding residues in CRF1R and CRF2R are different from other members of class B GPCRs, the residues involved in receptor activation and G protein coupling are conserved. In addition, both structures reveal bound cholesterol molecules to the receptor transmembrane helices. Our structures define the basis of ligand-binding specificity in the CRF receptor-hormone system, establish a common mechanism of class B GPCR activation and G protein coupling, and provide a paradigm for studying membrane protein-lipid interactions for class B GPCRs.

Protein-Protein Modeling Using Cryo-EM Restraints.

Recent improvements in cryo-electron microscopy (cryo-EM) in the past few years are now allowing to observe molecular complexes at atomic resolution. As a consequence, numerous structures derived from cryo-EM are now available in the Protein Data Bank. However, if for some complexes atomic resolution is reached, this is not true for all. This is also the case in cryo-electron tomography where the achievable resolution is still limited. Furthermore the resolution in a cryo-EM map is not a constant, with often outer regions being of lower resolution, possibly linked to conformational variability. Although those low- to medium-resolution EM maps (or regions thereof) cannot directly provide atomic structure of large molecular complexes, they provide valuable information to model the individual components and their assembly into them. Most approaches for this kind of modeling are performing rigid fitting of the individual components into the EM density map. While this would appear an obvious option, they ignore key aspects of molecular recognition, the energetics and flexibility of the interfaces. Moreover, this often restricts the modeling to a unique source of data, the EM density map.In this chapter, we describe a protocol where an EM map is used as restraint in HADDOCK to guide the modeling process. In the first step, rigid-body fitting is performed with PowerFit in order to identify the most likely locations of the molecules into the map. These are then used as centroids to which distance restraints are defined from the center of mass of the components of the complex for the initial rigid-body docking. The EM density is then directly used as an additional restraint energy term, which can be combined with all the other types of data supported by HADDOCK. This protocol relies on the new version 2.4 of both the HADDOCK webserver and software. Preparation steps consisting of cropping the EM map and rigid-body fitting of the atomic structure are explained. Then, the EM-driven docking protocol using HADDOCK is illustrated.

164. Reporting the High-resolution Structure of the Enterococcal Ribosome: A New Template for Antibiotic Discovery

BackgroundThe ribosome is a rich target for antibiotic design and its structural secrets have been described at the atomic level over the past 2 decades. However, most bacterial ribosome structures come from nonpathogenic species of Archaea or thermophilic bacteria. To aid in the development of modern antibiotics against the enterococcus, we report the structure of the ribosome from Enterococcus faecalis at 3.5 Å resolution using cryo-electron microscopy.Methods E. faecalis strain OG1 was grown in liquid culture, collected and lysed using a French press. 70S ribosomes were purified using centrifugation through a sucrose cushion followed by column chromatography and sucrose gradient centrifugation. 70S particles were diluted in buffer and applied to a holey carbon grid and using an FEI vitrobot were flash-frozen in liquid ethane. Data were collected on an FEI Titan Krios operating at 300 kV acceleration voltage. The particles classified into 6 distinct structures based on their composition. Completed maps were utilized for structure modelling using Coot and were then refined using real space refinement within Phenix.ResultsHigh-quality maps of the 70S ribosome were obtained at up to 3.5 Å resolution in several distinct conformations. The 23S, 16S, and 5S RNA structures were almost completely built into maps with clear density. All but 2 ribosome proteins L25 and L33 have been placed in density. The A, P, and E sites were built into unambiguous density and found to be consistent with other bacterial structures. Notably, 1 EM density map contains an uncharged t-RNA molecule in the E site. The sites identified for current antibiotics are also well defined and interpretable. This 70S structural platform is suitable for structural analysis of antibiotic binding sites, especially for those antibiotics directed specifically against the enterococcal ribosome.ConclusionFor the first time, the structure of the ribosome from the important human pathogen Enterococcus faecalis has been determined. The maps were obtained at high resolution and found to be suitable for antibiotic design. It is anticipated that the continued determination of the structures of ribosomes from pathogens will aid in the discovery of new treatments for infectious diseases.Disclosures All authors: No reported disclosures.

PRISM-EM: template interface-based modelling of multi-protein complexes guided by cryo-electron microscopy density maps.

The structures of protein assemblies are important for elucidating cellular processes at the molecular level. Three-dimensional electron microscopy (3DEM) is a powerful method to identify the structures of assemblies, especially those that are challenging to study by crystallography. Here, a new approach, PRISM-EM, is reported to computationally generate plausible structural models using a procedure that combines crystallographic structures and density maps obtained from 3DEM. The predictions are validated against seven available structurally different crystallographic complexes. The models display mean deviations in the backbone of <5 Å. PRISM-EM was further tested on different benchmark sets; the accuracy was evaluated with respect to the structure of the complex, and the correlation with EM density maps and interface predictions were evaluated and compared with those obtained using other methods. PRISM-EM was then used to predict the structure of the ternary complex of the HIV-1 envelope glycoprotein trimer, the ligand CD4 and the neutralizing protein m36.

Individual-Particle Electron Tomography: A Method for Monitoring the 3D Structure of Low-Dimensional DNA-Based Nanodevice

DNA base-pairing has been used for many years to direct the arrangement of inorganic nanocrystals into small groupings and arrays with tailored optical and electrical properties. The control of DNA-mediated assembly depends crucially on a better understanding of the three-dimensional (3D) structure of the DNA-nanocrystal hybridized building blocks. Existing techniques do not allow for the structural determination of these flexible and heterogeneous samples. Here, we employed cryo-electron microscopy (cryo-EM) and negative-staining (NS) techniques to investigate the morphology of DNA-nanogold conjugates that were self-assembled from a mixture of an 84-base-pair double-stranded DNA (dsDNA) conjugated with two 5-nm nanogold particles for potential substrates in plasmon coupling experiments. Using NS electron tomography and the individual-particle electron tomography (IPET) reconstruction method, we obtained 3D reconstructed 14 electron density maps at a resolution of ~2 nm from each individual dsDNA-nanogold particle. Using these 3D density maps as constraints, we derived 14 conformations of dsDNA by projected a standard flexible dsDNA model onto the observed EM density maps using Molecular Dynamics (MD) simulations. The variation of the conformations was largely consistent with the variation from liquid solution. Moreover, the IPET approach provides the most complete experimental determination of the flexibility and fluctuation range of these directed nanocrystal assemblies to date. The general features revealed by these experiments can be expected to occur in a broad range of DNA-assembled nanostructures.

Denoising of high-resolution single-particle electron-microscopy density maps by their approximation using three-dimensional Gaussian functions.

Cryo-electron microscopy (cryo-EM) of frozen-hydrated preparations of isolated macromolecular complexes is the method of choice to obtain the structure of complexes that cannot be easily studied by other experimental methods due to their flexibility or large size. An increasing number of macromolecular structures are currently being obtained at subnanometer resolution but the interpretation of structural details in such EM-derived maps is often difficult because of noise at these high-frequency signal components that reduces their contrast. In this paper, we show that the method for EM density-map approximation using Gaussian functions can be used for denoising of single-particle EM maps of high (typically subnanometer) resolution. We show its denoising performance using simulated and experimental EM density maps of several complexes.

Three Dimensional Dynamics and Fluctuations of DNA-Nanogold Dimers by Individual-Particle Electron Tomography

Organic-inorganic hybridized nanocrystals have been proven as an important additional tool for addressing many emerging challenges biological and material sciences. Nanogold and quantum dot conjugates have been used extensively as biomolecular markers, while DNA base pairing has directed the self-assembly of discrete groupings and arrays of organic and inorganic nanocrystals in formation of a network solid for electronic devices and memory components. The control of DNA-mediated assembly depends crucially upon the ability to exposure three-dimensional (3D) structure of DNA-nanocrystals hybridized building blocks. Current techniques limit in structure determination of flexible and heterogeneous sample. Here, we employed the electron tomography and negative-staining techniques to investigate the 3D structure of DNA-nanogold dimers that were self-assembled from a mixture of an 84-base-pair double-stranded DNA (dsDNA) conjugated with two 5-nm nanogold particles. We reconstructed a total of 14 electron density maps at a resolution of ∼2 nm from each individual dsDNA-nanogold particle using the individual-particle electron tomography (IPET) reconstruction method. Using these 3D density maps as constrain, we projected a standard flexible DNA model onto the observed EM density maps and derived 14 conformations of dsDNA by molecular dynamics simulations. By aligning these structures from one distal end nanogold, the distribution showed the DNA fluctuation and flexibility between nanogolds, which providing a blueprint research for studying of DNA-mediated assembly of complex inorganic nanostructures.

Cryo-EM of the pathogenic VCP variant R155P reveals long-range conformational changes in the D2 ATPase ring

Single amino acid mutations in valosin containing protein (VCP/p97), a highly conserved member of the ATPases associated with diverse cellular activities (AAA) family of ATPases has been linked to a severe degenerative disease affecting brain, muscle and bone tissue. Previous studies have demonstrated the role of VCP mutations in altering the ATPase activity of the D2 ring; however the structural consequences of these mutations remain unclear. In this study, we report the three-dimensional (3D) map of the pathogenic VCP variant, R155P, as revealed by single-particle Cryo-Electron Microscopy (EM) analysis at 14 Å resolution. We show that the N-terminal R155P mutation induces a large structural reorganisation of the D2 ATPase ring. Results from docking studies using crystal structure data of available wild-type VCP in the EM density maps indicate that the major difference is localized at the interface between two protomers within the D2 ring. Consistent with a conformational change, the VCP R155P variant shifted the isoelectric point of the protein and reduced its interaction with its well-characterized cofactor, nuclear protein localization-4 (Npl4). Together, our results demonstrate that a single amino acid substitution in the N-terminal domain can relay long-range conformational changes to the distal D2 ATPase ring. Our results provide the first structural clues of how VCP mutations may influence the activity and function of the D2 ATPase ring.

Fast Shape-Based Global and Local Electron Density Map Search

The Electron Microscopy DataBank (EMDB) is growing rapidly, accumulating biological structural data obtained cryo-electron microscopy (cryo-EM). Cryo-EM is an emerging technique for determining large biomolecular complexes and subcellular structures. Together with the Protein Data Bank (PDB), EMDB is becoming a fundamental resource of the tertiary structures of biological macromolecules. To take full advantage of this indispensable resource, the ability to search the database by structural similarity is essential. However, unlike high-resolution structures stored in PDB, methods for comparing low-resolution EM density maps are not well established. Here, we developed a novel computational method for efficiently searching EM maps. The method uses a compact fingerprint representation of EM maps based on the 3D Zernike descriptor, which is a mathematical series expansion for representing isosurface shape of EM maps. The method was implemented in a web server, named EM-SURFER (http://kiharalab.org/em-surfer/), which allows users to search against the entire EMDB with over 2400 entries in a few seconds. By combing with map segmentation, the method can also identify corresponding local regions in EM maps. Examples of search results from different types of query structures are discussed. The unique capability of EM-SURFER to detect 3D shape similarity of low-resolution EM maps should prove invaluable in structural biology.

Small-angle X-ray Solution Scattering Study of the Multi-aminoacyl-tRNA Synthetase Complex Reveals an Elongated and Multi-armed particle

In animal cells, nine aminoacyl-tRNA synthetases are associated with the three auxiliary proteins p18, p38, and p43 to form a stable and conserved large multi-aminoacyl-tRNA synthetase complex (MARS), whose molecular mass has been proposed to be between 1.0 and 1.5 MDa. The complex acts as a molecular hub for coordinating protein synthesis and diverse regulatory signal pathways. Electron microscopy studies defined its low resolution molecular envelope as an overall rather compact, asymmetric triangular shape. Here, we have analyzed the composition and homogeneity of the native mammalian MARS isolated from rabbit liver and characterized its overall internal structure, size, and shape at low resolution by hydrodynamic methods and small-angle x-ray scattering in solution. Our data reveal that the MARS exhibits a much more elongated and multi-armed shape than expected from previous reports. The hydrodynamic and structural features of the MARS are large compared with other supramolecular assemblies involved in translation, including ribosome. The large dimensions and non-compact structural organization of MARS favor a large protein surface accessibility for all its components. This may be essential to allow structural rearrangements between the catalytic and cis-acting tRNA binding domains of the synthetases required for binding the bulky tRNA substrates. This non-compact architecture may also contribute to the spatiotemporal controlled release of some of its components, which participate in non-canonical functions after dissociation from the complex.

Computational methods for constructing protein structure models from 3D electron microscopy maps

Protein structure determination by cryo-electron microscopy (EM) has made significant progress in the past decades. Resolutions of EM maps have been improving as evidenced by recently reported structures that are solved at high resolutions close to 3Å. Computational methods play a key role in interpreting EM data. Among many computational procedures applied to an EM map to obtain protein structure information, in this article we focus on reviewing computational methods that model protein three-dimensional (3D) structures from a 3D EM density map that is constructed from two-dimensional (2D) maps. The computational methods we discuss range from de novo methods, which identify structural elements in an EM map, to structure fitting methods, where known high resolution structures are fit into a low-resolution EM map. A list of available computational tools is also provided.

Electron Crystallography of Euglenoid Four-Transmembrane Protein Revealed the Linear Polymerization by a Combination of Three-Ways of Intermolecular Interaction

In Euglena gracilis, cell membranes on the ridge regions of the striped surface structure are covered with paracrystalline arrays, mainly composed of the integral membrane protein, called IP39. It has recently been cloned as a putative four-membrane-spanning protein with a conserved sequence motif of PMP-22/EMP/MP20/Claudin superfamily. Here, we report the three-dimensional structure of Euglena IP39 at 7 A resolution determined by image-based electron crystallography of two-dimensional (2D) crystals. The 2D crystal array of IP39 appears to be a striated pattern of the antiparallel double-rows, in which the trimeric units of IP39 are longitudinally polymerized, resulting in continuously extending zigzag lines. Each protomer shows an overall cylindrical density (approximately 20 A in diameter) with ‘nock’ and ‘pointed’ shapes at the respective sides of the membrane. Structural analysis of another 2D crystal bound with anti-phosphotyrosine Fab fragments reveals that the ‘nock’-like protruded structure is facing to the cytoplasmic side. A four-helical bundle model is consistently fitted to the EM density map and shows that the transmembrane regions are mainly involved in the intermolecular interactions to form the linear strands. The polymerization pattern of IP39 in the 2D crystal, however, exhibits unexpected interactions between respective protomers. In the trimeric unit of the single strand, one of the three protomers is likely to be rotated at approximately 180 degree in the opposite direction to the others, indicating that there are at least three different ways of possible intermolecular interactions in the transmembrane regions between neighboring protomers. A combination of such multiple interactions would be important for the continuous linear polymerization, thus providing important implications in the strand formation of other four-transmembrane proteins of this family in the lipid environment.

A unified convention for biological assemblies with helical symmetry.

Assemblies with helical symmetry can be conveniently formulated in many distinct ways. Here, a new convention is presented which unifies the two most commonly used helical systems for generating helical assemblies from asymmetric units determined by X-ray fibre diffraction and EM imaging. A helical assembly is viewed as being composed of identical repetitive units in a one- or two-dimensional lattice, named 1-D and 2-D helical systems, respectively. The unification suggests that a new helical description with only four parameters [n(1), n(2), twist, rise], which is called the augmented 1-D helical system, can generate the complete set of helical arrangements, including coverage of helical discontinuities (seams). A unified four-parameter characterization implies similar parameters for similar assemblies, can eliminate errors in reproducing structures of helical assemblies and facilitates the generation of polymorphic ensembles from helical atomic models or EM density maps. Further, guidelines are provided for such a unique description that reflects the structural signature of an assembly, as well as rules for manipulating the helical symmetry presentation.