Assembly of macromolecular complexes by satisfaction of spatial restraints from electron microscopy images

From integrative structural biology to cell biology

Integrative modeling is an increasingly important tool in structural biology, providing structures by combining data from varied experimental methods and prior information. As a result, molecular architectures of large, heterogeneous, and dynamic systems, such as the ∼52-MDa Nuclear Pore Complex, can be mapped with useful accuracy, precision, and completeness. Key challenges in improving integrative modeling include expanding model representations, increasing the variety of input data and prior information, quantifying a match between input information and a model in a Bayesian fashion, inventing more efficient structural sampling, as well as developing better model validation, analysis, and visualization. In addition, two community-level challenges in integrative modeling are being addressed under the auspices of the Worldwide Protein Data Bank (wwPDB). First, the impact of integrative structures is maximized by PDB-Development, a prototype wwPDB repository for archiving, validating, visualizing, and disseminating integrative structures. Second, the scope of structural biology is expanded by linking the wwPDB resource for integrative structures with archives of data that have not been generally used for structure determination but are increasingly important for computing integrative structures, such as data from various types of mass spectrometry, spectroscopy, optical microscopy, proteomics, and genetics. To address the largest of modeling problems, a type of integrative modeling called metamodeling is being developed; metamodeling combines different types of input models as opposed to different types of data to compute an output model. Collectively, these developments will facilitate the structural biology mindset in cell biology and underpin spatiotemporal mapping of the entire cell.

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.

SSEThread: Integrative threading of the DNA-PKcs sequence based on data from chemical cross-linking and hydrogen deuterium exchange

X-ray crystallography and electron microscopy maps resolved to 3-8 Å are generally sufficient for tracing the path of the polypeptide chain in space, while often insufficient for unambiguously registering the sequence on the path (i.e., threading). Frequently, however, additional information is available from other biophysical experiments, physical principles, statistical analyses, and other prior models. Here, we formulate an integrative approach for sequence assignment to a partial backbone model as an optimization problem, which requires three main components: the representation of the system, the scoring function, and the optimization method. The method is implemented in the open source Integrative Modeling Platform (IMP) (https://integrativemodeling.org), allowing a number of different terms in the scoring function. We apply this method to localizing the sequence assignment within a 199-residue disordered region of three structured and sequence unassigned helices in the DNA-PKcs crystallographic structure, using chemical crosslinks, hydrogen deuterium exchange, and sequence connectivity. The resulting ensemble of threading models provides two major solutions, one of which suggests that the crucial ABCDE cluster of phosphorylation sites cannot undergo intra-molecular autophosphorylation without a conformational rearrangement. The ensemble of solutions embodies the most accurate and precise sequence threading given the available information.

Archiving and disseminating integrative structure models

Limitations in the applicability, accuracy, and precision of individual structure characterization methods can sometimes be overcome via an integrative modeling approach that relies on information from all available sources, including all available experimental data and prior models. The open-source Integrative Modeling Platform (IMP) is one piece of software that implements all computational aspects of integrative modeling. To maximize the impact of integrative structures, the coordinates should be made publicly available, as is already the case for structures based on X-ray crystallography, NMR spectroscopy, and electron microscopy. Moreover, the associated experimental data and modeling protocols should also be archived, such that the original results can easily be reproduced. Finally, it is essential that the integrative structures are validated as part of their publication and deposition. A number of research groups have already developed software to implement integrative modeling and have generated a number of structures, prompting the formation of an Integrative/Hybrid Methods Task Force. Following the recommendations of this task force, the existing PDBx/mmCIF data representation used for atomic PDB structures has been extended to address the requirements for archiving integrative structural models. This IHM-dictionary adds a flexible model representation, including coarse graining, models in multiple states and/or related by time or other order, and multiple input experimental information sources. A prototype archiving system called PDB-Dev ( https://pdb-dev.wwpdb.org ) has also been created to archive integrative structural models, together with a Python library to facilitate handling of integrative models in PDBx/mmCIF format.

Single-Molecule 3D Images of "Hole-Hole" IgG1 Homodimers by Individual-Particle Electron Tomography

The engineering of immunoglobulin-G molecules (IgGs) is of wide interest for improving therapeutics, for example by modulating the activity or multiplexing the specificity of IgGs to recognize more than one antigen. Optimization of engineered IgG requires knowledge of three-dimensional (3D) structure of synthetic IgG. However, due to flexible nature of the molecules, their structural characterization is challenging. Here, we use our reported individual-particle electron tomography (IPET) method with optimized negative-staining (OpNS) for direct 3D reconstruction of individual IgG hole-hole homodimer molecules. The hole-hole homodimer is an undesired variant generated during the production of a bispecific antibody using the knob-into-hole heterodimer technology. A total of 64 IPET 3D density maps at ~15 Å resolutions were reconstructed from 64 individual molecules, revealing 64 unique conformations. In addition to the known Y-shaped conformation, we also observed an unusual X-shaped conformation. The 3D structure of the X-shaped conformation contributes to our understanding of the structural details of the interaction between two heavy chains in the Fc domain. The IPET approach, as an orthogonal technique to characterize the 3D structure of therapeutic antibodies, provides insight into the 3D structural variety and dynamics of heterogeneous IgG molecules.

Integrative Structure Modeling: Overview and Assessment

Integrative structure modeling computationally combines data from multiple sources of information with the aim of obtaining structural insights that are not revealed by any single approach alone. In the first part of this review, we survey the commonly used sources of structural information and the computational aspects of model building. Throughout the past decade, integrative modeling was applied to various biological systems, with a focus on large protein complexes. Recent progress in the field of cryo-electron microscopy (cryo-EM) has resolved many of these complexes to near-atomic resolution. In the second part of this review, we compare a range of published integrative models with their higher-resolution counterparts with the aim of critically assessing their accuracy. This comparison gives a favorable view of integrative modeling and demonstrates its ability to yield accurate and informative results. We discuss possible roles of integrative modeling in the new era of cryo-EM and highlight future challenges and directions.

Computational modeling of RNA 3D structure based on experimental data

RNA molecules are master regulators of cells. They are involved in a variety of molecular processes: they transmit genetic information, sense cellular signals and communicate responses, and even catalyze chemical reactions. As in the case of proteins, RNA function is dictated by its structure and by its ability to adopt different conformations, which in turn is encoded in the sequence. Experimental determination of high-resolution RNA structures is both laborious and difficult, and therefore the majority of known RNAs remain structurally uncharacterized. To address this problem, predictive computational methods were developed based on the accumulated knowledge of RNA structures determined so far, the physical basis of the RNA folding, and taking into account evolutionary considerations, such as conservation of functionally important motifs. However, all theoretical methods suffer from various limitations, and they are generally unable to accurately predict structures for RNA sequences longer than 100-nt residues unless aided by additional experimental data. In this article, we review experimental methods that can generate data usable by computational methods, as well as computational approaches for RNA structure prediction that can utilize data from experimental analyses. We outline methods and data types that can be potentially useful for RNA 3D structure modeling but are not commonly used by the existing software, suggesting directions for future development.

A Hybrid Hamiltonian for the Accelerated Sampling along Experimental Restraints

In this article, we present an enhanced sampling method based on a hybrid Hamiltonian which combines experimental distance restraints with a bias dependent from multiple path-dependent variables. This simulation method determines the bias-coordinates on the fly and does not require a priori knowledge about reaction coordinates. The hybrid Hamiltonian accelerates the sampling of proteins, and, combined with experimental distance information, the technique considers the restraints adaptively and in dependency of the system's intrinsic dynamics. We validate the methodology on the dipole relaxation of two water models and the conformational landscape of dialanine. Using experimental NMR-restraint data, we explore the folding landscape of the TrpCage mini-protein and in a second example apply distance restraints from chemical crosslinking/mass spectrometry experiments for the sampling of the conformation space of the Killer Cell Lectin-like Receptor Subfamily B Member 1A (NKR-P1A). The new methodology has the potential to adaptively introduce experimental restraints without affecting the conformational space of the system along an ergodic trajectory. Since only a limited number of input- and no-order parameters are required for the setup of the simulation, the method is broadly applicable and has the potential to be combined with coarse-graining methods.

Bayesian inference of rotor ring stoichiometry from electron microscopy images of archaeal ATP synthase

The 'Bayesian inference of electron microscopy' (BioEM) framework makes it possible to determine the stoichiometry of protein complexes using 3D coarse-grained models and a relatively small number of cryo-electron microscopy images as input. We applied the method to determine the most probable rotor ring stoichiometry of the archaeal Na+ ATP synthase from Pyrococcus furiosus, a multisubunit complex able to produce ATP under extreme conditions. Archaeal ATP synthases consist of a catalytic A1 part and a membrane-embedded AO portion. The AO portion is composed of a rotor ring and the a-subunit. The rotor ring of P. furiosus ATP synthase is composed of 16-kDa c-subunits, each consisting of four helices forming a bundle, with only one Na+-binding site per bundle. This ring was proposed to be decameric from LILBID-MS analysis of the entire ATP synthase. By contrast, the BioEM posterior favors a c9 ring stoichiometry. With BioEM, we ranked coarse-grained models of the whole complex with different ring geometry, using 6400 unprocessed particle images of the A1AO complex collected in vitreous ice. BioEM makes it possible to probabilistically establish the domain stoichiometry using low-resolution information and comparably few particle images.

Advances in image processing for single-particle analysis by electron cryomicroscopy and challenges ahead

Electron cryomicroscopy (cryoEM) is essential for the study and functional understanding of non-crystalline macromolecules such as proteins. These molecules cannot be imaged using X-ray crystallography or other popular methods. CryoEM has been successfully used to visualize macromolecular complexes such as ribosomes, viruses, and ion channels. Determination of structural models of these at various conformational states leads to insight on how these molecules function. Recent advances in imaging technology have given cryoEM a scientific rebirth. As a result of these technological advances image processing and analysis have yielded molecular structures at atomic resolution. Nevertheless there continue to be challenges in image processing, and in this article we will touch on the most essential in order to derive an accurate three-dimensional model from noisy projection images. Traditional approaches, such as k-means clustering for class averaging, will be provided as background. We will then highlight new approaches for each image processing subproblem, including a 3D reconstruction method for asymmetric molecules using just two projection images and deep learning algorithms for automated particle picking.

Bayesian Weighing of Electron Cryo-Microscopy Data for Integrative Structural Modeling

Cryo-electron microscopy (cryo-EM) has become a mainstream technique for determining the structures of complex biological systems. However, accurate integrative structural modeling has been hampered by the challenges in objectively weighing cryo-EM data against other sources of information due to the presence of random and systematic errors, as well as correlations, in the data. To address these challenges, we introduce a Bayesian scoring function that efficiently and accurately ranks alternative structural models of a macromolecular system based on their consistency with a cryo-EM density map as well as other experimental and prior information. The accuracy of this approach is benchmarked using complexes of known structure and illustrated in three applications: the structural determination of the GroEL/GroES, RNA polymerase II, and exosome complexes. The approach is implemented in the open-source Integrative Modeling Platform (http://integrativemodeling.org), thus enabling integrative structure determination by combining cryo-EM data with other sources of information.

Likelihood-based structural analysis of electron microscopy images

Likelihood-based analysis of single-particle electron microscopy images has contributed much to the recent improvements in resolution. By treating particle orientations and classes probabilistically, uncertainties in the reconstruction process are explicitly accounted for, and the risk of bias towards the initial model is diminished. As a result, the quality and reliability of the reconstructions have greatly improved at manageable computational cost. Likelihood-based analysis of electron microscopy images also offers a route to direct coordinate refinement for dynamic systems, as an alternative to 3D density reconstruction. Here, we review recent developments in the algorithms used for reconstructions of high-resolution maps, and in the integrative framework of combining likelihood methods with simulations to address conformational variability in cryo-electron microscopy.

The complex structure and function of Mediator

In eukaryotes, RNA polymerase II (pol II) transcribes all protein-coding genes and many noncoding RNAs. Whereas many factors contribute to the regulation of pol II activity, the Mediator complex is required for expression of most, if not all, pol II transcripts. Structural characterization of Mediator is challenging due to its large size (∼20 subunits in yeast and 26 subunits in humans) and conformational flexibility. However, recent studies have revealed structural details at higher resolution. Here, we summarize recent findings and place in context with previous results, highlighting regions within Mediator that are important for regulating its structure and function.

The molecular architecture of the yeast spindle pole body core determined by Bayesian integrative modeling

A model of the core of the yeast spindle pole body (SPB) was created by a Bayesian modeling approach that integrated a diverse data set of biophysical, biochemical, and genetic information. The model led to a proposed pathway for the assembly of Spc110, a protein related to pericentrin, and a mechanism for how calmodulin strengthens the SPB during mitosis.

Microtubule-organizing centers (MTOCs) form, anchor, and stabilize the polarized network of microtubules in a cell. The central MTOC is the centrosome that duplicates during the cell cycle and assembles a bipolar spindle during mitosis to capture and segregate sister chromatids. Yet, despite their importance in cell biology, the physical structure of MTOCs is poorly understood. Here we determine the molecular architecture of the core of the yeast spindle pole body (SPB) by Bayesian integrative structure modeling based on in vivo fluorescence resonance energy transfer (FRET), small-angle x-ray scattering (SAXS), x-ray crystallography, electron microscopy, and two-hybrid analysis. The model is validated by several methods that include a genetic analysis of the conserved PACT domain that recruits Spc110, a protein related to pericentrin, to the SPB. The model suggests that calmodulin can act as a protein cross-linker and Spc29 is an extended, flexible protein. The model led to the identification of a single, essential heptad in the coiled-coil of Spc110 and a minimal PACT domain. It also led to a proposed pathway for the integration of Spc110 into the SPB.

Reconstruction of 3D structures of MET antibodies from electron microscopy 2D class averages

Dynamics of three MET antibody constructs (IgG1, IgG2, and IgG4) and the IgG4-MET antigen complex was investigated by creating their atomic models with an integrative experimental and computational approach. In particular, we used two-dimensional (2D) Electron Microscopy (EM) images, image class averaging, homology modeling, Rapidly exploring Random Tree (RRT) structure sampling, and fitting of models to images, to find the relative orientations of antibody domains that are consistent with the EM images. We revealed that the conformational preferences of the constructs depend on the extent of the hinge flexibility. We also quantified how the MET antigen impacts on the conformational dynamics of IgG4. These observations allow to create testable hypothesis to investigate MET biology. Our protocol may also help describe structural diversity of other antigen systems at approximately 5 Å precision, as quantified by Root-Mean-Square Deviation (RMSD) among good-scoring models.

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.

Computational Refinement and Validation Protocol for Proteins with Large Variable Regions Applied to Model HIV Env Spike in CD4 and 17b Bound State

Envelope glycoprotein gp120 of HIV-1 possesses several variable regions; their precise structure has been difficult to establish. We report a new model of gp120, in complex with antibodies CD4 and 17b, complete with its variable regions. The model was produced by a computational protocol that uses cryo-electron microscopy (EM) maps, atomic-resolution structures of the core, and information on binding interactions. Our model has excellent fit with EMD: 5020, is stereochemically and energetically favorable, and has the expected binding interfaces. Comparison of the ternary arrangement of the loops in this model with those bound to PGT122 and PGV04 suggested a possible motion of the V1V2 away from the CCR5 binding site and toward CD4. Our study also revealed that the CD4-bound state of the V1V2 loop is not optimal for gp120 bound with several neutralizing antibodies.

The molecular architecture of the Dam1 kinetochore complex is defined by cross-linking based structural modelling

Accurate segregation of chromosomes during cell division is essential. The Dam1 complex binds kinetochores to microtubules and its oligomerization is required to form strong attachments. It is a key target of Aurora B kinase, which destabilizes erroneous attachments allowing subsequent correction. Understanding the roles and regulation of the Dam1 complex requires structural information. Here we apply cross-linking/mass spectrometry and structural modelling to determine the molecular architecture of the Dam1 complex. We find microtubule attachment is accompanied by substantial conformational changes, with direct binding mediated by the carboxy termini of Dam1p and Duo1p. Aurora B phosphorylation of Dam1p C terminus weakens direct interaction with the microtubule. Furthermore, the Dam1p amino terminus forms an interaction interface between Dam1 complexes, which is also disrupted by phosphorylation. Our results demonstrate that Aurora B inhibits both direct interaction with the microtubule and oligomerization of the Dam1 complex to drive error correction during mitosis.

Dissecting Fission Yeast Shelterin Interactions via MICro-MS Links Disruption of Shelterin Bridge to Tumorigenesis

Shelterin, a six-member complex, protects telomeres from nucleolytic attack and regulates their elongation by telomerase. Here, we have developed a strategy, called MICro-MS (Mapping Interfaces via Crosslinking-Mass Spectrometry), that combines crosslinking-mass spectrometry and phylogenetic analysis to identify contact sites within the complex. This strategy allowed identification of separation-of-function mutants of fission yeast Ccq1, Poz1, and Pot1 that selectively disrupt their respective interactions with Tpz1. The various telomere dysregulation phenotypes observed in these mutants further emphasize the critical regulatory roles of Tpz1-centered shelterin interactions in telomere homeostasis. Furthermore, the conservation between fission yeast Tpz1-Pot1 and human TPP1-POT1 interactions led us to map a human melanoma-associated POT1 mutation (A532P) to the TPP1-POT1 interface. Diminished TPP1-POT1 interaction caused by hPOT1-A532P may enable unregulated telomere extension, which, in turn, helps cancer cells to achieve replicative immortality. Therefore, our study reveals a connection between shelterin connectivity and tumorigenicity.

Determining protein complex structures based on a Bayesian model of in vivo Förster resonance energy transfer (FRET) data

The use of in vivo Förster resonance energy transfer (FRET) data to determine the molecular architecture of a protein complex in living cells is challenging due to data sparseness, sample heterogeneity, signal contributions from multiple donors and acceptors, unequal fluorophore brightness, photobleaching, flexibility of the linker connecting the fluorophore to the tagged protein, and spectral cross-talk. We addressed these challenges by using a Bayesian approach that produces the posterior probability of a model, given the input data. The posterior probability is defined as a function of the dependence of our FRET metric FRETR on a structure (forward model), a model of noise in the data, as well as prior information about the structure, relative populations of distinct states in the sample, forward model parameters, and data noise. The forward model was validated against kinetic Monte Carlo simulations and in vivo experimental data collected on nine systems of known structure. In addition, our Bayesian approach was validated by a benchmark of 16 protein complexes of known structure. Given the structures of each subunit of the complexes, models were computed from synthetic FRETR data with a distance root-mean-squared deviation error of 14 to 17 Å. The approach is implemented in the open-source Integrative Modeling Platform, allowing us to determine macromolecular structures through a combination of in vivo FRETR data and data from other sources, such as electron microscopy and chemical cross-linking.

Sequence co-evolution gives 3D contacts and structures of protein complexes

Protein-protein interactions are fundamental to many biological processes. Experimental screens have identified tens of thousands of interactions, and structural biology has provided detailed functional insight for select 3D protein complexes. An alternative rich source of information about protein interactions is the evolutionary sequence record. Building on earlier work, we show that analysis of correlated evolutionary sequence changes across proteins identifies residues that are close in space with sufficient accuracy to determine the three-dimensional structure of the protein complexes. We evaluate prediction performance in blinded tests on 76 complexes of known 3D structure, predict protein-protein contacts in 32 complexes of unknown structure, and demonstrate how evolutionary couplings can be used to distinguish between interacting and non-interacting protein pairs in a large complex. With the current growth of sequences, we expect that the method can be generalized to genome-wide elucidation of protein-protein interaction networks and used for interaction predictions at residue resolution.

Structural characterization by cross-linking reveals the detailed architecture of a coatomer-related heptameric module from the nuclear pore complex

Most cellular processes are orchestrated by macromolecular complexes. However, structural elucidation of these endogenous complexes can be challenging because they frequently contain large numbers of proteins, are compositionally and morphologically heterogeneous, can be dynamic, and are often of low abundance in the cell. Here, we present a strategy for the structural characterization of such complexes that has at its center chemical cross-linking with mass spectrometric readout. In this strategy, we isolate the endogenous complexes using a highly optimized sample preparation protocol and generate a comprehensive, high-quality cross-linking dataset using two complementary cross-linking reagents. We then determine the structure of the complex using a refined integrative method that combines the cross-linking data with information generated from other sources, including electron microscopy, X-ray crystallography, and comparative protein structure modeling. We applied this integrative strategy to determine the structure of the native Nup84 complex, a stable hetero-heptameric assembly (∼ 600 kDa), 16 copies of which form the outer rings of the 50-MDa nuclear pore complex (NPC) in budding yeast. The unprecedented detail of the Nup84 complex structure reveals previously unseen features in its pentameric structural hub and provides information on the conformational flexibility of the assembly. These additional details further support and augment the protocoatomer hypothesis, which proposes an evolutionary relationship between vesicle coating complexes and the NPC, and indicates a conserved mechanism by which the NPC is anchored in the nuclear envelope.

Integrative structure-function mapping of the nucleoporin Nup133 suggests a conserved mechanism for membrane anchoring of the nuclear pore complex

The nuclear pore complex (NPC) is the sole passageway for the transport of macromolecules across the nuclear envelope. Nup133, a major component in the essential Y-shaped Nup84 complex, is a large scaffold protein of the NPC's outer ring structure. Here, we describe an integrative modeling approach that produces atomic models for multiple states of Saccharomyces cerevisiae (Sc) Nup133, based on the crystal structures of the sequence segments and their homologs, including the related Vanderwaltozyma polyspora (Vp) Nup133 residues 55 to 502 (VpNup133(55-502)) determined in this study, small angle X-ray scattering profiles for 18 constructs of ScNup133 and one construct of VpNup133, and 23 negative-stain electron microscopy class averages of ScNup133(2-1157). Using our integrative approach, we then computed a multi-state structural model of the full-length ScNup133 and validated it with mutational studies and 45 chemical cross-links determined via mass spectrometry. Finally, the model of ScNup133 allowed us to annotate a potential ArfGAP1 lipid packing sensor (ALPS) motif in Sc and VpNup133 and discuss its potential significance in the context of the whole NPC; we suggest that ALPS motifs are scattered throughout the NPC's scaffold in all eukaryotes and play a major role in the assembly and membrane anchoring of the NPC in the nuclear envelope. Our results are consistent with a common evolutionary origin of Nup133 with membrane coating complexes (the protocoatomer hypothesis); the presence of the ALPS motifs in coatomer-like nucleoporins suggests an ancestral mechanism for membrane recognition present in early membrane coating complexes.

Protein modeling: what happened to the "protein structure gap"?

Computational modeling of three-dimensional macromolecular structures and complexes from their sequence has been a long-standing vision in structural biology. Over the last 2 decades, a paradigm shift has occurred: starting from a large "structure knowledge gap" between the huge number of protein sequences and small number of known structures, today, some form of structural information, either experimental or template-based models, is available for the majority of amino acids encoded by common model organism genomes. With the scientific focus of interest moving toward larger macromolecular complexes and dynamic networks of interactions, the integration of computational modeling methods with low-resolution experimental techniques allows the study of large and complex molecular machines. One of the open challenges for computational modeling and prediction techniques is to convey the underlying assumptions, as well as the expected accuracy and structural variability of a specific model, which is crucial to understanding its limitations.

Conformational States of macromolecular assemblies explored by integrative structure calculation

A detailed description of macromolecular assemblies in multiple conformational states can be very valuable for understanding cellular processes. At present, structural determination of most assemblies in different biologically relevant conformations cannot be achieved by a single technique and thus requires an integrative approach that combines information from multiple sources. Different techniques require different computational methods to allow efficient and accurate data processing and analysis. Here, we summarize the latest advances and future challenges in computational methods that help the interpretation of data from two techniques—mass spectrometry and three-dimensional cryo-electron microscopy (with focus on alignment and classification of heterogeneous subtomograms from cryo-electron tomography). We evaluate how new developments in these two broad fields will lead to further integration with atomic structures to broaden our picture of the dynamic behavior of assemblies in their native environment.

Modeling of proteins and their assemblies with the Integrative Modeling Platform

To understand the workings of the living cell, we need to characterize protein assemblies that constitute the cell (for example, the ribosome, 26S proteasome, and the nuclear pore complex). A reliable high-resolution structural characterization of these assemblies is frequently beyond the reach of current experimental methods, such as X-ray crystallography, NMR spectroscopy, electron microscopy, footprinting, chemical cross-linking, FRET spectroscopy, small angle X-ray scattering, and proteomics. However, the information garnered from different methods can be combined and used to build models of the assembly structures that are consistent with all of the available datasets, and therefore more accurate, precise, and complete. Here, we describe a protocol for this integration, whereby the information is converted to a set of spatial restraints and a variety of optimization procedures can be used to generate models that satisfy the restraints as well as possible. These generated models can then potentially inform about the precision and accuracy of structure determination, the accuracy of the input datasets, and further data generation. We also demonstrate the Integrative Modeling Platform (IMP) software, which provides the necessary computational framework to implement this protocol, and several applications for specific use cases.

Bayesian analysis of individual electron microscopy images: towards structures of dynamic and heterogeneous biomolecular assemblies

We develop a method to extract structural information from electron microscopy (EM) images of dynamic and heterogeneous molecular assemblies. To overcome the challenge of disorder in the imaged structures, we analyze each image individually, avoiding information loss through clustering or averaging. The Bayesian inference of EM (BioEM) method uses a likelihood-based probabilistic measure to quantify the consistency between each EM image and given structural models. The likelihood function accounts for uncertainties in the molecular position and orientation, variations in the relative intensities and noise in the experimental images. The BioEM formalism is physically intuitive and mathematically simple. We show that for experimental GroEL images, BioEM correctly identifies structures according to the functional state. The top-ranked structure is the corresponding X-ray crystal structure, followed by an EM structure generated previously from a superset of the EM images used here. To analyze EM images of highly flexible molecules, we propose an ensemble refinement procedure, and validate it with synthetic EM maps of the ESCRT-I-II supercomplex. Both the size of the ensemble and its structural members are identified correctly. BioEM offers an alternative to 3D-reconstruction methods, extracting accurate population distributions for highly flexible structures and their assemblies. We discuss limitations of the method, and possible applications beyond ensemble refinement, including the cross-validation and unbiased post-assessment of model structures, and the structural characterization of systems where traditional approaches fail. Overall, our results suggest that the BioEM framework can be used to analyze EM images of both ordered and disordered molecular systems.

Structural templates for modeling homodimers

Oligomeric proteins are more abundant in nature than monomeric proteins, and involved in all biological processes. In the absence of an experimental structure, their subunits can be modeled from their sequence like monomeric proteins, but reliable procedures to build the oligomeric assembly are scarce. Template-based methods, which start from known protein structures, are commonly applied to model subunits. We present a method to model homodimers that relies on a structural alignment of the subunits, and test it on a set of 511 target structures recently released by the Protein Data Bank, taking as templates the earlier released structures of 3108 homodimeric proteins (H-set), and 2691 monomeric proteins that form dimer-like assemblies in crystals (M-set). The structural alignment identifies a H-set template for 97% of the targets, and in half of the cases, it yields a correct model of the dimer geometry and residue-residue contacts in the target. It also identifies a M-set template for most of the targets, and some of the crystal dimers are very similar to the target homodimers. The procedure efficiently detects homology at low levels of sequence identities, and points to erroneous quaternary structures in the Protein Data Bank. The high coverage of the target set suggests that the content of the Protein Data Bank already approaches the structural diversity of protein assemblies in nature, and that template-based methods should become the choice method for modeling oligomeric as well as monomeric proteins.

Characterization of the free-energy landscapes of proteins by NMR-guided metadynamics

The use of free-energy landscapes rationalizes a wide range of aspects of protein behavior by providing a clear illustration of the different states accessible to these molecules, as well as of their populations and pathways of interconversion. The determination of the free-energy landscapes of proteins by computational methods is, however, very challenging as it requires an extensive sampling of their conformational spaces. We describe here a technique to achieve this goal with relatively limited computational resources by incorporating nuclear magnetic resonance (NMR) chemical shifts as collective variables in metadynamics simulations. As in this approach the chemical shifts are not used as structural restraints, the resulting free-energy landscapes correspond to the force fields used in the simulations. We illustrate this approach in the case of the third Ig-binding domain of protein G from streptococcal bacteria (GB3). Our calculations reveal the existence of a folding intermediate of GB3 with nonnative structural elements. Furthermore, the availability of the free-energy landscape enables the folding mechanism of GB3 to be elucidated by analyzing the conformational ensembles corresponding to the native, intermediate, and unfolded states, as well as the transition states between them. Taken together, these results show that, by incorporating experimental data as collective variables in metadynamics simulations, it is possible to enhance the sampling efficiency by two or more orders of magnitude with respect to standard molecular dynamics simulations, and thus to estimate free-energy differences among the different states of a protein with a k(B)T accuracy by generating trajectories of just a few microseconds.

Structure, dynamics, evolution, and function of a major scaffold component in the nuclear pore complex

The nuclear pore complex, composed of proteins termed nucleoporins (Nups), is responsible for nucleocytoplasmic transport in eukaryotes. Nuclear pore complexes (NPCs) form an annular structure composed of the nuclear ring, cytoplasmic ring, a membrane ring, and two inner rings. Nup192 is a major component of the NPC's inner ring. We report the crystal structure of Saccharomyces cerevisiae Nup192 residues 2-960 [ScNup192(2-960)], which adopts an α-helical fold with three domains (i.e., D1, D2, and D3). Small angle X-ray scattering and electron microscopy (EM) studies reveal that ScNup192(2-960) could undergo long-range transition between "open" and "closed" conformations. We obtained a structural model of full-length ScNup192 based on EM, the structure of ScNup192(2-960), and homology modeling. Evolutionary analyses using the ScNup192(2-960) structure suggest that NPCs and vesicle-coating complexes are descended from a common membrane-coating ancestral complex. We show that suppression of Nup192 expression leads to compromised nuclear transport and hypothesize a role for Nup192 in modulating the permeability of the NPC central channel.