Assessing Exhaustiveness of Stochastic Sampling for Integrative Modeling of Macromolecular Structures

Extensive structural rearrangement of intraflagellar transport trains underpins bidirectional cargo transport

Bidirectional transport in cilia is carried out by polymers of the IFTA and IFTB protein complexes, called anterograde and retrograde intraflagellar transport (IFT) trains. Anterograde trains deliver cargoes from the cell to the cilium tip, then convert into retrograde trains for cargo export. We set out to understand how the IFT complexes can perform these two directly opposing roles before and after conversion. We use cryoelectron tomography and in situ cross-linking mass spectrometry to determine the structure of retrograde IFT trains and compare it with the known structure of anterograde trains. The retrograde train is a 2-fold symmetric polymer organized around a central thread of IFTA complexes. We conclude that anterograde-to-retrograde remodeling involves global rearrangements of the IFTA/B complexes and requires complete disassembly of the anterograde train. Finally, we describe how conformational changes to cargo-binding sites facilitate unidirectional cargo transport in a bidirectional system.

The Molecular Architecture of the Nuclear Basket

The nuclear pore complex (NPC) is the sole mediator of nucleocytoplasmic transport. Despite great advances in understanding its conserved core architecture, the peripheral regions can exhibit considerable variation within and between species. One such structure is the cage-like nuclear basket. Despite its crucial roles in mRNA surveillance and chromatin organization, an architectural understanding has remained elusive. Using in-cell cryo-electron tomography and subtomogram analysis, we explored the NPC's structural variations and the nuclear basket across fungi (yeast; S. cerevisiae), mammals (mouse; M. musculus), and protozoa (T. gondii). Using integrative structural modeling, we computed a model of the basket in yeast and mammals that revealed how a hub of Nups in the nuclear ring binds to basket-forming Mlp/Tpr proteins: the coiled-coil domains of Mlp/Tpr form the struts of the basket, while their unstructured termini constitute the basket distal densities, which potentially serve as a docking site for mRNA preprocessing before nucleocytoplasmic transport.

The molecular architecture of the desmosomal outer dense plaque by integrative structural modeling

Desmosomes mediate cell-cell adhesion and are prevalent in tissues under mechanical stress. However, their detailed structural characterization is not available. Here, we characterized the molecular architecture of the desmosomal outer dense plaque (ODP) using Bayesian integrative structural modeling via the Integrative Modeling Platform. Starting principally from the structural interpretation of an electron cryo-tomogram, we integrated information from X-ray crystallography, an immuno-electron microscopy study, biochemical assays, in-silico predictions of transmembrane and disordered regions, homology modeling, and stereochemistry information. The integrative structure was validated by information from imaging, tomography, and biochemical studies that were not used in modeling. The ODP resembles a densely packed cylinder with a PKP layer and a PG layer; the desmosomal cadherins and PKP span these two layers. Our integrative approach allowed us to localize disordered regions, such as N-PKP and PG-C. We refined previous protein-protein interactions between desmosomal proteins and provided possible structural hypotheses for defective cell-cell adhesion in several diseases by mapping disease-related mutations on the structure. Finally, we point to features of the structure that could confer resilience to mechanical stress. Our model provides a basis for generating experimentally verifiable hypotheses on the structure and function of desmosomal proteins in normal and disease states.

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.

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.

Integrated modeling of the Nexin-dynein regulatory complex reveals its regulatory mechanism

Cilia are hairlike protrusions that project from the surface of eukaryotic cells and play key roles in cell signaling and motility. Ciliary motility is regulated by the conserved nexin-dynein regulatory complex (N-DRC), which links adjacent doublet microtubules and regulates and coordinates the activity of outer doublet complexes. Despite its critical role in cilia motility, the assembly and molecular basis of the regulatory mechanism are poorly understood. Here, using cryo-electron microscopy in conjunction with biochemical cross-linking and integrative modeling, we localize 12 DRC subunits in the N-DRC structure of Tetrahymena thermophila. We also find that the CCDC96/113 complex is in close contact with the DRC9/10 in the linker region. In addition, we reveal that the N-DRC is associated with a network of coiled-coil proteins that most likely mediates N-DRC regulatory activity.

Bayesian methods in integrative structure modeling

There is a growing interest in characterizing the structure and dynamics of large biomolecular assemblies and their interactions within the cellular environment. A diverse array of experimental techniques allows us to study biomolecular systems on a variety of length and time scales. These techniques range from imaging with light, X-rays or electrons, to spectroscopic methods, cross-linking mass spectrometry and functional genomics approaches, and are complemented by AI-assisted protein structure prediction methods. A challenge is to integrate all of these data into a model of the system and its functional dynamics. This review focuses on Bayesian approaches to integrative structure modeling. We sketch the principles of Bayesian inference, highlight recent applications to integrative modeling and conclude with a discussion of current challenges and future perspectives.

Integrated modeling of the Nexin-dynein regulatory complex reveals its regulatory mechanism

Cilia are hairlike protrusions that project from the surface of eukaryotic cells and play key roles in cell signaling and motility. Ciliary motility is regulated by the conserved nexin-dynein regulatory complex (N-DRC), which links adjacent doublet microtubules and regulates and coordinates the activity of outer doublet complexes. Despite its critical role in cilia motility, the assembly and molecular basis of the regulatory mechanism are poorly understood. Here, utilizing cryo-electron microscopy in conjunction with biochemical cross-linking and integrative modeling, we localized 12 DRC subunits in the N-DRC structure of Tetrahymena thermophila . We also found that the CCDC96/113 complex is in close contact with the N-DRC. In addition, we revealed that the N-DRC is associated with a network of coiled-coil proteins that most likely mediates N-DRC regulatory activity.

Does AlphaFold2 model proteins' intracellular conformations? An experimental test using cross-linking mass spectrometry of endogenous ciliary proteins

A major goal in structural biology is to understand protein assemblies in their biologically relevant states. Here, we investigate whether AlphaFold2 structure predictions match native protein conformations. We chemically cross-linked proteins in situ within intact Tetrahymena thermophila cilia and native ciliary extracts, identifying 1,225 intramolecular cross-links within the 100 best-sampled proteins, providing a benchmark of distance restraints obeyed by proteins in their native assemblies. The corresponding structure predictions were highly concordant, positioning 86.2% of cross-linked residues within Cɑ-to-Cɑ distances of 30 Å, consistent with the cross-linker length. 43% of proteins showed no violations. Most inconsistencies occurred in low-confidence regions or between domains. Overall, AlphaFold2 predictions with lower predicted aligned error corresponded to more correct native structures. However, we observe cases where rigid body domains are oriented incorrectly, as for ciliary protein BBC118, suggesting that combining structure prediction with experimental information will better reveal biologically relevant conformations.

Integrative structural analysis of the type III secretion system needle complex from Shigella flexneri

The type III secretion system (T3SS) is a large, transmembrane protein machinery used by various pathogenic gram-negative bacteria to transport virulence factors into the host cell during infection. Understanding the structure of T3SSs is crucial for future developments of therapeutics that could target this system. However, much of the knowledge about the structure of T3SS is available only for Salmonella, and it is unclear how this large assembly is conserved across species. Here, we combined cryo-electron microscopy, cross-linking mass spectrometry, and integrative modeling to determine the structure of the T3SS needle complex from Shigella flexneri. We show that the Shigella T3SS exhibits unique features distinguishing it from other structurally characterized T3SSs. The secretin pore complex adopts a new fold of its C-terminal S domain and the pilotin MxiM[SctG] locates around the outer surface of the pore. The export apparatus structure exhibits a conserved pseudohelical arrangement but includes the N-terminal domain of the SpaS[SctU] subunit, which was not present in any of the previously published virulence-related T3SS structures. Similar to other T3SSs, however, the apparatus is anchored within the needle complex by a network of flexible linkers that either adjust conformation to connect to equivalent patches on the secretin oligomer or bind distinct surface patches at the same height of the export apparatus. The conserved and unique features delineated by our analysis highlight the necessity to analyze T3SS in a species-specific manner, in order to fully understand the underlying molecular mechanisms of these systems. The structure of the type III secretion system from Shigella flexneri delineates conserved and unique features, which could be used for the development of broad-range therapeutics.

Structure and dynamics of the essential endogenous mycobacterial polyketide synthase Pks13

The mycolic acid layer of the Mycobacterium tuberculosis cell wall is essential for viability and virulence, and the enzymes responsible for its synthesis are targets for antimycobacterial drug development. Polyketide synthase 13 (Pks13) is a module encoding several enzymatic and transport functions that carries out the condensation of two different long-chain fatty acids to produce mycolic acids. We determined structures by cryogenic-electron microscopy of dimeric multi-enzyme Pks13 purified from mycobacteria under normal growth conditions, captured with native substrates. Structures define the ketosynthase (KS), linker and acyl transferase (AT) domains at 1.8 Å resolution and two alternative locations of the N-terminal acyl carrier protein. These structures suggest intermediate states on the pathway for substrate delivery to the KS domain. Other domains, visible at lower resolution, are flexible relative to the KS-AT core. The chemical structures of three bound endogenous long-chain fatty acid substrates were determined by electrospray ionization mass spectrometry.

Structural basis of a transcription pre-initiation complex on a divergent promoter

Most eukaryotic promoter regions are divergently transcribed. As the RNA polymerase II pre-initiation complex (PIC) is intrinsically asymmetric and responsible for transcription in a single direction, it is unknown how divergent transcription arises. Here, the Saccharomyces cerevisiae Mediator complexed with a PIC (Med-PIC) was assembled on a divergent promoter and analyzed by cryoelectron microscopy. The structure reveals two distinct Med-PICs forming a dimer through the Mediator tail module, induced by a homodimeric activator protein localized near the dimerization interface. The tail dimer is associated with ∼80-bp upstream DNA, such that two flanking core promoter regions are positioned and oriented in a suitable form for PIC assembly in opposite directions. Also, cryoelectron tomography visualized the progress of the PIC assembly on the two core promoter regions, providing direct evidence for the role of the Med-PIC dimer in divergent transcription.

New opportunities in integrative structural modeling

Integrative structural modeling enables structure determination of macromolecules and their complexes by integrating data from multiple sources. It has been successfully used to characterize macromolecular structures when a single structural biology technique was insufficient. Recent developments in cellular structural biology, including in-cell cryo-electron tomography and artificial intelligence-based structure prediction, have created new opportunities for integrative structural modeling. Here, we will review these opportunities along with the latest developments in integrative modeling methods and their applications. We also highlight open challenges and directions for further development.

Mapping of the plant SnRK1 kinase signalling network reveals a key regulatory role for the class II T6P synthase-like proteins

The central metabolic regulator SnRK1 controls plant growth and survival upon activation by energy depletion, but detailed molecular insight into its regulation and downstream targets is limited. Here we used phosphoproteomics to infer the sucrose-dependent processes targeted upon starvation by kinases as SnRK1, corroborating the relation of SnRK1 with metabolic enzymes and transcriptional regulators, while also pointing to SnRK1 control of intracellular trafficking. Next, we integrated affinity purification, proximity labelling and crosslinking mass spectrometry to map the protein interaction landscape, composition and structure of the SnRK1 heterotrimer, providing insight in its plant-specific regulation. At the intersection of this multi-dimensional interactome, we discovered a strong association of SnRK1 with class II T6P synthase (TPS)-like proteins. Biochemical and cellular assays show that TPS-like proteins function as negative regulators of SnRK1. Next to stable interactions with the TPS-like proteins, similar intricate connections were found with known regulators, suggesting that plants utilize an extended kinase complex to fine-tune SnRK1 activity for optimal responses to metabolic stress.

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.

The protein organization of a red blood cell

Red blood cells (RBCs) (erythrocytes) are the simplest primary human cells, lacking nuclei and major organelles and instead employing about a thousand proteins to dynamically control cellular function and morphology in response to physiological cues. In this study, we define a canonical RBC proteome and interactome using quantitative mass spectrometry and machine learning. Our data reveal an RBC interactome dominated by protein homeostasis, redox biology, cytoskeletal dynamics, and carbon metabolism. We validate protein complexes through electron microscopy and chemical crosslinking and, with these data, build 3D structural models of the ankyrin/Band 3/Band 4.2 complex that bridges the spectrin cytoskeleton to the RBC membrane. The model suggests spring-like compression of ankyrin may contribute to the characteristic RBC cell shape and flexibility. Taken together, our study provides an in-depth view of the global protein organization of human RBCs and serves as a comprehensive resource for future research.

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.

From systems to structure - using genetic data to model protein structures

Understanding the effects of genetic variation is a fundamental problem in biology that requires methods to analyse both physical and functional consequences of sequence changes at systems-wide and mechanistic scales. To achieve a systems view, protein interaction networks map which proteins physically interact, while genetic interaction networks inform on the phenotypic consequences of perturbing these protein interactions. Until recently, understanding the molecular mechanisms that underlie these interactions often required biophysical methods to determine the structures of the proteins involved. The past decade has seen the emergence of new approaches based on coevolution, deep mutational scanning and genome-scale genetic or chemical-genetic interaction mapping that enable modelling of the structures of individual proteins or protein complexes. Here, we review the emerging use of large-scale genetic datasets and deep learning approaches to model protein structures and their interactions, and discuss the integration of structural data from different sources.

Doublecortin engages the microtubule lattice through a cooperative binding mode involving its C-terminal domain

Doublecortin (DCX) is a microtubule (MT)-associated protein that regulates MT structure and function during neuronal development and mutations in DCX lead to a spectrum of neurological disorders. The structural properties of MT-bound DCX that explain these disorders are incompletely determined. Here, we describe the molecular architecture of the DCX–MT complex through an integrative modeling approach that combines data from X-ray crystallography, cryo-electron microscopy, and a high-fidelity chemical crosslinking method. We demonstrate that DCX interacts with MTs through its N-terminal domain and induces a lattice-dependent self-association involving the C-terminal structured domain and its disordered tail, in a conformation that favors an open, domain-swapped state. The networked state can accommodate multiple different attachment points on the MT lattice, all of which orient the C-terminal tails away from the lattice. As numerous disease mutations cluster in the C-terminus, and regulatory phosphorylations cluster in its tail, our study shows that lattice-driven self-assembly is an important property of DCX.

Integrative structure determination of histones H3 and H4 using genetic interactions

Integrative structure modeling is increasingly used for determining the architectures of biological assemblies, especially those that are structurally heterogeneous. Recently, we reported on how to convert in vivo genetic interaction measurements into spatial restraints for structural modeling: first, phenotypic profiles are generated for each point mutation and thousands of gene deletions or environmental perturbations. Following, the phenotypic profile similarities are converted into distance restraints on the pairs of mutated residues. We illustrate the approach by determining the structure of the histone H3-H4 complex. The method is implemented in our open-source IMP program, expanding the structural biology toolbox by allowing structural characterization based on in vivo data without the need to purify the target system. We compare genetic interaction measurements to other sources of structural information, such as residue coevolution and deep-learning structure prediction of complex subunits. We also suggest that determining genetic interactions could benefit from new technologies, such as CRISPR-Cas9 approaches to gene editing, especially for mammalian cells. Finally, we highlight the opportunity for using genetic interactions to determine recalcitrant biomolecular structures, such as those of disordered proteins, transient protein assemblies, and host-pathogen protein complexes.

Soluble TREM2 inhibits secondary nucleation of Aβ fibrillization and enhances cellular uptake of fibrillar Aβ

Triggering receptor expressed on myeloid cells 2 (TREM2) is a single-pass transmembrane receptor of the immunoglobulin superfamily that is secreted in a soluble (sTREM2) form. Mutations in TREM2 have been linked to increased risk of Alzheimer's disease (AD). A prominent neuropathological component of AD is deposition of the amyloid-β (Aβ) into plaques, particularly Aβ40 and Aβ42. While the membrane-bound form of TREM2 is known to facilitate uptake of Aβ fibrils and the polarization of microglial processes toward amyloid plaques, the role of its soluble ectodomain, particularly in interactions with monomeric or fibrillar Aβ, has been less clear. Our results demonstrate that sTREM2 does not bind to monomeric Aβ40 and Aβ42, even at a high micromolar concentration, while it does bind to fibrillar Aβ42 and Aβ40 with equal affinities (2.6 ± 0.3 µM and 2.3 ± 0.4 µM). Kinetic analysis shows that sTREM2 inhibits the secondary nucleation step in the fibrillization of Aβ, while having little effect on the primary nucleation pathway. Furthermore, binding of sTREM2 to fibrils markedly enhanced uptake of fibrils into human microglial and neuroglioma derived cell lines. The disease-associated sTREM2 mutant, R47H, displayed little to no effect on fibril nucleation and binding, but it decreased uptake and functional responses markedly. We also probed the structure of the WT sTREM2-Aβ fibril complex using integrative molecular modeling based primarily on the cross-linking mass spectrometry data. The model shows that sTREM2 binds fibrils along one face of the structure, leaving a second, mutation-sensitive site free to mediate cellular binding and uptake.

Comprehensive structure and functional adaptations of the yeast nuclear pore complex

Nuclear pore complexes (NPCs) mediate the nucleocytoplasmic transport of macromolecules. Here we provide a structure of the isolated yeast NPC in which the inner ring is resolved by cryo-EM at sub-nanometer resolution to show how flexible connectors tie together different structural and functional layers. These connectors may be targets for phosphorylation and regulated disassembly in cells with an open mitosis. Moreover, some nucleoporin pairs and transport factors have similar interaction motifs, which suggests an evolutionary and mechanistic link between assembly and transport. We provide evidence for three major NPC variants that may foreshadow functional specializations at the nuclear periphery. Cryo-electron tomography extended these studies, providing a model of the in situ NPC with a radially expanded inner ring. Our comprehensive model reveals features of the nuclear basket and central transporter, suggests a role for the lumenal Pom152 ring in restricting dilation, and highlights structural plasticity that may be required for transport.

Nuclear pores dilate and constrict in cellulo

[Figure: see text].

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.

Integrative structural modeling of macromolecular complexes using Assembline

Integrative modeling enables structure determination of macromolecular complexes by combining data from multiple experimental sources such as X-ray crystallography, electron microscopy or cross-linking mass spectrometry. It is particularly useful for complexes not amenable to high-resolution electron microscopy-complexes that are flexible, heterogeneous or imaged in cells with cryo-electron tomography. We have recently developed an integrative modeling protocol that allowed us to model multi-megadalton complexes as large as the nuclear pore complex. Here, we describe the Assembline software package, which combines multiple programs and libraries with our own algorithms in a streamlined modeling pipeline. Assembline builds ensembles of models satisfying data from atomic structures or homology models, electron microscopy maps and other experimental data, and provides tools for their analysis. Compared with other methods, Assembline enables efficient sampling of conformational space through a multistep procedure, provides new modeling restraints and includes a unique configuration system for setting up the modeling project. Our protocol achieves exhaustive sampling in less than 100-1,000 CPU-hours even for complexes in the megadalton range. For larger complexes, resources available in institutional or public computer clusters are needed and sufficient to run the protocol. We also provide step-by-step instructions for preparing the input, running the core modeling steps and assessing modeling performance at any stage.

A Framework for Stochastic Optimization of Parameters for Integrative Modeling of Macromolecular Assemblies

Integrative modeling of macromolecular assemblies requires stochastic sampling, for example, via MCMC (Markov Chain Monte Carlo), since exhaustively enumerating all structural degrees of freedom is infeasible. MCMC-based methods usually require tuning several parameters, such as the move sizes for coarse-grained beads and rigid bodies, for sampling to be efficient and accurate. Currently, these parameters are tuned manually. To automate this process, we developed a general heuristic for derivative-free, global, stochastic, parallel, multiobjective optimization, termed StOP (Stochastic Optimization of Parameters) and applied it to optimize sampling-related parameters for the Integrative Modeling Platform (IMP). Given an integrative modeling setup, list of parameters to optimize, their domains, metrics that they influence, and the target ranges of these metrics, StOP produces the optimal values of these parameters. StOP is adaptable to the available computing capacity and converges quickly, allowing for the simultaneous optimization of a large number of parameters. However, it is not efficient at high dimensions and not guaranteed to find optima in complex landscapes. We demonstrate its performance on several examples of random functions, as well as on two integrative modeling examples, showing that StOP enhances the efficiency of sampling the posterior distribution, resulting in more good-scoring models and better sampling precision.

IMProv: A Resource for Cross-link-Driven Structure Modeling that Accommodates Protein Dynamics

Proteomics methodology has expanded to include protein structural analysis, primarily through cross-linking mass spectrometry (XL-MS) and hydrogen-deuterium exchange mass spectrometry (HX-MS). However, while the structural proteomics community has effective tools for primary data analysis, there is a need for structure modeling pipelines that are accessible to the proteomics specialist. Integrative structural biology requires the aggregation of multiple distinct types of data to generate models that satisfy all inputs. Here, we describe IMProv, an app in the Mass Spec Studio that combines XL-MS data with other structural data, such as cryo-EM densities and crystallographic structures, for integrative structure modeling on high-performance computing platforms. The resource provides an easily deployed bundle that includes the open-source Integrative Modeling Platform program (IMP) and its dependencies. IMProv also provides functionality to adjust cross-link distance restraints according to the underlying dynamics of cross-linked sites, as characterized by HX-MS. A dynamics-driven conditioning of restraint values can improve structure modeling precision, as illustrated by an integrative structure of the five-membered Polycomb Repressive Complex 2. IMProv is extensible to additional types of data.

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.

Structural insight on assembly-line catalysis in terpene biosynthesis

Fusicoccadiene synthase from Phomopsis amygdali (PaFS) is a unique bifunctional terpenoid synthase that catalyzes the first two steps in the biosynthesis of the diterpene glycoside Fusicoccin A, a mediator of 14-3-3 protein interactions. The prenyltransferase domain of PaFS generates geranylgeranyl diphosphate, which the cyclase domain then utilizes to generate fusicoccadiene, the tricyclic hydrocarbon skeleton of Fusicoccin A. Here, we use cryo-electron microscopy to show that the structure of full-length PaFS consists of a central octameric core of prenyltransferase domains, with the eight cyclase domains radiating outward via flexible linker segments in variable splayed-out positions. Cryo-electron microscopy and chemical crosslinking experiments additionally show that compact conformations can be achieved in which cyclase domains are more closely associated with the prenyltransferase core. This structural analysis provides a framework for understanding substrate channeling, since most of the geranylgeranyl diphosphate generated by the prenyltransferase domains remains on the enzyme for cyclization to form fusicoccadiene.

Structure of the mycobacterial ESX-5 type VII secretion system pore complex

The ESX-5 type VII secretion system is a membrane-spanning protein complex key to the virulence of mycobacterial pathogens. However, the overall architecture of the fully assembled translocation machinery and the composition of the central secretion pore have remained unknown. Here, we present the high-resolution structure of the 2.1-megadalton ESX-5 core complex. Our structure captured a dynamic, secretion-competent conformation of the pore within a well-defined transmembrane section, sandwiched between two flexible protein layers at the cytosolic entrance and the periplasmic exit. We propose that this flexibility endows the ESX-5 machinery with large conformational plasticity required to accommodate targeted protein secretion. Compared to known secretion systems, a highly dynamic state of the pore may represent a fundamental principle of bacterial secretion machineries.

CM1-driven assembly and activation of yeast γ-tubulin small complex underlies microtubule nucleation

Microtubule (MT) nucleation is regulated by the γ-tubulin ring complex (γTuRC), conserved from yeast to humans. In Saccharomyces cerevisiae, γTuRC is composed of seven identical γ-tubulin small complex (γTuSC) sub-assemblies, which associate helically to template MT growth. γTuRC assembly provides a key point of regulation for the MT cytoskeleton. Here, we combine crosslinking mass spectrometry, X-ray crystallography, and cryo-EM structures of both monomeric and dimeric γTuSCs, and open and closed helical γTuRC assemblies in complex with Spc110p to elucidate the mechanisms of γTuRC assembly. γTuRC assembly is substantially aided by the evolutionarily conserved CM1 motif in Spc110p spanning a pair of adjacent γTuSCs. By providing the highest resolution and most complete views of any γTuSC assembly, our structures allow phosphorylation sites to be mapped, surprisingly suggesting that they are mostly inhibitory. A comparison of our structures with the CM1 binding site in the human γTuRC structure at the interface between GCP2 and GCP6 allows for the interpretation of significant structural changes arising from CM1 helix binding to metazoan γTuRC.

Improving integrative 3D modeling into low- to medium-resolution electron microscopy structures with evolutionary couplings

Electron microscopy (EM) continues to provide near-atomic resolution structures for well-behaved proteins and protein complexes. Unfortunately, structures of some complexes are limited to low- to medium-resolution due to biochemical or conformational heterogeneity. Thus, the application of unbiased systematic methods for fitting individual structures into EM maps is important. A method that employs co-evolutionary information obtained solely from sequence data could prove invaluable for quick, confident localization of subunits within these structures. Here, we incorporate the co-evolution of intermolecular amino acids as a new type of distance restraint in the integrative modeling platform in order to build three-dimensional models of atomic structures into EM maps ranging from 10-14 Å in resolution. We validate this method using four complexes of known structure, where we highlight the conservation of intermolecular couplings despite dynamic conformational changes using the BAM complex. Finally, we use this method to assemble the subunits of the bacterial holo-translocon into a model that agrees with previous biochemical data. The use of evolutionary couplings in integrative modeling improves systematic, unbiased fitting of atomic models into medium- to low-resolution EM maps, providing additional information to integrative models lacking in spatial data.

Nanobody Repertoires for Exposing Vulnerabilities of SARS-CoV-2

Despite the great promise of vaccines, the COVID-19 pandemic is ongoing and future serious outbreaks are highly likely, so that multi-pronged containment strategies will be required for many years. Nanobodies are the smallest naturally occurring single domain antigen binding proteins identified to date, possessing numerous properties advantageous to their production and use. We present a large repertoire of high affinity nanobodies against SARS-CoV-2 Spike protein with excellent kinetic and viral neutralization properties, which can be strongly enhanced with oligomerization. This repertoire samples the epitope landscape of the Spike ectodomain inside and outside the receptor binding domain, recognizing a multitude of distinct epitopes and revealing multiple neutralization targets of pseudoviruses and authentic SARS-CoV-2, including in primary human airway epithelial cells. Combinatorial nanobody mixtures show highly synergistic activities, and are resistant to mutational escape and emerging viral variants of concern. These nanobodies establish an exceptional resource for superior COVID-19 prophylactics and therapeutics.

Molecular architecture of the endocytic TPLATE complex

Eukaryotic cells rely on endocytosis to regulate their plasma membrane proteome and lipidome. Most eukaryotic groups, except fungi and animals, have retained the evolutionary ancient TSET complex as an endocytic regulator. Unlike other coatomer complexes, structural insight into TSET is lacking. Here, we reveal the molecular architecture of plant TSET [TPLATE complex (TPC)] using an integrative structural approach. We identify crucial roles for specific TSET subunits in complex assembly and membrane interaction. Our data therefore generate fresh insight into the differences between the hexameric TSET in Dictyostelium and the octameric TPC in plants. Structural elucidation of this ancient adaptor complex represents the missing piece in the coatomer puzzle and vastly advances our functional as well as evolutionary insight into the process of endocytosis.

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.

MMM: Integrative ensemble modeling and ensemble analysis

Proteins and their complexes can be heterogeneously disordered. In ensemble modeling of such systems with restraints from several experimental techniques the following problems arise: (a) integration of diverse restraints obtained on different samples under different conditions; (b) estimation of a realistic ensemble width; (c) sufficient sampling of conformational space; (d) representation of the ensemble by an interpretable number of conformers; (e) recognition of weak order with site resolution. Here, I introduce several tools that address these problems, focusing on utilization of distance distribution information for estimating ensemble width. The RigiFlex approach integrates such information with high-resolution structures of ordered domains and small-angle scattering data. The EnsembleFit module provides moderately sized ensembles by fitting conformer populations and discarding conformers with low population. EnsembleFit balances the loss in fit quality upon combining restraint subsets from different techniques. Pair correlation analysis for residues and local compaction analysis help in feature detection. The RigiFlex pipeline is tested on data simulated from the structure 70 kDa protein-RNA complex RsmE/RsmZ. It recovers this structure with ensemble width and difference from ground truth both being on the order of 4.2 Å. EnsembleFit reduces the ensemble of the proliferating-cell-nuclear-antigen-associated factor p15PAF from 4,939 to 75 conformers while maintaining good fit quality of restraints. Local compaction analysis for the PaaA2 antitoxin from E. coli O157 revealed correlations between compactness and enhanced residual dipolar couplings in the original NMR restraint set.

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).

In-cell architecture of an actively transcribing-translating expressome

Structural biology studies performed inside cells can capture molecular machines in action within their native context. In this work, we developed an integrative in-cell structural approach using the genome-reduced human pathogen Mycoplasma pneumoniae We combined whole-cell cross-linking mass spectrometry, cellular cryo-electron tomography, and integrative modeling to determine an in-cell architecture of a transcribing and translating expressome at subnanometer resolution. The expressome comprises RNA polymerase (RNAP), the ribosome, and the transcription elongation factors NusG and NusA. We pinpointed NusA at the interface between a NusG-bound elongating RNAP and the ribosome and propose that it can mediate transcription-translation coupling. Translation inhibition dissociated the expressome, whereas transcription inhibition stalled and rearranged it. Thus, the active expressome architecture requires both translation and transcription elongation within the cell.

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.

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.

High-density chemical cross-linking for modeling protein interactions

Detailed mechanistic understanding of protein complex function is greatly enhanced by insights from its 3-dimensional structure. Traditional methods of protein structure elucidation remain expensive and labor-intensive and require highly purified starting material. Chemical cross-linking coupled with mass spectrometry offers an alternative that has seen increased use, especially in combination with other experimental approaches like cryo-electron microscopy. Here we report advances in method development, combining several orthogonal cross-linking chemistries as well as improvements in search algorithms, statistical analysis, and computational cost to achieve coverage of 1 unique cross-linked position pair for every 7 amino acids at a 1% false discovery rate. This is accomplished without any peptide-level fractionation or enrichment. We apply our methods to model the complex between a carbonic anhydrase (CA) and its protein inhibitor, showing that the cross-links are self-consistent and define the interaction interface at high resolution. The resulting model suggests a scaffold for development of a class of protein-based inhibitors of the CA family of enzymes. We next cross-link the yeast proteasome, identifying 3,893 unique cross-linked peptides in 3 mass spectrometry runs. The dataset includes 1,704 unique cross-linked position pairs for the proteasome subunits, more than half of them intersubunit. Using multiple recently solved cryo-EM structures, we show that observed cross-links reflect the conformational dynamics and disorder of some proteasome subunits. We further demonstrate that this level of cross-linking density is sufficient to model the architecture of the 19-subunit regulatory particle de novo.

Federating Structural Models and Data: Outcomes from A Workshop on Archiving Integrative Structures

Structures of biomolecular systems are increasingly computed by integrative modeling. In this approach, a structural model is constructed by combining information from multiple sources, including varied experimental methods and prior models. In 2019, a Workshop was held as a Biophysical Society Satellite Meeting to assess progress and discuss further requirements for archiving integrative structures. The primary goal of the Workshop was to build consensus for addressing the challenges involved in creating common data standards, building methods for federated data exchange, and developing mechanisms for validating integrative structures. The summary of the Workshop and the recommendations that emerged are presented here.

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.

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.

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.