Modeling protein assemblies in the proteome

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.

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.

Modeling the Dynamics of Protein–Protein Interfaces, How and Why?

Protein–protein assemblies act as a key component in numerous cellular processes. Their accurate modeling at the atomic level remains a challenge for structural biology. To address this challenge, several docking and a handful of deep learning methodologies focus on modeling protein–protein interfaces. Although the outcome of these methods has been assessed using static reference structures, more and more data point to the fact that the interaction stability and specificity is encoded in the dynamics of these interfaces. Therefore, this dynamics information must be taken into account when modeling and assessing protein interactions at the atomistic scale. Expanding on this, our review initially focuses on the recent computational strategies aiming at investigating protein–protein interfaces in a dynamic fashion using enhanced sampling, multi-scale modeling, and experimental data integration. Then, we discuss how interface dynamics report on the function of protein assemblies in globular complexes, in fuzzy complexes containing intrinsically disordered proteins, as well as in active complexes, where chemical reactions take place across the protein–protein interface.

Structural analysis and functional study of phosphofructokinase B (PfkB) from Mycobacterium marinum

Phosphofructokinase B (PfkB) belongs to the ribokinase family, which uses the phosphorylated sugar as substrate, and catalyzes fructose-6-phosphate into fructose-1,6-diphosphate. However, the structural basis of Mycobacterium marinum PfkB is not clear. Here, we found that the PfkB protein was monomeric in solution, which was different from most enzymes in this family. The crystal structure of PfkB protein from M. marinum was solved at a resolution of 2.21 Å. The PfkB structure consists of two domains, a major three-layered α/β/α sandwich-like domain characteristic of the ribokinase-like superfamily, and a second domain composed of four-stranded β sheets. Structural comparison analysis suggested that residues G236, A237, G238, and D239 could be critical for ATP catalysis and substrate binding of PfkB. Our current work provides new insights into understanding the mechanism of the glycolysis in M. marinum.

Moving pictures: Reassessing docking experiments with a dynamic view of protein interfaces

The modeling of protein assemblies at the atomic level remains a central issue in structural biology, as protein interactions play a key role in numerous cellular processes. This problem is traditionally addressed using docking tools, where the quality of the models is based on their similarity to a single reference experimental structure. However, using a static reference does not take into account the dynamic quality of the protein interface. Here, we used all-atom classical Molecular Dynamics simulations to investigate the stability of the reference interface for three complexes that previously served as targets in the CAPRI competition. For each one of these targets, we also ran MD simulations for ten models that are distributed over the High, Medium and Acceptable accuracy categories. To assess the quality of these models from a dynamic perspective, we set up new criteria which take into account the stability of the reference experimental protein interface. We show that, when the protein interfaces are allowed to evolve along time, the original ranking based on the static CAPRI criteria no longer holds as over 50% of the docking models undergo a category change (which can be either toward a better or a lower accuracy group) when reassessing their quality using dynamic information.

Modeling of Multimolecular Complexes

Macromolecular complexes play a key role in cellular function. Predicting the structure and dynamics of these complexes is one of the key challenges in structural biology. Docking applications have traditionally been used to predict pairwise interactions between proteins. However, few methods exist for modeling multi-protein assemblies. Here we present two methods, CombDock and DockStar, that can predict multi-protein assemblies starting from subunit structural models. CombDock can assemble subunits without any assumptions about the pairwise interactions between subunits, while DockStar relies on the interaction graph or, alternatively, a homology model or a cryo-electron microscopy (EM) density map of the entire complex. We demonstrate the two methods using RNA polymerase II with 12 subunits and TRiC/CCT chaperonin with 16 subunits.

OncoBinder facilitates interpretation of proteomic interaction data by capturing coactivation pairs in cancer

High-throughput methods such as co-immunoprecipitationmass spectrometry (coIP-MS) and yeast 2 hybridization (Y2H) have suggested a broad range of unannotated protein-protein interactions (PPIs), and interpretation of these PPIs remains a challenging task. The advancements in cancer genomic researches allow for the inference of "coactivation pairs" in cancer, which may facilitate the identification of PPIs involved in cancer. Here we present OncoBinder as a tool for the assessment of proteomic interaction data based on the functional synergy of oncoproteins in cancer. This decision tree-based method combines gene mutation, copy number and mRNA expression information to infer the functional status of protein-coding genes. We applied OncoBinder to evaluate the potential binders of EGFR and ERK2 proteins based on the gastric cancer dataset of The Cancer Genome Atlas (TCGA). As a result, OncoBinder identified high confidence interactions (annotated by Kyoto Encyclopedia of Genes and Genomes (KEGG) or validated by low-throughput assays) more efficiently than co-expression based method. Taken together, our results suggest that evaluation of gene functional synergy in cancer may facilitate the interpretation of proteomic interaction data. The OncoBinder toolbox for Matlab is freely accessible online.

PRISM-EM: template interface-based modelling of multi-protein complexes guided by cryo-electron microscopy density maps

The structures of protein assemblies are important for elucidating cellular processes at the molecular level. Three-dimensional electron microscopy (3DEM) is a powerful method to identify the structures of assemblies, especially those that are challenging to study by crystallography. Here, a new approach, PRISM-EM, is reported to computationally generate plausible structural models using a procedure that combines crystallographic structures and density maps obtained from 3DEM. The predictions are validated against seven available structurally different crystallographic complexes. The models display mean deviations in the backbone of <5 Å. PRISM-EM was further tested on different benchmark sets; the accuracy was evaluated with respect to the structure of the complex, and the correlation with EM density maps and interface predictions were evaluated and compared with those obtained using other methods. PRISM-EM was then used to predict the structure of the ternary complex of the HIV-1 envelope glycoprotein trimer, the ligand CD4 and the neutralizing protein m36.

Large oligomeric complex structures can be computationally assembled by efficiently combining docked interfaces

Macromolecular oligomeric assemblies are involved in many biochemical processes of living organisms. The benefits of such assemblies in crowded cellular environments include increased reaction rates, efficient feedback regulation, cooperativity and protective functions. However, an atom-level structural determination of large assemblies is challenging due to the size of the complex and the difference in binding affinities of the involved proteins. In this study, we propose a novel combinatorial greedy algorithm for assembling large oligomeric complexes from information on the approximate position of interaction interfaces of pairs of monomers in the complex. Prior information on complex symmetry is not required but rather the symmetry is inferred during assembly. We implement an efficient geometric score, the transformation match score, that bypasses the model ranking problems of state-of-the-art scoring functions by scoring the similarity between the inferred dimers of the same monomer simultaneously with different binding partners in a (sub)complex with a set of pregenerated docking poses. We compiled a diverse benchmark set of 308 homo and heteromeric complexes containing 6 to 60 monomers. To explore the applicability of the method, we considered 48 sets of parameters and selected those three sets of parameters, for which the algorithm can correctly reconstruct the maximum number, namely 252 complexes (81.8%) in, at least one of the respective three runs. The crossvalidation coverage, that is, the mean fraction of correctly reconstructed benchmark complexes during crossvalidation, was 78.1%, which demonstrates the ability of the presented method to correctly reconstruct topology of a large variety of biological complexes.

DockStar: a novel ILP-based integrative method for structural modeling of multimolecular protein complexes

MOTIVATION

Atomic resolution modeling of large multimolecular assemblies is a key task in Structural Cell Biology. Experimental techniques can provide atomic resolution structures of single proteins and small complexes, or low resolution data of large multimolecular complexes.

RESULTS

We present a novel integrative computational modeling method, which integrates both low and high resolution experimental data. The algorithm accepts as input atomic resolution structures of the individual subunits obtained from X-ray, NMR or homology modeling, and interaction data between the subunits obtained from mass spectrometry. The optimal assembly of the individual subunits is formulated as an Integer Linear Programming task. The method was tested on several representative complexes, both in the bound and unbound cases. It placed correctly most of the subunits of multimolecular complexes of up to 16 subunits and significantly outperformed the CombDock and Haddock multimolecular docking methods.

AVAILABILITY AND IMPLEMENTATION

http://bioinfo3d.cs.tau.ac.il/DockStar

CONTACT

naamaamir@mail.tau.ac.il or wolfson@tau.ac.il

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The structural basis for cancer treatment decisions

Cancer treatment decisions rely on genetics, large data screens and clinical pharmacology. Here we point out that genetic analysis and treatment decisions may overlook critical elements in cancer development, progression and drug resistance. Two critical structural elements are missing in genetics-based decision-making: the mechanisms of oncogenic mutations and the cellular network which is rewired in cancer. These lay the foundation for the structural basis for cancer treatment decisions, which is rooted in the physical principles of the molecular conformational behavior of single molecules and their interactions. Improved tumor mutational analysis platforms and knowledge of the redundant pathways which can take over in cancer, may not only supplement known actionable findings, but forecast possible cancer progression and resistance. Such forward-looking can be powerful, endowing the oncologist with mechanistic insight and cancer prognosis, and consequently more informed treatment options. Examples include redundant pathways taking over after inhibition of EGFR constitutive activation, mutations in PIK3CA p110α and p85, and the non-hotspot AKT1 mutants conferring constitutive membrane localization.

ProBiS-ligands: a web server for prediction of ligands by examination of protein binding sites

The ProBiS-ligands web server predicts binding of ligands to a protein structure. Starting with a protein structure or binding site, ProBiS-ligands first identifies template proteins in the Protein Data Bank that share similar binding sites. Based on the superimpositions of the query protein and the similar binding sites found, the server then transposes the ligand structures from those sites to the query protein. Such ligand prediction supports many activities, e.g. drug repurposing. The ProBiS-ligands web server, an extension of the ProBiS web server, is open and free to all users at http://probis.cmm.ki.si/ligands.

Structural pathways of cytokines may illuminate their roles in regulation of cancer development and immunotherapy

Cytokines are messengers between tissues and the immune system. They play essential roles in cancer initiation, promotion, metastasis, and immunotherapy. Structural pathways of cytokine signaling which contain their interactions can help understand their action in the tumor microenvironment. Here, our aim is to provide an overview of the role of cytokines in tumor development from a structural perspective. Atomic details of protein-protein interactions can help in understanding how an upstream signal is transduced; how higher-order oligomerization modes of proteins can influence their function; how mutations, inhibitors or antagonists can change cellular consequences; why the same protein can lead to distinct outcomes, and which alternative parallel pathways can take over. They also help to design drugs/inhibitors against proteins de novo or by mimicking natural antagonists as in the case of interferon-γ. Since the structural database (PDB) is limited, structural pathways are largely built from a series of predicted binary protein-protein interactions. Below, to illustrate how protein-protein interactions can help illuminate roles played by cytokines, we model some cytokine interaction complexes exploiting a powerful algorithm (PRotein Interactions by Structural Matching-PRISM).