Uncertainty in integrative structural modeling

Small Angle Scattering and Structural Biology: Data Quality and Model Validation

This chapter provides a brief review of the current state-of-the-art in small-angle scattering (SAS) from biomolecules in solution in regard to: (1) sample preparation and instrumentation, (2) data reduction and analysis, and (3) three-dimensional structural modelling and validation. In this context, areas of ongoing research in regard to the interpretation of SAS data will be discussed with a particular focus on structural modelling using computational methods and data from different experimental techniques, including SAS (hybrid methods). Finally, progress made in establishing community accepted publication guidelines and a standard reporting framework that includes SAS data deposition in a public data bank will be described. Importantly, SAS data with associated meta-data can now be held in a format that supports exchange between data archives and seamless interoperability with the world-wide Protein Data Bank (wwPDB). Biomolecular SAS is thus well positioned to contribute to an envisioned federation of data archives in support of hybrid structural biology.

Integrative structure modeling with the Integrative Modeling Platform

Building models of a biological system that are consistent with the myriad data available is one of the key challenges in biology. Modeling the structure and dynamics of macromolecular assemblies, for example, can give insights into how biological systems work, evolved, might be controlled, and even designed. Integrative structure modeling casts the building of structural models as a computational optimization problem, for which information about the assembly is encoded into a scoring function that evaluates candidate models. Here, we describe our open source software suite for integrative structure modeling, Integrative Modeling Platform (https://integrativemodeling.org), and demonstrate its use.

Coming to Grips with Ambiguity: Ion Mobility-Mass Spectrometry for Protein Quaternary Structure Assignment

Multiprotein complexes are central to our understanding of cellular biology, as they play critical roles in nearly every biological process. Despite many impressive advances associated with structural characterization techniques, large and highly-dynamic protein complexes are too often refractory to analysis by conventional, high-resolution approaches. To fill this gap, ion mobility-mass spectrometry (IM-MS) methods have emerged as a promising approach for characterizing the structures of challenging assemblies due in large part to the ability of these methods to characterize the composition, connectivity, and topology of large, labile complexes. In this Critical Insight, we present a series of bioinformatics studies aimed at assessing the information content of IM-MS datasets for building models of multiprotein structure. Our computational data highlights the limits of current coarse-graining approaches, and compelled us to develop an improved workflow for multiprotein topology modeling, which we benchmark against a subset of the multiprotein complexes within the PDB. This improved workflow has allowed us to ascertain both the minimal experimental restraint sets required for generation of high-confidence multiprotein topologies, and quantify the ambiguity in models where insufficient IM-MS information is available. We conclude by projecting the future of IM-MS in the context of protein quaternary structure assignment, where we predict that a more complete knowledge of the ultimate information content and ambiguity within such models will undoubtedly lead to applications for a broader array of challenging biomolecular assemblies. Graphical Abstract ᅟ.

Bayesian refinement of protein structures and ensembles against SAXS data using molecular dynamics

Small-angle X-ray scattering is an increasingly popular technique used to detect protein structures and ensembles in solution. However, the refinement of structures and ensembles against SAXS data is often ambiguous due to the low information content of SAXS data, unknown systematic errors, and unknown scattering contributions from the solvent. We offer a solution to such problems by combining Bayesian inference with all-atom molecular dynamics simulations and explicit-solvent SAXS calculations. The Bayesian formulation correctly weights the SAXS data versus prior physical knowledge, it quantifies the precision or ambiguity of fitted structures and ensembles, and it accounts for unknown systematic errors due to poor buffer matching. The method further provides a probabilistic criterion for identifying the number of states required to explain the SAXS data. The method is validated by refining ensembles of a periplasmic binding protein against calculated SAXS curves. Subsequently, we derive the solution ensembles of the eukaryotic chaperone heat shock protein 90 (Hsp90) against experimental SAXS data. We find that the SAXS data of the apo state of Hsp90 is compatible with a single wide-open conformation, whereas the SAXS data of Hsp90 bound to ATP or to an ATP-analogue strongly suggest heterogenous ensembles of a closed and a wide-open state.

Methods of probing the interactions between small molecules and disordered proteins

It is generally recognized that a large fraction of the human proteome is made up of proteins that remain disordered in their native states. Despite the fact that such proteins play key biological roles and are involved in many major human diseases, they still represent challenging targets for drug discovery. A major bottleneck for the identification of compounds capable of interacting with these proteins and modulating their disease-promoting behaviour is the development of effective techniques to probe such interactions. The difficulties in carrying out binding measurements have resulted in a poor understanding of the mechanisms underlying these interactions. In order to facilitate further methodological advances, here we review the most commonly used techniques to probe three types of interactions involving small molecules: (1) those that disrupt functional interactions between disordered proteins; (2) those that inhibit the aberrant aggregation of disordered proteins, and (3) those that lead to binding disordered proteins in their monomeric states. In discussing these techniques, we also point out directions for future developments.

2017 publication guidelines for structural modelling of small-angle scattering data from biomolecules in solution: an update

In 2012, preliminary guidelines were published addressing sample quality, data acquisition and reduction, presentation of scattering data and validation, and modelling for biomolecular small-angle scattering (SAS) experiments. Biomolecular SAS has since continued to grow and authors have increasingly adopted the preliminary guidelines. In parallel, integrative/hybrid determination of biomolecular structures is a rapidly growing field that is expanding the scope of structural biology. For SAS to contribute maximally to this field, it is essential to ensure open access to the information required for evaluation of the quality of SAS samples and data, as well as the validity of SAS-based structural models. To this end, the preliminary guidelines for data presentation in a publication are reviewed and updated, and the deposition of data and associated models in a public archive is recommended. These guidelines and recommendations have been prepared in consultation with the members of the International Union of Crystallography (IUCr) Small-Angle Scattering and Journals Commissions, the Worldwide Protein Data Bank (wwPDB) Small-Angle Scattering Validation Task Force and additional experts in the field.

The In Vivo Architecture of the Exocyst Provides Structural Basis for Exocytosis

The structural characterization of protein complexes in their native environment is challenging but crucial for understanding the mechanisms that mediate cellular processes. We developed an integrative approach to reconstruct the 3D architecture of protein complexes in vivo. We applied this approach to the exocyst, a hetero-octameric complex of unknown structure that is thought to tether secretory vesicles during exocytosis with a poorly understood mechanism. We engineered yeast cells to anchor the exocyst on defined landmarks and determined the position of its subunit termini at nanometer precision using fluorescence microscopy. We then integrated these positions with the structural properties of the subunits to reconstruct the exocyst together with a vesicle bound to it. The exocyst has an open hand conformation made of rod-shaped subunits that are interlaced in the core. The exocyst architecture explains how the complex can tether secretory vesicles, placing them in direct contact with the plasma membrane.

Accelerating physical simulations of proteins by leveraging external knowledge

It is challenging to compute structure-function relationships of proteins using molecular physics. The problem arises from the exponential scaling of the computational searching and sampling of large conformational spaces. This scaling challenge is not met by today's methods, such as Monte Carlo, simulated annealing, genetic algorithms, or molecular dynamics (MD) or its variants such as replica exchange. Such methods of searching for optimal states on complex probabalistic landscapes are referred to more broadly as Explore-and-Exploit (EE), including in contexts such as computational learning, games, industrial planning and modeling military strategies. Here we describe a Bayesian method, called MELD, that 'melds' together explore-and-exploit approaches with externally added information that can be vague, combinatoric, noisy, intuitive, heuristic, or from experimental data. MELD is shown to accelerate physical MD simulations when using experimental data to determine protein structures; for predicting protein structures by using heuristic directives; and when predicting binding affinities of proteins from limited information about the binding site. Such Guided Explore-and-Exploit approaches might also be useful beyond proteins and beyond molecular science.

Principles of protein structural ensemble determination

The biological functions of protein molecules are intimately dependent on their conformational dynamics. This aspect is particularly evident for disordered proteins, which constitute perhaps one-third of the human proteome. Therefore, structural ensembles often offer more useful representations of proteins than individual conformations. Here, we describe how the well-established principles of protein structure determination should be extended to the case of protein structural ensembles determination. These principles concern primarily how to deal with conformationally heterogeneous states, and with experimental measurements that are averaged over such states and affected by a variety of errors. We first review the growing literature of recent methods that combine experimental and computational information to model structural ensembles, highlighting their similarities and differences. We then address some conceptual problems in the determination of structural ensembles and define future goals towards the establishment of objective criteria for the comparison, validation, visualization and dissemination of such ensembles.

Quantitative FRET studies and integrative modeling unravel the structure and dynamics of biomolecular systems

Förster Resonance Energy Transfer (FRET) combined with single-molecule spectroscopy probes macromolecular structure and dynamics and identifies coexisting conformational states. We review recent methodological developments in integrative structural modeling by satisfying spatial restraints on networks of FRET pairs (hybrid-FRET). We discuss procedures to incorporate prior structural knowledge and to obtain optimal distance networks. Finally, a workflow for hybrid-FRET is presented that automates integrative structural modeling and experiment planning to put hybrid-FRET on rails. To test this workflow, we simulate realistic single-molecule experiments and resolve three protein conformers, exchanging at 30μs and 10ms, with accuracies of 1-3Å RMSD versus the target structure. Incorporation of data from other spectroscopies and imaging is also discussed.

Structure and Function of the Nuclear Pore Complex Cytoplasmic mRNA Export Platform

The last steps in mRNA export and remodeling are performed by the Nup82 complex, a large conserved assembly at the cytoplasmic face of the nuclear pore complex (NPC). By integrating diverse structural data, we have determined the molecular architecture of the native Nup82 complex at subnanometer precision. The complex consists of two compositionally identical multiprotein subunits that adopt different configurations. The Nup82 complex fits into the NPC through the outer ring Nup84 complex. Our map shows that this entire 14-MDa Nup82-Nup84 complex assembly positions the cytoplasmic mRNA export factor docking sites and messenger ribonucleoprotein (mRNP) remodeling machinery right over the NPC's central channel rather than on distal cytoplasmic filaments, as previously supposed. We suggest that this configuration efficiently captures and remodels exporting mRNP particles immediately upon reaching the cytoplasmic side of the NPC.

Insights into HIV-1 proviral transcription from integrative structure and dynamics of the Tat:AFF4:P-TEFb:TAR complex

HIV-1 Tat hijacks the human superelongation complex (SEC) to promote proviral transcription. Here we report the 5.9 Å structure of HIV-1 TAR in complex with HIV-1 Tat and human AFF4, CDK9, and CycT1. The TAR central loop contacts the CycT1 Tat-TAR recognition motif (TRM) and the second Tat Zn²⁺-binding loop. Hydrogen-deuterium exchange (HDX) shows that AFF4 helix 2 is stabilized in the TAR complex despite not touching the RNA, explaining how it enhances TAR binding to the SEC 50-fold. RNA SHAPE and SAXS data were used to help model the extended (Tat Arginine-Rich Motif) ARM, which enters the TAR major groove between the bulge and the central loop. The structure and functional assays collectively support an integrative structure and a bipartite binding model, wherein the TAR central loop engages the CycT1 TRM and compact core of Tat, while the TAR major groove interacts with the extended Tat ARM.

DOI:http://dx.doi.org/10.7554/eLife.15910.001

Bayesian ensemble refinement by replica simulations and reweighting

We describe different Bayesian ensemble refinement methods, examine their interrelation, and discuss their practical application. With ensemble refinement, the properties of dynamic and partially disordered (bio)molecular structures can be characterized by integrating a wide range of experimental data, including measurements of ensemble-averaged observables. We start from a Bayesian formulation in which the posterior is a functional that ranks different configuration space distributions. By maximizing this posterior, we derive an optimal Bayesian ensemble distribution. For discrete configurations, this optimal distribution is identical to that obtained by the maximum entropy "ensemble refinement of SAXS" (EROS) formulation. Bayesian replica ensemble refinement enhances the sampling of relevant configurations by imposing restraints on averages of observables in coupled replica molecular dynamics simulations. We show that the strength of the restraints should scale linearly with the number of replicas to ensure convergence to the optimal Bayesian result in the limit of infinitely many replicas. In the "Bayesian inference of ensembles" method, we combine the replica and EROS approaches to accelerate the convergence. An adaptive algorithm can be used to sample directly from the optimal ensemble, without replicas. We discuss the incorporation of single-molecule measurements and dynamic observables such as relaxation parameters. The theoretical analysis of different Bayesian ensemble refinement approaches provides a basis for practical applications and a starting point for further investigations.

Identifying and Visualizing Macromolecular Flexibility in Structural Biology

Structural biology comprises a variety of tools to obtain atomic resolution data for the investigation of macromolecules. Conventional structural methodologies including crystallography, NMR and electron microscopy often do not provide sufficient details concerning flexibility and dynamics, even though these aspects are critical for the physiological functions of the systems under investigation. However, the increasing complexity of the molecules studied by structural biology (including large macromolecular assemblies, integral membrane proteins, intrinsically disordered systems, and folding intermediates) continuously demands in-depth analyses of the roles of flexibility and conformational specificity involved in interactions with ligands and inhibitors. The intrinsic difficulties in capturing often subtle but critical molecular motions in biological systems have restrained the investigation of flexible molecules into a small niche of structural biology. Introduction of massive technological developments over the recent years, which include time-resolved studies, solution X-ray scattering, and new detectors for cryo-electron microscopy, have pushed the limits of structural investigation of flexible systems far beyond traditional approaches of NMR analysis. By integrating these modern methods with powerful biophysical and computational approaches such as generation of ensembles of molecular models and selective particle picking in electron microscopy, more feasible investigations of dynamic systems are now possible. Using some prominent examples from recent literature, we review how current structural biology methods can contribute useful data to accurately visualize flexibility in macromolecular structures and understand its important roles in regulation of biological processes.

Metadynamic metainference: Enhanced sampling of the metainference ensemble using metadynamics

Accurate and precise structural ensembles of proteins and macromolecular complexes can be obtained with metainference, a recently proposed Bayesian inference method that integrates experimental information with prior knowledge and deals with all sources of errors in the data as well as with sample heterogeneity. The study of complex macromolecular systems, however, requires an extensive conformational sampling, which represents a separate challenge. To address such challenge and to exhaustively and efficiently generate structural ensembles we combine metainference with metadynamics and illustrate its application to the calculation of the free energy landscape of the alanine dipeptide.

Unraveling low-resolution structural data of large biomolecules by constructing atomic models with experiment-targeted parallel cascade selection simulations

Various low-resolution experimental techniques have gained more and more popularity in obtaining structural information of large biomolecules. In order to interpret the low-resolution structural data properly, one may need to construct an atomic model of the biomolecule by fitting the data using computer simulations. Here we develop, to our knowledge, a new computational tool for such integrative modeling by taking the advantage of an efficient sampling technique called parallel cascade selection (PaCS) simulation. For given low-resolution structural data, this PaCS-Fit method converts it into a scoring function. After an initial simulation starting from a known structure of the biomolecule, the scoring function is used to pick conformations for next cycle of multiple independent simulations. By this iterative screening-after-sampling strategy, the biomolecule may be driven towards a conformation that fits well with the low-resolution data. Our method has been validated using three proteins with small-angle X-ray scattering data and two proteins with electron microscopy data. In all benchmark tests, high-quality atomic models, with generally 1-3 Å from the target structures, are obtained. Since our tool does not need to add any biasing potential in the simulations to deform the structure, any type of low-resolution data can be implemented conveniently.

Principles and Overview of Sampling Methods for Modeling Macromolecular Structure and Dynamics

Investigation of macromolecular structure and dynamics is fundamental to understanding how macromolecules carry out their functions in the cell. Significant advances have been made toward this end in silico, with a growing number of computational methods proposed yearly to study and simulate various aspects of macromolecular structure and dynamics. This review aims to provide an overview of recent advances, focusing primarily on methods proposed for exploring the structure space of macromolecules in isolation and in assemblies for the purpose of characterizing equilibrium structure and dynamics. In addition to surveying recent applications that showcase current capabilities of computational methods, this review highlights state-of-the-art algorithmic techniques proposed to overcome challenges posed in silico by the disparate spatial and time scales accessed by dynamic macromolecules. This review is not meant to be exhaustive, as such an endeavor is impossible, but rather aims to balance breadth and depth of strategies for modeling macromolecular structure and dynamics for a broad audience of novices and experts.

A Hybrid Approach to Structure and Function Modeling of G Protein-Coupled Receptors

The recent GPCR Dock 2013 assessment of serotonin receptor 5-HT1B and 5-HT2B, and smoothened receptor SMO targets, exposed the strengths and weaknesses of the currently used computational approaches. The test cases of 5-HT1B and 5-HT2B demonstrated that both the receptor structure and the ligand binding mode can be predicted with the atomic-detail accuracy, as long as the target-template sequence similarity is relatively high. On the other hand, the observation of a low target-template sequence similarity, e.g., between SMO from the frizzled GPCR family and members of the rhodopsin family, hampers the GPCR structure prediction and ligand docking. Indeed, in GPCR Dock 2013, accurate prediction of the SMO target was still beyond the capabilities of most research groups. Another bottleneck in the current GPCR research, as demonstrated by the 5-HT2B target, is the reliable prediction of global conformational changes induced by activation of GPCRs. In this work, we report details of our protocol used during GPCR Dock 2013. Our structure prediction and ligand docking protocol was especially successful in the case of 5-HT1B and 5-HT2B-ergotamine complexes for which we provide one of the most accurate predictions. In addition to a description of the GPCR Dock 2013 results, we propose a novel hybrid computational methodology to improve GPCR structure and function prediction. This computational methodology employs two separate rankings for filtering GPCR models. The first ranking is ligand-based while the second is based on the scoring scheme of the recently published BCL method. In this work, we prove that the use of knowledge-based potentials implemented in BCL is an efficient way to cope with major bottlenecks in the GPCR structure prediction. Thereby, we also demonstrate that the knowledge-based potentials for membrane proteins were significantly improved, because of the recent surge in available experimental structures.

Structure of γ-tubulin small complex based on a cryo-EM map, chemical cross-links, and a remotely related structure

Modeling protein complex structures based on distantly related homologues can be challenging due to poor sequence and structure conservation. Therefore, utilizing even low-resolution experimental data can significantly increase model precision and accuracy. Here, we present models of the two key functional states of the yeast γ-tubulin small complex (γTuSC): one for the low-activity "open" state and another for the higher-activity "closed" state. Both models were computed based on remotely related template structures and cryo-EM density maps at 6.9Å and 8.0Å resolution, respectively. For each state, extensive sampling of alignments and conformations was guided by the fit to the corresponding cryo-EM density map. The resulting good-scoring models formed a tightly clustered ensemble of conformations in most regions. We found significant structural differences between the two states, primarily in the γ-tubulin subunit regions where the microtubule binds. We also report a set of chemical cross-links that were found to be consistent with equilibrium between the open and closed states. The protocols developed here have been incorporated into our open-source Integrative Modeling Platform (IMP) software package (http://integrativemodeling.org), and can therefore be applied to many other systems.

Cryo-EM Data Are Superior to Contact and Interface Information in Integrative Modeling

Protein-protein interactions carry out a large variety of essential cellular processes. Cryo-electron microscopy (cryo-EM) is a powerful technique for the modeling of protein-protein interactions at a wide range of resolutions, and recent developments have caused a revolution in the field. At low resolution, cryo-EM maps can drive integrative modeling of the interaction, assembling existing structures into the map. Other experimental techniques can provide information on the interface or on the contacts between the monomers in the complex. This inevitably raises the question regarding which type of data is best suited to drive integrative modeling approaches. Systematic comparison of the prediction accuracy and specificity of the different integrative modeling paradigms is unavailable to date. Here, we compare EM-driven, interface-driven, and contact-driven integrative modeling paradigms. Models were generated for the protein docking benchmark using the ATTRACT docking engine and evaluated using the CAPRI two-star criterion. At 20 Å resolution, EM-driven modeling achieved a success rate of 100%, outperforming the other paradigms even with perfect interface and contact information. Therefore, even very low resolution cryo-EM data is superior in predicting heterodimeric and heterotrimeric protein assemblies. Our study demonstrates that a force field is not necessary, cryo-EM data alone is sufficient to accurately guide the monomers into place. The resulting rigid models successfully identify regions of conformational change, opening up perspectives for targeted flexible remodeling.

Protein Structural Analysis via Mass Spectrometry-Based Proteomics

Modern mass spectrometry (MS) technologies have provided a versatile platform that can be combined with a large number of techniques to analyze protein structure and dynamics. These techniques include the three detailed in this chapter: (1) hydrogen/deuterium exchange (HDX), (2) limited proteolysis, and (3) chemical crosslinking (CX). HDX relies on the change in mass of a protein upon its dilution into deuterated buffer, which results in varied deuterium content within its backbone amides. Structural information on surface exposed, flexible or disordered linker regions of proteins can be achieved through limited proteolysis, using a variety of proteases and only small extents of digestion. CX refers to the covalent coupling of distinct chemical species and has been used to analyze the structure, function and interactions of proteins by identifying crosslinking sites that are formed by small multi-functional reagents, termed crosslinkers. Each of these MS applications is capable of revealing structural information for proteins when used either with or without other typical high resolution techniques, including NMR and X-ray crystallography.

Assembly and Molecular Architecture of the Phosphoinositide 3-Kinase p85α Homodimer

Phosphoinositide 3-kinases (PI3Ks) are a family of lipid kinases that are activated by growth factor and G-protein-coupled receptors and propagate intracellular signals for growth, survival, proliferation, and metabolism. p85α, a modular protein consisting of five domains, binds and inhibits the enzymatic activity of class IA PI3K catalytic subunits. Here, we describe the structural states of the p85α dimer, based on data from in vivo and in vitro solution characterization. Our in vitro assembly and structural analyses have been enabled by the creation of cysteine-free p85α that is functionally equivalent to native p85α. Analytical ultracentrifugation studies showed that p85α undergoes rapidly reversible monomer-dimer assembly that is highly exothermic in nature. In addition to the documented SH3-PR1 dimerization interaction, we identified a second intermolecular interaction mediated by cSH2 domains at the C-terminal end of the polypeptide. We have demonstrated in vivo concentration-dependent dimerization of p85α using fluorescence fluctuation spectroscopy. Finally, we have defined solution conditions under which the protein is predominantly monomeric or dimeric, providing the basis for small angle x-ray scattering and chemical cross-linking structural analysis of the discrete dimer. These experimental data have been used for the integrative structure determination of the p85α dimer. Our study provides new insight into the structure and assembly of the p85α homodimer and suggests that this protein is a highly dynamic molecule whose conformational flexibility allows it to transiently associate with multiple binding proteins.

Molecular architecture of the yeast Mediator complex

The 21-subunit Mediator complex transduces regulatory information from enhancers to promoters, and performs an essential role in the initiation of transcription in all eukaryotes. Structural information on two-thirds of the complex has been limited to coarse subunit mapping onto 2-D images from electron micrographs. We have performed chemical cross-linking and mass spectrometry, and combined the results with information from X-ray crystallography, homology modeling, and cryo-electron microscopy by an integrative modeling approach to determine a 3-D model of the entire Mediator complex. The approach is validated by the use of X-ray crystal structures as internal controls and by consistency with previous results from electron microscopy and yeast two-hybrid screens. The model shows the locations and orientations of all Mediator subunits, as well as subunit interfaces and some secondary structural elements. Segments of 20–40 amino acid residues are placed with an average precision of 20 Å. The model reveals roles of individual subunits in the organization of the complex.

DOI:http://dx.doi.org/10.7554/eLife.08719.001

Inside a cell, proteins are made from instructions encoded by DNA. To produce a particular protein, a section of DNA within a gene is copied into a molecule of messenger ribonucleic acid (or mRNA). This process is called transcription and is carried out by an enzyme known as RNA polymerase.

Transcription begins in a region of DNA called a promoter, which is found at the start of the gene. RNA polymerase is brought to the DNA by many proteins, including the so-called Mediator complex. Mediator receives signals from within the cell and from the environment, processes the information, and instructs RNA polymerase whether to transcribe the gene or not. Mediator performs this important role in all organisms from yeast to humans, but it is not clear how it works. A crucial step towards the solution of this problem is to understand the three-dimensional structure of the complex.

Previous research using a technique called ‘electron microscopy’ showed that Mediator is composed of three modules, referred to as Head, Middle and Tail. The images from electron microscopy were not sufficiently detailed to reveal the organization of the proteins within these modules. An open-source Integrative Modeling Platform (IMP for short) was recently developed to arrive at structural models of large protein complexes from a combination of experimental data and computer models. Now, Robinson, Trnka, Pellarin et al. have used this platform to study the Mediator complex.

First, Robinson, Trnka, Pellarin et al. collected experimental data on the structure of the Mediator complex using two approaches called ‘chemical cross-linking’ and ‘mass spectrometry’. This data was combined with biochemical and structural information from previous studies to generate a three-dimensional model of the structure of the entire Mediator using IMP. The model is detailed enough to show the location and orientation of all the proteins in the complex. For example, a protein called Med17 connects the Head and Middle modules, while another subunit—known as Med14—spans the entire complex and makes extensive contacts with other proteins in all three modules.

DOI:http://dx.doi.org/10.7554/eLife.08719.002

Outcome of the First wwPDB Hybrid/Integrative Methods Task Force Workshop

Structures of biomolecular systems are increasingly computed by integrative modeling that relies on varied types of experimental data and theoretical information. We describe here the proceedings and conclusions from the first wwPDB Hybrid/Integrative Methods Task Force Workshop held at the European Bioinformatics Institute in Hinxton, UK, on October 6 and 7, 2014. At the workshop, experts in various experimental fields of structural biology, experts in integrative modeling and visualization, and experts in data archiving addressed a series of questions central to the future of structural biology. How should integrative models be represented? How should the data and integrative models be validated? What data should be archived? How should the data and models be archived? What information should accompany the publication of integrative models?

TEMPy: a Python library for assessment of three-dimensional electron microscopy density fits

TEMPy is an object-oriented Python library that provides the means to validate density fits in electron microscopy reconstructions. This article highlights several features of particular interest for this purpose and includes some customized examples.

Three-dimensional electron microscopy is currently one of the most promising techniques used to study macromolecular assemblies. Rigid and flexible fitting of atomic models into density maps is often essential to gain further insights into the assemblies they represent. Currently, tools that facilitate the assessment of fitted atomic models and maps are needed. TEMPy (template and electron microscopy comparison using Python) is a toolkit designed for this purpose. The library includes a set of methods to assess density fits in intermediate-to-low resolution maps, both globally and locally. It also provides procedures for single-fit assessment, ensemble generation of fits, clustering, and multiple and consensus scoring, as well as plots and output files for visualization purposes to help the user in analysing rigid and flexible fits. The modular nature of TEMPy helps the integration of scoring and assessment of fits into large pipelines, making it a tool suitable for both novice and expert structural biologists.