cryo-EM Density Research Articles

Acceleration of cryo-EM Flexible Fitting for Large Biomolecular Systems by Efficient Space Partitioning

Flexible fitting is a powerful technique to build the 3D structures of biomolecules from cryoelectron microscopy (cryo-EM) density maps. One popular method is a cross-correlation coefficient-based approach, where the molecular dynamics (MD) simulation is carried out with the biasing potential that includes the cross-correlation coefficient between the experimental and simulated density maps. Here, we propose efficient parallelization schemes for the calculation of the cross-correlation coefficient to accelerate flexible fitting. Our schemes are tested for small, medium, and large biomolecules using CPU and hybrid CPU+ GPU architectures. The scheme for the atomic decomposition MD is suitable for small proteins such as Ca2+-ATPase with the all-atom Go model, while that for the domain decomposition MD is better for larger systems such as ribosome with the all-atom Go or the all-atom explicit solvent models. Our methods allow flexible fitting for various biomolecules with reasonable computational cost. This approach also connects high-resolution structure refinements with investigation of protein structure-function relationship.

Automation and assessment of de novo modeling with Pathwalking in near atomic resolution cryoEM density maps.

With the rapidly growing number of macromolecular structures solved to near-atomic resolution using electron cryomicroscopy (cryoEM), map interpretation and model building directly from the density without the use of structural templates has become increasingly important. As part of the 2015/2016 Map and Model Challenge, we attempted to assess our latest de novo modeling tool, Pathwalking, in terms of performance and usability, as well as identify areas for future improvements. In total, we applied Pathwalking to six density maps between 3 and 4.5 Å resolution selected from the challenge data sets. In five of the six cases, Pathwalking was able to accurately determine the protein fold and in three of these cases, the final all atom model had less than 1.6 Å RMSD when compared to the known structure. Model building and refinement was nearly completely automated, used default parameters and took less than 30 min to complete a refined all atom model. A direct outgrowth of this work was a more streamlined automated command line Pathwalking utility, as well as a novel sequence assignment and optimization routine, which attempts to register sidechain density with expected side chain volume. In total, Pathwalking offers a nearly complete, robust and efficient method for constructing atomistic protein structures directly from a density map without the aid of a template.

Assessment of structural features in Cryo-EM density maps using SSE and side chain Z-scores

We introduce a new method for assessing resolvability of structural features in density maps from Cryo-Electron Microscopy (Cryo-EM) using fitted or derived models. It calculates Z-scores for secondary structure elements (SSEs) and side chains. Z-scores capture how much larger the cross-correlation score (CCS) is for atoms in such features at their placed locations compared to the CCS at displaced positions. Z-scores are larger when the structural features are well-resolved, as confirmed by visual analysis. This method was applied to all 66 maps submitted to the 2015/2016 EMDB map challenge. For each map, the fitted model provided by the map committee was used in this assessment. The average Z-scores for each map and fitted model correlate moderately well with reported map resolutions (r2 = 0.45 for SSE Z-scores and r2 = 0.56 for side chain Z-scores). Rankings of the submitted maps based on average Z-scores seem to more closely agree with visual analysis. Z-scores can also be used to pinpoint which parts of a model are well-resolved in a map, and which parts of the model may need further fitting or refinement to make the model better match the density.

Reference-based restoration of local contrast for cryo-EM density interpretation

Reference-based restoration of local contrast for cryo-EM density interpretation

Exploratory Studies Detecting Secondary Structures in Medium Resolution 3D Cryo-EM Images Using Deep Convolutional Neural Networks.

Cryo-electron microscopy (cryo-EM) is an emerging biophysical technique for structural determination of protein complexes. However, accurate detection of secondary structures is still challenging when cryo-EM density maps are at medium resolutions (5-10 Å). Most of existing methods are image processing methods that do not fully utilize available images in the cryo-EM database. In this paper, we present a deep learning approach to segment secondary structure elements as helices and β-sheets from medium-resolution density maps. The proposed 3D convolutional neural network is shown to detect secondary structure locations with an F1 score between 0.79 and 0.88 for six simulated test cases. The architecture was also applied to an experimentally-derived cryo-EM density map with good accuracy.

Using Combined Features to Analyze Atomic Structures Derived from Cryo-EM Density Maps.

Cryo-electron microscopy (cryo-EM) has become a major technique for protein structure determination. Many atomic structures have been derived from cryo-EM density maps of about 3Å resolution. Side-chain conformations are well determined in density maps with super-resolutions such as 1-2Å. It is desirable to have a statistical method to detect anomalous side-chains without a super-resolution density map. In this study, we analyzed structures derived from X-ray density maps with higher than 1.5Å resolution and those from cryo-EM density maps with 2-4 Å and 4-6 Å resolutions respectively. We introduce a histogram-based outlier score (HBOS) for anomaly detection in protein models built from cryo-EM density maps. This method uses the statistics derived from X-ray dataset (<1.5Å) as the reference and combines five features involving the distal block distance, side-chain length, phi, psi, and first chi angle of the residue. Higher percentages of anomalies were observed in the cryo-EM models than in the super-resolution X-ray models. Lower percentages of anomalies were observed in cryo-EM models derived after January 2017 than those derived before 2017.

Constructing atomic structural models into cryo-EM densities using molecular dynamics - Pros and cons.

Accurate structure determination from electron density maps at 3-5 Å resolution necessitates a balance between extensive global and local sampling of atomistic models, yet with the stereochemical correctness of backbone and sidechain geometries. Molecular Dynamics Flexible Fitting (MDFF), particularly through a resolution-exchange scheme, ReMDFF, provides a robust way of achieving this balance for hybrid structure determination. Employing two high-resolution density maps, namely that of β-galactosidase at 3.2 Å and TRPV1 at 3.4 Å, we showcase the quality of ReMDFF-generated models, comparing them against ones submitted by independent research groups for the 2015-2016 Cryo-EM Model Challenge. This comparison offers a clear evaluation of ReMDFF's strengths and shortcomings, and those of data-guided real-space refinements in general. ReMDFF results scored highly on the various metric for judging the quality-of-fit and quality-of-model. However, some systematic discrepancies are also noted employing a Molprobity analysis, that are reproducible across multiple competition entries. A space of key refinement parameters is explored within ReMDFF to observe their impact within the final model. Choice of force field parameters and initial model seem to have the most significant impact on ReMDFF model-quality. To this end, very recently developed CHARMM36m force field parameters provide now more refined ReMDFF models than the ones originally submitted to the Cryo-EM challenge. Finally, a set of good-practices is prescribed for the community to benefit from the MDFF developments.

Amyloid fibril structure of α-synuclein determined by cryo-electron microscopy.

α-Synuclein (α-syn) amyloid fibrils are the major component of Lewy bodies, which are the pathological hallmark of Parkinson's disease (PD) and other synucleinopathies. High-resolution structure of α-syn fibril is important for understanding its assembly and pathological mechanism. Here, we determined a fibril structure of full-length α-syn (1-140) at the resolution of 3.07 Å by cryo-electron microscopy (cryo-EM). The fibrils are cytotoxic, and transmissible to induce endogenous α-syn aggregation in primary neurons. Based on the reconstructed cryo-EM density map, we were able to unambiguously build the fibril structure comprising residues 37-99. The α-syn amyloid fibril structure shows two protofilaments intertwining along an approximate 21 screw axis into a left-handed helix. Each protofilament features a Greek key-like topology. Remarkably, five out of the six early-onset PD familial mutations are located at the dimer interface of the fibril (H50Q, G51D, and A53T/E) or involved in the stabilization of the protofilament (E46K). Furthermore, these PD mutations lead to the formation of fibrils with polymorphic structures distinct from that of the wild-type. Our study provides molecular insight into the fibrillar assembly of α-syn at the atomic level and sheds light on the molecular pathogenesis caused by familial PD mutations of α-syn.

A methodology using Gaussian-based density map approximation to assess sets of cryo-electron microscopy density maps

This article presents a methodology to assess a set of density maps, as used in the Blind Assessment Phase of the 2015/2016 Map Challenge (EMDataBank Validation Challenges). The synthetic and experimental cryo-electron microscopy (cryo-EM) density maps obtained by different single particle analysis protocols and by different participants, submitted in the Challenge Phase for assessment, were analyzed with this methodology and the obtained results are presented and discussed here. The goal of using such a methodology was to blindly identify the density maps with globally similar structural information, meaning the maps with the structural information mostly "reproduced" by different protocols. To this end, the density maps are "coarsened" using Gaussian-based approximations, with the same input approximation parameters for all maps of the target biological complex. The approximated maps are then represented in a common reduced-dimension (here, 3D) space of their correlation-coefficient-based distances, in which close maps mean similar maps. The distance analysis allows identifying maps with the most "reproduced" structural information by different protocols. The obtained results are also discussed taking into account the detailed information about the protocols that has been released after the end of the Blind Assessment Phase.

FragFit: a web-application for interactive modeling of protein segments into cryo-EM density maps.

Cryo-electron microscopy (cryo-EM) is a standard method to determine the three-dimensional structures of molecular complexes. However, easy to use tools for modeling of protein segments into cryo-EM maps are sparse. Here, we present the FragFit web-application, a web server for interactive modeling of segments of up to 35 amino acids length into cryo-EM density maps. The fragments are provided by a regularly updated database containing at the moment about 1 billion entries extracted from PDB structures and can be readily integrated into a protein structure. Fragments are selected based on geometric criteria, sequence similarity and fit into a given cryo-EM density map. Web-based molecular visualization with the NGL Viewer allows interactive selection of fragments. The FragFit web-application, accessible at http://proteinformatics.de/FragFit, is free and open to all users, without any login requirements.

Resolution Measurement from a Single Reconstructed Cryo-EM Density Map with Multiscale Spectral Analysis.

As a relatively new technology to solve the three-dimensional (3D) structure of a protein or protein complex, single-particle reconstruction (SPR) of cryogenic electron microscopy (cryo-EM) images shows much superiority and is in a rapidly developing stage. Resolution measurement in SPR, which evaluates the quality of a reconstructed 3D density map, plays a critical role in promoting methodology development of SPR and structural biology. Because there is no benchmark map in the generation of a new structure, how to realize the resolution estimation of a new map is still an open problem. Existing approaches try to generate a hypothetical benchmark map by reconstructing two 3D models from two halves of the original 2D images for cross-reference, which may result in a premature estimation with a half-data model. In this paper, we report a new self-reference-based resolution estimation protocol, called SRes, that requires only a single reconstructed 3D map. The core idea of SRes is to perform a multiscale spectral analysis (MSSA) on the map through multiple size-variable masks segmenting the map. The MSSA-derived multiscale spectral signal-to-noise ratios (mSSNRs) reveal that their corresponding estimated resolutions will show a cliff jump phenomenon, indicating a significant change in the SSNR properties. The critical point on the cliff borderline is demonstrated to be the right estimator for the resolution of the map.

Atomic Resolution Cryo-EM Structure of β-Galactosidase

The advent of direct electron detectors has enabled the routine use of single-particle cryo-electron microscopy (EM) approaches to determine structures of a variety of protein complexes at near-atomic resolution. Here, we report the development of methods to account for local variations in defocus and beam-induced drift, and the implementation of a data-driven dose compensation scheme that significantly improves the extraction of high-resolution information recorded during exposure of the specimen to the electron beam. These advances enable determination of a cryo-EM density map for β-galactosidase bound to the inhibitor phenylethyl β-D-thiogalactopyranoside where the ordered regions are resolved at a level of detail seen in X-ray maps at ∼ 1.5Å resolution. Using this density map in conjunction with constrained molecular dynamics simulations provides a measure of the local flexibility of the non-covalently bound inhibitor and offers further opportunities for structure-guided inhibitor design.

Cryo-EM structure of the adenosine A2A receptor coupled to an engineered heterotrimeric G protein.

The adenosine A2A receptor (A2AR) is a prototypical G protein-coupled receptor (GPCR) that couples to the heterotrimeric G protein GS. Here, we determine the structure by electron cryo-microscopy (cryo-EM) of A2AR at pH 7.5 bound to the small molecule agonist NECA and coupled to an engineered heterotrimeric G protein, which contains mini-GS, the βγ subunits and nanobody Nb35. Most regions of the complex have a resolution of ~3.8 Å or better. Comparison with the 3.4 Å resolution crystal structure shows that the receptor and mini-GS are virtually identical and that the density of the side chains and ligand are of comparable quality. However, the cryo-EM density map also indicates regions that are flexible in comparison to the crystal structures, which unexpectedly includes regions in the ligand binding pocket. In addition, an interaction between intracellular loop 1 of the receptor and the β subunit of the G protein was observed.

An Investigation of Atomic Structures Derived from X-ray Crystallography and Cryo-Electron Microscopy Using Distal Blocks of Side-Chains.

Cryo-electron microscopy (cryo-EM) is a structure determination method for large molecular complexes. As more and more atomic structures are determined using this technique, it is becoming possible to perform statistical characterization of side-chain conformations. Two data sets were involved to characterize block lengths for each of the 18 types of amino acids. One set contains 9131 structures resolved using X-ray crystallography from density maps with better than or equal to 1.5 Å resolutions, and the other contains 237 protein structures derived from cryo-EM density maps with 2–4 Å resolutions. The results show that the normalized probability density function of block lengths is similar between the X-ray data set and the cryo-EM data set for most of the residue types, but differences were observed for ARG, GLU, ILE, LYS, PHE, TRP, and TYR for which conformations with certain shorter block lengths are more likely to be observed in the cryo-EM set with 2–4 Å resolutions.

Iterative Molecular Dynamics−Rosetta Membrane Protein Structure Refinement Guided by Cryo-EM Densities

Knowing atomistic details of proteins is essential not only for the understanding of protein function but also for the development of drugs. Experimental methods such as X-ray crystallography, NMR, and cryo-electron microscopy (cryo-EM) are the preferred forms of protein structure determination and have achieved great success over the most recent decades. Computational methods may be an alternative when experimental techniques fail. However, computational methods are severely limited when it comes to predicting larger macromolecule structures with little sequence similarity to known structures. The incorporation of experimental restraints in computational methods is becoming increasingly important to more reliably predict protein structure. One such experimental input used in structure prediction and refinement is cryo-EM densities. Recent advances in cryo-EM have arguably revolutionized the field of structural biology. Our previously developed cryo-EM-guided Rosetta–MD protocol has shown great promise in...

Binding and Regulation of the Cardiac Ryanodine Receptor by PKA and CaMKII

Ryanodine Receptors (RyRs) are the largest known (2.2 MDa) ion channels, controlling the release of Ca2+ from the sarcoplasmic and endoplasmic reticula. A previously solved high-resolution crystal structure of a tandem repeat domain contains multiple phosphorylation sites for kinases PKA and CaMKII, but the mechanism through which phosphorylation alters channel activity is not understood. Aberrant phosphorylation has been linked to a multitude of disorders, affecting channel gating in both skeletal muscle (RyR1) and cardiac (RyR2) isoforms. Within the RyR2 phosphorylation domain, residues Ser2808 and Ser2814 have received extensive attention as targets for PKA and CaMKII. Docking into cryo-electron microscopy maps suggests a location in the turret region, implying that mutations and phosphorylation may affect the allosteric motions within this area. However, the cryo-EM density in this area is poorly resolved, so other methods are required to look at the precise impact of phosphorylation. Here we present NMR, ITC, and X-ray crystallographic data on the influence of PKA and CaMKII binding to RyR2 as well as the effects of phosphomimetics (S2808D, S2814D) on the structure and dynamics of RyR2 repeat domain. In agreement with previous mass-spectrometry investigation, multiple sites can be phosphorylated by both kinases. Most of these are located directly within a flexible loop, but another site is located on a different face of the domain. Importantly, both S2808D and S2814D show different fingerprints on the spectra, suggesting that they may have different functional impacts on the channel. Finally, we show the preference of PKA for S2808 via ITC studies and are currently investigating the interplay between PKA and CaMKII where we test the impact of prior incubation with one kinase on the activity/affinity of the other.

Fully Automated Correlation-Based Refinement of Atomic Models into High Resolution Cryo-EM Density Maps

Fully Automated Correlation-Based Refinement of Atomic Models into High Resolution Cryo-EM Density Maps

Structural and functional analysis of ribosome assembly factor Efg1.

Ribosome biogenesis in eukaryotes is a complicated process that involves association and dissociation of numerous assembly factors and snoRNAs. The yeast small ribosomal subunit is first assembled into 90S pre-ribosomes in an ordered and dynamic manner. Efg1 is a protein with no recognizable domain that is associated with early 90S particles. Here, we determine the crystal structure of Efg1 from Chaetomium thermophilum at 3.3 Å resolution, revealing a novel elongated all-helical structure. Efg1 is not located in recently determined cryo-EM densities of 90S likely due to its low abundance in mature 90S. Genetic analysis in Saccharomyces cerevisiae shows that the functional core of Efg1 contains two helical hairpins composed of highly conserved residues. Depletion of Efg1 blocks 18S rRNA processing at sites A1 and A2, but not at site A0, and production of small ribosomal subunits. Efg1 is initially recruited by the 5′ domain of 18S rRNA. Its absence disturbs the assembly of the 5′ domain and inhibits release of U14 snoRNA from 90S. Our study shows that Efg1 is required for early assembly and reorganization of the 5′ domain of 18S rRNA.

Allosteric effects in bacteriophage HK97 procapsids revealed directly from covariance analysis of cryo EM data.

The information content of cryo EM data sets exceeds that of the electron scattering potential (cryo EM) density initially derived for structure determination. Previously we demonstrated the power of data variance analysis for characterizing regions of cryo EM density that displayed functionally important variance anomalies associated with maturation cleavage events in Nudaurelia Omega Capensis Virus and the presence or absence of a maturation protease in bacteriophage HK97 procapsids. Here we extend the analysis in two ways. First, instead of imposing icosahedral symmetry on every particle in the data set during the variance analysis, we only assume that the data set as a whole has icosahedral symmetry. This change removes artifacts of high variance along icosahedral symmetry axes, but retains all of the features previously reported in the HK97 data set. Second we present a covariance analysis that reveals correlations in structural dynamics (variance) between the interior of the HK97 procapsid with the protease and regions of the exterior (not seen in the absence of the protease). The latter analysis corresponds well with hydrogen deuterium exchange studies previously published that reveal the same correlation.

A Pattern Recognition Tool for Medium-resolution Cryo-EM Density Maps and Low-resolution Cryo-ET Density maps.

Cryo-electron microscopy (Cryo-EM) and cryo-electron tomography (cryo-ET) produce 3-D density maps of biological molecules at a range of resolution levels. Pattern recognition tools are important in distinguishing biological components from volumetric maps with the available resolutions. One of the most distinct characters in density maps at medium (5-10 Å) resolution is the visibility of protein secondary structures. Although computational methods have been developed, the accurate detection of helices and β-strands from cryo-EM density maps is still an active research area. We have developed a tool for protein secondary structure detection and evaluation of medium resolution 3-D cryo-EM density maps which combines three computational methods (SSETracer, StrandTwister, and AxisComparison). The program was integrated in UCSF Chimera, a popular visualization software in the cryo-EM community. In related work, we have developed BundleTrac, a computational method to trace filaments in a bundle from lower resolution cryo-ET density maps. It has been applied to actin filament tracing in stereocilia with good accuracy and can be potentially added as a tool in Chimera.