Specificity in protein-protein interactions: the structural basis for dual recognition in endonuclease colicin-immunity protein complexes

Bacteria producing endonuclease colicins are protected against their cytotoxic activity by virtue of a small immunity protein that binds with high affinity and specificity to inactivate the endonuclease. DNase binding by the immunity protein occurs through a "dual recognition" mechanism in which conserved residues from helix III act as the binding-site anchor, while variable residues from helix II define specificity. We now report the 1.7 A crystal structure of the 24.5 kDa complex formed between the endonuclease domain of colicin E9 and its cognate immunity protein Im9, which provides a molecular rationale for this mechanism. Conserved residues of Im9 form a binding-energy hotspot through a combination of backbone hydrogen bonds to the endonuclease, many via buried solvent molecules, and hydrophobic interactions at the core of the interface, while the specificity-determining residues interact with corresponding specificity side-chains on the enzyme. Comparison between the present structure and that reported recently for the colicin E7 endonuclease domain in complex with Im7 highlights how specificity is achieved by very different interactions in the two complexes, predominantly hydrophobic in nature in the E9-Im9 complex but charged in the E7-Im7 complex. A key feature of both complexes is the contact between a conserved tyrosine residue from the immunity proteins (Im9 Tyr54) with a specificity residue on the endonuclease directing it toward the specificity sites of the immunity protein. Remarkably, this tyrosine residue and its neighbour (Im9 Tyr55) are the pivots of a 19 degrees rigid-body rotation that relates the positions of Im7 and Im9 in the two complexes. This rotation does not affect conserved immunity protein interactions with the endonuclease but results in different regions of the specificity helix being presented to the enzyme.     

PRI-Modeler: extracting RNA structural elements from PDB files of protein-RNA complexes

A complete understanding of protein and RNA structures and their interactions is important for determining the binding sites in protein-RNA complexes. Computational approaches exist for identifying secondary structural elements in proteins from atomic coordinates. However, similar methods have not been developed for RNA, due in part to the very limited structural data so far available. We have developed a set of algorithms for extracting and visualizing secondary and tertiary structures of RNA and for analyzing protein-RNA complexes. These algorithms have been implemented in a web-based program called PRI-Modeler (protein-RNA interaction modeler). Given one or more protein data bank files of protein-RNA complexes, PRI-Modeler analyzes the conformation of the RNA, calculates the hydrogen bond (H bond) and van der Waals interactions between amino acids and nucleotides, extracts secondary and tertiary RNA structure elements, and identifies the patterns of interactions between the proteins and RNAs. This paper presents PRI-Modeler and its application to the hydrogen bond and van der Waals interactions in the most representative set of protein-RNA complexes. The analysis reveals several interesting interaction patterns at various levels. The information provided by PRI-Modeler should prove useful for determining the binding sites in protein-RNA complexes. PRI-Modeler is accessible at http://wilab.inha.ac.kr/primodeler/, and supplementary materials are available in the analysis results section at http://wilab.inha.ac.kr/primodeler/.     

First observation of solution structures of bradykinin–penta-O-galloyl-d-glucopyranose complexes as determined by NMR and simulated annealing

Polyphenols (tannins) are known for their high propensity to precipitate proteins. They bind most strongly to proteins with a high proline content. Understanding the mechanism of this association is of prime interest because this interaction might induce protein conformational changes that may modify their biological activity. To investigate the interaction, an NMR study was carried out on the binding of a representative polyphenol, penta-O-galloyl-d-glucopyranose, to a nonapeptide hormone, bradykinin (BDK), where proline accounts for 30% of residues. Series of 1D and 2D-NMR experiments were performed. For the first time, a three-dimensional structure of complexes was determined using 2D-NMR experiments and molecular modeling. These structure calculations are a potent tool to understand how the association arises. They clearly show that the interaction is a complex phenomenon where several parameters are involved. The PGG/BDK complexes are formed by multiple weak interactions between peptide side chains and galloyl rings. Proline and arginine are good anchoring points and the glycine gives a certain flexibility in the peptide backbone that allows the polyphenol to approach and interact. Therefore, it is not only the hydrophobic stackings between galloyl rings and proline and hydrogen bonding involving arginine and aromatic rings which are important. The residue sequence and the side chain steric bulk also intervene.     

Coulombic and noncoulombic effect of polyanions on cytochrome c structure

The properties of the complexes of ferricytochrome c with two different polyanions—poly(vinylsulfate) and poly(4-styrene-sulfonate)—with a comparable charge density but with the different size of the uncharged part of their molecules have been studied by means of optical spectroscopy, differential scanning calorimetry, and gel chromatography. Ferriccytochrome c formed a complex with the former one through coulombic interactions and remained in a native-like state. The addition of the second polyanion to a solution of ferric cytochrome c at a low ionic strength, pH 7.0, resulted in profound conformational change in the hydrophobic core of protein (opening of the heme crevice with a perturbation of the methionine 80-heme iron bond and the hydrophobic core of the protein). These may be understood as an involvement of noncoulombic (hydrophobic, H-bonding) interactions of the uncharged part of the polyanion molecule. Conformational changes and the observed shift in acidic transition from low spin to high spin state of ferric cytochrome c detected in the presence of the polyanions may have biological implication in understanding the origin of conformational changes in proteins induced in the course of their interaction with membrane lipids and membrane proteins. © 1998 John Wiley & Sons, Inc. Biopoly 46: 145–154, 1998     

First observation of solution structures of bradykinin-penta-O-galloyl-D-glucopyranose complexes as determined by NMR and simulated annealing

Polyphenols (tannins) are known for their high propensity to precipitate proteins. They bind most strongly to proteins with a high proline content. Understanding the mechanism of this association is of prime interest because this interaction might induce protein conformational changes that may modify their biological activity. To investigate the interaction, an NMR study was carried out on the binding of a representative polyphenol, penta-O-galloyl-D-glucopyranose, to a nonapeptide hormone, bradykinin (BDK), where proline accounts for 30% of residues. Series of 1D and 2D-NMR experiments were performed. For the first time, a three-dimensional structure of complexes was determined using 2D-NMR experiments and molecular modeling. These structure calculations are a potent tool to understand how the association arises. They clearly show that the interaction is a complex phenomenon where several parameters are involved. The PGG/BDK complexes are formed by multiple weak interactions between peptide side chains and galloyl rings. Proline and arginine are good anchoring points and the glycine gives a certain flexibility in the peptide backbone that allows the polyphenol to approach and interact. Therefore, it is not only the hydrophobic stackings between galloyl rings and proline and hydrogen bonding involving arginine and aromatic rings which are important. The residue sequence and the side chain steric bulk also intervene.     

Protein docking using a genetic algorithm

A genetic algorithm (GA) for protein-protein docking is described, in which the proteins are represented by dot surfaces calculated using the Connolly program. The GA is used to move the surface of one protein relative to the other to locate the area of greatest surface complementarity between the two. Surface dots are deemed complementary if their normals are opposed, their Connolly shape type is complementary, and their hydrogen bonding or hydrophobic potential is fulfilled. Overlap of the protein interiors is penalized. The GA is tested on 34 large protein-protein complexes where one or both proteins has been crystallized separately. Parameters are established for which 30 of the complexes have at least one near-native solution ranked in the top 100. We have also successfully reassembled a 1,400-residue heptamer based on the top-ranking GA solution obtained when docking two bound subunits.     

Structural analysis of protein‐ligand interactions: the binding of endogenous compounds and of synthetic drugs

The large number of macromolecular structures deposited with the Protein Data Bank (PDB) describing complexes between proteins and either physiological compounds or synthetic drugs made it possible a systematic analysis of the interactions occurring between proteins and their ligands. In this work, the binding pockets of about 4000 PDB protein‐ligand complexes were investigated and amino acid and interaction types were analyzed. The residues observed with lowest frequency in protein sequences, Trp, His, Met, Tyr, and Phe, turned out to be the most abundant in binding pockets. Significant differences between drug‐like and physiological compounds were found. On average, physiological compounds establish with respect to drugs about twice as many hydrogen bonds with protein atoms, whereas drugs rely more on hydrophobic interactions to establish target selectivity. The large number of PDB structures describing homologous proteins in complex with the same ligand made it possible to analyze the conservation of binding pocket residues among homologous protein structures bound to the same ligand, showing that Gly, Glu, Arg, Asp, His, and Thr are more conserved than other amino acids. Also in the cases in which the same ligand is bound to unrelated proteins, the binding pockets showed significant conservation in the residue types. In this case, the probability of co‐occurrence of the same amino acid type in the binding pockets could be up to thirteen times higher than that expected on a random basis. The trends identified in this study may provide an useful guideline in the process of drug design and lead optimization. Copyright © 2014 John Wiley & Sons, Ltd.     

Protein domain interfaces: characterization and comparison with oligomeric protein interfaces

The physical and chemical properties of domain-domain interactions have been analysed in two-domain proteins selected from the protein classification, CATH. The two-domain structures were divided into those derived from (i) monomeric proteins, or (ii) oligomeric or complexed proteins. The size, polarity, hydrogen bonding and packing of the intra-chain domain interface were calculated for both sets of two-domain structures. The results were compared with inter-chain interface parameters from permanent and non-obligate protein-protein complexes. In general, the intra-chain domain and inter-chain interfaces were remarkably similar. Many of the intra-chain interface properties are intermediate between those calculated for permanent and non-obligate inter-chain complexes. Residue interface propensities were also found to be very similar, with hydrophobic residues playing a major role, together with positively charged arginine residues. In addition, the residue composition of the domain interfaces were found to be more comparable with domain surfaces than domain cores. The implications of these results for domain swapping and protein folding are discussed.     

Potential of mean force for protein-protein interaction studies.

Calculating protein-protein interaction energies is crucial for understanding protein-protein associations. On the basis of the methodology of mean-field potential, we have developed an empirical approach to estimate binding free energy for protein-protein interactions. This knowledge-based approach has been used to derive distance-dependent free energies of protein complexes from a nonredundant training set in the Protein Data Bank (PDB), with a careful treatment of homology. We calculate atom pair potentials for 16 pair interactions, which can reflect the importance of hydrophobic interactions and specific hydrogen-bonding interactions. The derived potentials for hydrogen-bonding interactions show a valley of favorable interactions at a distance of approximately 3 A, corresponding to that of an established hydrogen bond. For the test set of 28 protein complexes, the calculated energies have a correlation coefficient of 0.75 compared with experimental binding free energies. The performance of the method in ranking the binding energies of different protein-protein complexes shows that the energy estimation can be applied to value binding free energies for protein-protein associations.     

Effects of low urea concentrations on protein-water interactions

Solvent properties of aqueous media (dipolarity/polarizability, hydrogen bond donor acidity, and hydrogen bond acceptor basicity) were measured in the coexisting phases of Dextran–PEG aqueous two-phase systems (ATPSs) containing .5 and 2.0 M urea. The differences between the electrostatic and hydrophobic properties of the phases in the ATPSs were quantified by analysis of partitioning of the homologous series of sodium salts of dinitrophenylated amino acids with aliphatic alkyl side chains. Furthermore, partitioning of eleven different proteins in the ATPSs was studied. The analysis of protein partition behavior in a set of ATPSs with protective osmolytes (sorbitol, sucrose, trehalose, and TMAO) at the concentration of .5 M, in osmolyte-free ATPS, and in ATPSs with .5 or 2.0 M urea in terms of the solvent properties of the phases was performed. The results show unambiguously that even at the urea concentration of .5 M, this denaturant affects partitioning of all proteins (except concanavalin A) through direct urea–protein interactions and via its effect on the solvent properties of the media. The direct urea–protein interactions seem to prevail over the urea effects on the solvent properties of water at the concentration of .5 M urea and appear to be completely dominant at 2.0 M urea concentration.     

Charge centers and formation of the protein folding core

Electrostatic interactions are important for protein folding. At low resolution, the electrostatic field of the whole molecule can be described in terms of charge center(s). To study electrostatic effects, the centers of positive and negative charge were calculated for 20 small proteins of known structure, for which hydrogen exchange cores had been determined experimentally. Two observations seem to be important. First, in all 20 proteins studied 30-100% of the residues forming hydrogen exchange core(s) were clustered around the charge centers. Moreover, in each protein more than half of the core sequences are located near the centers of charge. Second, the general architecture of all-alpha proteins from the set seems to be stabilized by interactions of residues surrounding the charge centers. In most of the alpha-beta proteins, either or both of the centers are located near a pair of consecutive strands, and this is even more characteristic for alpha/Beta and all-beta structures. Consecutive strands are very probable sites of early folding events. These two points lead to the conclusion that charge centers, defined solely from the structure of the folded protein may indicate the location of a protein's hydrogen exchange/folding core. In addition, almost all the proteins contain well-conserved continuous hydrophobic sequences of three or more residues located in the vicinity of the charge centers. These hydrophobic sequences may be primary nucleation sites for protein folding. The results suggest the following scheme for the order of events in folding: local hydrophobic nucleation, electrostatic collapse of the core, global hydrophobic collapse, and slow annealing to the native state. This analysis emphasizes the importance of treating electrostatics during protein-folding simulations.     

Inversion of the balance between hydrophobic and hydrogen bonding interactions in protein folding and aggregation

Identifying the forces that drive proteins to misfold and aggregate, rather than to fold into their functional states, is fundamental to our understanding of living systems and to our ability to combat protein deposition disorders such as Alzheimer's disease and the spongiform encephalopathies. We report here the finding that the balance between hydrophobic and hydrogen bonding interactions is different for proteins in the processes of folding to their native states and misfolding to the alternative amyloid structures. We find that the minima of the protein free energy landscape for folding and misfolding tend to be respectively dominated by hydrophobic and by hydrogen bonding interactions. These results characterise the nature of the interactions that determine the competition between folding and misfolding of proteins by revealing that the stability of native proteins is primarily determined by hydrophobic interactions between side-chains, while the stability of amyloid fibrils depends more on backbone intermolecular hydrogen bonding interactions.     

Structure-based analysis of protein-RNA interactions using the program ENTANGLE

Until recently, drawing general conclusions about RNA recognition by proteins has been hindered by the paucity of high-resolution structures. We have analyzed 45 PDB entries of protein-RNA complexes to explore the underlying chemical principles governing both specific and non-sequence specific binding. To facilitate the analysis, we have constructed a database of interactions using ENTANGLE, a JAVA-based program that uses available structural models in their PDB format and searches for appropriate hydrogen bonding, stacking, electrostatic, hydrophobic and van der Waals interactions. The resulting database of interactions reveals correlations that suggest the basis for the discrimination of RNA from DNA and for base-specific recognition. The data illustrate both major and minor interaction strategies employed by families of proteins such as tRNA synthetases, ribosomal proteins, or RNA recognition motifs with their RNA targets. Perhaps most surprisingly, specific RNA recognition appears to be mediated largely by interactions of amide and carbonyl groups in the protein backbone with the edge of the RNA base. In cases where a base accepts a proton, the dominant amino acid donor is arginine, whereas in cases where the base donates a proton, the predominant acceptor is the backbone carbonyl group, not a side-chain group. This is in marked contrast to DNA-protein interactions, which are governed predominantly by amino acid side-chain interactions with functional groups that are presented in the accessible major groove. RNA recognition often proceeds through loops, bulges, kinks and other irregular structures that permit use of all the RNA functional groups and this is seen throughout the protein-RNA interaction database.     

Protein‐spanning water networks and implications for prediction of protein–protein interactions mediated through hydrophobic effects

Hydrophobic effects, often conflated with hydrophobic forces, are implicated as major determinants in biological association and self‐assembly processes. Protein–protein interactions involved in signaling pathways in living systems are a prime example where hydrophobic effects have profound implications. In the context of protein–protein interactions, a priori knowledge of relevant binding interfaces (i.e., clusters of residues involved directly with binding interactions) is difficult. In the case of hydrophobically mediated interactions, use of hydropathy‐based methods relying on single residue hydrophobicity properties are routinely and widely used to predict propensities for such residues to be present in hydrophobic interfaces. However, recent studies suggest that consideration of hydrophobicity for single residues on a protein surface require accounting of the local environment dictated by neighboring residues and local water. In this study, we use a method derived from percolation theory to evaluate spanning water networks in the first hydration shells of a series of small proteins. We use residue‐based water density and single‐linkage clustering methods to predict hydrophobic regions of proteins; these regions are putatively involved in binding interactions. We find that this simple method is able to predict with sufficient accuracy and coverage the binding interface residues of a series of proteins. The approach is competitive with automated servers. The results of this study highlight the importance of accounting of local environment in determining the hydrophobic nature of individual residues on protein surfaces. Proteins 2014; 82:3312–3326. © 2014 Wiley Periodicals, Inc.     

Aggregating Ability of Seed Storage Proteins from Cereals Differing in Gluten Quality

The effects of medium pH, ionic strength, and composition on the formation of macrocomplexes of seed storage proteins from wheat, rye, and barley have been studied. Various noncovalent interactions (electrostatic and hydrophobic interactions and hydrogen bonds) are involved in protein aggregation. Their combined action depends significantly on the biochemical nature of the storage proteins and on the medium.     

Electron paramagnetic resonance studies of radical species of proanthocyanidins and gallate esters

The polyphenols present in green tea or red wine comprise both regular flavon(ol)s and proanthocyanidins, i.e., derivatives of flavan-3-ols, whose distinct antioxidative potential is of great importance for explaining the beneficial effects of these nutrient beverages. Using EPR spectroscopy, we investigated radical structures obtained after oxidation of the parent compounds either by horseradish peroxidase/hydrogen peroxide or after autoxidation in moderately alkaline solutions. Both proanthocyanidins (monomers of condensed tannins, e.g., (+)-catechin, (-)-epicatechin, (-)-epicatechin gallate, (-)-epigallocatechin, (-)-epigallocatechin gallate, Pycnogenol) and gallate esters (hydrolyzable tannins, e.g., propylgallate, beta-glucogallin, pentagalloyl glucose and tannic acid) yielded predominantly semiquinone structures derived from the parent catechol or pyrogallol moieties. Evidence for a time-dependent oligomerization was obtained for (-)-epigallocatechin gallate, supporting our hypothesis that o-quinones formed from the initial semiquinone disproportionation are prone to nucleophilic addition reactions. Such phenolic coupling reactions would retain the numbers of reactive catechol/pyrogallol structures and thus the antioxidative potential. We therefore propose that proanthocyanidins are superior antioxidants as compared to flavon(ol)s proper, whose quinones are more likely to redox-cycle and act as prooxidants.     

Hydrophobic environment is a key factor for the stability of thermophilic proteins

The stability of thermophilic proteins has been viewed from different perspectives and there is yet no unified principle to understand this stability. It would be valuable to reveal the most important interactions for designing thermostable proteins for such applications as industrial protein engineering. In this work, we have systematically analyzed the importance of various interactions by computing different parameters such as surrounding hydrophobicity, inter‐residue interactions, ion‐pairs and hydrogen bonds. The importance of each interaction has been determined by its predicted relative contribution in thermophiles versus the same contribution in mesophilic homologues based on a dataset of 373 protein families. We predict that hydrophobic environment is the major factor for the stability of thermophilic proteins and found that 80% of thermophilic proteins analyzed showed higher hydrophobicity than their mesophilic counterparts. Ion pairs, hydrogen bonds, and interaction energy are also important and favored in 68%, 50%, and 62% of thermophilic proteins, respectively. Interestingly, thermophilic proteins with decreased hydrophobic environments display a greater number of hydrogen bonds and/or ion pairs. The systematic elimination of mesophilic proteins based on surrounding hydrophobicity, interaction energy, and ion pairs/hydrogen bonds, led to correctly identifying 95% of the thermophilic proteins in our analyses. Our analysis was also applied to another, more refined set of 102 thermophilic–mesophilic pairs, which again identified hydrophobicity as a dominant property in 71% of the thermophilic proteins. Further, the notion of surrounding hydrophobicity, which characterizes the hydrophobic behavior of residues in a protein environment, has been applied to the three‐dimensional structures of elongation factor‐Tu proteins and we found that the thermophilic proteins are enriched with a hydrophobic environment. The results obtained in this work highlight the importance of hydrophobicity as the dominating characteristic in the stability of thermophilic proteins, and we anticipate this will be useful in our attempts to engineering thermostable proteins. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Does structural and chemical divergence play a role in precluding undesirable protein interactions?

To understand the evolutionary forces establishing, maintaining, breaking, or precluding protein-protein interactions, a comprehensive data set of protein complexes has been analyzed to examine the overlap between protein interfaces and the most conserved or divergent protein surface areas. The most divergent areas tend to be found predominantly away from protein interfaces, although when found at interfaces, they are associated with specific lack of cross-reactivity between close homologues, like in antibody-antigen complexes. Moreover, the amino acid composition of highly variable regions is significantly different from any other protein surfaces. The variable regions present higher structural plasticity as a result of insertions and deletions, and favor charged over hydrophobic residues, a known strategy to minimize aggregation. This suggests that (1) a rapid rate of mutations at these regions might be continuously altering their properties, making difficult the coadaptation, in shape and chemical complementarity, to potential interacting partners; and (2) the existence of some form of selective pressure for variable areas away from interfaces to accumulate charged residues, perhaps as an evolutionary mechanism to increase solubility and minimize undesirable interactions within the crowded cellular environment. Finally, these results are placed into the context of the aberrant oligomerization of sickle-cell anemia hemoglobin and prion proteins.     

Aggregating ability of seed storage proteins from cereals differing in gluten quality

The effects of pH, ionic strength, and medium composition on formation of macrocomplexes of seed storage proteins from wheat, rye, and barley have been studied. It has been found that various noncovalent interactions (electrostatic and hydrophobic interactions and hydrogen bonds) are involved in protein aggregation. Their combined action depends significantly on the biochemical nature of storage proteins and on the medium.     

C‐H…O hydrogen bonds in FK506‐binding protein–ligand interactions

Hydrogen bonds are important interaction forces observed in protein structures. They can be classified as stronger or weaker depending on their energy, thereby reflecting on the type of donor. The contribution of weak hydrogen bonds is deemed as an important factor toward structure stability along with the stronger bonds. One such bond, the C‐H…O type hydrogen bond, is shown to make a contribution in maintaining three dimensional structures of proteins. Apart from their presence within protein structures, the role of these bonds in protein–ligand interactions is also noteworthy. In this study, we present a statistical analysis on the presence of C‐H…O hydrogen bonds observed between FKBPs and their cognate ligands. The FK506‐binding proteins (FKBPs) carry peptidyl cis–trans isomerase activity apart from the immunosuppressive property by binding to the immunosuppressive drugs FK506 or rapamycin. Because the active site of FKBPs is lined up by many hydrophobic residues, we speculated that the prevalence of C‐H…O hydrogen bonds will be considerable. In a total of 25 structures analyzed, a higher frequency of C‐H…O hydrogen bonds is observed in comparison with the stronger hydrogen bonds. These C‐H…O hydrogen bonds are dominated by a highly conserved donor, the C^α/β of Val55 and an acceptor, the backbone oxygen of Glu54. Both these residues are positioned in the β4‐α1 loop, whereas the other residues Tyr26, Phe36 and Phe99 with higher frequencies are lined up at the opposite face of the active site. These preferences could be implicated in FKBP pharmacophore models toward enhancing the ligand affinity. This study could be a prelude to studying other proteins with hydrophobic pockets to gain better insights into ligand recognition. Copyright © 2013 John Wiley & Sons, Ltd.