Protein-protein interactions: hot spots and structurally conserved residues often locate in complemented pockets that pre-organized in the unbound states: implications for docking

Energetic hot spots account for a significant portion of the total binding free energy and correlate with structurally conserved interface residues. Here, we map experimentally determined hot spots and structurally conserved residues to investigate their geometrical organization. Unfilled pockets are pockets that remain unfilled after protein-protein complexation, while complemented pockets are pockets that disappear upon binding, representing tightly fit regions. We find that structurally conserved residues and energetic hot spots are strongly favored to be located in complemented pockets, and are disfavored in unfilled pockets. For the three available protein-protein complexes with complemented pockets where both members of the complex were alanine-scanned, 62% of all hot spots (DeltaDeltaG>2kcal/mol) are within these pockets, and 60% of the residues in the complemented pockets are hot spots. 93% of all red-hot residues (DeltaDeltaG>/=4kcal/mol) either protrude into or are located in complemented pockets. The occurrence of hot spots and conserved residues in complemented pockets highlights the role of local tight packing in protein associations, and rationalizes their energetic contribution and conservation. Complemented pockets and their corresponding protruding residues emerge among the most important geometric features in protein-protein interactions. By screening the solvent, this organization shields backbone hydrogen bonds and charge-charge interactions. Complemented pockets often pre-exist binding. For 18 protein-protein complexes with complemented pockets whose unbound structures are available, in 16 the pockets are identified to pre-exist in the unbound structures. The root-mean-squared deviations of the atoms lining the pockets between the bound and unbound states is as small as 0.9A, suggesting that such pockets constitute features of the populated native state that may be used in docking.     

Conservation of polar residues as hot spots at protein interfaces

A number of studies have addressed the question of which are the critical residues at protein-binding sites. These studies examined either a single or a few protein-protein interfaces. The most extensive study to date has been an analysis of alanine-scanning mutagenesis. However, although the total number of mutations was large, the number of protein interfaces was small, with some of the interfaces closely related. Here we show that although overall binding sites are hydrophobic, they are studded with specific, conserved polar residues at specific locations, possibly serving as energy "hot spots." Our results confirm and generalize the alanine-scanning data analysis, despite its limited size. Previously Trp, Arg, and Tyr were shown to constitute energetic hot spots. These were rationalized by their polar interactions and by their surrounding rings of hydrophobic residues. However, there was no compelling reason as to why specifically these residues were conserved. Here we show that other polar residues are similarly conserved. These conserved residues have been detected consistently in all interface families that we have examined. Our results are based on an extensive examination of residues which are in contact across protein interfaces. We utilize all clustered interface families with at least five members and with sequence similarity between the members in the range of 20-90%. There are 11 such clustered interface families, comprising a total of 97 crystal structures. Our three-dimensional superpositioning analysis of the occurrences of matched residues in each of the families identifies conserved residues at spatially similar environments. Additionally, in enzyme inhibitors, we observe that residues are more conserved at the interfaces than at other locations. On the other hand, antibody-protein interfaces have similar surface conservation as compared to their corresponding linear sequence alignment, consistent with the suggestion that evolution has optimized protein interfaces for function.     

Hot regions in protein--protein interactions: the organization and contribution of structurally conserved hot spot residues

Structurally conserved residues at protein-protein interfaces correlate with the experimental alanine-scanning hot spots. Here, we investigate the organization of these conserved, computational hot spots and their contribution to the stability of protein associations. We find that computational hot spots are not homogeneously distributed along the protein interfaces; rather they are clustered within locally tightly packed regions. Within the dense clusters, they form a network of interactions and consequently their contributions to the stability of the complex are cooperative; however the contributions of independent clusters are additive. This suggests that the binding free energy is not a simple summation of the single hot spot residue contributions. As expected, around the hot spot residues we observe moderately conserved residues, further highlighting the crucial role of the conserved interactions in the local densely packed environment. The conserved occurrence of these organizations suggests that they are advantageous for protein-protein associations. Interestingly, the total number of hydrogen bonds and salt bridges contributed by hot spots is as expected. Thus, H-bond forming residues may use a "hot spot for water exclusion" mechanism. Since conserved residues are located within highly packed regions, water molecules are easily removed upon binding, strengthening electrostatic contributions of charge-charge interactions. Hence, the picture that emerges is that protein-protein associations are optimized locally, with the clustered, networked, highly packed structurally conserved residues contributing dominantly and cooperatively to the stability of the complex. When addressing the crucial question of "what are the preferred ways of proteins to associate", these findings point toward a critical involvement of hot regions in protein-protein interactions.     

Protein-protein interactions within an E2-RING finger complex. Implications for ubiquitin-dependent DNA damage repair

The RING finger protein RAD5 interacts and cooperates with the UBC13-MMS2 ubiquitin-conjugating enzyme in postreplication DNA damage repair in yeast. Previous observations implied that the function of UBC13 and MMS2 is dependent on the presence of RAD5, suggesting that the RING finger protein might act as a ubiquitin-protein ligase specific for the UBC13-MMS2 complex. In support of this notion it is shown here that the contact surfaces between the RAD5 RING domain and UBC13 correspond to those found in other pairs of ubiquitin-conjugating enzymes and ubiquitin-protein ligases. Mutations that compromise the protein-protein interactions either between the RING domain and UBC13 or within the UBC13-MMS2 dimer were found to have variable effects on repair activity in vivo that strongly depended on the expression levels of the corresponding mutants. Quantitative analysis of the affinity and kinetics of the UBC13-MMS2 interaction suggests a highly dynamic association model in which compromised mutual interactions result in phenotypic effects only under conditions where protein levels become limiting. Finally, this study demonstrates that beyond its cooperation with the UBC13-MMS2 dimer, RAD5 must have an additional role in DNA damage repair independent of its RING finger domain.     

蛋白质-蛋白质界面热点残基预测及其在线工具

蛋白质-蛋白质结合热点是界面中对结合自由能有着显著贡献的一小簇残基。捕捉和揭示这类热点残基可以加深对蛋白质间相互作用机制的理解,为蛋白质工程和药物设计提供指导。但实验技术费时费力且代价昂贵。计算工具可用于辅助和补充实验上的尝试。该文较详细、系统地介绍了蛋白质界面热点的特性、计算预测的策略与技术,并应用实例进一步说明这些方法学的特征;还介绍了界面热点的数据库和一些主要的在线预测工具,旨在为设计、挑选和应用这类工具解决特定问题的研究人员提供指南。     

Prediction of hot spot residues at protein-protein interfaces by combining machine learning and energy-based methods

Background  

Alanine scanning mutagenesis is a powerful experimental methodology for investigating the structural and energetic characteristics of protein complexes. Individual amino-acids are systematically mutated to alanine and changes in free energy of binding (ΔΔG) measured. Several experiments have shown that protein-protein interactions are critically dependent on just a few residues ("hot spots") at the interface. Hot spots make a dominant contribution to the free energy of binding and if mutated they can disrupt the interaction. As mutagenesis studies require significant experimental efforts, there is a need for accurate and reliable computational methods. Such methods would also add to our understanding of the determinants of affinity and specificity in protein-protein recognition.     

Mutagenesis of histidinol dehydrogenase reveals roles for conserved histidine residues.

The dimeric zinc metalloenzyme L-histidinol dehydrogenase (HDH) catalyzes an unusual four-electron oxidation of the amino alcohol histidinol via the histidinaldehyde intermediate to the acid product histidine with the reduction of two molecules of NAD. An essential base, with pKa about 8, is involved in catalysis. Here we report site-directed mutagenesis studies to replace each of the five histidine residues (His-98, His-261, His-326, His-366, and His-418) in Salmonella typhimurium with either asparagine or glutamine. In all cases, the overexpressed enzymes were readily purified and behaved as dimers. Substitution of His-261 and His-326 by asparagine caused about 7000- and 500-fold decreases in kcat, respectively, with little change in KM values. Similar loss of activity was also reported for a H261N mutant Brassica HDH [Nagai, A., and Ohta, D. (1994) J. Biochem. 115, 22-25]. Kinetic isotope effects, pH profiles, substrate rescue, and stopped-flow experiments suggested that His-261 and His-326 are involved in proton transfers during catalysis. Sensitivity to metal ion chelator and decreased affinities for metal ions with substitutions at His-261 and His-418 suggested that these two residues are candidates for zinc ion ligands.     

An analysis of the periodicity of conserved residues in sequence alignments of G-protein coupled receptors. Implications for the three-dimensional structure

Twenty-three sequences from the family of G-protein coupled receptors have been aligned according to the 'historical alignment' procedure of Feng and Doolittle. Fourier transform analysis of this reveals that parts of five of the seven putative membrane-spanning regions exhibit a periodicity of conserved/nonconserved residues which is compatible with the periodicity of the alpha-helix. This would place the conserved residues on one side of the helix, which may face the inside of the proposed seven membered helical bundle.     

Evidence for structurally conserved recognition of the major carbohydrate xenoantigen by natural antibodies.

Natural or preformed antibodies that react with oligosaccharides bearing terminal galactose-alpha(1,3)-galactose [Gal alpha(1,3)Gal] stuctures are present in the sera of all humans. Antibodies against Gal alpha(1,3)Gal epitopes initiate hyperacute rejection of xenografts of porcine organs in human recipients. Despite the enormous clinical potential for xenotransplantation, very little is known about the 3D structural basis for natural antibody recognition of the major xenoantigen (i.e. Gal alpha(1,3)Gal). In this review, we discuss general binding patterns that have been repeatedly identified in antibody complexes with small molecules (haptens), carbohydrate and peptide ligands because similar mechanisms will almost certainly mediate recognition of the major xenoantigen by natural antibodies.     

Protein-DNA interactions: structural, thermodynamic and clustering patterns of conserved residues in DNA-binding proteins

Amino acid residues, which play important roles in protein function, are often conserved. Here, we analyze thermodynamic and structural data of protein-DNA interactions to explore a relationship between free energy, sequence conservation and structural cooperativity. We observe that the most stabilizing residues or putative hotspots are those which occur as clusters of conserved residues. The higher packing density of the clusters and available experimental thermodynamic data of mutations suggest cooperativity between conserved residues in the clusters. Conserved singlets contribute to the stability of protein-DNA complexes to a lesser extent. We also analyze structural features of conserved residues and their clusters and examine their role in identifying DNA-binding sites. We show that about half of the observed conserved residue clusters are in the interface with the DNA, which could be identified from their amino acid composition; whereas the remaining clusters are at the protein-protein or protein-ligand interface, or embedded in the structural scaffolds. In protein-protein interfaces, conserved residues are highly correlated with experimental residue hotspots, contributing dominantly and often cooperatively to the stability of protein-protein complexes. Overall, the conservation patterns of the stabilizing residues in DNA-binding proteins also highlight the significance of clustering as compared to single residue conservation.     

A catalytic function for the structurally conserved residue Phe 100 of ribonuclease T1.

The function of the conserved Phe 100 residue of RNase T1 (EC 3.1.27.3) has been investigated by site-directed mutagenesis and X-ray crystallography. Replacement of Phe 100 by alanine results in a mutant enzyme with kcat reduced 75-fold and a small increase in Km for the dinucleoside phosphate substrate GpC. The Phe 100 Ala substitution has similar effects on the turnover rates of GpC and its minimal analogue GpOMe, in which the leaving cytidine is replaced by methanol. The contribution to catalysis is independent of the nature of the leaving group, indicating that Phe 100 belongs to the primary site. The contribution of Phe 100 to catalysis may result from a direct van der Waals contact between its aromatic ring and the phosphate moiety of the substrate. Phe 100 may also contribute to the positioning of the pentacovalent phosphorus of the transition state, relative to other catalytic residues. If compared to the corresponding wild-type data, the structural implications of the mutation in the present crystal structure of Phe 100 Ala RNase T1 complexed with the specific inhibitor 2'-GMP are restricted to the active site. Repositioning of 2'-GMP, caused by the Phe 100 Ala mutation, generates new or improved contacts of the phosphate moiety with Arg 77 and His 92. In contrast, interactions with the Glu 58 carboxylate appear to be weakened. The effects of the His 92 Gln and Phe 100 Ala mutations on GpC turnover are additive in the corresponding double mutant, indicating that the contribution of Phe 100 to catalysis is independent of the catalytic acid His 92. The present results lead to the conclusion that apolar residues may contribute considerably to catalyze conversions of charged molecules to charged products, involving even more polar transition states.     

Protein-protein interactions of the light-harvesting chlorophyll protein. II. Evidence for two stages of cation independent association

In a previous paper, we observed a two-stage cation-independent association of the light-harvesting chlorophyll protein from spinach chloroplasts based on concentration-dependent changes in the sedimentation coefficient. The two stages of association occurred between (2–4) and (4–7) μg/ml chlorophyll. In this paper, we provide further evidence for this association.This includes: (1) A decrease in the number of divalent cation binding sites in the second stage of association. (2) A corresponding decrease in the extent of the cation-dependent association. (3) A positive deviation from Beer's law for chlorophyll b for both stages of the cation-independent association and a positive deviation for chlorophyll a for the second stage of association only. (4) A change in the fluorescence emission of both chlorophyll a and b. The change for chlorophyll b was observed for both steps of association whereas that for chlorophyll a was observed for the second step of association only. Therefore, the first stage of association affects only chlorophyll b whereas the second stage alters the environment of both chlorophyll a and b. (5) In addition, divalent cations quenched chlorophyll fluorescence. However, the quenching which required 200–300 μM divalent cation for half-maximal effects was related neither to divalent cation binding nor to the divalent cation-induced association of the protein.     

Protein-protein docking on molecular models of Aspergillus niger RNase and human actin: novel target for anticancer therapeutics

The 3D models of human actin protein and A.niger RNase were designed using the templates ACTBIND (PDB ID: 3D3Z) and crystalline profilin-beta-actin (PDB ID: 2BTF), respectively in Modeller9v5. These models are testified using several validation methods including PROCHECK, ERRAT, WHAT-IF, PROSA2003 and VERIFY-3D. The stereo-chemical quality of the models was judged by Ramachandran plot with PROCHECK. The total quality G-factor −0.2, shows a good quality model. The ERRAT score for the human actin and A.niger RNase models are 86.104 and 84.615, respectively, fit well within the range of a high quality model. The ERRAT score for the templates 2BTF and 3D3Z are 91.111 and 97.391, respectively. The WHAT-IF evaluation justifies a reasonable homology model structure as none of the scores for each residue in the homology model is lower than −5.0. The energy-minimized model of human actin with PROSA reveals the Z-score value −10.52 between native conformations of the crystal structures. The VERIFY 3D average score is 0.36. All evidence suggests that the geometric quality of the backbone conformation, the residue interaction, the residue contact and the energy profile of the structures were well within the limits of reliable structures. The interaction energy of docking was calculated using the HEX server. The Etotal, lowest docked energy, and calculated RMSD values were −1.608 kcal mol⁻¹, -8.369 kcal mol⁻¹ and 0.617 Å, respectively. The study presented in the current project may be useful to design molecules that may have anticancer activity.

     

Novel aromatic inhibitors of influenza virus neuraminidase make selective interactions with conserved residues and water molecules in the active site.

The active site of type A or B influenza virus neuraminidase is composed of 11 conserved residues that directly interact with the substrate, sialic acid. An aromatic benzene ring has been used to replace the pyranose of sialic acid in our design of novel neuraminidase inhibitors. A bis(hydroxymethyl)pyrrolidinone ring was constructed in place of the N-acetyl group on the sialic acid. The hydroxymethyl groups replace two active site water molecules, which resulted in the high affinity of the nanomolar inhibitors. However, these inhibitors have greater potency for type A influenza virus than for type B influenza virus. To resolve the differences, we determined the X-ray crystal structure of three benzoic acid substituted inhibitors bound to the active site of B/Lee/40 neuraminidase. The investigation of a hydrophobic aliphatic group and a hydrophilic guanidino group on the aromatic inhibitors shows changes in the interaction with the active site residue Glu275. The results provide an explanation for the difference in efficacy of these inhibitors against types A and B viruses, even though the 11 active site residues of the neuraminidase are conserved.     

Electromechanical coupling between skeletal and cardiac muscle. Implications for infarct repair

Skeletal myoblasts form grafts of mature muscle in injured hearts, and these grafts contract when exogenously stimulated. It is not known, however, whether cardiac muscle can form electromechanical junctions with skeletal muscle and induce its synchronous contraction. Here, we report that undifferentiated rat skeletal myoblasts expressed N-cadherin and connexin43, major adhesion and gap junction proteins of the intercalated disk, yet both proteins were markedly downregulated after differentiation into myo-tubes. Similarly, differentiated skeletal muscle grafts in injured hearts had no detectable N-cadherin or connexin43; hence, electromechanical coupling did not occur after in vivo grafting. In contrast, when neonatal or adult cardiomyocytes were cocultured with skeletal muscle, approximately 10% of the skeletal myotubes contracted in synchrony with adjacent cardiomyocytes. Isoproterenol increased myotube contraction rates by 25% in coculture without affecting myotubes in monoculture, indicating the cardiomyocytes were the pacemakers. The gap junction inhibitor heptanol aborted myotube contractions but left spontaneous contractions of individual cardiomyocytes intact, suggesting myotubes were activated via gap junctions. Confocal microscopy revealed the expression of cadherin and connexin43 at junctions between myotubes and neonatal or adult cardiomyocytes in vitro. After microinjection, myotubes transferred dye to neonatal cardiomyocytes via gap junctions. Calcium imaging revealed synchronous calcium transients in cardiomyocytes and myotubes. Thus, cardiomyocytes can form electromechanical junctions with some skeletal myotubes in coculture and induce their synchronous contraction via gap junctions. Although the mechanism remains to be determined, if similar junctions could be induced in vivo, they might be sufficient to make skeletal muscle grafts beat synchronously with host myocardium.     

Site-specific integration of genes into hot spots for recombination flanking aadA in Tn21 transposons.

Summary Tn21-related transposons are widespread among bacteria and carry various resistance determinants at preferential sites, hs1 and hs2. In an in vivo integrative recombination assay it was demonstrated that these hot spots direct the integration of aminoglycoside resistance genes like aadB from Klebsiella pneumoniae and aacAI from Serratia marcescens, in a recA ^– background. The maximum required recognition sequence which must be present in both the donor and recipient plasmids is 5 CTAAAACAAAGTTA 3 (hs2). The double-site-specific recombination occurred with a frequency of 10^–5–10^–6. The resulting structures include not only replicon fusion products but also more complex structures carrying two copies of the donor plasmid or simply the donor gene flanked by hs elements. hs1 and hs2 are thought to act as recognition sites for a trans-acting site-specific recombinase. By the use of Tn21 deletion derivatives, it has been shown that the recombinase is not encoded by Tn21. This new integrative recombination system is involved in the acquisition of new genes by Tn21-related transposons and their spread among bacterial populations.     

Patchiness in the dispersion of nectar resources: Evidence for hot and cold spots

Summary The dispersion pattern of resources can have a significant effect on foraging behavior. We examined the dispersion pattern of standing crop of nectar in a population of Delphinium nelsonii. The nectar content of flowers was measured for subject inflorescences and their two nearest neighbors. Inflorescences were divided into those which had some nectar (hot plants) and those which had none (cold plants). Tests for independence showed that subject plants and their neighbors were likely to have the same temperature, indicating that resources were patchy. The implications of this pattern for the foraging movements of bumblebees are discussed.The order of authorship was determined by the flip of a coin     

Multiple origins of the human glycophorin Sta gene. Identification of hot spots for independent unequal homologous recombinations.

Human glycophorin Sta (HGpSta), one of the structural variants of erythrocyte membrane sialoglycoproteins, is encoded by a delta-alpha hybrid gene that arose from a single unequal crossover between the parent HGpB(delta) and HGpA(alpha) genes. We report here the identification of two new HGpSta genes (type A and type B) in four unrelated Sta heterozygotes from two ethnic groups. These Sta genes represent distinct genetic isoforms that differ from the previously reported Sta gene (type C) in the location of crossing-over sites. Comparison of nucleotide sequences among HGpB(delta), HGpA(alpha), and HGpSta type A genes revealed that the delta-alpha unequal crossover for the Sta type A gene occurred 110-246 base pairs downstream from pseudoexon III. In the crossing-over site of this Sta gene, an AT-rich sequence lying 3' to a nonameric palindrome was found to be highly similar to the lambda phage attachment site, att B, in inverted orientation. In the Sta type B gene, the delta-alpha crossing-over point was localized to an AG-rich sequence that is 302-490 base pairs downstream from pseudoexon III. Multiple lambda chi-like elements were identified at the crossover boundaries and within the breakpoint of this Sta gene. These results suggest strongly that recurrent and independent unequal recombination events have occurred in the formation of multiple Sta genes and that particular genomic sequences are important in defining the recombination sites for these homology-driven processes.