Homology inference of protein-protein interactions via conserved binding sites

The coverage and reliability of protein-protein interactions determined by high-throughput experiments still needs to be improved, especially for higher organisms, therefore the question persists, how interactions can be verified and predicted by computational approaches using available data on protein structural complexes. Recently we developed an approach called IBIS (Inferred Biomolecular Interaction Server) to predict and annotate protein-protein binding sites and interaction partners, which is based on the assumption that the structural location and sequence patterns of protein-protein binding sites are conserved between close homologs. In this study first we confirmed high accuracy of our method and found that its accuracy depends critically on the usage of all available data on structures of homologous complexes, compared to the approaches where only a non-redundant set of complexes is employed. Second we showed that there exists a trade-off between specificity and sensitivity if we employ in the prediction only evolutionarily conserved binding site clusters or clusters supported by only one observation (singletons). Finally we addressed the question of identifying the biologically relevant interactions using the homology inference approach and demonstrated that a large majority of crystal packing interactions can be correctly identified and filtered by our algorithm. At the same time, about half of biological interfaces that are not present in the protein crystallographic asymmetric unit can be reconstructed by IBIS from homologous complexes without the prior knowledge of crystal parameters of the query protein.     

Modulating protein-protein interactions with small molecules: the importance of binding hotspots

The modulation of protein-protein interactions (PPIs) by small drug-like molecules is a relatively new area of research and has opened up new opportunities in drug discovery. However, the progress made in this area is limited to a handful of known cases of small molecules that target specific diseases. With the increasing availability of protein structure complexes, it is highly important to devise strategies exploiting homologous structure space on a large scale for discovering putative PPIs that could be attractive drug targets. Here, we propose a scheme that allows performing large-scale screening of all protein complexes and finding putative small-molecule and/or peptide binding sites overlapping with protein-protein binding sites (so-called "multibinding sites"). We find more than 600 nonredundant proteins from 60 protein families with multibinding sites. Moreover, we show that the multibinding sites are mostly observed in transient complexes, largely overlap with the binding hotspots and are more evolutionarily conserved than other interface sites. We investigate possible mechanisms of how small molecules may modulate protein-protein binding and discuss examples of new candidates for drug design.     

There are more small amino acids and fewer aromatic rings in proteins of ionizing radiation-resistant bacteria

The identification of specific amino acids (AAs) or groups of functionally important AA residues in ionizing radiation-resistant bacteria (IRRB) is an important challenge in understanding the biological basis of resistance to ionizing radiation (IR; X-rays and gamma-rays). To address this problem, we compared homologous sites in multiple alignments of proteins of IRRB and IR-sensitive bacteria (IRSB) using the DeltaProt Toolbox. Substitution patterns were used as evidence for selection of certain AAs over others. Our results show that, in contrast to aromatic AAs, small/tiny AAs tend to be preferred in IRRB compared to IRSB. In agreement with previous experimental data showing that oxidation of AA residues is causative in the killing of irradiated cells and that IR resistance is correlated with the accumulation of divalent manganese ions (Mn²⁺)–peptide–orthophosphate (Pi) complexes, we proposed a chemical interpretation based on the Hard and Soft (Lewis) Acids and Bases (HSAB) concept. These findings should assist future efforts in selecting mutations for rational design of proteins with enhanced IR tolerance properties.     

MAGIIC-PRO: detecting functional signatures by efficient discovery of long patterns in protein sequences

This paper presents a web service named MAGIICPRO,which aims to discover functional signatures of a query protein by sequential pattern mining. Automatic discovery of patterns from unaligned biological sequences is an important problem in molecular biology. MAGIIC-PRO is different from several previously established methods performing similar tasks in two major ways. The first remarkable feature of MAGIIC-PRO is its efficiency in delivering long patterns. With incorporating a new type of gap constraints and some of the state-of-theart data mining techniques, MAGIIC-PRO usually identifies satisfied patterns within an acceptable response time. The efficiency of MAGIIC-PRO enables the users to quickly discover functional signatures of which the residues are not from only one region of the protein sequences or are only conserved in few members of a protein family. The second remarkable feature of MAGIIC-PRO is its effort in refining the mining results. Considering large flexible gaps improves the completeness of the derived functional signatures. The users can be directly guided to the patterns with as many blocks as that are conserved simultaneously. In this paper,we show by experiments that MAGIIC-PRO is efficient and effective in identifying ligand-binding sites and hot regions in protein-protein interactions directly from sequences. The web service is availableat http://biominer.bime.ntu.edu.tw/magiicproand a mirror site at http://biominer.cse.yzu.edu.tw/magiicpro.     

Molecular principles of the interactions of disordered proteins

Thorough knowledge of the molecular principles of protein-protein recognition is essential to our understanding of protein function at the cellular level. Whereas interactions of ordered proteins have been analyzed in great detail, complexes of intrinsically unstructured/disordered proteins (IUPs) have hardly been addressed so far. Here, we have collected a database of 39 complexes of experimentally verified IUPs, and compared their interfaces with those of 72 complexes of ordered, globular proteins. The characteristic differences found between the two types of complexes suggest that IUPs represent a distinct molecular implementation of the principles of protein-protein recognition. The interfaces do not differ in size, but those of IUPs cover a much larger part of the surface of the protein than for their ordered counterparts. Moreover, IUP interfaces are significantly more hydrophobic relative to their overall amino acid composition, but also in absolute terms. They rely more on hydrophobic-hydrophobic than on polar-polar interactions. Their amino acids in the interface realize more intermolecular contacts, which suggests a better fit with the partner due to induced folding upon binding that results in a better adaptation to the partner. The two modes of interaction also differ in that IUPs usually use only a single continuous segment for partner binding, whereas the binding sites of ordered proteins are more segmented. Probably, all these features contribute to the increased evolutionary conservation of IUP interface residues. These noted molecular differences are also manifested in the interaction energies of IUPs. Our approximation of these by low-resolution force-fields shows that IUPs gain much more stabilization energy from intermolecular contacts, than from folding, i.e. they use their binding energy for folding. Overall, our findings provide a structural rationale to the prior suggestions that many IUPs are specialized for functions realized by protein-protein interactions.     

A novel method for protein-protein interaction site prediction using phylogenetic substitution models

Protein-protein binding events mediate many critical biological functions in the cell. Typically, functionally important sites in proteins can be well identified by considering sequence conservation. However, protein-protein interaction sites exhibit higher sequence variation than other functional regions, such as catalytic sites of enzymes. Consequently, the mutational behavior leading to weak sequence conservation poses significant challenges to the protein-protein interaction site prediction. Here, we present a phylogenetic framework to capture critical sequence variations that favor the selection of residues essential for protein-protein binding. Through the comprehensive analysis of diverse protein families, we show that protein binding interfaces exhibit distinct amino acid substitution as compared with other surface residues. On the basis of this analysis, we have developed a novel method, BindML, which utilizes the substitution models to predict protein-protein binding sites of protein with unknown interacting partners. BindML estimates the likelihood that a phylogenetic tree of a local surface region in a query protein structure follows the substitution patterns of protein binding interface and nonbinding surfaces. BindML is shown to perform well compared to alternative methods for protein binding interface prediction. The methodology developed in this study is very versatile in the sense that it can be generally applied for predicting other types of functional sites, such as DNA, RNA, and membrane binding sites in proteins.     

Rapid evolution exposes the boundaries of domain structure and function in natively unfolded FG nucleoporins

Nucleoporins with phenylalanine-glycine repeats (FG Nups) function at the nuclear pore complex (NPC) to facilitate nucleocytoplasmic transport. In Saccharomyces cerevisiae, each FG Nup contains a large natively unfolded domain that is punctuated by FG repeats. These FG repeats are surrounded by hydrophilic amino acids (AAs) common to disordered protein domains. Here we show that the FG domain of Nups from human, fly, worm, and other yeast species is also enriched in these disorder-associated AAs, indicating that structural disorder is a conserved feature of FG Nups and likely serves an important role in NPC function. Despite the conservation of AA composition, FG Nup sequences from different species show extensive divergence. A comparison of the AA substitution rates of proteins with syntenic orthologs in four Saccharomyces species revealed that FG Nups have evolved at twice the rate of average yeast proteins with most substitutions occurring in sequences between FG repeats. The rapid evolution of FG Nups is poorly explained by parameters known to influence AA substitution rate, such as protein expression level, interactivity, and essentiality; instead their rapid evolution may reflect an intrinsic permissiveness of natively unfolded structures to AA substitutions. The overall lack of AA sequence conservation in FG Nups is sharply contrasted by discrete stretches of conserved sequences. These conserved sequences highlight known karyopherin and nucleoporin binding sites as well as other uncharacterized sites that may have important structural and functional properties.     

Highly discriminating protein-protein interaction specificities in the context of a conserved binding energy hotspot

We explore the thermodynamic basis for high affinity binding and specificity in conserved protein complexes using colicin endonuclease-immunity protein complexes as our model system. We investigated the ability of each colicin-specific immunity protein (Im2, Im7, Im8 and Im9) to bind the endonuclease (DNase) domains of colicins E2, E7 and E8 in vitro and compared these to the previously studied colicin E9. We find that high affinity binding (Kd < or = 10(-14) M) is a common feature of cognate colicin DNase-Im protein complexes as are non-cognate protein-protein associations, which are generally 10(6)-10(8)-fold weaker. Comparative alanine scanning of Im2 and Im9 residues involved in binding the E2 DNase revealed similar behaviour to that of the two proteins binding the E9 DNase; helix III forms a conserved binding energy hotspot with specificity residues from helix II only contributing favourably in a cognate interaction, a combination we have termed as "dual recognition". Significant differences are seen, however, in the number and side-chain chemistries of specificity sites that contribute to cognate binding. In Im2, Asp33 from helix II dominates colicin E2 specificity, whereas in Im9 several hydrophobic residues, including position 33 (leucine), help define its colicin specificity. A similar distribution of specificity sites was seen using phage display where, with Im2 as the template, a library of randomised sequences was generated in helix II and the library panned against either the E2 or E9 DNase. Position 33 was the dominant specificity site recovered in all E2 DNase-selected clones, whereas a number of Im9 specificity sites were recovered in E9 DNase-selected clones, including position 33. In order to probe the relationship between biological specificity and in vitro binding affinity we compared the degree of protection afforded to bacteria against colicin E9 toxicity by a set of immunity proteins whose affinities for the E9 DNase differed by up to ten orders of magnitude. This analysis indicated that the Kd required for complete biological protection is <10(-10)M and that the "affinity window" over which the selection of novel immunity protein specificities likely evolves is 10(-6)-10(-10)M. This comprehensive survey of colicin DNase-immunity protein complexes illustrates how high affinity protein-protein interactions can be very discriminating even though binding is dominated by a conserved hotspot, with single or multiple specificity sites modulating the overall binding free energy. We discuss these results in the context of other conserved protein complexes and suggest that they point to a generic specificity mechanism in divergently evolved protein-protein interactions.     

Chromophore/protein and chromophore/anion interactions in halorhodopsin.

Halorhodopsin (HR), the light-driven chloride transport pigment of Halobacterium halobium, was bleached and reconstituted with retinal analogues with the pi electron system interrupted at different locations (dihydroretinals). The absorption maxima of the artificial pigments formed with the dihydroretinals are found to be very similar to those of the corresponding pigments formed by reconstitution of bacteriorhodopsin (BR) and sensory rhodopsin (SR). This strongly suggests that the distribution of charges around the retinal is similar in all three bacterial rhodopsins. Comparison of the primary, and proposed secondary, structures for HR and BR reveal conserved asparagine (asp) and arginine (arg) residues, which are likely candidates for the ionizable amino acids that interact with the retinal. In a second set of experiments absorption shifts due to the binding of anions to Sites I and II in HR, reconstituted with different retinal analogues, were used to estimate the locations of these binding sites relative to the retinal. Site I is localized near the Schiff base, and Site II near the ionone ring. On the basis of these results a structural model for HR is proposed, which accounts for the spectroscopic properties of HR in terms of the three buried arg residues and two of the buried asp residues in the protein.     

DNA sequence dependent and independent conformational changes in multipartite operator recognition by lambda-repressor

Binding of regulatory proteins to multipartite DNA binding sites often occurs with protein-protein interaction, resulting in cooperative binding. The operators of bacteriophage lambda have several pairs of repressor binding sites (O(R)1-O(R)2, O(R)2-O(R)3, O(L)1-O(L)2, and O(L)2-O(L)3) separated by a variable number of base pairs, and thus, bacteriophage lambda is a model system for studying multipartite operator recognition by DNA-binding proteins. Near-UV circular dichroism spectra show that the DNA is distorted in O(R)1-O(R)2 and O(L)2-O(L)3 but much less so in O(R)2-O(R)3. Upon titration of lambda-repressor with single-operator sites O(R)1, O(R)2, and O(R)3, it was observed that the tryptophan fluorescence quenches to different degrees, suggesting different conformations of the protein in the three DNA-protein complexes. Acrylamide quenching of tryptophan fluorescence of lambda-repressor bound to these single operators also shows different Stern-Volmer constants, supporting the above conclusions. Titration of lambda-repressor with oligonucleotides containing pairs of operator sites also causes different degrees of fluorescence quenching. In particular, fluorescence quenching induced by O(R)1-O(R)2 binding is less than the quenching induced by either of the single operators alone, suggesting additional conformational changes upon establishment of protein-protein contact. Stern-Volmer constants obtained from acrylamide quenching of tryptophan fluorescence of lambda-repressor bound cooperatively to pairs of operator sites are different from those of the single-operator-site-bound repressors. For example, O(R)2-O(R)3-bound repressor has significantly higher acrylamide quenchable components than either of the O(R)2- or O(R)3-bound proteins, again suggesting additional conformational changes upon establishment of protein-protein contact. We conclude that the strategy of recognition of multipartite operator by lambda-repressor is complex and varied, involving conformational changes in both DNA and protein that are determined by the separation of the binding sites as well as the nucleic acid sequence.     

Purification of RNA and RNA-protein complexes by an R17 coat protein affinity method.

We describe an affinity chromatography method to isolate specific RNAs and RNA-protein complexes formed in vivo or in vitro. It exploits the highly selective binding of the coat protein of bacteriophage R17 to a short hairpin in its genomic RNA. RNA containing that hairpin binds to coat protein that has been covalently bound to a solid support. Bound RNA-protein complexes can be eluted with excess R17 recognition sites. Using purified RNA, we demonstrate that binding to immobilized coat protein is highly specific and enables one to separate an RNA of interest from a large excess of other RNAs in a single step. Surprisingly, binding of an RNA containing non-R17 sequences to the support requires two recognition sites in tandem; a single site is insufficient. We determine optimal conditions for purification of specific RNAs by comparing specific binding (retention of RNAs with recognition sites) to non-specific binding (retention of RNAs without recognition sites) over a range of experimental conditions. These results suggest that binding of immobilized coat protein to RNAs containing two sites is cooperative. We illustrate the potential utility of the approach in purifying RNA-protein complexes by demonstrating that a U1 snRNP formed in vivo on an RNA containing tandem recognition sites is selectively retained by the coat protein support.     

In vivo assay for protein-protein interactions using Drosophila chromosomes

The ability of a chimeric HP1-Polycomb (Pc) protein to bind both to heterochromatin and to euchromatic sites of Pc protein binding was exploited to detect stable protein-protein interactions in vivo. Previously, we showed that endogenous Pc protein was recruited to ectopic heterochromatic binding sites by the chimeric protein. Here, we examine the association of other Pc group (Pc-G) proteins. We show that Posterior sex combs (Psc) protein also is recruited to heterochromatin by the chimeric protein, demonstrating that Psc protein participates in direct protein-protein interaction with Pc protein or Pc-associated protein. In flies carrying temperature-sensitive alleles of Enhancer of zeste[E(z)] the general decondensation of polytene chromosomes that occurs at the restrictive temperature is associated with loss of binding of endogenous Pc and chimeric HP1-Polycomb protein to euchromatin, but binding of HP1 and chimeric HP1-Polycomb protein to the heterochromatin is maintained. The E(z) mutation also results in the loss of chimera-dependent binding to heterochromatin by endogenous Pc and Psc proteins at the restrictive temperature, suggesting that interaction of these proteins is mediated by E(z) protein. A myc-tagged full-length Suppressor 2 of zeste [Su(z)2] protein interacts poorly or not at all with ectopic Pc-G complexes, but a truncated Su(z)2 protein is strongly recruited to all sites of chimeric protein binding. Trithorax protein is not recruited to the heterochromatin by the chimeric HP1-Polycomb protein, suggesting either that this protein does not interact directly with Pc-G complexes or that such interactions are regulated. Ectopic binding of chimeric chromosomal proteins provides a useful tool for distinguishing specific protein-protein interactions from specific protein-DNA interactions important for complex assembly in vivo.     

Computational design of a new hydrogen bond network and at least a 300-fold specificity switch at a protein-protein interface

The redesign of protein-protein interactions is a stringent test of our understanding of molecular recognition and specificity. Previously we engineered a modest specificity switch into the colicin E7 DNase-Im7 immunity protein complex by identifying mutations that are disruptive in the native complex, but can be compensated by mutations on the interacting partner. Here we extend the approach by systematically sampling alternate rigid body orientations to optimize the interactions in a binding mode specific manner. Using this protocol we designed a de novo hydrogen bond network at the DNase-immunity protein interface and confirmed the design with X-ray crystallographic analysis. Subsequent design of the second shell of interactions guided by insights from the crystal structure on tightly bound water molecules, conformational strain, and packing defects yielded new binding partners that exhibited specificities of at least 300-fold between the cognate and the non-cognate complexes. This multi-step approach should be applicable to the design of polar protein-protein interactions and contribute to the re-engineering of regulatory networks mediated by protein-protein interactions.     

A quantitative analysis of interfacial amino acid conservation in protein-protein hetero complexes

A long-standing question in molecular biology is whether interfaces of protein-protein complexes are more conserved than the rest of the protein surfaces. Although it has been reported that conservation can be used as an indicator for predicting interaction sites on proteins, there are recent reports stating that the interface regions are only slightly more conserved than the rest of the protein surfaces, with conservation signals not being statistically significant enough for predicting protein-protein binding sites. In order to properly address these controversial reports we have studied a set of 28 well resolved hetero complex structures of proteins that consists of transient and non-transient complexes. The surface positions were classified into four conservation classes and the conservation index of the surface positions was quantitatively analyzed. The results indicate that the surface density of highly conserved positions is significantly higher in the protein-protein interface regions compared with the other regions of the protein surface. However, the average conservation index of the patches in the interface region is not significantly higher compared with other surface regions of the protein structures. This finding demonstrates that the number of conserved residue positions is a more appropriate indicator for predicting protein-protein binding sites than the average conservation index in the interacting region. We have further validated our findings on a set of 59 benchmark complex structures. Furthermore, an analysis of 19 complexes of antigen-antibody interactions shows that there is no conservation of amino acid positions in the interacting regions of these complexes, as expected, with the variable region of the immunoglobulins interacting mostly with the antigens. Interestingly, antigen interacting regions also have a higher number of non-conserved residue positions in the interacting region than the rest of the protein surface.     

Hydrophilicity Matching – A Potential Prerequisite for the Formation of Protein-Protein Complexes in the Cell

A binding event between two proteins typically consists of a diffusional search of binding partners for one another, followed by a specific recognition of the compatible binding sites resulting in the formation of the complex. However, it is unclear how binding partners find each other in the context of the crowded, constantly fluctuating, and interaction-rich cellular environment. Here we examine the non-specific component of protein-protein interactions, which refers to those physicochemical properties of the binding partners that are independent of the exact details of their binding sites, but which can affect their localization or diffusional search for one another. We show that, for a large set of high-resolution experimental 3D structures of binary, transient protein complexes taken from the DOCKGROUND database, the binding partners display a surprising, statistically significant similarity in terms of their total hydration free energies normalized by a size-dependent variable. We hypothesize that colocalization of binding partners, even within individual cellular compartments such as the cytoplasm, may be influenced by their relative hydrophilicity, potentially in response to local hydrophilic gradients.     

Deep learning-based method for predicting and classifying the binding affinity of protein-protein complexes

Protein-protein interactions (PPIs) play a critical role in various biological processes. Accurately estimating the binding affinity of PPIs is essential for understanding the underlying molecular recognition mechanisms. In this study, we employed a deep learning approach to predict the binding affinity (ΔG) of protein-protein complexes. To this end, we compiled a dataset of 903 protein-protein complexes, each with its corresponding experimental binding affinity, which belong to six functional classes. We extracted 8 to 20 non-redundant features from the sequence information as well as the predicted three-dimensional structures using feature selection methods for each protein functional class. Our method showed an overall mean absolute error of 1.05 kcal/mol and a correlation of 0.79 between experimental and predicted ΔG values. Additionally, we evaluated our model for discriminating high and low affinity protein-protein complexes and it achieved an accuracy of 87% with an F1 score of 0.86 using 10-fold cross-validation on the selected features. Our approach presents an efficient tool for studying PPIs and provides crucial insights into the underlying mechanisms of the molecular recognition process. The web server can be freely accessed at https://web.iitm.ac.in/bioinfo2/DeepPPAPred/index.html     

Ecotin: lessons on survival in a protease-filled world.

Ecotin, an Escherichia coli periplasmic protein of 142 amino acids, has been shown to be a potent inhibitor of a group of homologous serine proteases with widely differing substrate recognition. It is highly effective against a number of enzymes, including both pancreatic and neutrophil-derived elastases, chymotrypsin, trypsin, factor Xa, and kallikrein. Recent structural and functional studies on ecotin and its interactions with different serine proteases have clarified these initial observations and revealed the remarkable features of this protein in inhibiting a strikingly large subset of the chymotrypsin family of serine proteases. The structures of the ecotin:serine protease complexes provide the first examples of protein-protein recognition where the concept of specificity of interactions needs to be reexamined. The binding sites show a fluidity of protein contacts derived from ecotin's innate flexibility in fitting itself to proteases while strongly interfering with their function.