Cross-validation in cryo-EM-based structural modeling

Structures of human γδ T cell receptor-CD3 complex

Gamma delta (γδ) T cells, a unique T cell subgroup, are crucial in various immune responses and immunopathology1-3. The γδ T cell receptor (TCR), which is generated by γδ T cells, recognizes a diverse range of antigens independently of the major histocompatibility complex2. The γδ TCR associates with CD3 subunits, initiating T cell activation and holding great potential in immunotherapy4. Here we report the structures of two prototypical human Vγ9Vδ2 and Vγ5Vδ1 TCR-CD3 complexes5,6, revealing two distinct assembly mechanisms that depend on Vγ usage. The Vγ9Vδ2 TCR-CD3 complex is monomeric, with considerable conformational flexibility in the TCRγ-TCRδ extracellular domain and connecting peptides. The length of the connecting peptides regulates the ligand association and T cell activation. A cholesterol-like molecule wedges into the transmembrane region, exerting an inhibitory role in TCR signalling. The Vγ5Vδ1 TCR-CD3 complex displays a dimeric architecture, whereby two protomers nestle back to back through the Vγ5 domains of the TCR extracellular domains. Our biochemical and biophysical assays further corroborate the dimeric structure. Importantly, the dimeric form of the Vγ5Vδ1 TCR is essential for T cell activation. These findings reveal organizing principles of the γδ TCR-CD3 complex, providing insights into the unique properties of γδ TCR and facilitating immunotherapeutic interventions.

Cryo-EM of Aβ fibrils from mouse models find tg-APPArcSwe fibrils resemble those found in patients with sporadic Alzheimer's disease

The use of transgenic mice displaying amyloid-β (Aβ) brain pathology has been essential for the preclinical assessment of new treatment strategies for Alzheimer's disease. However, the properties of Aβ in such mice have not been systematically compared to Aβ in the brains of patients with Alzheimer's disease. Here, we determined the structures of nine ex vivo Aβ fibrils from six different mouse models by cryogenic-electron microscopy. We found novel Aβ fibril structures in the APP/PS1, ARTE10 and tg-SwDI models, whereas the human type II filament fold was found in the ARTE10, tg-APPSwe and APP23 models. The tg-APPArcSwe mice showed an Aβ fibril whose structure resembles the human type I filament found in patients with sporadic Alzheimer's disease. A detailed assessment of the Aβ fibril structure is key to the selection of adequate mouse models for the preclinical development of novel plaque-targeting therapeutics and positron emission tomography imaging tracers in Alzheimer's disease.

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-EM structure of the human IgM B cell receptor

The B cell receptor (BCR) initiates immune responses through antigen recognition. We report a 3.3-angstrom cryo-electron microscopy structure of human immunoglobulin M (IgM)-BCR in the resting state. IgM-BCR comprises two heavy chains, two light chains, and the Igα/Igβ heterodimer. The ectodomains of the heavy chains closely stack against those of Igα/Igβ, with one heavy chain locked between Igα and Igβ in the juxtamembrane region. Extracellular interactions may determine isotype specificity of the BCR. The transmembrane helices of IgM-BCR form a four-helix bundle that appears to be conserved among all BCR isotypes. This structure contains 14 glycosylation sites on the IgM-BCR ectodomains and reveals three potential surface binding sites. Our work reveals the organizational principles of the BCR and may facilitate the design of antibody-based therapeutics.

Atomic model validation using the CCP-EM software suite

Recently, there has been a dramatic improvement in the quality and quantity of data derived using cryogenic electron microscopy (cryo-EM). This is also associated with a large increase in the number of atomic models built. Although the best resolutions that are achievable are improving, often the local resolution is variable, and a significant majority of data are still resolved at resolutions worse than 3 Å. Model building and refinement is often challenging at these resolutions, and hence atomic model validation becomes even more crucial to identify less reliable regions of the model. Here, a graphical user interface for atomic model validation, implemented in the CCP-EM software suite, is presented. It is aimed to develop this into a platform where users can access multiple complementary validation metrics that work across a range of resolutions and obtain a summary of evaluations. Based on the validation estimates from atomic models associated with cryo-EM structures from SARS-CoV-2, it was observed that models typically favor adopting the most common conformations over fitting the observations when compared with the model agreement with data. At low resolutions, the stereochemical quality may be favored over data fit, but care should be taken to ensure that the model agrees with the data in terms of resolvable features. It is demonstrated that further re-refinement can lead to improvement of the agreement with data without the loss of geometric quality. This also highlights the need for improved resolution-dependent weight optimization in model refinement and an effective test for overfitting that would help to guide the refinement process.

Predicting RNA Scaffolds with a Hybrid Method of Vfold3D and VfoldLA

The ever-increasing discoveries of noncoding RNA functions draw a strong demand for RNA structure determination from the sequence. In recently years, computational studies for RNA structures, at both the two-dimensional and the three-dimensional levels, led to several highly promising new developments. In this chapter, we describe a hybrid method, which combines the motif template-based Vfold3D model and the loop template-based VfoldLA model, to predict RNA 3D structures. The main emphasis is placed on the definition of motifs and loops, the treatment of no-template motifs, and the 3D structure assembly from templates of motifs and loops. For illustration, we use the ZIKV xrRNA1 as an example to show the template-based prediction of RNA 3D structures from the 2D structure. The web server for the hybrid model is freely accessible at http://rna.physics.missouri.edu/vfold3D2 .

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.

Structural interpretation of cryo-EM image reconstructions

The productivity of single-particle cryo-EM as a structure determination method has rapidly increased as many novel biological structures are being elucidated. The ultimate result of the cryo-EM experiment is an atomic model that should faithfully represent the computed image reconstruction. Although the principal approach of atomic model building and refinement from maps resembles that of the X-ray crystallographic methods, there are important differences due to the unique properties resulting from the 3D image reconstructions. In this review, we discuss the practiced work-flow from the cryo-EM image reconstruction to the atomic model. We give an overview of (i) resolution determination methods in cryo-EM including local and directional resolution variation, (ii) cryo-EM map contrast optimization including complementary map types that can help in identifying ambiguous density features, (iii) atomic model building and (iv) refinement in various resolution regimes including (v) their validation and (vi) discuss differences between X-ray and cryo-EM maps. Based on the methods originally developed for X-ray crystallography, the path from 3D image reconstruction to atomic coordinates has become an integral and important part of the cryo-EM structure determination work-flow that routinely delivers atomic models.

Validation tests for cryo-EM maps using an independent particle set

Cryo-electron microscopy (cryo-EM) has revolutionized structural biology by providing 3D density maps of biomolecules at near-atomic resolution. However, map validation is still an open issue. Despite several efforts from the community, it is possible to overfit 3D maps to noisy data. Here, we develop a novel methodology that uses a small independent particle set (not used during the 3D refinement) to validate the maps. The main idea is to monitor how the map probability evolves over the control set during the 3D refinement. The method is complementary to the gold-standard procedure, which generates two reconstructions at each iteration. We low-pass filter the two reconstructions for different frequency cutoffs, and we calculate the probability of each filtered map given the control set. For high-quality maps, the probability should increase as a function of the frequency cutoff and the refinement iteration. We also compute the similarity between the densities of probability distributions of the two reconstructions. As higher frequencies are included, the distributions become more dissimilar. We optimized the BioEM package to perform these calculations, and tested it over systems ranging from quality data to pure noise. Our results show that with our methodology, it possible to discriminate datasets that are constructed from noise particles. We conclude that validation against a control particle set provides a powerful tool to assess the quality of cryo-EM maps.

Cryo-EM structure of islet amyloid polypeptide fibrils reveals similarities with amyloid-β fibrils

Amyloid deposits consisting of fibrillar islet amyloid polypeptide (IAPP) in pancreatic islets are associated with beta-cell loss and have been implicated in type 2 diabetes (T2D). Here, we applied cryo-EM to reconstruct densities of three dominant IAPP fibril polymorphs, formed in vitro from synthetic human IAPP. An atomic model of the main polymorph, built from a density map of 4.2-Å resolution, reveals two S-shaped, intertwined protofilaments. The segment 21-NNFGAIL-27, essential for IAPP amyloidogenicity, forms the protofilament interface together with Tyr37 and the amidated C terminus. The S-fold resembles polymorphs of Alzheimer's disease (AD)-associated amyloid-β (Aβ) fibrils, which might account for the epidemiological link between T2D and AD and reports on IAPP-Aβ cross-seeding in vivo. The results structurally link the early-onset T2D IAPP genetic polymorphism (encoding Ser20Gly) with the AD Arctic mutation (Glu22Gly) of Aβ and support the design of inhibitors and imaging probes for IAPP fibrils.

Cryo-EM structure of islet amyloid polypeptide fibrils reveals similarities with amyloid-β fibrils

Amyloid deposits consisting of fibrillar islet amyloid polypeptide (IAPP) in pancreatic islets are associated with beta-cell loss and have been implicated in type 2 diabetes (T2D). Here, we applied cryo-EM to reconstruct densities of three dominant IAPP fibril polymorphs, formed in vitro from synthetic human IAPP. An atomic model of the main polymorph, built from a density map of 4.2-Å resolution, reveals two S-shaped, intertwined protofilaments. The segment 21-NNFGAIL-27, essential for IAPP amyloidogenicity, forms the protofilament interface together with Tyr37 and the amidated C terminus. The S-fold resembles polymorphs of Alzheimer's disease (AD)-associated amyloid-β (Aβ) fibrils, which might account for the epidemiological link between T2D and AD and reports on IAPP-Aβ cross-seeding in vivo. The results structurally link the early-onset T2D IAPP genetic polymorphism (encoding Ser20Gly) with the AD Arctic mutation (Glu22Gly) of Aβ and support the design of inhibitors and imaging probes for IAPP fibrils.

Simulation-Based Methods for Model Building and Refinement in Cryoelectron Microscopy

Advances in cryoelectron microscopy (cryo-EM) have revolutionized the structural investigation of large macromolecular assemblies. In this review, we first provide a broad overview of modeling methods used for flexible fitting of molecular models into cryo-EM density maps. We give special attention to approaches rooted in molecular simulations-atomistic molecular dynamics and Monte Carlo. Concise descriptions of the methods are given along with discussion of their advantages, limitations, and most popular alternatives. We also describe recent extensions of the widely used molecular dynamics flexible fitting (MDFF) method and discuss how different model-building techniques could be incorporated into new hybrid modeling schemes and simulation workflows. Finally, we provide two illustrative examples of model-building and refinement strategies employing MDFF, cascade MDFF, and RosettaCM. These examples come from recent cryo-EM studies that elucidated transcription preinitiation complexes and shed light on the functional roles of these assemblies in gene expression and gene regulation.

Need for Cross-Validation of Single Particle Cryo-EM

Cross-validation is used to determine the validity of a model on unseen data by assessing if the model is overfitted to noise. It is widely used in many fields, from artificial intelligence to structural biology in X-ray crystallography and nuclear magnetic resonance. Although there are concerns of map overfitting in cryo-electron microscopy (cryo-EM), cross-validation is rarely used. The problem is that establishing a performance metric of the maps over unseen data (given by 2D-projection images) is difficult due to the low signal-to-noise ratios in the individual particles. Here, I present recent advances for cryo-EM map reconstruction. I highlight that the gold-standard procedure can fail to detect map overfitting in certain cases, showing the necessity of assessing the map quality on unbiased data. Finally, I describe the challenges and advantages of developing a robust cross-validation methodology for cryo-EM.

Salt Dependence of A-Form RNA Duplexes: Structures and Implications

The biological functions of RNA range from gene regulation through catalysis and depend critically on its structure and flexibility. Conformational variations of flexible, non-base-paired components, including RNA hinges, bulges, or single-stranded tails, are well documented. Recent work has also identified variations in the structure of ubiquitous, base-paired duplexes found in almost all functional RNAs. Duplexes anchor the structures of folded RNAs, and their surface features are recognized by partner molecules. To date, no consistent picture has been obtained that describes the range of conformations assumed by RNA duplexes. Here, we apply wide angle, solution X-ray scattering (WAXS) to quantify these variations, by sampling length scales characteristic of helical geometries under different solution conditions. To identify the radius, helical rise, twist, and length of dsRNA helices, we exploit molecular dynamics generated structures, explicit solvent models, and ensemble optimization methods. Our results quantify the substantial and salt-dependent deviations of double-stranded (ds) RNA duplexes from the assumed canonical A-form conformation. Recent experiments underscore the need to properly describe the structures of RNA duplexes when interpreting the salt dependence of RNA conformations.

CryoEM-based hybrid modeling approaches for structure determination

Recent advances in cryo-electron microscopy (cryoEM) have dramatically improved the resolutions at which vitrified biological specimens can be studied, revealing new structural and mechanistic insights over a broad range of spatial scales. Bolstered by these advances, much effort has been directed toward the development of hybrid modeling methodologies for the construction and refinement of high-fidelity atomistic models from cryoEM data. In this brief review, we will survey the key elements of cryoEM-based hybrid modeling, providing an overview of available computational tools and strategies as well as several recent applications.

Variability of Protein Structure Models from Electron Microscopy

An increasing number of biomolecular structures are solved by electron microscopy (EM). However, the quality of structure models determined from EM maps vary substantially. To understand to what extent structure models are supported by information embedded in EM maps, we used two computational structure refinement methods to examine how much structures can be refined using a dataset of 49 maps with accompanying structure models. The extent of structure modification as well as the disagreement between refinement models produced by the two computational methods scaled inversely with the global and the local map resolutions. A general quantitative estimation of deviations of structures for particular map resolutions are provided. Our results indicate that the observed discrepancy between the deposited map and the refined models is due to the lack of structural information present in EM maps and thus these annotations must be used with caution for further applications.

Archaeal flagellin combines a bacterial type IV pilin domain with an Ig-like domain

The bacterial flagellar apparatus, which involves ∼40 different proteins, has been a model system for understanding motility and chemotaxis. The bacterial flagellar filament, largely composed of a single protein, flagellin, has been a model for understanding protein assembly. This system has no homology to the eukaryotic flagellum, in which the filament alone, composed of a microtubule-based axoneme, contains more than 400 different proteins. The archaeal flagellar system is simpler still, in some cases having ∼13 different proteins with a single flagellar filament protein. The archaeal flagellar system has no homology to the bacterial one and must have arisen by convergent evolution. However, it has been understood that the N-terminal domain of the archaeal flagellin is a homolog of the N-terminal domain of bacterial type IV pilin, showing once again how proteins can be repurposed in evolution for different functions. Using cryo-EM, we have been able to generate a nearly complete atomic model for a flagellar-like filament of the archaeon Ignicoccus hospitalis from a reconstruction at ∼4-Å resolution. We can now show that the archaeal flagellar filament contains a β-sandwich, previously seen in the FlaF protein that forms the anchor for the archaeal flagellar filament. In contrast to the bacterial flagellar filament, where the outer globular domains make no contact with each other and are not necessary for either assembly or motility, the archaeal flagellin outer domains make extensive contacts with each other that largely determine the interesting mechanical properties of these filaments, allowing these filaments to flex.

Tools for Model Building and Optimization into Near-Atomic Resolution Electron Cryo-Microscopy Density Maps

Electron cryo-microscopy (cryoEM) has advanced dramatically to become a viable tool for high-resolution structural biology research. The ultimate outcome of a cryoEM study is an atomic model of a macromolecule or its complex with interacting partners. This chapter describes a variety of algorithms and software to build a de novo model based on the cryoEM 3D density map, to optimize the model with the best stereochemistry restraints and finally to validate the model with proper protocols. The full process of atomic structure determination from a cryoEM map is described. The tools outlined in this chapter should prove extremely valuable in revealing atomic interactions guided by cryoEM data.

Refinement of atomic models in high resolution EM reconstructions using Flex-EM and local assessment

As the resolutions of Three Dimensional Electron Microscopic reconstructions of biological macromolecules are being improved, there is a need for better fitting and refinement methods at high resolutions and robust approaches for model assessment. Flex-EM/MODELLER has been used for flexible fitting of atomic models in intermediate-to-low resolution density maps of different biological systems. Here, we demonstrate the suitability of the method to successfully refine structures at higher resolutions (2.5-4.5Å) using both simulated and experimental data, including a newly processed map of Apo-GroEL. A hierarchical refinement protocol was adopted where the rigid body definitions are relaxed and atom displacement steps are reduced progressively at successive stages of refinement. For the assessment of local fit, we used the SMOC (segment-based Manders' overlap coefficient) score, while the model quality was checked using the Qmean score. Comparison of SMOC profiles at different stages of refinement helped in detecting regions that are poorly fitted. We also show how initial model errors can have significant impact on the goodness-of-fit. Finally, we discuss the implementation of Flex-EM in the CCP-EM software suite.

Multifrequency Force Microscopy of Helical Protein Assembly on a Virus

High-resolution microscopy techniques have been extensively used to investigate the structure of soft, biological matter at the nanoscale, from very thin membranes to small objects, like viruses. Electron microscopy techniques allow for obtaining extraordinary resolution by averaging signals from multiple identical structures. In contrast, atomic force microscopy (AFM) collects data from single entities. Here, it is possible to finely modulate the interaction with the samples, in order to be sensitive to their top surface, avoiding mechanical deformations. However, most biological surfaces are highly curved, such as fibers or tubes, and ultimate details of their surface are in the vicinity of steep height variations. This limits lateral resolution, even when sharp probes are used. We overcome this problem by using multifrequency force microscopy on a textbook example, the Tobacco Mosaic Virus (TMV). We achieved unprecedented resolution in local maps of amplitude and phase shift of the second excited mode, recorded together with sample topography. Our data, which combine multifrequency imaging and Fourier analysis, confirm the structure deduced from averaging techniques (XRD, cryoEM) for surface features of single virus particles, down to the helical pitch of the coat protein subunits, 2.3 nm. Remarkably, multifrequency AFM images do not require any image postprocessing.

Resolution and Probabilistic Models of Components in CryoEM Maps of Mature P22 Bacteriophage

CryoEM continues to produce density maps of larger and more complex assemblies with multiple protein components of mixed symmetries. Resolution is not always uniform throughout a cryoEM map, and it can be useful to estimate the resolution in specific molecular components of a large assembly. In this study, we present procedures to 1) estimate the resolution in subcomponents by gold-standard Fourier shell correlation (FSC); 2) validate modeling procedures, particularly at medium resolutions, which can include loop modeling and flexible fitting; and 3) build probabilistic models that combine high-accuracy priors (such as crystallographic structures) with medium-resolution cryoEM densities. As an example, we apply these methods to new cryoEM maps of the mature bacteriophage P22, reconstructed without imposing icosahedral symmetry. Resolution estimates based on gold-standard FSC show the highest resolution in the coat region (7.6 Å), whereas other components are at slightly lower resolutions: portal (9.2 Å), hub (8.5 Å), tailspike (10.9 Å), and needle (10.5 Å). These differences are indicative of inherent structural heterogeneity and/or reconstruction accuracy in different subcomponents of the map. Probabilistic models for these subcomponents provide new insights, to our knowledge, and structural information when taking into account uncertainty given the limitations of the observed density.

Ligand-induced structural changes in the cyclic nucleotide-modulated potassium channel MloK1

Cyclic nucleotide-modulated ion channels are important for signal transduction and pacemaking in eukaryotes. The molecular determinants of ligand gating in these channels are still unknown, mainly because of a lack of direct structural information. Here we report ligand-induced conformational changes in full-length MloK1, a cyclic nucleotide-modulated potassium channel from the bacterium Mesorhizobium loti, analysed by electron crystallography and atomic force microscopy. Upon cAMP binding, the cyclic nucleotide-binding domains move vertically towards the membrane, and directly contact the S1–S4 voltage sensor domains. This is accompanied by a significant shift and tilt of the voltage sensor domain helices. In both states, the inner pore-lining helices are in an ‘open’ conformation. We propose a mechanism in which ligand binding can favour pore opening via a direct interaction between the cyclic nucleotide-binding domains and voltage sensors. This offers a simple mechanistic hypothesis for the coupling between ligand gating and voltage sensing in eukaryotic HCN channels.

The molecular determinants underlying ligand gating of cyclic nucleotide-modulated ion channels remain unclear. Kowal et al. determine the conformational changes underlying cAMP binding to the bacterial channel MloK1, and propose a mechanism for coupling of ligand gating and voltage sensing in eukaryotic HCN channels.

EMRinger: side chain-directed model and map validation for 3D cryo-electron microscopy

Advances in high resolution electron cryomicroscopy (cryo-EM) have been accompanied by the development of validation metrics to independently assess map quality and model geometry. EMRinger assesses the precise fitting of an atomic model into the map during refinement and shows how radiation damage alters scattering from negatively charged amino acids. EMRinger will be useful for monitoring progress in resolving and modeling high-resolution features in cryo-EM.

Validation methods for low-resolution fitting of atomic structures to electron microscopy data

Fitting of atomic-resolution structures into reconstructions from electron cryo-microscopy is routinely used to understand the structure and function of macromolecular machines. Despite the fact that a plethora of fitting methods has been developed over recent years, standard protocols for quality assessment and validation of these fits have not been established. Here, we present the general concepts underlying current validation ideas as they relate to fitting of atomic-resolution models into electron cryo-microscopy reconstructions, with an emphasis on reconstructions with resolutions below the sub-nanometer range.

Validation of 3D EM Reconstructions: The Phantom in the Noise

Validation is a necessity to trust the structures solved by electron microscopy by single particle techniques. The impressive achievements in single particle reconstruction fuel its expansion beyond a small community of image processing experts. This poses the risk of inappropriate data processing with dubious results. Nowhere is it more clearly illustrated than in the recovery of a reference density map from pure noise aligned to that map—a phantom in the noise. Appropriate use of existing validating methods such as resolution-limited alignment and the processing of independent data sets (“gold standard”) avoid this pitfall. However, these methods can be undermined by biases introduced in various subtle ways. How can we test that a map is a coherent structure present in the images selected from the micrographs? In stead of viewing the phantom emerging from noise as a cautionary tale, it should be used as a defining baseline. Any map is always recoverable from noise images, provided a sufficient number of images are aligned and used in reconstruction. However, with smaller numbers of images, the expected coherence in the real particle images should yield better reconstructions than equivalent numbers of noise or background images, even without masking or imposing resolution limits as potential biases. The validation test proposed is therefore a simple alignment of a limited number of micrograph and noise images against the final reconstruction as reference, demonstrating that the micrograph images yield a better reconstruction. I examine synthetic cases to relate the resolution of a reconstruction to the alignment error as a function of the signal-to-noise ratio. I also administered the test to real cases of publicly available data. Adopting such a test can aid the microscopist in assessing the usefulness of the micrographs taken before committing to lengthy processing with questionable outcomes.

Tools for macromolecular model building and refinement into electron cryo-microscopy reconstructions

The recent rapid development of single-particle electron cryo-microscopy (cryo-EM) now allows structures to be solved by this method at resolutions close to 3 Å. Here, a number of tools to facilitate the interpretation of EM reconstructions with stereochemically reasonable all-atom models are described. The BALBES database has been repurposed as a tool for identifying protein folds from density maps. Modifications to Coot, including new Jiggle Fit and morphing tools and improved handling of nucleic acids, enhance its functionality for interpreting EM maps. REFMAC has been modified for optimal fitting of atomic models into EM maps. As external structural information can enhance the reliability of the derived atomic models, stabilize refinement and reduce overfitting, ProSMART has been extended to generate interatomic distance restraints from nucleic acid reference structures, and a new tool, LIBG, has been developed to generate nucleic acid base-pair and parallel-plane restraints. Furthermore, restraint generation has been integrated with visualization and editing in Coot, and these restraints have been applied to both real-space refinement in Coot and reciprocal-space refinement in REFMAC.

High-resolution modeling of protein structures based on flexible fitting of low-resolution structural data

To circumvent the difficulty of directly solving high-resolution biomolecular structures, low-resolution structural data from Cryo-electron microscopy (EM) and small angle solution X-ray scattering (SAXS) are increasingly used to explore multiple conformational states of biomolecular assemblies. One promising venue to obtain high-resolution structural models from low-resolution data is via data-constrained flexible fitting. To this end, we have developed a new method based on a coarse-grained Cα-only protein representation, and a modified form of the elastic network model (ENM) that allows large-scale conformational changes while maintaining the integrity of local structures including pseudo-bonds and secondary structures. Our method minimizes a pseudo-energy which linearly combines various terms of the modified ENM energy with an EM/SAXS-fitting score and a collision energy that penalizes steric collisions. Unlike some previous flexible fitting efforts using the lowest few normal modes, our method effectively utilizes all normal modes so that both global and local structural changes can be fully modeled with accuracy. This method is also highly efficient in computing time. We have demonstrated our method using adenylate kinase as a test case which undergoes a large open-to-close conformational change. The EM-fitting method is available at a web server (http://enm.lobos.nih.gov), and the SAXS-fitting method is available as a pre-compiled executable upon request.

Unified polymerization mechanism for the assembly of ASC-dependent inflammasomes

Inflammasomes elicit host defense inside cells by activating caspase-1 for cytokine maturation and cell death. AIM2 and NLRP3 are representative sensor proteins in two major families of inflammasomes. The adaptor protein ASC bridges the sensor proteins and caspase-1 to form ternary inflammasome complexes, achieved through pyrin domain (PYD) interactions between sensors and ASC and through caspase activation and recruitment domain (CARD) interactions between ASC and caspase-1. We found that PYD and CARD both form filaments. Activated AIM2 and NLRP3 nucleate PYD filaments of ASC, which, in turn, cluster the CARD of ASC. ASC thus nucleates CARD filaments of caspase-1, leading to proximity-induced activation. Endogenous NLRP3 inflammasome is also filamentous. The cryoelectron microscopy structure of ASC(PYD) filament at near-atomic resolution provides a template for homo- and hetero-PYD/PYD associations, as confirmed by structure-guided mutagenesis. We propose that ASC-dependent inflammasomes in both families share a unified assembly mechanism that involves two successive steps of nucleation-induced polymerization. PAPERFLICK:

M-free: scoring the reference bias in sub-tomogram averaging and template matching

Cryo-electron tomography provides a snapshot of the cellular proteome. With template matching, the spatial positions of various macromolecular complexes within their native cellular context can be detected. However, the growing awareness of the reference bias introduced by the cross-correlation based approaches, and more importantly the lack of a reliable confidence measurement in the selection of these macromolecular complexes, has restricted the use of these applications. Here we propose a heuristic, in which the reference bias is measured in real space in an analogous way to the R-free value in X-ray crystallography. We measure the reference bias within the mask used to outline the area of the template, and do not modify the template itself. The heuristic works by splitting the mask into a working and a testing area in a volume ratio of 9:1. While the working area is used during the calculation of the cross-correlation function, the information from both areas is explored to calculate the M-free score. We show using artificial data, that the M-free score gives a reliable measure for the reference bias. The heuristic can be applied in template matching and in sub-tomogram averaging. We further test the applicability of the heuristic in tomograms of purified macromolecules, and tomograms of whole Mycoplasma cells.

Finding the right fit: chiseling structures out of cryo-electron microscopy maps

Cryo-electron microscopy is a central tool for studying the architecture of macromolecular complexes at subnanometer resolution. Interpretation of an electron microscopy map requires its computational integration with data about the structure's components from all available sources, notably atomic models. Selecting a protocol for EM density-guided integrative structural modeling depends on the resolution and quality of the EM map as well as the available complimentary datasets. Here, we review rigid, flexible, and de novo integrative fitting into EM maps and provide guidelines and considerations for the design of modeling experiments. Finally, we discuss efforts towards establishing unified criteria for map and model assessment and validation.

Asymmetric perturbations of signalling oligomers

This review focuses on rapid and reversible noncovalent interactions for symmetric oligomers of signalling proteins. Symmetry mismatch, transient symmetry breaking and asymmetric perturbations via chemical (ligand binding) and physical (electric or mechanic) effects can initiate the signalling events. Advanced biophysical methods can reveal not only structural symmetries of stable membrane-bound signalling proteins but also asymmetric functional transition states. Relevant techniques amenable to distinguish between symmetric and asymmetric architectures are discussed including those with the capability of capturing low-populated transient conformational states. Typical examples of signalling proteins are overviewed for symmetry breaking in dimers (GPCRs, growth factor receptors, transcription factors); trimers (acid-sensing ion channels); tetramers (voltage-gated cation channels, ionotropic glutamate receptor, CNG and CHN channels); pentameric ligand-gated and mechanosensitive channels; higher order oligomers (gap junction channel, chaperonins, proteasome, virus capsid); as well as primary and secondary transporters. In conclusion, asymmetric perturbations seem to play important functional roles in a broad range of communicating networks.

The joys and perils of flexible fitting

While performing their functions, biological macromolecules often form large, dynamically changing macromolecular assemblies. Only a relatively small number of such assemblies have been accessible to the atomic-resolution techniques X-ray crystallography and NMR. Electron microscopy in conjunction with image reconstruction has become the preferred alternative for revealing the structures of such macromolecular complexes. However, for most assemblies the achievable resolution is too low to allow accurate atomic modeling directly from the data. Yet, useful models often can be obtained by fitting atomic models of individual components into a low-resolution reconstruction of the entire assembly. Several algorithms for achieving optimal fits in this context were developed recently, many allowing considerable degrees of flexibility to account for binding-induced conformational changes of the assembly components. This chapter describes the advantages and potential pitfalls of these methods and puts them into perspective with alternative approaches such as iterative modular fitting of rigid-body domains.

Building a pseudo-atomic model of the anaphase-promoting complex

The anaphase-promoting complex (APC/C) is a large E3 ubiquitin ligase that regulates progression through specific stages of the cell cycle by coordinating the ubiquitin-dependent degradation of cell-cycle regulatory proteins. Depending on the species, the active form of the APC/C consists of 14-15 different proteins that assemble into a 20-subunit complex with a mass of approximately 1.3 MDa. A hybrid approach of single-particle electron microscopy and protein crystallography of individual APC/C subunits has been applied to generate pseudo-atomic models of various functional states of the complex. Three approaches for assigning regions of the EM-derived APC/C density map to specific APC/C subunits are described. This information was used to dock atomic models of APC/C subunits, determined either by protein crystallography or homology modelling, to specific regions of the APC/C EM map, allowing the generation of a pseudo-atomic model corresponding to 80% of the entire complex.