Template-based protein–protein docking exploiting pairwise interfacial residue restraints

Although many advanced and sophisticated ab initio approaches for modeling protein–protein complexes have been proposed in past decades, template-based modeling (TBM) remains the most accurate and widely used approach, given a reliable template is available. However, there are many different ways to exploit template information in the modeling process. Here, we systematically evaluate and benchmark a TBM method that uses conserved interfacial residue pairs as docking distance restraints [referred to as alpha carbon–alpha carbon (CA-CA)-guided docking]. We compare it with two other template-based protein–protein modeling approaches, including a conserved non-pairwise interfacial residue restrained docking approach [referred to as the ambiguous interaction restraint (AIR)-guided docking] and a simple superposition-based modeling approach. Our results show that, for most cases, the CA-CA-guided docking method outperforms both superposition with refinement and the AIR-guided docking method. We emphasize the superiority of the CA-CA-guided docking on cases with medium to large conformational changes, and interactions mediated through loops, tails or disordered regions. Our results also underscore the importance of a proper refinement of superimposition models to reduce steric clashes. In summary, we provide a benchmarked TBM protocol that uses conserved pairwise interface distance as restraints in generating realistic 3D protein–protein interaction models, when reliable templates are available. The described CA-CA-guided docking protocol is based on the HADDOCK platform, which allows users to incorporate additional prior knowledge of the target system to further improve the quality of the resulting models

Exploring the potential of 3D Zernike descriptors and SVM for protein\u2013protein interface prediction

Abstract Background The correct determination of protein–protein interaction interfaces is important for understanding disease mechanisms and for rational drug design. To date, several computational methods for the prediction of protein interfaces have been developed, but the interface prediction problem is still not fully understood. Experimental evidence suggests that the location of binding sites is imprinted in the protein structure, but there are major differences among the interfaces of the various protein types: the characterising properties can vary a lot depending on the interaction type and function. The selection of an optimal set of features characterising the protein interface and the development of an effective method to represent and capture the complex protein recognition patterns are of paramount importance for this task. Results In this work we investigate the potential of a novel local surface descriptor based on 3D Zernike moments for the interface prediction task. Descriptors invariant to roto-translations are extracted from circular patches of the protein surface enriched with physico-chemical properties from the HQI8 amino acid index set, and are used as samples for a binary classification problem. Support Vector Machines are used as a classifier to distinguish interface local surface patches from non-interface ones. The proposed method was validated on 16 classes of proteins extracted from the Protein–Protein Docking Benchmark 5.0 and compared to other state-of-the-art protein interface predictors (SPPIDER, PrISE and NPS-HomPPI). Conclusions The 3D Zernike descriptors are able to capture the similarity among patterns of physico-chemical and biochemical properties mapped on the protein surface arising from the various spatial arrangements of the underlying residues, and their usage can be easily extended to other sets of amino acid properties. The results suggest that the choice of a proper set of features characterising the protein interface is crucial for the interface prediction task, and that optimality strongly depends on the class of proteins whose interface we want to characterise. We postulate that different protein classes should be treated separately and that it is necessary to identify an optimal set of features for each protein class

Methods for Molecular Modelling of Protein Complexes.

Biological processes are often mediated by complexes formed between proteins and various biomolecules. The 3D structures of such protein-biomolecule complexes provide insights into the molecular mechanism of their action. The structure of these complexes can be predicted by various computational methods. Choosing an appropriate method for modelling depends on the category of biomolecule that a protein interacts with and the availability of structural information about the protein and its interacting partner. We intend for the contents of this chapter to serve as a guide as to what software would be the most appropriate for the type of data at hand and the kind of 3D complex structure required. Particularly, we have dealt with protein-small molecule ligand, protein-peptide, protein-protein, and protein-nucleic acid interactions.Most, if not all, model building protocols perform some sampling and scoring. Typically, several alternate conformations and configurations of the interactors are sampled. Each such sample is then scored for optimization. To boost the confidence in these predicted models, their assessment using other independent scoring schemes besides the inbuilt/default ones would prove to be helpful. This chapter also lists such software and serves as a guide to gauge the fidelity of modelled structures of biomolecular complexes

Template-based protein–protein docking exploiting pairwise interfacial residue restraints

Although many advanced and sophisticated ab initio approaches for modeling protein–protein complexes have been proposed in past decades, template-based modeling (TBM) remains the most accurate and widely used approach, given a reliable template is available. However, there are many different ways to exploit template information in the modeling process. Here, we systematically evaluate and benchmark a TBM method that uses conserved interfacial residue pairs as docking distance restraints [referred to as alpha carbon–alpha carbon (CA-CA)-guided docking]. We compare it with two other template-based protein–protein modeling approaches, including a conserved non-pairwise interfacial residue restrained docking approach [referred to as the ambiguous interaction restraint (AIR)-guided docking] and a simple superposition-based modeling approach. Our results show that, for most cases, the CA-CA-guided docking method outperforms both superposition with refinement and the AIR-guided docking method. We emphasize the superiority of the CA-CA-guided docking on cases with medium to large conformational changes, and interactions mediated through loops, tails or disordered regions. Our results also underscore the importance of a proper refinement of superimposition models to reduce steric clashes. In summary, we provide a benchmarked TBM protocol that uses conserved pairwise interface distance as restraints in generating realistic 3D protein–protein interaction models, when reliable templates are available. The described CA-CA-guided docking protocol is based on the HADDOCK platform, which allows users to incorporate additional prior knowledge of the target system to further improve the quality of the resulting models.This article is from Briefings in Bioinformatics (2016), doi: 10.1093/bib/bbw027. Posted with permission.</p

Template-based protein-protein docking exploiting pairwise interfacial residue restraints

Although many advanced and sophisticatedab initioapproaches for modeling protein-protein complexes have been proposed in past decades, template-based modeling (TBM) remains the most accurate and widely used approach, given a reliable template is available. However, there are many different ways to exploit template information in the modeling process. Here, we systematically evaluate and benchmark a TBM method that uses conserved interfacial residue pairs as docking distance restraints [referred to as alpha carbon-alpha carbon (CA-CA)-guided docking]. We compare it with two other template-based protein-protein modeling approaches, including a conserved non-pairwise interfacial residue restrained docking approach [referred to as the ambiguous interaction restraint (AIR)-guided docking] and a simple superposition-based modeling approach. Our results show that, for most cases, the CA-CA-guided docking method outperforms both superposition with refinement and the AIR-guided docking method. We emphasize the superiority of the CA-CA-guided docking on cases with medium to large conformational changes, and interactions mediated through loops, tails or disordered regions. Our results also underscore the importance of a proper refinement of superimposition models to reduce steric clashes. In summary, we provide a benchmarked TBM protocol that uses conserved pairwise interface distance as restraints in generating realistic 3D protein-protein interaction models, when reliable templates are available. The described CA-CA-guided docking protocol is based on the HADDOCK platform, which allows users to incorporate additional prior knowledge of the target system to further improve the quality of the resulting models

Modeling Interactions of Flexible Proteins

Proteins are dynamic molecules that mediate most biological processes through interactions with other proteins and biomolecules. A fundamental understanding of the mechanisms governing protein interactions requires intricate knowledge of the three-dimensional structures of biomolecular complexes. Despite advances in experimental structure determination, we have structural insights into only a small fraction of known complexes. Computational modeling provides an invaluable complementary tool to explore protein interactions in a rapid and high-throughput manner. A principal challenge limiting the accuracy of current computational methods is the ability to predict binding-induced conformational changes during protein–protein association. In this dissertation, I address this challenge by creating new tools to predict atomistic models of flexible protein complexes. First, I develop a heterodimer docking protocol that incorporates flexibility by efficiently simulating conformational selection from hundreds of pre-generated backbone conformations and identifies the near-native models with a novel, coarse-grained score function called Motif Dock Score (MDS). On a benchmark of 88 complexes with different degrees of flexibility, this protocol, RosettaDock 4.0, is the first method to successfully dock approximately 50% of complexes with conformational change of up to 2.2 Å. Next, I present the results of our participation in the community-wide blind experiment, Critical Assessment of PRedicted Interactions (CAPRI) rounds 37–45, where I use various docking methods to predict the structures of protein homomer, heteromer and oligosaccharide complexes. In the process, I identify inadequacies in these methods and propose enhancements. Based on the shortcomings identified in CAPRI, I develop a protocol to predict the structure of symmetric homomers from monomeric inputs with a focus on tightly-packed complexes. This method, Rosetta SymDock2, leverages MDS in the coarse-grained phase and simulates subunit flexibility through induced fit by all-atom flexible-backbone refinement. It outperforms competing algorithms by docking 61% of cyclic complexes and 42% of dihedral complexes in a diverse benchmark of 43 homomers. In the course of developing these algorithms, I also discover that the binding energy wells of homomers are narrower, steeper and deeper than those of heterodimers, thus explaining their increased stability. Finally, I present preliminary results to propose data-driven strategies that can overcome current barriers to accurate modeling