Uncertainty in integrative structural modeling

Advancing Molecular Dynamics: Toward Standardization, Integration, and Data Accessibility in Structural Biology

Molecular dynamics (MD) simulations have become a valuable tool in structural biology, offering insights into complex biological systems that are difficult to obtain through experimental techniques alone. The lack of available data sets and structures in most published computational work has limited other researchers' use of these models. In recent years, the emergence of online sharing platforms and MD database initiatives favor the deposition of ensembles and structures to accompany publications, favoring reuse of the data sets. However, the lack of uniform metadata collection, formats, and what data are deposited limits the impact and its use by different communities that are not necessarily experts in MD. This Perspective highlights the need for standardization and better resource sharing for processing and interpreting MD simulation results, akin to efforts in other areas of structural biology. As the field moves forward, we will see an increase in popularity and benefits of MD-based integrative approaches combining experimental data and simulations through probabilistic reasoning, but these too are limited by uniformity in experimental data availability and choices on how the data are modeled that are not trivial to decipher from papers. Other fields have addressed similar challenges comprehensively by establishing task forces with different degrees of success. The large scope and number of communities to represent the breadth of types of MD simulations complicates a parallel approach that would fit all. Thus, each group typically decides what data and which format to upload on servers like Zenodo. Uploading data with FAIR (findable, accessible, interoperable, reusable) principles in mind including optimal metadata collection will make the data more accessible and actionable by the community. Such a wealth of simulation data will foster method development and infrastructure advancements, thus propelling the field forward.

Optimizing representations for integrative structural modeling using Bayesian model selection

MOTIVATION

Integrative structural modeling combines data from experiments, physical principles, statistics of previous structures, and prior models to obtain structures of macromolecular assemblies that are challenging to characterize experimentally. The choice of model representation is a key decision in integrative modeling, as it dictates the accuracy of scoring, efficiency of sampling, and resolution of analysis. But currently, the choice is usually made ad hoc, manually.

RESULTS

Here, we report NestOR (Nested Sampling for Optimizing Representation), a fully automated, statistically rigorous method based on Bayesian model selection to identify the optimal coarse-grained representation for a given integrative modeling setup. Given an integrative modeling setup, it determines the optimal representations from given candidate representations based on their model evidence and sampling efficiency. The performance of NestOR was evaluated on a benchmark of four macromolecular assemblies.

AVAILABILITY AND IMPLEMENTATION

NestOR is implemented in the Integrative Modeling Platform (https://integrativemodeling.org) and is available at https://github.com/isblab/nestor. Data for the benchmark is at https://www.doi.org/10.5281/zenodo.10360718.

CombFold: predicting structures of large protein assemblies using a combinatorial assembly algorithm and AlphaFold2

Deep learning models, such as AlphaFold2 and RosettaFold, enable high-accuracy protein structure prediction. However, large protein complexes are still challenging to predict due to their size and the complexity of interactions between multiple subunits. Here we present CombFold, a combinatorial and hierarchical assembly algorithm for predicting structures of large protein complexes utilizing pairwise interactions between subunits predicted by AlphaFold2. CombFold accurately predicted (TM-score >0.7) 72% of the complexes among the top-10 predictions in two datasets of 60 large, asymmetric assemblies. Moreover, the structural coverage of predicted complexes was 20% higher compared to corresponding Protein Data Bank entries. We applied the method on complexes from Complex Portal with known stoichiometry but without known structure and obtained high-confidence predictions. CombFold supports the integration of distance restraints based on crosslinking mass spectrometry and fast enumeration of possible complex stoichiometries. CombFold's high accuracy makes it a promising tool for expanding structural coverage beyond monomeric proteins.

Integrative spatiotemporal map of nucleocytoplasmic transport

The Nuclear Pore Complex (NPC) facilitates rapid and selective nucleocytoplasmic transport of molecules as large as ribosomal subunits and viral capsids. It is not clear how key emergent properties of this transport arise from the system components and their interactions. To address this question, we constructed an integrative coarse-grained Brownian dynamics model of transport through a single NPC, followed by coupling it with a kinetic model of Ran-dependent transport in an entire cell. The microscopic model parameters were fitted to reflect experimental data and theoretical information regarding the transport, without making any assumptions about its emergent properties. The resulting reductionist model is validated by reproducing several features of transport not used for its construction, such as the morphology of the central transporter, rates of passive and facilitated diffusion as a function of size and valency, in situ radial distributions of pre-ribosomal subunits, and active transport rates for viral capsids. The model suggests that the NPC functions essentially as a virtual gate whose flexible phenylalanine-glycine (FG) repeat proteins raise an entropy barrier to diffusion through the pore. Importantly, this core functionality is greatly enhanced by several key design features, including 'fuzzy' and transient interactions, multivalency, redundancy in the copy number of FG nucleoporins, exponential coupling of transport kinetics and thermodynamics in accordance with the transition state theory, and coupling to the energy-reliant RanGTP concentration gradient. These design features result in the robust and resilient rate and selectivity of transport for a wide array of cargo ranging from a few kilodaltons to megadaltons in size. By dissecting these features, our model provides a quantitative starting point for rationally modulating the transport system and its artificial mimics.

Optimizing representations for integrative structural modeling using Bayesian model selection

Motivation

Integrative structural modeling combines data from experiments, physical principles, statistics of previous structures, and prior models to obtain structures of macromolecular assemblies that are challenging to characterize experimentally. The choice of model representation is a key decision in integrative modeling, as it dictates the accuracy of scoring, efficiency of sampling, and resolution of analysis. But currently, the choice is usually made ad hoc, manually.

Results

Here, we report NestOR (Nested Sampling for Optimizing Representation), a fully automated, statistically rigorous method based on Bayesian model selection to identify the optimal coarse-grained representation for a given integrative modeling setup. Given an integrative modeling setup, it determines the optimal representations from given candidate representations based on their model evidence and sampling efficiency. The performance of NestOR was evaluated on a benchmark of four macromolecular assemblies.

Availability

NestOR is implemented in the Integrative Modeling Platform (https://integrativemodeling.org) and is available at https://github.com/isblab/nestor.

Modern Alchemical Free Energy Methods for Drug Discovery Explained

This Perspective provides a contextual explanation of the current state-of-the-art alchemical free energy methods and their role in drug discovery as well as highlights select emerging technologies. The narrative attempts to answer basic questions about what goes on "under the hood" in free energy simulations and provide general guidelines for how to run simulations and analyze the results. It is the hope that this work will provide a valuable introduction to students and scientists in the field.

Hybrid computational methods combining experimental information with molecular dynamics

A goal of structural biology is to understand how macromolecules carry out their biological roles by identifying their metastable states, mechanisms of action, pathways leading to conformational changes, and the thermodynamic and kinetic relationships between those states. Integrative modeling brings structural insights into systems where traditional structure determination approaches cannot help. We focus on the synergies and challenges of integrative modeling combining experimental data with molecular dynamics simulations.

Epistemic expression in the determination of biomolecular structure

Scientific research is constrained by limited resources, so it is imperative that it be conducted efficiently. This paper introduces the notion of epistemic expression, a kind of representation that expedites the solution of research problems. Epistemic expressions are representations that (i) contain information in a way that enables more reliable information to place the most stringent constraints on possible solutions and (ii) make new information readily extractible by biasing the search through that space. I illustrate these conditions using historical and contemporary examples of biomolecular structure determination. Then, I argue that the notion of epistemic expression parts ways with pragmatic accounts of scientific representation and an understanding of models as artifacts, neither of which require models to accurately represent. Explicating epistemic expression thus fills a gap in our understanding of scientific practice, extending Morrison and Morgan's (1999) conception of models as investigative instruments.

Predicting structures of large protein assemblies using combinatorial assembly algorithm and AlphaFold2

Deep learning models, such as AlphaFold2 and RosettaFold, enable high-accuracy protein structure prediction. However, large protein complexes are still challenging to predict due to their size and the complexity of interactions between multiple subunits. Here we present CombFold, a combinatorial and hierarchical assembly algorithm for predicting structures of large protein complexes utilizing pairwise interactions between subunits predicted by AlphaFold2. CombFold accurately predicted (TM-score > 0.7) 72% of the complexes among the Top-10 predictions in two datasets of 60 large, asymmetric assemblies. Moreover, the structural coverage of predicted complexes was 20% higher compared to corresponding PDB entries. We applied the method on complexes from Complex Portal with known stoichiometry but without known structure and obtained high-confidence predictions. CombFold supports the integration of distance restraints based on crosslinking mass spectrometry and fast enumeration of possible complex stoichiometries. CombFold's high accuracy makes it a promising tool for expanding structural coverage beyond monomeric proteins.

Multi-state modeling of antibody-antigen complexes with SAXS profiles and deep-learning models

Antibodies are an established class of human therapeutics. Epitope characterization is an important part of therapeutic antibody discovery. However, structural characterization of antibody-antigen complexes remains challenging. On the one hand, X-ray crystallography or cryo-electron microscopy provide atomic resolution characterization of the epitope, but the data collection process is typically long and the success rate is low. On the other hand, computational methods for modeling antibody-antigen structures from the individual components frequently suffer from a high false positive rate, rarely resulting in a unique solution. Recent deep learning models for structure prediction are also successful in predicting protein-protein complexes. However, they do not perform well for antibody-antigen complexes. Small Angle X-ray Scattering (SAXS) is a reliable technique for rapid structural characterization of protein samples in solution albeit at low resolution. Here, we present an integrative approach for modeling antigen-antibody complexes using the antibody sequence, antigen structure, and experimentally determined SAXS profiles of the antibody, antigen, and the complex. The method models antibody structures using a novel deep-learning approach, NanoNet. The structures of the antibodies and antigens are represented using multiple 3D conformations to account for compositional and conformational heterogeneity of the protein samples that are used to collect the SAXS data. The complexes are predicted by integrating the SAXS profiles with scoring functions for protein-protein interfaces that are based on statistical potentials and antibody-specific deep-learning models. We validated the method via application to four Fab:EGFR and one Fab:PCSK9 antibody:antigen complexes with experimentally available SAXS datasets. The integrative approach returns accurate predictions (interface RMSD<4Å) in the top five predictions for four out of five complexes (respective interface RMSD values of 1.95, 2.18, 2.66 and 3.87Å), providing support for the utility of such a computational pipeline for epitope characterization during therapeutic antibody discovery.

Molecular architecture of nucleosome remodeling and deacetylase sub-complexes by integrative structure determination

The nucleosome remodeling and deacetylase (NuRD) complex is a chromatin-modifying assembly that regulates gene expression and DNA damage repair. Despite its importance, limited structural information describing the complete NuRD complex is available and a detailed understanding of its mechanism is therefore lacking. Drawing on information from SEC-MALLS, DIA-MS, XLMS, negative-stain EM, X-ray crystallography, NMR spectroscopy, secondary structure predictions, and homology models, we applied Bayesian integrative structure determination to investigate the molecular architecture of three NuRD sub-complexes: MTA1-HDAC1-RBBP4, MTA1N -HDAC1-MBD3GATAD2CC , and MTA1-HDAC1-RBBP4-MBD3-GATAD2A [nucleosome deacetylase (NuDe)]. The integrative structures were corroborated by examining independent crosslinks, cryo-EM maps, biochemical assays, known cancer-associated mutations, and structure predictions from AlphaFold. The robustness of the models was assessed by jack-knifing. Localization of the full-length MBD3, which connects the deacetylase and chromatin remodeling modules in NuRD, has not previously been possible; our models indicate two different locations for MBD3, suggesting a mechanism by which MBD3 in the presence of GATAD2A asymmetrically bridges the two modules in NuRD. Further, our models uncovered three previously unrecognized subunit interfaces in NuDe: HDAC1C -MTA1BAH , MTA1BAH -MBD3MBD , and HDAC160-100 -MBD3MBD . Our approach also allowed us to localize regions of unknown structure, such as HDAC1C and MBD3IDR , thereby resulting in the most complete and robustly cross-validated structural characterization of these NuRD sub-complexes so far.

Integrative modeling of the cell

A whole-cell model represents certain aspects of the cell structure and/or function. Due to the high complexity of the cell, an integrative modeling approach is often taken to utilize all available information including experimental data, prior knowledge and prior models. In this review, we summarize an emerging workflow of whole-cell modeling into five steps: (i) gather information; (ii) represent the modeled system into modules; (iii) translate input information into scoring function; (iv) sample the whole-cell model; (v) validate and interpret the model. In particular, we propose the integrative modeling of the cell by combining available (whole-cell) models to maximize the accuracy, precision, and completeness. In addition, we list quantitative predictions of various aspects of cell biology from existing whole-cell models. Moreover, we discuss the remaining challenges and future directions, and highlight the opportunity to establish an integrative spatiotemporal multi-scale whole-cell model based on a community approach.

PrISM: Precision for Integrative Structural Models

MOTIVATION

A single precision value is currently reported for an integrative model. However, precision may vary for different regions of an integrative model owing to varying amounts of input information.

RESULTS

We develop PrISM (Precision for Integrative Structural Models), to efficiently identify high and low-precision regions for integrative models.

AVAILABILITY

PrISM is written in Python and available under the GNU General Public License v3.0 at https://github.com/isblab/prism; benchmark data used in this paper is available at doi:10.5281/zenodo.6241200.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Leveraging crosslinking mass spectrometry in structural and cell biology

Crosslinking mass spectrometry (crosslinking-MS) is a versatile tool providing structural insights into protein conformation and protein-protein interactions. Its medium-resolution residue-residue distance restraints have been used to validate protein structures proposed by other methods and have helped derive models of protein complexes by integrative structural biology approaches. The use of crosslinking-MS in integrative approaches is underpinned by progress in estimating error rates in crosslinking-MS data and in combining these data with other information. The flexible and high-throughput nature of crosslinking-MS has allowed it to complement the ongoing resolution revolution in electron microscopy by providing system-wide residue-residue distance restraints, especially for flexible regions or systems. Here, we review how crosslinking-MS information has been leveraged in structural model validation and integrative modeling. Crosslinking-MS has also been a key technology for cell biology studies and structural systems biology where, in conjunction with cryoelectron tomography, it can provide structural and mechanistic insights directly in situ.

Characterization of an A3G-VifHIV-1-CRL5-CBFβ Structure Using a Cross-linking Mass Spectrometry Pipeline for Integrative Modeling of Host-Pathogen Complexes

Structural analysis of host-pathogen protein complexes remains challenging, largely due to their structural heterogeneity. Here, we describe a pipeline for the structural characterization of these complexes using integrative structure modeling based on chemical cross-links and residue-protein contacts inferred from mutagenesis studies. We used this approach on the HIV-1 Vif protein bound to restriction factor APOBEC3G (A3G), the Cullin-5 E3 ring ligase (CRL5), and the cellular transcription factor Core Binding Factor Beta (CBFβ) to determine the structure of the (A3G-Vif-CRL5-CBFβ) complex. Using the MS-cleavable DSSO cross-linker to obtain a set of 132 cross-links within this reconstituted complex along with the atomic structures of the subunits and mutagenesis data, we computed an integrative structure model of the heptameric A3G-Vif-CRL5-CBFβ complex. The structure, which was validated using a series of tests, reveals that A3G is bound to Vif mostly through its N-terminal domain. Moreover, the model ensemble quantifies the dynamic heterogeneity of the A3G C-terminal domain and Cul5 positions. Finally, the model was used to rationalize previous structural, mutagenesis and functional data not used for modeling, including information related to the A3G-bound and unbound structures as well as mapping functional mutations to the A3G-Vif interface. The experimental and computational approach described here is generally applicable to other challenging host-pathogen protein complexes.

High resolution ensemble description of metamorphic and intrinsically disordered proteins using an efficient hybrid parallel tempering scheme

Mapping free energy landscapes of complex multi-funneled metamorphic proteins and weakly-funneled intrinsically disordered proteins (IDPs) remains challenging. While rare-event sampling molecular dynamics simulations can be useful, they often need to either impose restraints or reweigh the generated data to match experiments. Here, we present a parallel-tempering method that takes advantage of accelerated water dynamics and allows efficient and accurate conformational sampling across a wide variety of proteins. We demonstrate the improved sampling efficiency by benchmarking against standard model systems such as alanine di-peptide, TRP-cage and β-hairpin. The method successfully scales to large metamorphic proteins such as RFA-H and to highly disordered IDPs such as Histatin-5. Across the diverse proteins, the calculated ensemble averages match well with the NMR, SAXS and other biophysical experiments without the need to reweigh. By allowing accurate sampling across different landscapes, the method opens doors for sampling free energy landscape of complex uncharted proteins.

The active DNA-PK holoenzyme occupies a tensed state in a staggered synaptic complex

In the non-homologous end-joining (NHEJ) of a DNA double-strand break, DNA ends are bound and protected by DNA-PK, which synapses across the break to tether the broken ends and initiate repair. There is little clarity surrounding the nature of the synaptic complex and the mechanism governing the transition to repair. We report an integrative structure of the synaptic complex at a precision of 13.5 Å, revealing a symmetric head-to-head arrangement with a large offset in the DNA ends and an extensive end-protection mechanism involving a previously uncharacterized plug domain. Hydrogen/deuterium exchange mass spectrometry identifies an allosteric pathway connecting DNA end-binding with the kinase domain that places DNA-PK under tension in the kinase-active state. We present a model for the transition from end-protection to repair, where the synaptic complex supports hierarchical processing of the ends and scaffold assembly, requiring displacement of the catalytic subunit and tension release through kinase activity.

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.

Genetic interaction mapping informs integrative structure determination of protein complexes

Determining structures of protein complexes is crucial for understanding cellular functions. Here, we describe an integrative structure determination approach that relies on in vivo measurements of genetic interactions. We construct phenotypic profiles for point mutations crossed against gene deletions or exposed to environmental perturbations, followed by converting similarities between two profiles into an upper bound on the distance between the mutated residues. We determine the structure of the yeast histone H3-H4 complex based on ~500,000 genetic interactions of 350 mutants. We then apply the method to subunits Rpb1-Rpb2 of yeast RNA polymerase II and subunits RpoB-RpoC of bacterial RNA polymerase. The accuracy is comparable to that based on chemical cross-links; using restraints from both genetic interactions and cross-links further improves model accuracy and precision. The approach provides an efficient means to augment integrative structure determination with in vivo observations.

Using Integrative Modeling Platform to compute, validate, and archive a model of a protein complex structure

Biology is advanced by producing structural models of biological systems, such as protein complexes. Some systems are recalcitrant to traditional structure determination methods. In such cases, it may still be possible to produce useful models by integrative structure determination that depends on simultaneous use of multiple types of data. An ensemble of models that are sufficiently consistent with the data is produced by a structural sampling method guided by a data-dependent scoring function. The variation in the ensemble of models quantified the uncertainty of the structure, generally resulting from the uncertainty in the input information and actual structural heterogeneity in the samples used to produce the data. Here, we describe how to generate, assess, and interpret ensembles of integrative structural models using our open source Integrative Modeling Platform program (https://integrativemodeling.org).

Structural basis of CD4 downregulation by HIV-1 Nef

The HIV-1 Nef protein suppresses multiple immune surveillance mechanisms to promote viral pathogenesis and is an attractive target for the development of novel therapeutics. A key function of Nef is to remove the CD4 receptor from the cell surface by hijacking clathrin- and adaptor protein complex 2 (AP2)-dependent endocytosis. However, exactly how Nef does this has been elusive. Here, we describe the underlying mechanism as revealed by a 3.0-Å crystal structure of a fusion protein comprising Nef and the cytoplasmic domain of CD4 bound to the tetrameric AP2 complex. An intricate combination of conformational changes occurs in both Nef and AP2 to enable CD4 binding and downregulation. A pocket on Nef previously identified as crucial for recruiting class I MHC is also responsible for recruiting CD4, revealing a potential approach to inhibit two of Nef's activities and sensitize the virus to immune clearance.

Integrative structure and function of the yeast exocyst complex

Exocyst is an evolutionarily conserved hetero-octameric tethering complex that plays a variety of roles in membrane trafficking, including exocytosis, endocytosis, autophagy, cell polarization, cytokinesis, pathogen invasion, and metastasis. Exocyst serves as a platform for interactions between the Rab, Rho, and Ral small GTPases, SNARE proteins, and Sec1/Munc18 regulators that coordinate spatial and temporal fidelity of membrane fusion. However, its mechanism is poorly described at the molecular level. Here, we determine the molecular architecture of the yeast exocyst complex by an integrative approach, based on a 3D density map from negative-stain electron microscopy (EM) at ~16 Å resolution, 434 disuccinimidyl suberate and 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide hydrochloride cross-links from chemical-crosslinking mass spectrometry, and partial atomic models of the eight subunits. The integrative structure is validated by a previously determined cryo-EM structure, cross-links, and distances from in vivo fluorescence microscopy. Our subunit configuration is consistent with the cryo-EM structure, except for Sec5. While not observed in the cryo-EM map, the integrative model localizes the N-terminal half of Sec3 near the Sec6 subunit. Limited proteolysis experiments suggest that the conformation of Exo70 is dynamic, which may have functional implications for SNARE and membrane interactions. This study illustrates how integrative modeling based on varied low-resolution structural data can inform biologically relevant hypotheses, even in the absence of high-resolution data.

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.

A Practical Guide to the Simultaneous Determination of Protein Structure and Dynamics Using Metainference

Accurate protein structural ensembles can be determined with metainference, a Bayesian inference method that integrates experimental information with prior knowledge of the system and deals with all sources of uncertainty and errors as well as with system heterogeneity. Furthermore, metainference can be implemented using the metadynamics approach, which enables the computational study of complex biological systems requiring extensive conformational sampling. In this chapter, we provide a step-by-step guide to perform and analyse metadynamic metainference simulations using the ISDB module of the open-source PLUMED library, as well as a series of practical tips to avoid common mistakes. Specifically, we will guide the reader in the process of learning how to model the structural ensemble of a small disordered peptide by combining state-of-the-art molecular mechanics force fields with nuclear magnetic resonance data, including chemical shifts, scalar couplings and residual dipolar couplings.

Structural proteomics, electron cryo-microscopy and structural modeling approaches in bacteria-human protein interactions

A central challenge in infection medicine is to determine the structure and function of host-pathogen protein-protein interactions to understand how these interactions facilitate bacterial adhesion, dissemination and survival. In this review, we focus on proteomics, electron cryo-microscopy and structural modeling to showcase instances where affinity-purification (AP) and cross-linking (XL) mass spectrometry (MS) has advanced our understanding of host-pathogen interactions. We highlight cases where XL-MS in combination with structural modeling has provided insight into the quaternary structure of interspecies protein complexes. We further exemplify how electron cryo-tomography has been used to visualize bacterial-human interactions during attachment and infection. Lastly, we discuss how AP-MS, XL-MS and electron cryo-microscopy and -tomography together with structural modeling approaches can be used in future studies to broaden our knowledge regarding the function, dynamics and evolution of such interactions. This knowledge will be of relevance for future drug and vaccine development programs.

Structural dynamics of the human COP9 signalosome revealed by cross-linking mass spectrometry and integrative modeling

The COP9 signalosome (CSN) is an evolutionarily conserved eight-subunit (CSN1-8) protein complex that controls protein ubiquitination by deneddylating Cullin-RING E3 ligases (CRLs). The activation and function of CSN hinges on its structural dynamics, which has been challenging to decipher by conventional tools. Here, we have developed a multichemistry cross-linking mass spectrometry approach enabled by three mass spectometry-cleavable cross-linkers to generate highly reliable cross-link data. We applied this approach with integrative structure modeling to determine the interaction and structural dynamics of CSN with the recently discovered ninth subunit, CSN9, in solution. Our results determined the localization of CSN9 binding sites and revealed CSN9-dependent structural changes of CSN. Together with biochemical analysis, we propose a structural model in which CSN9 binding triggers CSN to adopt a configuration that facilitates CSN-CRL interactions, thereby augmenting CSN deneddylase activity. Our integrative structure analysis workflow can be generalized to define in-solution architectures of dynamic protein complexes that remain inaccessible to other approaches.

Integrating cryo-EM and NMR data

Single-particle cryo-electron microscopy (cryo-EM) is increasingly used as a technique to determine the atomic structure of challenging biological systems. Recent advances in microscope engineering, electron detection, and image processing have allowed the structural determination of bigger and more flexible targets than possible with the complementary techniques X-ray crystallography and NMR spectroscopy. However, there exist many biological targets for which atomic resolution cannot be currently achieved with cryo-EM, making unambiguous determination of the protein structure impossible. Although determining the structure of large biological systems using solely NMR is often difficult, highly complementary experimental atomic-level data for each molecule can be derived from the spectra, and used in combination with cryo-EM data. We review here strategies with which both techniques can be synergistically combined, in order to reach detail and understanding unattainable by each technique acting alone; and the types of biological systems for which such an approach would be desirable.

Bayesian-Maximum-Entropy Reweighting of IDP Ensembles Based on NMR Chemical Shifts

Bayesian and Maximum Entropy approaches allow for a statistically sound and systematic fitting of experimental and computational data. Unfortunately, assessing the relative confidence in these two types of data remains difficult as several steps add unknown error. Here we propose the use of a validation-set method to determine the balance, and thus the amount of fitting. We apply the method to synthetic NMR chemical shift data of an intrinsically disordered protein. We show that the method gives consistent results even when other methods to assess the amount of fitting cannot be applied. Finally, we also describe how the errors in the chemical shift predictor can lead to an incorrect fitting and how using secondary chemical shifts could alleviate this problem.

Structural mass spectrometry goes viral

Over the last 20 years, mass spectrometry (MS), with its ability to analyze small sample amounts with high speed and sensitivity, has more and more entered the field of structural virology, aiming to investigate the structure and dynamics of viral proteins as close to their native environment as possible. The use of non-perturbing labels in hydrogen-deuterium exchange MS allows for the analysis of interactions between viral proteins and host cell factors as well as their dynamic responses to the environment. Cross-linking MS, on the other hand, can analyze interactions in viral protein complexes and identify virus-host interactions in cells. Native MS allows transferring viral proteins, complexes and capsids into the gas phase and has broken boundaries to overcome size limitations, so that now even the analysis of intact virions is possible. Different MS approaches not only inform about size, stability, interactions and dynamics of virus assemblies, but also bridge the gap to other biophysical techniques, providing valuable constraints for integrative structural modeling of viral complex assemblies that are often inaccessible by single technique approaches. In this review, recent advances are highlighted, clearly showing that structural MS approaches in virology are moving towards systems biology and ever more experiments are performed on cellular level.

Role of integrative structural biology in understanding transcriptional initiation

Integrative structural biology combines data from multiple experimental techniques to generate complete structural models for the biological system of interest. Most commonly cross-linking data sets are employed alongside electron microscopy maps, crystallographic structures, and other data by computational methods that integrate all known information and produce structural models at a level of resolution that is appropriate to the input data. The precision of these modelled solutions is limited by the sparseness of cross-links observed, the length of the cross-linking reagent, the ambiguity arisen from the presence of multiple copies of the same protein, and structural and compositional heterogeneity. In recent years integrative structural biology approaches have been successfully applied to a range of RNA polymerase II complexes. Here we will provide a general background to integrative structural biology, a description of how it should be practically implemented and how it has furthered our understanding of the biology of large transcriptional assemblies. Finally, in the context of recent breakthroughs in microscope and direct electron detector technology, where increasingly EM is capable of resolving structural features directly without the aid of other structural techniques, we will discuss the future role of integrative structural techniques.

Simultaneous Determination of Protein Structure and Dynamics Using Cryo-Electron Microscopy

Cryo-electron microscopy is rapidly emerging as a powerful technique to determine the structures of complex macromolecular systems elusive to other techniques. Because many of these systems are highly dynamical, characterizing their movements is also a crucial step to unravel their biological functions. To achieve this goal, we report an integrative modeling approach to simultaneously determine structure and dynamics of macromolecular systems from cryo-electron microscopy density maps. By quantifying the level of noise in the data and dealing with their ensemble-averaged nature, this approach enables the integration of multiple sources of information to model ensembles of structures and infer their populations. We illustrate the method by characterizing structure and dynamics of the integral membrane receptor STRA6, thus providing insights into the mechanisms by which it interacts with retinol binding protein and translocates retinol across the membrane.

Weighing-up protein dynamics: the combination of native mass spectrometry and molecular dynamics simulations

Structural dynamics underpin biological function at the molecular level, yet many biophysical and structural biology approaches give only a static or averaged view of proteins. Native mass spectrometry yields spectra of the many states and interactions in the structural ensemble, but its spatial resolution is limited. Conversely, molecular dynamics simulations are innately high-resolution, but have a limited capacity for exploring all structural possibilities. The two techniques hence differ fundamentally in the information they provide, returning data that reflect different length scales and time scales, making them natural bedfellows. Here we discuss how the combination of native mass spectrometry with molecular dynamics simulations is enabling unprecedented insights into a range of biological questions by interrogating the motions of proteins, their assemblies, and interactions.

Optimizing model representation for integrative structure determination of macromolecular assemblies

Integrative structure determination of macromolecular assemblies requires specifying the representation of the modeled structure, a scoring function for ranking alternative models based on diverse types of data, and a sampling method for generating these models. Structures are often represented at atomic resolution, although ad hoc simplified representations based on generic guidelines and/or trial and error are also used. In contrast, we introduce here the concept of optimizing representation. To illustrate this concept, the optimal representation is selected from a set of candidate representations based on an objective criterion that depends on varying amounts of information available for different parts of the structure. Specifically, an optimal representation is defined as the highest-resolution representation for which sampling is exhaustive at a precision commensurate with the precision of the representation. Thus, the method does not require an input structure and is applicable to any input information. We consider a space of representations in which a representation is a set of nonoverlapping, variable-length segments (i.e., coarse-grained beads) for each component protein sequence. We also implement a method for efficiently finding an optimal representation in our open-source Integrative Modeling Platform (IMP) software (https://integrativemodeling.org/). The approach is illustrated by application to three complexes of two subunits and a large assembly of 10 subunits. The optimized representation facilitates exhaustive sampling and thus can produce a more accurate model and a more accurate estimate of its uncertainty for larger structures than were possible previously.

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.

Treating Biomacromolecular Conformational Variability

The function of a biomacromolecule is related not only to its structure but also to the different conformations that its structural elements can sample. It is therefore important to determine the extent of the structural fluctuations and to identify the states that are actually populated as a result of the rearrangement. However, this accomplishment is undermined by an intrinsic limitation: the amount of experimental data is by and large inferior to the number of the states that a biomacromolecule can actually sample. This means that additional, a priori information must be applied in order to derive the most from the available experimental data but not to run into overinterpretation. In this chapter we will give a summary of the experimental observables that can be used towards the reconstruction of structural ensembles, how the data can be profitably combined and to what extent the data are affected by error; finally we will give an overview of the computational methods that have been developed to model structural ensembles, highlighting their difference and similarities, advantages and disadvantages.

CryoEM-based hybrid modeling approaches for structure determination

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

Conformational ensembles of RNA oligonucleotides from integrating NMR and molecular simulations

RNA molecules are key players in numerous cellular processes and are characterized by a complex relationship between structure, dynamics, and function. Despite their apparent simplicity, RNA oligonucleotides are very flexible molecules, and understanding their internal dynamics is particularly challenging using experimental data alone. We show how to reconstruct the conformational ensemble of four RNA tetranucleotides by combining atomistic molecular dynamics simulations with nuclear magnetic resonance spectroscopy data. The goal is achieved by reweighting simulations using a maximum entropy/Bayesian approach. In this way, we overcome problems of current simulation methods, as well as in interpreting ensemble- and time-averaged experimental data. We determine the populations of different conformational states by considering several nuclear magnetic resonance parameters and point toward properties that are not captured by state-of-the-art molecular force fields. Although our approach is applied on a set of model systems, it is fully general and may be used to study the conformational dynamics of flexible biomolecules and to detect inaccuracies in molecular dynamics force fields.

Opportunities and Challenges in Building a Spatiotemporal Multi-scale Model of the Human Pancreatic β Cell

The construction of a predictive model of an entire eukaryotic cell that describes its dynamic structure from atomic to cellular scales is a grand challenge at the intersection of biology, chemistry, physics, and computer science. Having such a model will open new dimensions in biological research and accelerate healthcare advancements. Developing the necessary experimental and modeling methods presents abundant opportunities for a community effort to realize this goal. Here, we present a vision for creation of a spatiotemporal multi-scale model of the pancreatic β-cell, a relevant target for understanding and modulating the pathogenesis of diabetes.

Integrative structure and functional anatomy of a nuclear pore complex

Despite the central role of Nuclear Pore Complexes (NPCs) as gatekeepers of RNA and protein transport between the cytoplasm and nucleoplasm, their large size and dynamic nature have impeded a full structural and functional elucidation. Here, we have determined a subnanometer precision structure for the entire 552-protein yeast NPC by satisfying diverse data including stoichiometry, a cryo-electron tomography map, and chemical cross-links. The structure reveals the NPC’s functional elements in unprecedented detail. The NPC is built of sturdy diagonal columns to which are attached connector cables, imbuing both strength and flexibility, while tying together all other elements of the NPC, including membrane-interacting regions and RNA processing platforms. Inwardly-directed anchors create a high density of transport factor-docking Phe-Gly repeats in the central channel, organized in distinct functional units. Taken together, this integrative structure allows us to rationalize the architecture, transport mechanism, and evolutionary origins of the NPC.

Data-assisted protein structure modeling by global optimization in CASP12

In CASP12, 2 types of data-assisted protein structure modeling were experimented. Either SAXS experimental data or cross-linking experimental data was provided for a selected number of CASP12 targets that the CASP12 predictor could utilize for better protein structure modeling. We devised 2 separate energy terms for SAXS data and cross-linking data to drive the model structures into more native-like structures that satisfied the given experimental data as much as possible. In CASP11, we successfully performed protein structure modeling using simulated sparse and ambiguously assigned NOE data and/or correct residue-residue contact information, where the only energy term that folded the protein into its native structure was the term which was originated from the given experimental data. However, the 2 types of experimental data provided in CASP12 were far from being sufficient enough to fold the target protein into its native structure because SAXS data provides only the overall shape of the molecule and the cross-linking contact information provides only very low-resolution distance information. For this reason, we combined the SAXS or cross-linking energy term with our regular modeling energy function that includes both the template energy term and the de novo energy terms. By optimizing the newly formulated energy function, we obtained protein models that fit better with provided SAXS data than the X-ray structure of the target. However, the improvement of the model relative to the 1 modeled without the SAXS data, was not significant. Consistent structural improvement was achieved by incorporating cross-linking data into the protein structure modeling.

Structural Architecture of the Nucleosome Remodeler ISWI Determined from Cross-Linking, Mass Spectrometry, SAXS, and Modeling

Chromatin remodeling factors assume critical roles by regulating access to nucleosomal DNA. To determine the architecture of the Drosophila ISWI remodeling enzyme, we developed an integrative structural approach that combines protein cross-linking, mass spectrometry, small-angle X-ray scattering, and computational modeling. The resulting structural model shows the ATPase module in a resting state with both ATPase lobes twisted against each other, providing support for a conformation that was recently trapped by crystallography. The autoinhibiting NegC region does not protrude from the ATPase module as suggested previously. The regulatory NTR domain is located near both ATPase lobes. The full-length enzyme is flexible and can adopt a compact structure in solution with the C-terminal HSS domain packing against the ATPase module. Our data imply a series of conformational changes upon activation of the enzyme and illustrate how the NTR, NegC, and HSS domains contribute to regulation of the ATPase module.

Assessing Exhaustiveness of Stochastic Sampling for Integrative Modeling of Macromolecular Structures

Modeling of macromolecular structures involves structural sampling guided by a scoring function, resulting in an ensemble of good-scoring models. By necessity, the sampling is often stochastic, and must be exhaustive at a precision sufficient for accurate modeling and assessment of model uncertainty. Therefore, the very first step in analyzing the ensemble is an estimation of the highest precision at which the sampling is exhaustive. Here, we present an objective and automated method for this task. As a proxy for sampling exhaustiveness, we evaluate whether two independently and stochastically generated sets of models are sufficiently similar. The protocol includes testing 1) convergence of the model score, 2) whether model scores for the two samples were drawn from the same parent distribution, 3) whether each structural cluster includes models from each sample proportionally to its size, and 4) whether there is sufficient structural similarity between the two model samples in each cluster. The evaluation also provides the sampling precision, defined as the smallest clustering threshold that satisfies the third, most stringent test. We validate the protocol with the aid of enumerated good-scoring models for five illustrative cases of binary protein complexes. Passing the proposed four tests is necessary, but not sufficient for thorough sampling. The protocol is general in nature and can be applied to the stochastic sampling of any set of models, not just structural models. In addition, the tests can be used to stop stochastic sampling as soon as exhaustiveness at desired precision is reached, thereby improving sampling efficiency; they may also help in selecting a model representation that is sufficiently detailed to be informative, yet also sufficiently coarse for sampling to be exhaustive.

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.