Protein binding hot spots and the residue-residue pairing preference: a water exclusion perspective

Crystal structures of glycoprotein D of equine alphaherpesviruses reveal potential binding sites to the entry receptor MHC-I

Cell entry of most alphaherpesviruses is mediated by the binding of glycoprotein D (gD) to different cell surface receptors. Equine herpesvirus type 1 (EHV-1) and EHV-4 gDs interact with equine major histocompatibility complex I (MHC-I) to initiate entry into equine cells. We have characterized the gD-MHC-I interaction by solving the crystal structures of EHV-1 and EHV-4 gDs (gD1, gD4), performing protein-protein docking simulations, surface plasmon resonance (SPR) analysis, and biological assays. The structures of gD1 and gD4 revealed the existence of a common V-set immunoglobulin-like (IgV-like) core comparable to those of other gD homologs. Molecular modeling yielded plausible binding hypotheses and identified key residues (F213 and D261) that are important for virus binding. Altering the key residues resulted in impaired virus growth in cells, which highlights the important role of these residues in the gD-MHC-I interaction. Taken together, our results add to our understanding of the initial herpesvirus-cell interactions and will contribute to the targeted design of antiviral drugs and vaccine development.

Cold spots are universal in protein–protein interactions

Proteins interact with each other through binding interfaces that differ greatly in size and physico‐chemical properties. Within the binding interface, a few residues called hot spots contribute the majority of the binding free energy and are hence irreplaceable. In contrast, cold spots are occupied by suboptimal amino acids, providing possibility for affinity enhancement through mutations. In this study, we identify cold spots due to cavities and unfavorable charge interactions in multiple protein–protein interactions (PPIs). For our cold spot analysis, we first use a small affinity database of PPIs with known structures and affinities and then expand our search to nearly 4000 homo‐ and heterodimers in the Protein Data Bank (PDB). We observe that cold spots due to cavities are present in nearly all PPIs unrelated to their binding affinity, while unfavorable charge interactions are relatively rare. We also find that most cold spots are located in the periphery of the binding interface, with high‐affinity complexes showing fewer centrally located colds spots than low‐affinity complexes. A larger number of cold spots is also found in non‐cognate interactions compared to their cognate counterparts. Furthermore, our analysis reveals that cold spots are more frequent in homo‐dimeric complexes compared to hetero‐complexes, likely due to symmetry constraints imposed on sequences of homodimers. Finally, we find that glycines, glutamates, and arginines are the most frequent amino acids appearing at cold spot positions. Our analysis emphasizes the importance of cold spot positions to protein evolution and facilitates protein engineering studies directed at enhancing binding affinity and specificity in a wide range of applications.

Getting Insights into Structural and Energetic Properties of Reciprocal Peptide-Protein Interactions

Peptide-protein interactions play a key role for many cellular and metabolic processes involved in the onset of largely spread diseases such as cancer and neurodegenerative pathologies. Despite the progress in the structural characterization of peptide-protein interfaces, the in-depth knowledge of the molecular details behind their interactions is still a daunting task. Here, we present the first comprehensive in silico morphological and energetic study of peptide binding sites by focusing on both peptide and protein standpoints. Starting from the PixelDB database, a nonredundant benchmark collection of high-quality 3D crystallographic structures of peptide-protein complexes, a classification analysis of the most representative categories based on the nature of each cocrystallized peptide has been carried out. Several interpretable geometrical and energetic descriptors have been computed both from peptide and target protein sides in the attempt to unveil physicochemical and structural causative correlations. Finally, we investigated the most frequent peptide-protein residue pairs at the binding interface and made extensive energetic analyses, based on GRID MIFs, with the aim to study the peptide affinity-enhancing interactions to be further exploited in rational drug design strategies.

Positive-unlabeled learning for the prediction of conformational B-cell epitopes

Background

The incomplete ground truth of training data of B-cell epitopes is a demanding issue in computational epitope prediction. The challenge is that only a small fraction of the surface residues of an antigen are confirmed as antigenic residues (positive training data); the remaining residues are unlabeled. As some of these uncertain residues can possibly be grouped to form novel but currently unknown epitopes, it is misguided to unanimously classify all the unlabeled residues as negative training data following the traditional supervised learning scheme.

Results

We propose a positive-unlabeled learning algorithm to address this problem. The key idea is to distinguish between epitope-likely residues and reliable negative residues in unlabeled data. The method has two steps: (1) identify reliable negative residues using a weighted SVM with a high recall; and (2) construct a classification model on the positive residues and the reliable negative residues. Complex-based 10-fold cross-validation was conducted to show that this method outperforms those commonly used predictors DiscoTope 2.0, ElliPro and SEPPA 2.0 in every aspect. We conducted four case studies, in which the approach was tested on antigens of West Nile virus, dihydrofolate reductase, beta-lactamase, and two Ebola antigens whose epitopes are currently unknown. All the results were assessed on a newly-established data set of antigen structures not bound by antibodies, instead of on antibody-bound antigen structures. These bound structures may contain unfair binding information such as bound-state B-factors and protrusion index which could exaggerate the epitope prediction performance. Source codes are available on request.

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.

Protein painting reveals solvent-excluded drug targets hidden within native protein-protein interfaces

Identifying the contact regions between a protein and its binding partners is essential for creating therapies that block the interaction. Unfortunately, such contact regions are extremely difficult to characterize because they are hidden inside the binding interface. Here we introduce protein painting as a new tool that employs small molecules as molecular paints to tightly coat the surface of protein–protein complexes. The molecular paints, which block trypsin cleavage sites, are excluded from the binding interface. Following mass spectrometry, only peptides hidden in the interface emerge as positive hits, revealing the functional contact regions that are drug targets. We use protein painting to discover contact regions between the three-way interaction of IL1β ligand, the receptor IL1RI and the accessory protein IL1RAcP. We then use this information to create peptides and monoclonal antibodies that block the interaction and abolish IL1β cell signalling. The technology is broadly applicable to discover protein interaction drug targets.

Identifying the site where a protein binds to another molecule is an important factor for the design of therapeutics intended to prevent this interaction. Here, the authors coat protein–receptor complexes with surface-binding molecules, and determine their interacting regions using mass spectrometry.

Hot spot-based design of small-molecule inhibitors for protein-protein interactions

Protein-protein interactions (PPIs) are important targets for the development of chemical probes and therapeutic agents. From the initial discovery of the existence of hot spots at PPI interfaces, it has been proposed that hot spots might provide the key for developing small-molecule PPI inhibitors. However, there has been no review on the ways in which the knowledge of hot spots can be used to achieve inhibitor design, nor critical examination of successful examples. This Digest discusses the characteristics of hot spots and the identification of druggable hot spot pockets. An analysis of four examples of hot spot-based design reveals the importance of this strategy in discovering potent and selective PPI inhibitors. A general procedure for hot spot-based design of PPI inhibitors is outlined.

Solvent-accessible surface area: How well can be applied to hot-spot detection?

A detailed comprehension of protein-based interfaces is essential for the rational drug development. One of the key features of these interfaces is their solvent accessible surface area profile. With that in mind, we tested a group of 12 SASA-based features for their ability to correlate and differentiate hot- and null-spots. These were tested in three different data sets, explicit water MD, implicit water MD, and static PDB structure. We found no discernible improvement with the use of more comprehensive data sets obtained from molecular dynamics. The features tested were shown to be capable of discerning between hot- and null-spots, while presenting low correlations. Residue standardization such as rel SASAi or rel/res SASAi , improved the features as a tool to predict ΔΔGbinding values. A new method using support machine learning algorithms was developed: SBHD (Sasa-Based Hot-spot Detection). This method presents a precision, recall, and F1 score of 0.72, 0.81, and 0.76 for the training set and 0.91, 0.73, and 0.81 for an independent test set.

Computational Alanine Scanning Mutagenesis: MM-PBSA vs TI

Understanding protein-protein association and being able to determine the crucial residues responsible for their association (hot-spots) is a key issue with huge practical applications such as rational drug design and protein engineering. A variety of computational methods exist to detect hot-spots residues, but the development of a fast and accurate quantitative alanine scanning mutagenesis (ASM) continues to be crucial. Using four protein-protein complexes, we have compared a variation of the standard computational ASM protocol developed at our group, based on the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM-PBSA) approach, against Thermodynamic Integration (TI), a well-known and accurate but computationally expensive method. To compare the efficiency and the accuracy of the two methods, we have calculated the protein-protein binding free energy differences upon alanine mutation of interfacial residues (ΔΔGbind). In relation to the experimental ΔΔGbind values, the average error obtained with TI was 1.53 kcal/mol, while the ASM protocol resulted in an average error of 1.18 kcal/mol. The results demonstrate that the much faster ASM protocol gives results at the same level of accuracy as the TI method but at a fraction of the computational time required to run TI. This ASM protocol is therefore a strong and efficient alternative to the systematic evaluation of protein-protein interfaces, involving hundreds of amino acid residues in search of hot-spots.

A simple recipe for the non-expert bioinformaticist for building experimentally-testable hypotheses for proteins with no known homologs

The study of the protein-protein interactions (PPIs) of unique ORFs is a strategy for deciphering the biological roles of unique ORFs of interest. For uniform reference, we define unique ORFs as those for which no matching protein is found after PDB-BLAST search with default parameters. The uniqueness of the ORFs generally precludes the straightforward use of structure-based approaches in the design of experiments to explore PPIs. Many open-source bioinformatics tools, from the commonly-used to the relatively esoteric, have been built and validated to perform analyses and/or predictions of sorts on proteins. How can these available tools be combined into a protocol that helps the non-expert bioinformaticist researcher to design experiments to explore the PPIs of their unique ORF? Here we define a pragmatic protocol based on accessibility of software to achieve this and we make it concrete by applying it on two proteins-the ImuB and ImuA' proteins from Mycobacterium tuberculosis. The protocol is pragmatic in that decisions are made largely based on the availability of easy-to-use freeware. We define the following basic and user-friendly software pathway to build testable PPI hypotheses for a query protein sequence: PSI-PRED → MUSTER → metaPPISP → ASAView and ConSurf. Where possible, other analytical and/or predictive tools may be included. Our protocol combines the software predictions and analyses with general bioinformatics principles to arrive at consensus, prioritised and testable PPI hypotheses.

Z-score biological significance of binding hot spots of protein interfaces by using crystal packing as the reference state

Characterization of binding hot spots of protein interfaces is a fundamental study in molecular biology. Many computational methods have been proposed to identify binding hot spots. However, there are few studies to assess the biological significance of binding hot spots. We introduce the notion of biological significance of a contact residue for capturing the probability of the residue occurring in or contributing to protein binding interfaces. We take a statistical Z-score approach to the assessment of the biological significance. The method has three main steps. First, the potential score of a residue is defined by using a knowledge-based potential function with relative accessible surface area calculations. A null distribution of this potential score is then generated from artifact crystal packing contacts. Finally, the Z-score significance of a contact residue with a specific potential score is determined according to this null distribution. We hypothesize that residues at binding hot spots have big absolute values of Z-score as they contribute greatly to binding free energy. Thus, we propose to use Z-score to predict whether a contact residue is a hot spot residue. Comparison with previously reported methods on two benchmark datasets shows that this Z-score method is mostly superior to earlier methods. This article is part of a Special Issue entitled: Computational Methods for Protein Interaction and Structural Prediction.

Protein-protein interactions: principles, techniques, and their potential role in new drug development

A vast network of genes is inter-linked through protein-protein interactions and is critical component of almost every biological process under physiological conditions. Any disruption of the biologically essential network leads to pathological conditions resulting into related diseases. Therefore, proper understanding of biological functions warrants a comprehensive knowledge of protein-protein interactions and the molecular mechanisms that govern such processes. The importance of protein-protein interaction process is highlighted by the fact that a number of powerful techniques/methods have been developed to understand how such interactions take place under various physiological and pathological conditions. Many of the key protein-protein interactions are known to participate in disease-associated signaling pathways, and represent novel targets for therapeutic intervention. Thus, controlling protein-protein interactions offers a rich dividend for the discovery of new drug targets. Availability of various tools to study and the knowledge of human genome have put us in a unique position to understand highly complex biological network, and the mechanisms involved therein. In this review article, we have summarized protein-protein interaction networks, techniques/methods of their binding/kinetic parameters, and the role of these interactions in the development of potential tools for drug designing.

Architecture of a full-length retroviral integrase monomer and dimer, revealed by small angle X-ray scattering and chemical cross-linking

We determined the size and shape of full-length avian sarcoma virus (ASV) integrase (IN) monomers and dimers in solution using small angle x-ray scattering. The low resolution data obtained establish constraints for the relative arrangements of the three component domains in both forms. Domain organization within the small angle x-ray envelopes was determined by combining available atomic resolution data for individual domains with results from cross-linking coupled with mass spectrometry. The full-length dimer architecture so revealed is unequivocally different from that proposed from x-ray crystallographic analyses of two-domain fragments, in which interactions between the catalytic core domains play a prominent role. Core-core interactions are detected only in cross-linked IN tetramers and are required for concerted integration. The solution dimer is stabilized by C-terminal domain (CTD-CTD) interactions and by interactions of the N-terminal domain in one subunit with the core and CTD in the second subunit. These results suggest a pathway for formation of functional IN-DNA complexes that has not previously been considered and possible strategies for preventing such assembly.