Intuitive,but not simple: Including explicit water molecules in protein–protein docking simulations improves model quality

Characterizing the nature of interaction between proteins that have not been experimentally cocrystallized requires a computational docking approach that can successfully predict the spatial conformation adopted in the complex. In this work, the Hydropathic INTeractions (HINT) force field model was used for scoring docked models in a data set of 30 high‐resolution crystallographically characterized “dry” protein–protein complexes and was shown to reliably identify native‐like models. However, most current protein–protein docking algorithms fail to explicitly account for water molecules involved in bridging interactions that mediate and stabilize the association of the protein partners, so we used HINT to illuminate the physical and chemical properties of bridging waters and account for their energetic stabilizing contributions. The HINT water Relevance metric identified the “truly” bridging waters at the 30 protein–protein interfaces and we utilized them in “solvated” docking by manually inserting them into the input files for the rigid body ZDOCK program. By accounting for these interfacial waters, a statistically significant improvement of ～24% in the average hit‐count within the top‐10 predictions the protein–protein dataset was seen, compared to standard “dry” docking. The results also show scoring improvement, with medium and high accuracy models ranking much better than incorrect ones. These improvements can be attributed to the physical presence of water molecules that alter surface properties and better represent native shape and hydropathic complementarity between interacting partners, with concomitantly more accurate native‐like structure predictions. Proteins 2014; 82:916–932. © 2013 Wiley Periodicals, Inc.     

Interolog interfaces in protein–protein docking

Proteins are essential elements of biological systems, and their function typically relies on their ability to successfully bind to specific partners. Recently, an emphasis of study into protein interactions has been on hot spots, or residues in the binding interface that make a significant contribution to the binding energetics. In this study, we investigate how conservation of hot spots can be used to guide docking prediction. We show that the use of evolutionary data combined with hot spot prediction highlights near‐native structures across a range of benchmark examples. Our approach explores various strategies for using hot spots and evolutionary data to score protein complexes, using both absolute and chemical definitions of conservation along with refinements to these strategies that look at windowed conservation and filtering to ensure a minimum number of hot spots in each binding partner. Finally, structure‐based models of orthologs were generated for comparison with sequence‐based scoring. Using two data sets of 22 and 85 examples, a high rate of top 10 and top 1 predictions are observed, with up to 82% of examples returning a top 10 hit and 35% returning top 1 hit depending on the data set and strategy applied; upon inclusion of the native structure among the decoys, up to 55% of examples yielded a top 1 hit. The 20 common examples between data sets show that more carefully curated interolog data yields better predictions, particularly in achieving top 1 hits. Proteins 2015; 83:1940–1946. © 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

A non‐redundant protein–RNA docking benchmark version 2.0

We present an updated version of the protein–RNA docking benchmark, which we first published four years back. The non‐redundant protein–RNA docking benchmark version 2.0 consists of 126 test cases, a threefold increase in number compared to its previous version. The present version consists of 21 unbound–unbound cases, of which, in 12 cases, the unbound RNAs are taken from another complex. It also consists of 95 unbound–bound cases where only the protein is available in the unbound state. Besides, we introduce 10 new bound–unbound cases where only the RNA is found in the unbound state. Based on the degree of conformational change of the interface residues upon complex formation the benchmark is classified into 72 rigid‐body cases, 25 semiflexible cases and 19 full flexible cases. It also covers a wide range of conformational flexibility including small side chain movement to large domain swapping in protein structures as well as flipping and restacking in RNA bases. This benchmark should provide the docking community with more test cases for evaluating rigid‐body as well as flexible docking algorithms. Besides, it will also facilitate the development of new algorithms that require large number of training set. The protein–RNA docking benchmark version 2.0 can be freely downloaded from http://www.csb.iitkgp.ernet.in/applications/PRDBv2 . Proteins 2017; 85:256–267. © 2016 Wiley Periodicals, Inc.     

iATTRACT: Simultaneous global and local interface optimization for protein–protein docking refinement

Protein‐protein interactions are abundant in the cell but to date structural data for a large number of complexes is lacking. Computational docking methods can complement experiments by providing structural models of complexes based on structures of the individual partners. A major caveat for docking success is accounting for protein flexibility. Especially, interface residues undergo significant conformational changes upon binding. This limits the performance of docking methods that keep partner structures rigid or allow limited flexibility. A new docking refinement approach, iATTRACT, has been developed which combines simultaneous full interface flexibility and rigid body optimizations during docking energy minimization. It employs an atomistic molecular mechanics force field for intermolecular interface interactions and a structure‐based force field for intramolecular contributions. The approach was systematically evaluated on a large protein‐protein docking benchmark, starting from an enriched decoy set of rigidly docked protein–protein complexes deviating by up to 15 Å from the native structure at the interface. Large improvements in sampling and slight but significant improvements in scoring/discrimination of near native docking solutions were observed. Complexes with initial deviations at the interface of up to 5.5 Å were refined to significantly better agreement with the native structure. Improvements in the fraction of native contacts were especially favorable, yielding increases of up to 70%. Proteins 2015; 83:248–258. © 2014 Wiley Periodicals, Inc.     

Discovery of binding proteins for a protein target using protein–protein docking‐based virtual screening

Target structure‐based virtual screening, which employs protein‐small molecule docking to identify potential ligands, has been widely used in small‐molecule drug discovery. In the present study, we used a protein–protein docking program to identify proteins that bind to a specific target protein. In the testing phase, an all‐to‐all protein–protein docking run on a large dataset was performed. The three‐dimensional rigid docking program SDOCK was used to examine protein–protein docking on all protein pairs in the dataset. Both the binding affinity and features of the binding energy landscape were considered in the scoring function in order to distinguish positive binding pairs from negative binding pairs. Thus, the lowest docking score, the average Z‐score, and convergency of the low‐score solutions were incorporated in the analysis. The hybrid scoring function was optimized in the all‐to‐all docking test. The docking method and the hybrid scoring function were then used to screen for proteins that bind to tumor necrosis factor‐α (TNFα), which is a well‐known therapeutic target for rheumatoid arthritis and other autoimmune diseases. A protein library containing 677 proteins was used for the screen. Proteins with scores among the top 20% were further examined. Sixteen proteins from the top‐ranking 67 proteins were selected for experimental study. Two of these proteins showed significant binding to TNFα in an in vitro binding study. The results of the present study demonstrate the power and potential application of protein–protein docking for the discovery of novel binding proteins for specific protein targets. Proteins 2014; 82:2472–2482. © 2014 Wiley Periodicals, Inc.     

DynaDock: A new molecular dynamics‐based algorithm for protein–peptide docking including receptor flexibility

Molecular docking programs play an important role in drug development and many well‐established methods exist. However, there are two situations for which the performance of most approaches is still not satisfactory, namely inclusion of receptor flexibility and docking of large, flexible ligands like peptides. In this publication a new approach is presented for docking peptides into flexible receptors. For this purpose a two step procedure was developed: first, the protein–peptide conformational space is scanned and approximate ligand poses are identified and second, the identified ligand poses are refined by a new molecular dynamics‐based method, optimized potential molecular dynamics (OPMD). The OPMD approach uses soft‐core potentials for the protein–peptide interactions and applies a new optimization scheme to the soft‐core potential. Comparison with refinement results obtained by conventional molecular dynamics and a soft‐core scaling approach shows significant improvements in the sampling capability for the OPMD method. Thus, the number of starting poses needed for successful refinement is much lower than for the other methods. The algorithm was evaluated on 15 protein–peptide complexes with 2–16mer peptides. Docking poses with peptide RMSD values <2.10 Å from the equilibrated experimental structures were obtained in all cases. For four systems docking into the unbound receptor structures was performed, leading to peptide RMSD values <2.12 Å. Using a specifically fitted scoring function in 11 of 15 cases the best scoring poses featured a peptide RMSD ≤2.10 Å. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

DockRank: Ranking docked conformations using partner‐specific sequence homology‐based protein interface prediction

Selecting near‐native conformations from the immense number of conformations generated by docking programs remains a major challenge in molecular docking. We introduce DockRank, a novel approach to scoring docked conformations based on the degree to which the interface residues of the docked conformation match a set of predicted interface residues. DockRank uses interface residues predicted by partner‐specific sequence homology‐based protein–protein interface predictor (PS‐HomPPI), which predicts the interface residues of a query protein with a specific interaction partner. We compared the performance of DockRank with several state‐of‐the‐art docking scoring functions using Success Rate (the percentage of cases that have at least one near‐native conformation among the top m conformations) and Hit Rate (the percentage of near‐native conformations that are included among the top m conformations). In cases where it is possible to obtain partner‐specific (PS) interface predictions from PS‐HomPPI, DockRank consistently outperforms both (i) ZRank and IRAD, two state‐of‐the‐art energy‐based scoring functions (improving Success Rate by up to 4‐fold); and (ii) Variants of DockRank that use predicted interface residues obtained from several protein interface predictors that do not take into account the binding partner in making interface predictions (improving success rate by up to 39‐fold). The latter result underscores the importance of using partner‐specific interface residues in scoring docked conformations. We show that DockRank, when used to re‐rank the conformations returned by ClusPro, improves upon the original ClusPro rankings in terms of both Success Rate and Hit Rate. DockRank is available as a server at http://einstein.cs.iastate.edu/DockRank/ . Proteins 2014; 82:250–267. © 2013 Wiley Periodicals, Inc.     

Linkers in the structural biology of protein–protein interactions

Linkers or spacers are short amino acid sequences created in nature to separate multiple domains in a single protein. Most of them are rigid and function to prohibit unwanted interactions between the discrete domains. However, Gly‐rich linkers are flexible, connecting various domains in a single protein without interfering with the function of each domain. The advent of recombinant DNA technology made it possible to fuse two interacting partners with the introduction of artificial linkers. Often, independent proteins may not exist as stable or structured proteins until they interact with their binding partner, following which they gain stability and the essential structural elements. Gly‐rich linkers have been proven useful for these types of unstable interactions, particularly where the interaction is weak and transient, by creating a covalent link between the proteins to form a stable protein–protein complex. Gly‐rich linkers are also employed to form stable covalently linked dimers, and to connect two independent domains that create a ligand‐binding site or recognition sequence. The lengths of linkers vary from 2 to 31 amino acids, optimized for each condition so that the linker does not impose any constraints on the conformation or interactions of the linked partners. Various structures of covalently linked protein complexes have been described using X‐ray crystallography, nuclear magnetic resonance and cryo‐electron microscopy techniques. In this review, we evaluate several structural studies where linkers have been used to improve protein quality, to produce stable protein–protein complexes, and to obtain protein dimers.     

Electron‐deficient p‐benzoyl‐l‐phenylalanine derivatives increase covalent chemical capture yields for protein–protein interactions

The photoactivatable amino acid p‐benzoyl‐l ‐phenylalanine (pBpa) has been used for the covalent capture of protein–protein interactions (PPIs) in vitro and in living cells. However, this technique often suffers from poor photocrosslinking yields due to the low reactivity of the active species. Here we demonstrate that the incorporation of halogenated pBpa analogs into proteins leads to increased crosslinking yields for protein–protein interactions. The analogs can be incorporated into live yeast and upon irradiation capture endogenous PPIs. Halogenated pBpas will extend the scope of PPIs that can be captured and expand the toolbox for mapping PPIs in their native environment.     

Computational design of second‐site suppressor mutations at protein–protein interfaces

The importance of a protein–protein interaction to a signaling pathway can be established by showing that amino acid mutations that weaken the interaction disrupt signaling, and that additional mutations that rescue the interaction recover signaling. Identifying rescue mutations, often referred to as second‐site suppressor mutations, controls against scenarios in which the initial deleterious mutation inactivates the protein or disrupts alternative protein–protein interactions. Here, we test a structure‐based protocol for identifying second‐site suppressor mutations that is based on a strategy previously described by Kortemme and Baker. The molecular modeling software Rosetta is used to scan an interface for point mutations that are predicted to weaken binding but can be rescued by mutations on the partner protein. The protocol typically identifies three types of specificity switches: knob‐in‐to‐hole redesigns, switching hydrophobic interactions to hydrogen bond interactions, and replacing polar interactions with nonpolar interactions. Computational predictions were tested with two separate protein complexes; the G‐protein Gα_i1 bound to the RGS14 GoLoco motif, and UbcH7 bound to the ubiquitin ligase E6AP. Eight designs were experimentally tested. Swapping a buried hydrophobic residue with a polar residue dramatically weakened binding affinities. In none of these cases were we able to identify compensating mutations that returned binding to wild‐type affinity, highlighting the challenges inherent in designing buried hydrogen bond networks. The strongest specificity switches were a knob‐in‐to‐hole design (20‐fold) and the replacement of a charge–charge interaction with nonpolar interactions (55‐fold). In two cases, specificity was further tuned by including mutations distant from the initial design. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Predicting permanent and transient protein–protein interfaces

Protein–protein interactions (PPIs) are involved in diverse functions in a cell. To optimize functional roles of interactions, proteins interact with a spectrum of binding affinities. Interactions are conventionally classified into permanent and transient, where the former denotes tight binding between proteins that result in strong complexes, whereas the latter compose of relatively weak interactions that can dissociate after binding to regulate functional activity at specific time point. Knowing the type of interactions has significant implications for understanding the nature and function of PPIs. In this study, we constructed amino acid substitution models that capture mutation patterns at permanent and transient type of protein interfaces, which were found to be different with statistical significance. Using the substitution models, we developed a novel computational method that predicts permanent and transient protein binding interfaces (PBIs) in protein surfaces. Without knowledge of the interacting partner, the method uses a single query protein structure and a multiple sequence alignment of the sequence family. Using a large dataset of permanent and transient proteins, we show that our method, BindML+, performs very well in protein interface classification. A very high area under the curve (AUC) value of 0.957 was observed when predicted protein binding sites were classified. Remarkably, near prefect accuracy was achieved with an AUC of 0.991 when actual binding sites were classified. The developed method will be also useful for protein design of permanent and transient PBIs. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Examining post‐translational modification‐mediated protein–protein interactions using a chemical proteomics approach

Post‐translational modifications (PTM) of proteins can control complex and dynamic cellular processes via regulating interactions between key proteins. To understand these regulatory mechanisms, it is critical that we can profile the PTM‐dependent protein–protein interactions. However, identifying these interactions can be very difficult using available approaches, as PTMs can be dynamic and often mediate relatively weak protein–protein interactions. We have recently developed CLASPI (cross‐linking‐assisted and stable isotope labeling in cell culture‐based protein identification), a chemical proteomics approach to examine protein–protein interactions mediated by methylation in human cell lysates. Here, we report three extensions of the CLASPI approach. First, we show that CLASPI can be used to analyze methylation‐dependent protein–protein interactions in lysates of fission yeast, a genetically tractable model organism. For these studies, we examined trimethylated histone H3 lysine‐9 (H3K9Me₃)‐dependent protein–protein interactions. Second, we demonstrate that CLASPI can be used to examine phosphorylation‐dependent protein–protein interactions. In particular, we profile proteins recognizing phosphorylated histone H3 threonine‐3 (H3T3‐Phos), a mitotic histone “mark” appearing exclusively during cell division. Our approach identified survivin, the only known H3T3‐Phos‐binding protein, as well as other proteins, such as MCAK and KIF2A, that are likely to be involved in weak but selective interactions with this histone phosphorylation “mark”. Finally, we demonstrate that the CLASPI approach can be used to study the interplay between histone H3T3‐Phos and trimethylation on the adjacent residue lysine 4 (H3K4Me₃). Together, our findings indicate the CLASPI approach can be broadly applied to profile protein–protein interactions mediated by PTMs.     

Towards the prediction of protein interaction partners using physical docking

Deciphering the whole network of protein interactions for a given proteome (‘interactome’) is the goal of many experimental and computational efforts in Systems Biology. Separately the prediction of the structure of protein complexes by docking methods is a well‐established scientific area. To date, docking programs have not been used to predict interaction partners. We provide a proof of principle for such an approach. Using a set of protein complexes representing known interactors in their unbound form, we show that a standard docking program can distinguish the true interactors from a background of 922 non‐redundant potential interactors. We additionally show that true interactions can be distinguished from non‐likely interacting proteins within the same structural family. Our approach may be put in the context of the proposed ‘funnel‐energy model’; the docking algorithm may not find the native complex, but it distinguishes binding partners because of the higher probability of favourable models compared with a collection of non‐binders. The potential exists to develop this proof of principle into new approaches for predicting interaction partners and reconstructing biological networks.     

Importance of the interaction protein–protein of the CaM–PDE1A and CaM–MLCK complexes in the development of new anti‐CaM drugs

Protein–protein interactions play central roles in physiological and pathological processes. The bases of the mechanisms of drug action are relevant to the discovery of new therapeutic targets. This work focuses on understanding the interactions in protein–protein–ligands complexes, using proteins calmodulin (CaM), human calcium/calmodulin‐dependent 3′,5′‐cyclic nucleotide phosphodiesterase 1A active human (PDE1A), and myosin light chain kinase (MLCK) and ligands αII–spectrin peptide (αII–spec), and two inhibitors of CaM (chlorpromazine (CPZ) and malbrancheamide (MBC)). The interaction was monitored with a fluorescent biosensor of CaM (hCaM M124C–mBBr). The results showed changes in the affinity of CPZ and MBC depending on the CaM–protein complex under analysis. For the Ca²⁺–CaM, Ca²⁺–CaM–PDE1A, and Ca²⁺–CaM–MLCK complexes, CPZ apparent dissociation constants (K_ds) were 1.11, 0.28, and 0.55 μM, respectively; and for MBC K_ds were 1.43, 1.10, and 0.61 μM, respectively. In competition experiments the addition of calmodulin binding peptide 1 (αII–spec) to Ca²⁺–hCaM M124C–mBBr quenched the fluorescence (K_d = 2.55 ± 1.75 pM) and the later addition of MBC (up to 16 μM) did not affect the fluorescent signal. Instead, the additions of αII–spec to a preformed Ca²⁺–hCaM M124C–mBBr–MBC complex modified the fluorescent signal. However, MBC was able to displace the PDE1A and MLCK from its complex with Ca²⁺–CaM. In addition, docking studies were performed for all complexes with both ligands showing an excellent correlation with experimental data. These experiments may help to explain why in vivo many CaM drugs target prefer only a subset of the Ca²⁺–CaM regulated proteins and adds to the understanding of molecular interactions between protein complexes and small ligands. Copyright © 2013 John Wiley & Sons, Ltd.     

Monte Carlo replica‐exchange based ensemble docking of protein conformations

A replica‐exchange Monte Carlo (REMC) ensemble docking approach has been developed that allows efficient exploration of protein–protein docking geometries. In addition to Monte Carlo steps in translation and orientation of binding partners, possible conformational changes upon binding are included based on Monte Carlo selection of protein conformations stored as ordered pregenerated conformational ensembles. The conformational ensembles of each binding partner protein were generated by three different approaches starting from the unbound partner protein structure with a range spanning a root mean square deviation of 1–2.5 Å with respect to the unbound structure. Because MC sampling is performed to select appropriate partner conformations on the fly the approach is not limited by the number of conformations in the ensemble compared to ensemble docking of each conformer pair in ensemble cross docking. Although only a fraction of generated conformers was in closer agreement with the bound structure the REMC ensemble docking approach achieved improved docking results compared to REMC docking with only the unbound partner structures or using docking energy minimization methods. The approach has significant potential for further improvement in combination with more realistic structural ensembles and better docking scoring functions. Proteins 2017; 85:924–937. © 2016 Wiley Periodicals, Inc.     

LuTHy: a double‐readout bioluminescence‐based two‐hybrid technology for quantitative mapping of protein–protein interactions in mammalian cells

Information on protein–protein interactions (PPIs) is of critical importance for studying complex biological systems and developing therapeutic strategies. Here, we present a double‐readout bioluminescence‐based two‐hybrid technology, termed LuTHy, which provides two quantitative scores in one experimental procedure when testing binary interactions. PPIs are first monitored in cells by quantification of bioluminescence resonance energy transfer (BRET) and, following cell lysis, are again quantitatively assessed by luminescence‐based co‐precipitation (LuC). The double‐readout procedure detects interactions with higher sensitivity than traditional single‐readout methods and is broadly applicable, for example, for detecting the effects of small molecules or disease‐causing mutations on PPIs. Applying LuTHy in a focused screen, we identified 42 interactions for the presynaptic chaperone CSPα, causative to adult‐onset neuronal ceroid lipofuscinosis (ANCL), a progressive neurodegenerative disease. Nearly 50% of PPIs were found to be affected when studying the effect of the disease‐causing missense mutations L115R and ?L116 in CSPα with LuTHy. Our study presents a robust, sensitive research tool with high utility for investigating the molecular mechanisms by which disease‐associated mutations impair protein activity in biological systems.     

A combinatorial scoring function for protein–RNA docking

Protein–RNA docking is still an open question. One of the main challenges is to develop an effective scoring function that can discriminate near‐native structures from the incorrect ones. To solve the problem, we have constructed a knowledge‐based residue‐nucleotide pairwise potential with secondary structure information considered for nonribosomal protein–RNA docking. Here we developed a weighted combined scoring function RpveScore that consists of the pairwise potential and six physics‐based energy terms. The weights were optimized using the multiple linear regression method by fitting the scoring function to L_rmsd for the bound docking decoys from Benchmark II. The scoring functions were tested on 35 unbound docking cases. The results show that the scoring function RpveScore including all terms performs best. Also RpveScore was compared with the statistical mechanics‐based method derived potential ITScore‐PR, and the united atom‐based statistical potentials QUASI‐RNP and DARS‐RNP. The success rate of RpveScore is 71.6% for the top 1000 structures and the number of cases where a near‐native structure is ranked in top 30 is 25 out of 35 cases. For 32 systems (91.4%), RpveScore can find the binding mode in top 5 that has no lower than 50% native interface residues on protein and nucleotides on RNA. Additionally, it was found that the long‐range electrostatic attractive energy plays an important role in distinguishing near‐native structures from the incorrect ones. This work can be helpful for the development of protein–RNA docking methods and for the understanding of protein–RNA interactions. RpveScore program is available to the public at http://life.bjut.edu.cn/kxyj/kycg/2017116/14845362285362368_1.html Proteins 2017; 85:741–752. © 2016 Wiley Periodicals, Inc.     

Feature selection and classification of protein–protein complexes based on their binding affinities using machine learning approaches

Protein–protein interactions are intrinsic to virtually every cellular process. Predicting the binding affinity of protein–protein complexes is one of the challenging problems in computational and molecular biology. In this work, we related sequence features of protein–protein complexes with their binding affinities using machine learning approaches. We set up a database of 185 protein–protein complexes for which the interacting pairs are heterodimers and their experimental binding affinities are available. On the other hand, we have developed a set of 610 features from the sequences of protein complexes and utilized Ranker search method, which is the combination of Attribute evaluator and Ranker method for selecting specific features. We have analyzed several machine learning algorithms to discriminate protein‐protein complexes into high and low affinity groups based on their K_d values. Our results showed a 10‐fold cross‐validation accuracy of 76.1% with the combination of nine features using support vector machines. Further, we observed accuracy of 83.3% on an independent test set of 30 complexes. We suggest that our method would serve as an effective tool for identifying the interacting partners in protein–protein interaction networks and human–pathogen interactions based on the strength of interactions. Proteins 2014; 82:2088–2096. © 2014 Wiley Periodicals, Inc.