Putting the pieces together: integrative modeling platform software for structure determination of macromolecular assemblies

IHMCIF: An Extension of the PDBx/mmCIF Data Standard for Integrative Structure Determination Methods

IHMCIF (github.com/ihmwg/IHMCIF) is a data information framework that supports archiving and disseminating macromolecular structures determined by integrative or hybrid modeling (IHM), and making them Findable, Accessible, Interoperable, and Reusable (FAIR). IHMCIF is an extension of the Protein Data Bank Exchange/macromolecular Crystallographic Information Framework (PDBx/mmCIF) that serves as the framework for the Protein Data Bank (PDB) to archive experimentally determined atomic structures of biological macromolecules and their complexes with one another and small molecule ligands (e.g., enzyme cofactors and drugs). IHMCIF serves as the foundational data standard for the PDB-Dev prototype system, developed for archiving and disseminating integrative structures. It utilizes a flexible data representation to describe integrative structures that span multiple spatiotemporal scales and structural states with definitions for restraints from a variety of experimental methods contributing to integrative structural biology. The IHMCIF extension was created with the benefit of considerable community input and recommendations gathered by the Worldwide Protein Data Bank (wwPDB) Task Force for Integrative or Hybrid Methods (wwpdb.org/task/hybrid). Herein, we describe the development of IHMCIF to support evolving methodologies and ongoing advancements in integrative structural biology. Ultimately, IHMCIF will facilitate the unification of PDB-Dev data and tools with the PDB archive so that integrative structures can be archived and disseminated through PDB.

On the Power and Challenges of Atomistic Molecular Dynamics to Investigate RNA Molecules

RNA molecules play a vital role in biological processes within the cell, with significant implications for science and medicine. Notably, the biological functions exerted by specific RNA molecules are often linked to the RNA conformational ensemble. However, the experimental characterization of such three-dimensional RNA structures is challenged by the structural heterogeneity of RNA and by its multiple dynamic interactions with binding partners such as small molecules, proteins, and metal ions. Consequently, our current understanding of the structure-function relationship of RNA molecules is still limited. In this context, we highlight molecular dynamics (MD) simulations as a powerful tool to complement experimental efforts on RNAs. Despite the recognized limitations of current force fields for RNA MD simulations, examining the dynamics of selected RNAs has provided valuable functional insights into their structures.

An integrated structural model of the DNA damage-responsive H3K4me3 binding WDR76:SPIN1 complex with the nucleosome

Serial capture affinity purification (SCAP) is a powerful method to isolate a specific protein complex. When combined with cross-linking mass spectrometry and computational approaches, one can build an integrated structural model of the isolated complex. Here, we applied SCAP to dissect a subpopulation of WDR76 in complex with SPIN1, a histone reader that recognizes trimethylated histone H3 lysine4 (H3K4me3). In contrast to a previous SCAP analysis of the SPIN1:SPINDOC complex, histones and the H3K4me3 mark were enriched with the WDR76:SPIN1 complex. Next, interaction network analysis of copurifying proteins and microscopy analysis revealed a potential role of the WDR76:SPIN1 complex in the DNA damage response. Since we detected 149 pairs of cross-links between WDR76, SPIN1, and histones, we then built an integrated structural model of the complex where SPIN1 recognized the H3K4me3 epigenetic mark while interacting with WDR76. Finally, we used the powerful Bayesian Integrative Modeling approach as implemented in the Integrative Modeling Platform to build a model of WDR76 and SPIN1 bound to the nucleosome.

DomainFit: Identification of protein domains in cryo-EM maps at intermediate resolution using AlphaFold2-predicted models

Cryoelectron microscopy (cryo-EM) has revolutionized the structural determination of macromolecular complexes. With the paradigm shift to structure determination of highly complex endogenous macromolecular complexes ex vivo and in situ structural biology, there are an increasing number of structures of native complexes. These complexes often contain unidentified proteins, related to different cellular states or processes. Identifying proteins at resolutions lower than 4 Å remains challenging because side chains cannot be visualized reliably. Here, we present DomainFit, a program for semi-automated domain-level protein identification from cryo-EM maps, particularly at resolutions lower than 4 Å. By fitting domains from AlphaFold2-predicted models into cryo-EM maps, the program performs statistical analyses and attempts to identify the domains and protein candidates forming the density. Using DomainFit, we identified two microtubule inner proteins, one of which contains a CCDC81 domain and is exclusively localized in the proximal region of the doublet microtubule in Tetrahymena thermophila.

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 nucleoporins (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.

Structures, dynamics, complexes, and functions: From classic computation to artificial intelligence

Computational approaches can provide highly detailed insight into the molecular recognition processes that underlie drug binding, the assembly of protein complexes, and the regulation of biological functional processes. Classical simulation methods can bridge a wide range of length- and time-scales typically involved in such processes. Lately, automated learning and artificial intelligence methods have shown the potential to expand the reach of physics-based approaches, ushering in the possibility to model and even design complex protein architectures. The synergy between atomistic simulations and AI methods is an emerging frontier with a huge potential for advances in structural biology. Herein, we explore various examples and frameworks for these approaches, providing select instances and applications that illustrate their impact on fundamental biomolecular problems.

Integrative modeling meets deep learning: Recent advances in modeling protein assemblies

Recent progress in protein structure prediction based on deep learning revolutionized the field of Structural Biology. Beyond single proteins, it also enabled high-throughput prediction of structures of protein-protein interactions. Despite the success in predicting complex structures, large macromolecular assemblies still require specialized approaches. Here we describe recent advances in modeling macromolecular assemblies using integrative and hierarchical approaches. We highlight applications that predict protein-protein interactions and challenges in modeling complexes based on the interaction networks, including the prediction of complex stoichiometry and heterogeneity.

IDRWalker: A Random Walk Based Tool for Generating Intrinsically Disordered Regions in Large Protein Complexes

Intrinsically disordered regions (IDRs), which may be functionally important, are common in proteins. However, the structures of IDRs are often missing due to their highly dynamic nature. In the study of IDRs, integrative modeling combining computational simulations and experimental data is a common approach, for which initial structures of the IDRs need to be built. However, applying this method to large protein complexes is challenging because existing structure generation tools are sometimes unsuitable for IDRs in large systems. To facilitate convenient and rapid structure generation of IDRs in large protein complexes, we developed a computational tool named IDRWalker based on self-avoiding random walks. Three protein complexes were used to illustrate the efficiency of the tool, and it was found that IDRs in more than 800 chains of the nuclear pore complex could be generated in minutes. These structures of large protein complexes with added IDRs can be further used to run computational simulations for integrative modeling.

Cell Maps for Artificial Intelligence: AI-Ready Maps of Human Cell Architecture from Disease-Relevant Cell Lines

This article describes the Cell Maps for Artificial Intelligence (CM4AI) project and its goals, methods, standards, current datasets, software tools , status, and future directions. CM4AI is the Functional Genomics Data Generation Project in the U.S. National Institute of Health's (NIH) Bridge2AI program. Its overarching mission is to produce ethical, AI-ready datasets of cell architecture, inferred from multimodal data collected for human cell lines, to enable transformative biomedical AI research.

Connexin Gap Junction Channels and Hemichannels: Insights from High-Resolution Structures

Connexins (Cxs) are a family of integral membrane proteins, which function as both hexameric hemichannels (HCs) and dodecameric gap junction channels (GJCs), behaving as conduits for the electrical and molecular communication between cells and between cells and the extracellular environment, respectively. Their proper functioning is crucial for many processes, including development, physiology, and response to disease and trauma. Abnormal GJC and HC communication can lead to numerous pathological states including inflammation, skin diseases, deafness, nervous system disorders, and cardiac arrhythmias. Over the last 15 years, high-resolution X-ray and electron cryomicroscopy (cryoEM) structures for seven Cx isoforms have revealed conservation in the four-helix transmembrane (TM) bundle of each subunit; an αβ fold in the disulfide-bonded extracellular loops and inter-subunit hydrogen bonding across the extracellular gap that mediates end-to-end docking to form a tight seal between hexamers in the GJC. Tissue injury is associated with cellular Ca2+ overload. Surprisingly, the binding of 12 Ca2+ ions in the Cx26 GJC results in a novel electrostatic gating mechanism that blocks cation permeation. In contrast, acidic pH during tissue injury elicits association of the N-terminal (NT) domains that sterically blocks the pore in a "ball-and-chain" fashion. The NT domains under physiologic conditions display multiple conformational states, stabilized by protein-protein and protein-lipid interactions, which may relate to gating mechanisms. The cryoEM maps also revealed putative lipid densities within the pore, intercalated among transmembrane α-helices and between protomers, the functions of which are unknown. For the future, time-resolved cryoEM of isolated Cx channels as well as cryotomography of GJCs and HCs in cells and tissues will yield a deeper insight into the mechanisms for channel regulation. The cytoplasmic loop (CL) and C-terminal (CT) domains are divergent in sequence and length, are likely involved in channel regulation, but are not visualized in the high-resolution X-ray and cryoEM maps presumably due to conformational flexibility. We expect that the integrated use of synergistic physicochemical, spectroscopic, biophysical, and computational methods will reveal conformational dynamics relevant to functional states. We anticipate that such a wealth of results under different pathologic conditions will accelerate drug discovery related to Cx channel modulation.

Physical models of bacterial chromosomes

The interplay between bacterial chromosome organization and functions such as transcription and replication can be studied in increasing detail using novel experimental techniques. Interpreting the resulting quantitative data, however, can be theoretically challenging. In this minireview, we discuss how connecting experimental observations to biophysical theory and modeling can give rise to new insights on bacterial chromosome organization. We consider three flavors of models of increasing complexity: simple polymer models that explore how physical constraints, such as confinement or plectoneme branching, can affect bacterial chromosome organization; bottom-up mechanistic models that connect these constraints to their underlying causes, for instance, chromosome compaction to macromolecular crowding, or supercoiling to transcription; and finally, data-driven methods for inferring interpretable and quantitative models directly from complex experimental data. Using recent examples, we discuss how biophysical models can both deepen our understanding of how bacterial chromosomes are structured and give rise to novel predictions about bacterial chromosome organization.

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.

Leveraging Cross-Linking Mass Spectrometry for Modeling Antibody-Antigen Complexes

Elucidating antibody-antigen complexes at the atomic level is of utmost interest for understanding immune responses and designing better therapies. Cross-linking mass spectrometry (XL-MS) has emerged as a powerful tool for mapping protein-protein interactions, suggesting valuable structural insights. However, the use of XL-MS studies to enable epitope/paratope mapping of antibody-antigen complexes is still limited up to now. XL-MS data can be used to drive integrative modeling of antibody-antigen complexes, where cross-links information serves as distance restraints for the precise determination of binding interfaces. In this approach, XL-MS data are employed to identify connections between binding interfaces of the antibody and the antigen, thus informing molecular modeling. Current literature provides minimal input about the impact of XL-MS data on the integrative modeling of antibody-antigen complexes. Here, we applied XL-MS to retrieve information about binding interfaces of three antibody-antigen complexes. We leveraged XL-MS data to perform integrative modeling using HADDOCK (active-passive residues and distance restraints strategies) and AlphaLink2. We then compared these three approaches with initial predictions of investigated antibody-antigen complexes by AlphaFold Multimer. This work emphasizes the importance of cross-linking data in resolving conformational dynamics of antibody-antigen complexes, ultimately enhancing the design of better protein therapeutics and vaccines.

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.

Understanding the cell: Future views of structural biology

Determining the structure and mechanisms of all individual functional modules of cells at high molecular detail has often been seen as equal to understanding how cells work. Recent technical advances have led to a flush of high-resolution structures of various macromolecular machines, but despite this wealth of detailed information, our understanding of cellular function remains incomplete. Here, we discuss present-day limitations of structural biology and highlight novel technologies that may enable us to analyze molecular functions directly inside cells. We predict that the progression toward structural cell biology will involve a shift toward conceptualizing a 4D virtual reality of cells using digital twins. These will capture cellular segments in a highly enriched molecular detail, include dynamic changes, and facilitate simulations of molecular processes, leading to novel and experimentally testable predictions. Transferring biological questions into algorithms that learn from the existing wealth of data and explore novel solutions may ultimately unveil how cells work.

Bonds and bytes: The odyssey of structural biology

Characterizing structural and dynamic properties of proteins and large macromolecular assemblies is crucial to understand the molecular mechanisms underlying biological functions. In the field of structural biology, no single method comprehensively reveals the behavior of biological systems across various spatiotemporal scales. Instead, we have a versatile toolkit of techniques, each contributing a piece to the overall puzzle. Integrative structural biology combines different techniques to create accurate and precise multi-scale models that expand our understanding of complex biological systems. This review outlines recent advancements in computational and experimental methods in structural biology, with special focus on recent Artificial Intelligence techniques, emphasizes integrative approaches that combine different types of data for precise spatiotemporal modeling, and provides an outlook into future directions of this field.

Chromosome structure modeling tools and their evaluation in bacteria

The three-dimensional (3D) structure of bacterial chromosomes is crucial for understanding chromosome function. With the growing availability of high-throughput chromosome conformation capture (3C/Hi-C) data, the 3D structure reconstruction algorithms have become powerful tools to study bacterial chromosome structure and function. It is highly desired to have a recommendation on the chromosome structure reconstruction tools to facilitate the prokaryotic 3D genomics. In this work, we review existing chromosome 3D structure reconstruction algorithms and classify them based on their underlying computational models into two categories: constraint-based modeling and thermodynamics-based modeling. We briefly compare these algorithms utilizing 3C/Hi-C datasets and fluorescence microscopy data obtained from Escherichia coli and Caulobacter crescentus, as well as simulated datasets. We discuss current challenges in the 3D reconstruction algorithms for bacterial chromosomes, primarily focusing on software usability. Finally, we briefly prospect future research directions for bacterial chromosome structure reconstruction algorithms.

KARR-seq reveals cellular higher-order RNA structures and RNA-RNA interactions

RNA fate and function are affected by their structures and interactomes. However, how RNA and RNA-binding proteins (RBPs) assemble into higher-order structures and how RNA molecules may interact with each other to facilitate functions remain largely unknown. Here we present KARR-seq, which uses N3-kethoxal labeling and multifunctional chemical crosslinkers to covalently trap and determine RNA-RNA interactions and higher-order RNA structures inside cells, independent of local protein binding to RNA. KARR-seq depicts higher-order RNA structure and detects widespread intermolecular RNA-RNA interactions with high sensitivity and accuracy. Using KARR-seq, we show that translation represses mRNA compaction under native and stress conditions. We determined the higher-order RNA structures of respiratory syncytial virus (RSV) and vesicular stomatitis virus (VSV) and identified RNA-RNA interactions between the viruses and the host RNAs that potentially regulate viral replication.

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.

Integrative modeling of JCVI-Syn3A nucleoids with a modular approach

A lattice-based method is presented for creating 3D models of entire bacterial nucleoids integrating ultrastructural information cryoelectron tomography, genomic and proteomic data, and experimental atomic structures of biomolecules and assemblies. The method is used to generate models of the minimal genome bacterium JCVI-Syn3A, producing a series of models that test hypotheses about transcription, condensation, and overall distribution of the genome. Lattice models are also used to generate atomic models of an entire JCVI-Syn3A cell.

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.

Implications of a multiscale structure of the yeast nuclear pore complex

Nuclear pore complexes (NPCs) direct the nucleocytoplasmic transport of macromolecules. Here, we provide a composite multiscale structure of the yeast NPC, based on improved 3D density maps from cryogenic electron microscopy and AlphaFold2 models. Key features of the inner and outer rings were integrated into a comprehensive model. We resolved flexible connectors that tie together the core scaffold, along with equatorial transmembrane complexes and a lumenal ring that anchor this channel within the pore membrane. The organization of the nuclear double outer ring reveals an architecture that may be shared with ancestral NPCs. Additional connections between the core scaffold and the central transporter suggest that under certain conditions, a degree of local organization is present at the periphery of the transport machinery. These connectors may couple conformational changes in the scaffold to the central transporter to modulate transport. Collectively, this analysis provides insights into assembly, transport, and NPC evolution.

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.

A deep learning model for predicting optimal distance range in crosslinking mass spectrometry data

Macromolecular assemblies play an important role in all cellular processes. While there has recently been significant progress in protein structure prediction based on deep learning, large protein complexes cannot be predicted with these approaches. The integrative structure modeling approach characterizes multi-subunit complexes by computational integration of data from fast and accessible experimental techniques. Crosslinking mass spectrometry is one such technique that provides spatial information about the proximity of crosslinked residues. One of the challenges in interpreting crosslinking datasets is designing a scoring function that, given a structure, can quantify how well it fits the data. Most approaches set an upper bound on the distance between Cα atoms of crosslinked residues and calculate a fraction of satisfied crosslinks. However, the distance spanned by the crosslinker greatly depends on the neighborhood of the crosslinked residues. Here, we design a deep learning model for predicting the optimal distance range for a crosslinked residue pair based on the structures of their neighborhoods. We find that our model can predict the distance range with the area under the receiver-operator curve of 0.86 and 0.7 for intra- and inter-protein crosslinks, respectively. Our deep scoring function can be used in a range of structure modeling applications.

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.

Contrastive learning in protein language space predicts interactions between drugs and protein targets

Sequence-based prediction of drug-target interactions has the potential to accelerate drug discovery by complementing experimental screens. Such computational prediction needs to be generalizable and scalable while remaining sensitive to subtle variations in the inputs. However, current computational techniques fail to simultaneously meet these goals, often sacrificing performance of one to achieve the others. We develop a deep learning model, ConPLex, successfully leveraging the advances in pretrained protein language models ("PLex") and employing a protein-anchored contrastive coembedding ("Con") to outperform state-of-the-art approaches. ConPLex achieves high accuracy, broad adaptivity to unseen data, and specificity against decoy compounds. It makes predictions of binding based on the distance between learned representations, enabling predictions at the scale of massive compound libraries and the human proteome. Experimental testing of 19 kinase-drug interaction predictions validated 12 interactions, including four with subnanomolar affinity, plus a strongly binding EPHB1 inhibitor (KD = 1.3 nM). Furthermore, ConPLex embeddings are interpretable, which enables us to visualize the drug-target embedding space and use embeddings to characterize the function of human cell-surface proteins. We anticipate that ConPLex will facilitate efficient drug discovery by making highly sensitive in silico drug screening feasible at the genome scale. ConPLex is available open source at https://ConPLex.csail.mit.edu.

Spt6 directly interacts with Cdc73 and is required for Paf1 complex occupancy at active genes in Saccharomyces cerevisiae

The Paf1 complex (Paf1C) is a conserved transcription elongation factor that regulates transcription elongation efficiency, facilitates co-transcriptional histone modifications, and impacts molecular processes linked to RNA synthesis, such as polyA site selection. Coupling of the activities of Paf1C to transcription elongation requires its association with RNA polymerase II (Pol II). Mutational studies in yeast identified Paf1C subunits Cdc73 and Rtf1 as important mediators of Paf1C association with Pol II on active genes. While the interaction between Rtf1 and the general elongation factor Spt5 is relatively well-understood, the interactions involving Cdc73 have not been fully elucidated. Using a site-specific protein cross-linking strategy in yeast cells, we identified direct interactions between Cdc73 and two components of the Pol II elongation complex, the elongation factor Spt6 and the largest subunit of Pol II. Both of these interactions require the tandem SH2 domain of Spt6. We also show that Cdc73 and Spt6 can interact in vitro and that rapid depletion of Spt6 dissociates Paf1 from chromatin, altering patterns of Paf1C-dependent histone modifications genome-wide. These results reveal interactions between Cdc73 and the Pol II elongation complex and identify Spt6 as a key factor contributing to the occupancy of Paf1C at active genes in Saccharomyces cerevisiae.

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.

Molecular structure of soluble vimentin tetramers

Intermediate filaments (IFs) are essential constituents of the metazoan cytoskeleton. A vast family of cytoplasmic IF proteins are capable of self-assembly from soluble tetrameric species into typical 10-12 nm wide filaments. The primary structure of these proteins includes the signature central 'rod' domain of ~ 300 residues which forms a dimeric α-helical coiled coil composed of three segments (coil1A, coil1B and coil2) interconnected by non-helical, flexible linkers (L1 and L12). The rod is flanked by flexible terminal head and tail domains. At present, the molecular architecture of mature IFs is only poorly known, limiting our capacity to rationalize the effect of numerous disease-related mutations found in IF proteins. Here we addressed the molecular structure of soluble vimentin tetramers which are formed by two antiparallel, staggered dimers with coil1B domains aligned (A₁₁ tetramers). By examining a series of progressive truncations, we show that the presence of the coil1A domain is essential for the tetramer formation. In addition, we employed a novel chemical cross-linking pipeline including isotope labelling to identify intra- and interdimeric cross-links within the tetramer. We conclude that the tetramer is synergistically stabilized by the interactions of the aligned coil1B domains, the interactions between coil1A and the N-terminal portion of coil2, and the electrostatic attraction between the oppositely charged head and rod domains. Our cross-linking data indicate that, starting with a straight A₁₁ tetramer, flexibility of linkers L1 and L12 enables 'backfolding' of both the coil1A and coil2 domains onto the tetrameric core formed by the coil1B domains. Through additional small-angle X-ray scattering experiments we show that the elongated A₁₁ tetramers dominate in low ionic strength solutions, while there is also a significant structural flexibility especially in the terminal domains.

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.

CryoFold 2.0: Cryo-EM Structure Determination with MELD

Cryo-electron microscopy data are becoming more prevalent and accessible at higher resolution levels, leading to the development of new computational tools to determine the atomic structure of macromolecules. However, while existing tools adapted from X-ray crystallography are suitable for the highest-resolution maps, new tools are needed for lower-resolution levels and to account for map heterogeneity. In this article, we introduce CryoFold 2.0, an integrative physics-based approach that combines Bayesian inference and the ability to handle multiple data sources with the molecular dynamics flexible fitting (MDFF) approach to determine the structures of macromolecules by using cryo-EM data. CryoFold 2.0 is incorporated into the MELD (modeling employing limited data) plugin, resulting in a pipeline that is more computationally efficient and accurate than running MELD or MDFF alone. The approach requires fewer computational resources and shorter simulation times than the original CryoFold, and it minimizes manual intervention. We demonstrate the effectiveness of the approach on eight different systems, highlighting its various benefits.

Breaking free from the crystal lattice: Structural biology in solution to study protein degraders

Structural biology offers a versatile arsenal of techniques and methods to investigate the structure and conformational dynamics of proteins and their assemblies. The growing field of targeted protein degradation centres on the premise of developing small molecules, termed degraders, to induce proximity between an E3 ligase and a protein of interest to be signalled for degradation. This new drug modality brings with it new opportunities and challenges to structural biologists. Here we discuss how several structural biology techniques, including nuclear magnetic resonance, cryo-electron microscopy, structural mass spectrometry and small angle scattering, have been explored to complement X-ray crystallography in studying degraders and their ternary complexes. Together the studies covered in this review make a case for the invaluable perspectives that integrative structural biology techniques in solution can bring to understanding ternary complexes and designing degraders.

Cross-linking mass spectrometry for mapping protein complex topologies in situ

Cross-linking mass spectrometry has become an established technology to provide structural information on the topology and dynamics of protein complexes. Readily accessible workflows can provide detailed data on simplified systems, such as purified complexes. However, using this technology to study the structure of protein complexes in situ, such as in organelles, cells, and even tissues, is still a technological frontier. The complexity of these systems remains a considerable challenge, but there have been dramatic improvements in sample handling, data acquisition, and data processing. Here, we summarise these developments and describe the paths towards comprehensive and comparative structural interactomes by cross-linking mass spectrometry.

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.

Modeling of Protein Complexes

The recent advances in structural biology, combined with continuously increasing computational capabilities and development of advanced softwares, have drastically simplified the workflow for protein homology modeling. Modeling of individual proteins is nowadays quick and straightforward for a large variety of protein targets, thanks to guided pipelines relying on advanced computational tools and user-friendly interfaces, which have extended and promoted the use of modeling also to scientists not focusing on molecular structures of proteins. Nevertheless, construction of models of multi-protein complexes remains quite challenging for the non-experts, often due to the usage of specific procedures depending on the system under investigation and the need for experimental validation approaches to strengthen the generated output.In this chapter, we provide a brief overview of the approaches enabling generation of multi-protein complex models starting from homology models of individual protein components. Using real-life examples, we include two examples to guide the reader in the generation of homomeric and heteromeric protein models.

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.

Electron microscopy holdings of the Protein Data Bank: the impact of the resolution revolution, new validation tools, and implications for the future

As a discipline, structural biology has been transformed by the three-dimensional electron microscopy (3DEM) "Resolution Revolution" made possible by convergence of robust cryo-preservation of vitrified biological materials, sample handling systems, and measurement stages operating a liquid nitrogen temperature, improvements in electron optics that preserve phase information at the atomic level, direct electron detectors (DEDs), high-speed computing with graphics processing units, and rapid advances in data acquisition and processing software. 3DEM structure information (atomic coordinates and related metadata) are archived in the open-access Protein Data Bank (PDB), which currently holds more than 11,000 3DEM structures of proteins and nucleic acids, and their complexes with one another and small-molecule ligands (~ 6% of the archive). Underlying experimental data (3DEM density maps and related metadata) are stored in the Electron Microscopy Data Bank (EMDB), which currently holds more than 21,000 3DEM density maps. After describing the history of the PDB and the Worldwide Protein Data Bank (wwPDB) partnership, which jointly manages both the PDB and EMDB archives, this review examines the origins of the resolution revolution and analyzes its impact on structural biology viewed through the lens of PDB holdings. Six areas of focus exemplifying the impact of 3DEM across the biosciences are discussed in detail (icosahedral viruses, ribosomes, integral membrane proteins, SARS-CoV-2 spike proteins, cryogenic electron tomography, and integrative structure determination combining 3DEM with complementary biophysical measurement techniques), followed by a review of 3DEM structure validation by the wwPDB that underscores the importance of community engagement.

Integrative modeling reveals the molecular architecture of the intraflagellar transport A (IFT-A) complex

Intraflagellar transport (IFT) is a conserved process of cargo transport in cilia that is essential for development and homeostasis in organisms ranging from algae to vertebrates. In humans, variants in genes encoding subunits of the cargo-adapting IFT-A and IFT-B protein complexes are a common cause of genetic diseases known as ciliopathies. While recent progress has been made in determining the atomic structure of IFT-B, little is known of the structural biology of IFT-A. Here, we combined chemical cross-linking mass spectrometry and cryo-electron tomography with AlphaFold2-based prediction of both protein structures and interaction interfaces to model the overall architecture of the monomeric six-subunit IFT-A complex, as well as its polymeric assembly within cilia. We define monomer-monomer contacts and membrane-associated regions available for association with transported cargo, and we also use this model to provide insights into the pleiotropic nature of human ciliopathy-associated genetic variants in genes encoding IFT-A subunits. Our work demonstrates the power of integration of experimental and computational strategies both for multi-protein structure determination and for understanding the etiology of human genetic disease.

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.

Predicting the structure of large protein complexes using AlphaFold and Monte Carlo tree search

AlphaFold can predict the structure of single- and multiple-chain proteins with very high accuracy. However, the accuracy decreases with the number of chains, and the available GPU memory limits the size of protein complexes which can be predicted. Here we show that one can predict the structure of large complexes starting from predictions of subcomponents. We assemble 91 out of 175 complexes with 10-30 chains from predicted subcomponents using Monte Carlo tree search, with a median TM-score of 0.51. There are 30 highly accurate complexes (TM-score ≥0.8, 33% of complete assemblies). We create a scoring function, mpDockQ, that can distinguish if assemblies are complete and predict their accuracy. We find that complexes containing symmetry are accurately assembled, while asymmetrical complexes remain challenging. The method is freely available and accesible as a Colab notebook https://colab.research.google.com/github/patrickbryant1/MoLPC/blob/master/MoLPC.ipynb .

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.

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.

Exploring the Conformations and Binding Location of HMGA2·DNA Complexes Using Ion Mobility Spectrometry and 193 nm Ultraviolet Photodissociation Mass Spectrometry

Although it is widely accepted that protein function is largely dependent on its structure, intrinsically disordered proteins (IDPs) lack defined structure but are essential in proper cellular processes. Mammalian high mobility group proteins (HMGA) are one such example of IDPs that perform a number of crucial nuclear activities and have been highly studied due to their involvement in the proliferation of a variety of disease and cancers. Traditional structural characterization methods have had limited success in understanding HMGA proteins and their ability to coordinate to DNA. Ion mobility spectrometry and mass spectrometry provide insights into the diversity and heterogeneity of structures adopted by IDPs and are employed here to interrogate HMGA2 in its unbound states and bound to two DNA hairpins. The broad distribution of collision cross sections observed for the apo-protein are restricted when HMGA2 is bound to DNA, suggesting that increased protein organization is promoted in the holo-form. Ultraviolet photodissociation was utilized to probe the changes in structures for the compact and elongated structures of HMGA2 by analyzing backbone cleavage propensities and solvent accessibility based on charge-site analysis, which revealed a spectrum of conformational possibilities. Namely, preferential binding of the DNA hairpins with the second of three AT-hooks of HMGA2 is suggested based on the suppression of backbone fragmentation and distribution of DNA-containing protein fragments.

Recent advances in RNA structurome

RNA structures are essential to support RNA functions and regulation in various biological processes. Recently, a range of novel technologies have been developed to decode genome-wide RNA structures and novel modes of functionality across a wide range of species. In this review, we summarize key strategies for probing the RNA structurome and discuss the pros and cons of representative technologies. In particular, these new technologies have been applied to dissect the structural landscape of the SARS-CoV-2 RNA genome. We also summarize the functionalities of RNA structures discovered in different regulatory layers-including RNA processing, transport, localization, and mRNA translation-across viruses, bacteria, animals, and plants. We review many versatile RNA structural elements in the context of different physiological and pathological processes (e.g., cell differentiation, stress response, and viral replication). Finally, we discuss future prospects for RNA structural studies to map the RNA structurome at higher resolution and at the single-molecule and single-cell level, and to decipher novel modes of RNA structures and functions for innovative applications.

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.