Real-space refinement with DireX: from global fitting to side-chain improvements

One bead per residue can describe all-atom protein structures

Atomistic resolution is the standard for high-resolution biomolecular structures, but experimental structural data are often at lower resolution. Coarse-grained models are also used extensively in computational studies to reach biologically relevant spatial and temporal scales. This study explores the use of advanced machine learning networks for reconstructing atomistic models from reduced representations. The main finding is that a single bead per amino acid residue allows construction of accurate and stereochemically realistic all-atom structures with minimal loss of information. This suggests that lower resolution representations of proteins may be sufficient for many applications when combined with a machine learning framework that encodes knowledge from known structures. Practical applications include the rapid addition of atomistic detail to low-resolution structures from experiment or computational coarse-grained models. The application of rapid, deterministic all-atom reconstruction within multi-scale frameworks is further demonstrated with a rapid protocol for the generation of accurate models from cryo-EM densities close to experimental structures.

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.

Automated simulation-based membrane protein refinement into cryo-EM data

The resolution revolution has increasingly enabled single-particle cryogenic electron microscopy (cryo-EM) reconstructions of previously inaccessible systems, including membrane proteins-a category that constitutes a disproportionate share of drug targets. We present a protocol for using density-guided molecular dynamics simulations to automatically refine atomistic models into membrane protein cryo-EM maps. Using adaptive force density-guided simulations as implemented in the GROMACS molecular dynamics package, we show how automated model refinement of a membrane protein is achieved without the need to manually tune the fitting force ad hoc. We also present selection criteria to choose the best-fit model that balances stereochemistry and goodness of fit. The proposed protocol was used to refine models into a new cryo-EM density of the membrane protein maltoporin, either in a lipid bilayer or detergent micelle, and we found that results do not substantially differ from fitting in solution. Fitted structures satisfied classical model-quality metrics and improved the quality and the model-to-map correlation of the x-ray starting structure. Additionally, the density-guided fitting in combination with generalized orientation-dependent all-atom potential was used to correct the pixel-size estimation of the experimental cryo-EM density map. This work demonstrates the applicability of a straightforward automated approach to fitting membrane protein cryo-EM densities. Such computational approaches promise to facilitate rapid refinement of proteins under different conditions or with various ligands present, including targets in the highly relevant superfamily of membrane proteins.

Gentle and fast all-atom model refinement to cryo-EM densities via a maximum likelihood approach

Better detectors and automated data collection have generated a flood of high-resolution cryo-EM maps, which in turn has renewed interest in improving methods for determining structure models corresponding to these maps. However, automatically fitting atoms to densities becomes difficult as their resolution increases and the refinement potential has a vast number of local minima. In practice, the problem becomes even more complex when one also wants to achieve a balance between a good fit of atom positions to the map, while also establishing good stereochemistry or allowing protein secondary structure to change during fitting. Here, we present a solution to this challenge using a maximum likelihood approach by formulating the problem as identifying the structure most likely to have produced the observed density map. This allows us to derive new types of smooth refinement potential-based on relative entropy-in combination with a novel adaptive force scaling algorithm to allow balancing of force-field and density-based potentials. In a low-noise scenario, as expected from modern cryo-EM data, the relative-entropy based refinement potential outperforms alternatives, and the adaptive force scaling appears to aid all existing refinement potentials. The method is available as a component in the GROMACS molecular simulation toolkit.

MarkovFit: Structure Fitting for Protein Complexes in Electron Microscopy Maps Using Markov Random Field

An increasing number of protein complex structures are determined by cryo-electron microscopy (cryo-EM). When individual protein structures have been determined and are available, an important task in structure modeling is to fit the individual structures into the density map. Here, we designed a method that fits the atomic structures of proteins in cryo-EM maps of medium to low resolutions using Markov random fields, which allows probabilistic evaluation of fitted models. The accuracy of our method, MarkovFit, performed better than existing methods on datasets of 31 simulated cryo-EM maps of resolution 10 Å , nine experimentally determined cryo-EM maps of resolution less than 4 Å , and 28 experimentally determined cryo-EM maps of resolution 6 to 20 Å .

EMDA: A Python package for Electron Microscopy Data Analysis

An open-source Python library EMDA for cryo-EM map and model manipulation is presented with a specific focus on validation. The use of several functionalities in the library is presented through several examples. The utility of local correlation as a metric for identifying map-model differences and unmodeled regions in maps, and how it is used as a metric of map-model validation is demonstrated. The mapping of local correlation to individual atoms, and its use to draw insights on local signal variations are discussed. EMDA’s likelihood-based map overlay is demonstrated by carrying out a superposition of two domains in two related structures. The overlay is carried out first to bring both maps into the same coordinate frame and then to estimate the relative movement of domains. Finally, the map magnification refinement in EMDA is presented with an example to highlight the importance of adjusting the map magnification in structural comparison studies.

Redeployment of automated MrBUMP search-model identification for map fitting in cryo-EM

In crystallography, the phase problem can often be addressed by the careful preparation of molecular-replacement search models. This has led to the development of pipelines such as MrBUMP that can automatically identify homologous proteins from an input sequence and edit them to focus on the areas that are most conserved. Many of these approaches can be applied directly to cryo-EM to help discover, prepare and correctly place models (here called cryo-EM search models) into electrostatic potential maps. This can significantly reduce the amount of manual model building that is required for structure determination. Here, MrBUMP is repurposed to fit automatically obtained PDB-derived chains and domains into cryo-EM maps. MrBUMP was successfully able to identify and place cryo-EM search models across a range of resolutions. Methods such as map segmentation are also explored as potential routes to improved performance. Map segmentation was also found to improve the effectiveness of the pipeline for higher resolution (<8 Å) data sets.

Cryo-EM model validation recommendations based on outcomes of the 2019 EMDataResource challenge

This paper describes outcomes of the 2019 Cryo-EM Model Challenge. The goals were to (1) assess the quality of models that can be produced from cryogenic electron microscopy (cryo-EM) maps using current modeling software, (2) evaluate reproducibility of modeling results from different software developers and users and (3) compare performance of current metrics used for model evaluation, particularly Fit-to-Map metrics, with focus on near-atomic resolution. Our findings demonstrate the relatively high accuracy and reproducibility of cryo-EM models derived by 13 participating teams from four benchmark maps, including three forming a resolution series (1.8 to 3.1 Å). The results permit specific recommendations to be made about validating near-atomic cryo-EM structures both in the context of individual experiments and structure data archives such as the Protein Data Bank. We recommend the adoption of multiple scoring parameters to provide full and objective annotation and assessment of the model, reflective of the observed cryo-EM map density.

Architecture of the flexible tail tube of bacteriophage SPP1

Bacteriophage SPP1 is a double-stranded DNA virus of the Siphoviridae family that infects the bacterium Bacillus subtilis. This family of phages features a long, flexible, non-contractile tail that has been difficult to characterize structurally. Here, we present the atomic structure of the tail tube of phage SPP1. Our hybrid structure is based on the integration of structural restraints from solid-state nuclear magnetic resonance (NMR) and a density map from cryo-EM. We show that the tail tube protein gp17.1 organizes into hexameric rings that are stacked by flexible linker domains and, thus, form a hollow flexible tube with a negatively charged lumen suitable for the transport of DNA. Additionally, we assess the dynamics of the system by combining relaxation measurements with variances in density maps.

Bacteriophages of the Siphoviridae family have a long, flexible, non-contractile tail that has been difficult to characterize structurally. Here, the authors present the atomic structure of the tail tube of one of these phages, showing a hollow flexible tube formed by hexameric rings stacked by flexible linkers.

Atomistic structure and dynamics of the human MHC-I peptide-loading complex

The major histocompatibility complex class-I (MHC-I) peptide-loading complex (PLC) is a cornerstone of the human adaptive immune system, being responsible for processing antigens that allow killer T cells to distinguish between healthy and compromised cells. Based on a recent low-resolution cryo-electron microscopy (cryo-EM) structure of this large membrane-bound protein complex, we report an atomistic model of the PLC and study its conformational dynamics on the multimicrosecond time scale using all-atom molecular dynamics (MD) simulations in an explicit lipid bilayer and water environment (1.6 million atoms in total). The PLC has a layered structure, with two editing modules forming a flexible protein belt surrounding a stable, catalytically active core. Tapasin plays a central role in the PLC, stabilizing the MHC-I binding groove in a conformation reminiscent of antigen-loaded MHC-I. The MHC-I-linked glycan steers a tapasin loop involved in peptide editing toward the binding groove. Tapasin conformational dynamics are also affected by calreticulin through a conformational selection mechanism that facilitates MHC-I recruitment into the complex.

Cryo-EM structure of the complete and ligand-saturated insulin receptor ectodomain

Glucose homeostasis and growth essentially depend on the hormone insulin engaging its receptor. Despite biochemical and structural advances, a fundamental contradiction has persisted in the current understanding of insulin ligand-receptor interactions. While biochemistry predicts two distinct insulin binding sites, 1 and 2, recent structural analyses have resolved only site 1. Using a combined approach of cryo-EM and atomistic molecular dynamics simulation, we present the structure of the entire dimeric insulin receptor ectodomain saturated with four insulin molecules. Complementing the previously described insulin-site 1 interaction, we present the first view of insulin bound to the discrete insulin receptor site 2. Insulin binding stabilizes the receptor ectodomain in a T-shaped conformation wherein the membrane-proximal domains converge and contact each other. These findings expand the current models of insulin binding to its receptor and of its regulation. In summary, we provide the structural basis for a comprehensive description of ligand-receptor interactions that ultimately will inform new approaches to structure-based drug design.

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.

Principles for Integrative Structural Biology Studies

Integrative structure determination is a powerful approach to modeling the structures of biological systems based on data produced by multiple experimental and theoretical methods, with implications for our understanding of cellular biology and drug discovery. This Primer introduces the theory and methods of integrative approaches, emphasizing the kinds of data that can be effectively included in developing models and using the nuclear pore complex as an example to illustrate the practice and challenges involved. These guidelines are intended to aid the researcher in understanding and applying integrative structural methods to systems of their interest and thus take advantage of this rapidly evolving field.

Advances in Structure Modeling Methods for Cryo-Electron Microscopy Maps

Cryo-electron microscopy (cryo-EM) has now become a widely used technique for structure determination of macromolecular complexes. For modeling molecular structures from density maps of different resolutions, many algorithms have been developed. These algorithms can be categorized into rigid fitting, flexible fitting, and de novo modeling methods. It is also observed that machine learning (ML) techniques have been increasingly applied following the rapid progress of the ML field. Here, we review these different categories of macromolecule structure modeling methods and discuss their advances over time.

Structural mechanism of a Rag GTPase activation checkpoint by the lysosomal folliculin complex

The tumor suppressor folliculin (FLCN) enables nutrient-dependent activation of the mechanistic target of rapamycin complex 1 (mTORC1) protein kinase via its guanosine triphosphatase (GTPase) activating protein (GAP) activity toward the GTPase RagC. Concomitant with mTORC1 inactivation by starvation, FLCN relocalizes from the cytosol to lysosomes. To determine the lysosomal function of FLCN, we reconstituted the human lysosomal FLCN complex (LFC) containing FLCN, its partner FLCN-interacting protein 2 (FNIP2), and the RagA^GDP:RagC^GTP GTPases as they exist in the starved state with their lysosomal anchor Ragulator complex and determined its cryo-electron microscopy structure to 3.6 angstroms. The RagC-GAP activity of FLCN was inhibited within the LFC, owing to displacement of a catalytically required arginine in FLCN from the RagC nucleotide. Disassembly of the LFC and release of the RagC-GAP activity of FLCN enabled mTORC1-dependent regulation of the master regulator of lysosomal biogenesis, transcription factor E3, implicating the LFC as a checkpoint in mTORC1 signaling.

Modelling structures in cryo-EM maps

Recent advances in structure determination of sub-cellular structures using cryo-electron microscopy and tomography have enabled us to understand their architecture in a more detailed manner and gain insight into their function. The choice of approach to use for atomic model building, fitting, refinement and validation in the 3D map resulting from these experiments depends primarily on the resolution of the map and the prior information on the corresponding model. Here, we survey some of such methods and approaches and highlight their uses in specific recent examples.

De novo main-chain modeling for EM maps using MAINMAST

An increasing number of protein structures are determined by cryo-electron microscopy (cryo-EM) at near atomic resolution. However, tracing the main-chains and building full-atom models from EM maps of ~4-5 Å is still not trivial and remains a time-consuming task. Here, we introduce a fully automated de novo structure modeling method, MAINMAST, which builds three-dimensional models of a protein from a near-atomic resolution EM map. The method directly traces the protein's main-chain and identifies Cα positions as tree-graph structures in the EM map. MAINMAST performs significantly better than existing software in building global protein structure models on data sets of 40 simulated density maps at 5 Å resolution and 30 experimentally determined maps at 2.6-4.8 Å resolution. In another benchmark of building missing fragments in protein models for EM maps, MAINMAST builds fragments of 11-161 residues long with an average RMSD of 2.68 Å.

Tools for the cryo-EM gold rush: going from the cryo-EM map to the atomistic model

As cryo-electron microscopy (cryo-EM) enters mainstream structural biology, the demand for fitting methods is high. Here, we review existing flexible fitting methods for cryo-EM. We discuss their importance, potential concerns and assessment strategies. We aim to give readers concrete descriptions of cryo-EM flexible fitting methods with corresponding examples.

Cryo-EM structures reveal specialization at the myosin VI-actin interface and a mechanism of force sensitivity

Despite extensive scrutiny of the myosin superfamily, the lack of high-resolution structures of actin-bound states has prevented a complete description of its mechanochemical cycle and limited insight into how sequence and structural diversification of the motor domain gives rise to specialized functional properties. Here we present cryo-EM structures of the unique minus-end directed myosin VI motor domain in rigor (4.6 Å) and Mg-ADP (5.5 Å) states bound to F-actin. Comparison to the myosin IIC-F-actin rigor complex reveals an almost complete lack of conservation of residues at the actin-myosin interface despite preservation of the primary sequence regions composing it, suggesting an evolutionary path for motor specialization. Additionally, analysis of the transition from ADP to rigor provides a structural rationale for force sensitivity in this step of the mechanochemical cycle. Finally, we observe reciprocal rearrangements in actin and myosin accompanying the transition between these states, supporting a role for actin structural plasticity during force generation by myosin VI.

A fragment based method for modeling of protein segments into cryo-EM density maps

BACKGROUND

Single-particle analysis of electron cryo-microscopy (cryo-EM) is a key technology for elucidation of macromolecular structures. Recent technical advances in hardware and software developments significantly enhanced the resolution of cryo-EM density maps and broadened the applicability and the circle of users. To facilitate modeling of macromolecules into cryo-EM density maps, fast and easy to use methods for modeling are now demanded.

RESULTS

Here we investigated and benchmarked the suitability of a classical and well established fragment-based approach for modeling of segments into cryo-EM density maps (termed FragFit). FragFit uses a hierarchical strategy to select fragments from a pre-calculated set of billions of fragments derived from structures deposited in the Protein Data Bank, based on sequence similarly, fit of stem atoms and fit to a cryo-EM density map. The user only has to specify the sequence of the segment and the number of the N- and C-terminal stem-residues in the protein. Using a representative data set of protein structures, we show that protein segments can be accurately modeled into cryo-EM density maps of different resolution by FragFit. Prediction quality depends on segment length, the type of secondary structure of the segment and local quality of the map.

CONCLUSION

Fast and automated calculation of FragFit renders it applicable for implementation of interactive web-applications e.g. to model missing segments, flexible protein parts or hinge-regions into cryo-EM density maps.

Flexible Fitting of Atomic Models into Cryo-EM Density Maps Guided by Helix Correspondences

Although electron cryo-microscopy (cryo-EM) has recently achieved resolutions of better than 3 Å, at which point molecular modeling can be done directly from the density map, analysis and annotation of a cryo-EM density map still primarily rely on fitting atomic or homology models to the density map. In this article, we present, to our knowledge, a new method for flexible fitting of known or modeled protein structures into cryo-EM density maps. Unlike existing methods that are guided by local density gradients, our method is guided by correspondences between the α-helices in the density map and model, and does not require an initial rigid-body fitting step. Compared with current methods on both simulated and experimental density maps, our method not only achieves greater accuracy for proteins with large deformations but also runs as fast or faster than many of the other flexible fitting routines.

Improved metrics for comparing structures of macromolecular assemblies determined by 3D electron-microscopy

Recent developments in 3-dimensional electron microcopy (3D-EM) techniques and a concomitant drive to look at complex molecular structures, have led to a rapid increase in the amount of volume data available for biomolecules. This creates a demand for better methods to analyse the data, including improved scores for comparison, classification and integration of data at different resolutions. To this end, we developed and evaluated a set of scoring functions that compare 3D-EM volumes. To test our scores we used a benchmark set of volume alignments derived from the Electron Microscopy Data Bank. We find that the performance of different scores vary with the map-type, resolution and the extent of overlap between volumes. Importantly, adding the overlap information to the local scoring functions can significantly improve their precision and accuracy in a range of resolutions. A combined score involving the local mutual information and overlap (LMI_OV) performs best overall, irrespective of the map category, resolution or the extent of overlap, and we recommend this score for general use. The local mutual information score itself is found to be more discriminatory than cross-correlation coefficient for intermediate-to-low resolution maps or when the map size and density distribution differ significantly. For comparing map surfaces, we implemented two filters to detect the surface points, including one based on the 'extent of surface exposure'. We show that scores that compare surfaces are useful at low resolutions and for maps with evident surface features. All the scores discussed are implemented in TEMPy (http://tempy.ismb.lon.ac.uk/).

Electron Cryo-microscopy as a Tool for Structure-Based Drug Development

For decades, X-ray crystallography and NMR have been the most important techniques for studying the atomic structure of macromolecules. However, as a result of size, instability, low yield, and other factors, many macromolecules are difficult to crystallize or unsuitable for NMR studies. Electron cryo-microscopy (cryo-EM) does not depend on crystals and has therefore been the method of choice for many macromolecular complexes that cannot be crystallized, but atomic resolution has mostly been beyond its reach. A new generation of detectors that are capable of sensing directly the incident electrons has recently revolutionized the field, with structures of macromolecules now routinely being solved to near-atomic resolution. In this review, we summarize some of the most recent examples of high-resolution cryo-EM structures. We put particular emphasis on proteins with pharmacological relevance that have traditionally been inaccessible to crystallography. Furthermore, we discuss examples where interactions with small molecules have been fully characterized at atomic resolution. Finally, we stress the current limits of cryo-EM, and methodological issues related to its usage as a tool for drug development.

The Central domain of RyR1 is the transducer for long-range allosteric gating of channel opening

The ryanodine receptors (RyRs) are intracellular calcium channels responsible for rapid release of Ca²⁺ from the sarcoplasmic/endoplasmic reticulum (SR/ER) to the cytoplasm, which is essential for the excitation-contraction (E-C) coupling of cardiac and skeletal muscles. The near-atomic resolution structure of closed RyR1 revealed the molecular details of this colossal channel, while the long-range allosteric gating mechanism awaits elucidation. Here, we report the cryo-EM structures of rabbit RyR1 in three closed conformations at about 4 Å resolution and an open state at 5.7 Å. Comparison of the closed RyR1 structures shows a breathing motion of the cytoplasmic platform, while the channel domain and its contiguous Central domain remain nearly unchanged. Comparison of the open and closed structures shows a dilation of the S6 tetrahelical bundle at the cytoplasmic gate that leads to channel opening. During the pore opening, the cytoplasmic “O-ring” motif of the channel domain and the U-motif of the Central domain exhibit coupled motion, while the Central domain undergoes domain-wise displacement. These structural analyses provide important insight into the E-C coupling in skeletal muscles and identify the Central domain as the transducer that couples the conformational changes of the cytoplasmic platform to the gating of the central pore.

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.

Using Molecular Simulation to Model High-Resolution Cryo-EM Reconstructions

An explosion of new data from high-resolution cryo-electron microscopy (cryo-EM) studies has produced a large number of data sets for many species of ribosomes in various functional states over the past few years. While many methods exist to produce structural models for lower resolution cryo-EM reconstructions, high-resolution reconstructions are often modeled using crystallographic techniques and extensive manual intervention. Here, we present an automated fitting technique for high-resolution cryo-EM data sets that produces all-atom models highly consistent with the EM density. Using a molecular dynamics approach, atomic positions are optimized with a potential that includes the cross-correlation coefficient between the structural model and the cryo-EM electron density, as well as a biasing potential preserving the stereochemistry and secondary structure of the biomolecule. Specifically, we use a hybrid structure-based/ab initio molecular dynamics potential to extend molecular dynamics fitting. In addition, we find that simulated annealing integration, as opposed to straightforward molecular dynamics integration, significantly improves performance. We obtain atomistic models of the human ribosome consistent with high-resolution cryo-EM reconstructions of the human ribosome. Automated methods such as these have the potential to produce atomistic models for a large number of ribosome complexes simultaneously that can be subsequently refined manually.

The importance of dynamics in integrative modeling of supramolecular assemblies

Revealing the atomistic architecture of supramolecular complexes is a fundamental step toward a deeper understanding of cellular functioning. To date, this formidable task is facilitated by an emerging array of integrative modeling approaches that combine experimental data from different sources. One major challenge these methods have to face is the treatment of the dynamic rearrangements of the individual subunits upon assembly. While this flexibility can be sampled at different levels, integrating native dynamic determinants with available experimental inputs can provide an effective way to reveal the molecular recognition mechanisms at the basis of supramolecular assembly.

A flipped ion pair at the dynein-microtubule interface is critical for dynein motility and ATPase activation

Dynein is a motor protein that moves on microtubules (MTs) using the energy of adenosine triphosphate (ATP) hydrolysis. To understand its motility mechanism, it is crucial to know how the signal of MT binding is transmitted to the ATPase domain to enhance ATP hydrolysis. However, the molecular basis of signal transmission at the dynein-MT interface remains unclear. Scanning mutagenesis of tubulin identified two residues in α-tubulin, R403 and E416, that are critical for ATPase activation and directional movement of dynein. Electron cryomicroscopy and biochemical analyses revealed that these residues form salt bridges with the residues in the dynein MT-binding domain (MTBD) that work in concert to induce registry change in the stalk coiled coil and activate the ATPase. The R403-E3390 salt bridge functions as a switch for this mechanism because of its reversed charge relative to other residues at the interface. This study unveils the structural basis for coupling between MT binding and ATPase activation and implicates the MTBD in the control of directional movement.

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.

Structure of the F-actin-tropomyosin complex

Filamentous actin (F-actin) is the major protein of muscle thin filaments, and actin microfilaments are the main component of the eukaryotic cytoskeleton. Mutations in different actin isoforms lead to early-onset autosomal dominant non-syndromic hearing loss, familial thoracic aortic aneurysms and dissections, and multiple variations of myopathies. In striated muscle fibres, the binding of myosin motors to actin filaments is mainly regulated by tropomyosin and troponin. Tropomyosin also binds to F-actin in smooth muscle and in non-muscle cells and stabilizes and regulates the filaments there in the absence of troponin. Although crystal structures for monomeric actin (G-actin) are available, a high-resolution structure of F-actin is still missing, hampering our understanding of how disease-causing mutations affect the function of thin muscle filaments and microfilaments. Here we report the three-dimensional structure of F-actin at a resolution of 3.7 Å in complex with tropomyosin at a resolution of 6.5 Å, determined by electron cryomicroscopy. The structure reveals that the D-loop is ordered and acts as a central region for hydrophobic and electrostatic interactions that stabilize the F-actin filament. We clearly identify map density corresponding to ADP and Mg(2+) and explain the possible effect of prominent disease-causing mutants. A comparison of F-actin with G-actin reveals the conformational changes during filament formation and identifies the D-loop as their key mediator. We also confirm that negatively charged tropomyosin interacts with a positively charged groove on F-actin. Comparison of the position of tropomyosin in F-actin-tropomyosin with its position in our previously determined F-actin-tropomyosin-myosin structure reveals a myosin-induced transition of tropomyosin. Our results allow us to understand the role of individual mutations in the genesis of actin- and tropomyosin-related diseases and will serve as a strong foundation for the targeted development of drugs.

Deformable elastic network refinement for low-resolution macromolecular crystallography

Crystals of membrane proteins and protein complexes often diffract to low resolution owing to their intrinsic molecular flexibility, heterogeneity or the mosaic spread of micro-domains. At low resolution, the building and refinement of atomic models is a more challenging task. The deformable elastic network (DEN) refinement method developed previously has been instrumental in the determinion of several structures at low resolution. Here, DEN refinement is reviewed, recommendations for its optimal usage are provided and its limitations are discussed. Representative examples of the application of DEN refinement to challenging cases of refinement at low resolution are presented. These cases include soluble as well as membrane proteins determined at limiting resolutions ranging from 3 to 7 Å. Potential extensions of the DEN refinement technique and future perspectives for the interpretation of low-resolution crystal structures are also discussed.

Cryo-EM study of the chromatin fiber reveals a double helix twisted by tetranucleosomal units

The hierarchical packaging of eukaryotic chromatin plays a central role in transcriptional regulation and other DNA-related biological processes. Here, we report the 11-angstrom-resolution cryogenic electron microscopy (cryo-EM) structures of 30-nanometer chromatin fibers reconstituted in the presence of linker histone H1 and with different nucleosome repeat lengths. The structures show a histone H1-dependent left-handed twist of the repeating tetranucleosomal structural units, within which the four nucleosomes zigzag back and forth with a straight linker DNA. The asymmetric binding and the location of histone H1 in chromatin play a role in the formation of the 30-nanometer fiber. Our results provide mechanistic insights into how nucleosomes compact into higher-order chromatin fibers.

iMODFIT: efficient and robust flexible fitting based on vibrational analysis in internal coordinates

Here, we employed the collective motions extracted from Normal Mode Analysis (NMA) in internal coordinates (torsional space) for the flexible fitting of atomic-resolution structures into electron microscopy (EM) density maps. The proposed methodology was validated using a benchmark of simulated cases, highlighting its robustness over the full range of EM resolutions and even over coarse-grained representations. A systematic comparison with other methods further showcased the advantages of this proposed methodology, especially at medium to lower resolutions. Using this method, computational costs and potential overfitting problems are naturally reduced by constraining the search in low-frequency NMA space, where covalent geometry is implicitly maintained. This method also effectively captures the macromolecular changes of a representative set of experimental test cases. We believe that this novel approach will extend the currently available EM hybrid methods to the atomic-level interpretation of large conformational changes and their functional implications.

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.

ΔΔPT: a comprehensive toolbox for the analysis of protein motion

Background

Normal Mode Analysis is one of the most successful techniques for studying motions in proteins and macromolecules. It can provide information on the mechanism of protein functions, used to aid crystallography and NMR data reconstruction, and calculate protein free energies.

Results

ΔΔPT is a toolbox allowing calculation of elastic network models and principle component analysis. It allows the analysis of pdb files or trajectories taken from; Gromacs, Amber, and DL_POLY. As well as calculation of the normal modes it also allows comparison of the modes with experimental protein motion, variation of modes with mutation or ligand binding, and calculation of molecular dynamic entropies.

Conclusions

This toolbox makes the respective tools available to a wide community of potential NMA users, and allows them unrivalled ability to analyse normal modes using a variety of techniques and current software.

Small-angle X-ray scattering of apolipoprotein A-IV reveals the importance of its termini for structural stability

ApoA-IV is an amphipathic protein that can emulsify lipids and has been linked to protective roles against cardiovascular disease and obesity. We previously reported an x-ray crystal structure of apoA-IV that was truncated at its N and C termini. Here, we have extended this work by demonstrating that self-associated states of apoA-IV are stable and can be structurally studied using small-angle x-ray scattering. Both the full-length monomeric and dimeric forms of apoA-IV were examined, with the dimer showing an elongated rod core with two nodes at opposing ends. The monomer is roughly half the length of the dimer with a single node. Small-angle x-ray scattering visualization of several deletion mutants revealed that removal of both termini can have substantial conformational effects throughout the molecule. Additionally, the F334A point mutation, which we previously showed increases apoA-IV lipid binding, also exhibited large conformational effects on the entire dimer. Merging this study's low-resolution structural information with the crystal structure provides insight on the conformation of apoA-IV as a monomer and as a dimer and further defines that a clasp mechanism may control lipid binding and, ultimately, protein function.

Consensus among multiple approaches as a reliability measure for flexible fitting into cryo-EM data

Cryo-electron microscopy (cryo-EM) can provide low-resolution density maps of large macromolecular assemblies. As the number of structures deposited in the Protein Data Bank by fitting a high-resolution structure into a low-resolution cryo-EM map is increasing, there is a need to revise the protocols and improve the measures for fitting. A recent study suggested using a combination of multiple automated flexible fitting approaches to improve the interpretation of cryo-EM data. The current work further explores the use of multiple approaches by validating this "consensus" fitting approach and deriving a local reliability measure. Here four different flexible fitting approaches are applied for fitting an initial structure into a simulated density map of known target structure from a dataset of proteins. It is found that the models produced from different approaches often have a consensus in conformation and are also near to the target structure, whereas cases not showing consensus are away from the target. A high correlation is also observed between the RMSF profiles calculated with respect to the average and the target structures, which indicates that the relation between consensus and accuracy can also be extended to a per-residue level. Therefore, the RMSF among the fitted models is proposed as a local reliability measure, which can be used to assess the reliability of the fit at specific regions. Hence, we encourage the community to use consensus flexible fitting with different methods to report on local reliability of the resulting models and improve the interpretation of cryo-EM data.

Atomic modeling of cryo-electron microscopy reconstructions--joint refinement of model and imaging parameters

When refining the fit of component atomic structures into electron microscopic reconstructions, use of a resolution-dependent atomic density function makes it possible to jointly optimize the atomic model and imaging parameters of the microscope. Atomic density is calculated by one-dimensional Fourier transform of atomic form factors convoluted with a microscope envelope correction and a low-pass filter, allowing refinement of imaging parameters such as resolution, by optimizing the agreement of calculated and experimental maps. A similar approach allows refinement of atomic displacement parameters, providing indications of molecular flexibility even at low resolution. A modest improvement in atomic coordinates is possible following optimization of these additional parameters. Methods have been implemented in a Python program that can be used in stand-alone mode for rigid-group refinement, or embedded in other optimizers for flexible refinement with stereochemical restraints. The approach is demonstrated with refinements of virus and chaperonin structures at resolutions of 9 through 4.5 Å, representing regimes where rigid-group and fully flexible parameterizations are appropriate. Through comparisons to known crystal structures, flexible fitting by RSRef is shown to be an improvement relative to other methods and to generate models with all-atom rms accuracies of 1.5-2.5 Å at resolutions of 4.5-6 Å.