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.     

Effect of methionine oxidation of a recombinant monoclonal antibody on the binding affinity to protein A and protein G

Oxidation of methionine (Met) residues is one of the most common protein degradation pathways. Two Met residues, Met256 and Met432, of a recombinant fully human monoclonal IgG1 antibody have been shown to be susceptible to oxidation. Met256 and Met432 are located in the antibody CH2-CH3 interface and in close proximity to protein A and protein G binding sites. The effect of oxidation of these susceptible Met residues on the binding to protein A and protein G was investigated in the current study. Incubation of the antibody with 5% tert-butyl hydroperoxide (tBHP) resulted in a nearly complete oxidation of Met256 and Met432, while incubation with 1% tBHP resulted in mixed populations of the antibody with different degrees of Met oxidation. Oxidation of Met256 and Met432 resulted in earlier elution of the antibody from protein A and protein G columns when eluted with a gradient of decreasing pH. Analysis by ELISA and surface plasmon resonance (SPR) revealed decreased binding affinity of the oxidized antibody to protein A and protein G. It is therefore concluded that oxidation of the Met256 and Met432 residues of the recombinant monoclonal antibody altered its interaction with protein A and protein G resulting in a decrease in binding affinity.     

Non‐interacting surface solvation and dynamics in protein–protein interactions

Protein–protein interactions control a plethora of cellular processes, including cell proliferation, differentiation, apoptosis, and signal transduction. Understanding how and why proteins interact will inevitably lead to novel structure‐based drug design methods, as well as design of de novo binders with preferred interaction properties. At a structural and molecular level, interface and rim regions are not enough to fully account for the energetics of protein–protein binding, even for simple lock‐and‐key rigid binders. As we have recently shown, properties of the global surface might also play a role in protein–protein interactions. Here, we report on molecular dynamics simulations performed to understand solvent effects on protein–protein surfaces. We compare properties of the interface, rim, and non‐interacting surface regions for five different complexes and their free components. Interface and rim residues become, as expected, less mobile upon complexation. However, non‐interacting surface appears more flexible in the complex. Fluctuations of polar residues are always lower compared with charged ones, independent of the protein state. Further, stable water molecules are often observed around polar residues, in contrast to charged ones. Our analysis reveals that (a) upon complexation, the non‐interacting surface can have a direct entropic compensation for the lower interface and rim entropy and (b) the mobility of the first hydration layer, which is linked to the stability of the protein–protein complex, is influenced by the local chemical properties of the surface. These findings corroborate previous hypotheses on the role of the hydration layer in shielding protein–protein complexes from unintended protein–protein interactions. Proteins 2015; 83:445–458. © 2014 Wiley Periodicals, Inc.     

Crystal structures of five bovine chymotrypsin complexes with P1 BPTI variants

The bovine chymotrypsin-bovine pancreatic trypsin inhibitor (BPTI) interaction belongs to extensively studied models of protein-protein recognition. The accommodation of the inhibitor P1 residue in the S1 binding site of the enzyme forms the hot spot of this interaction. Mutations introduced at the P1 position of BPTI result in a more than five orders of magnitude difference of the association constant values with the protease. To elucidate the structural aspects of the discrimination between different P1 residues, crystal structures of five bovine chymotrypsin-P1 BPTI variant complexes have been determined at pH 7.8 to a resolution below 2 A. The set includes polar (Thr), ionizable (Glu, His), medium-sized aliphatic (Met) and large aromatic (Trp) P1 residues and complements our earlier studies of the interaction of different P1 side-chains with the S1 pocket of chymotrypsin. The structures have been compared to the complexes of proteases with similar and dissimilar P1 preferences, including Streptomyces griseus proteases B and E, human neutrophil elastase, crab collagenase, bovine trypsin and human thrombin. The S1 sites of these enzymes share a common general shape of significant rigidity. Large and branched P1 residues adapt in their complexes similar conformations regardless of the polarity and size differences between their S1 pockets. Conversely, long and flexible residues such as P1 Met are present in the disordered form and display a conformational diversity despite similar inhibitory properties with respect to most enzymes studied. Thus, the S1 specificity profiles of the serine proteases appear to result from the precise complementarity of the P1-S1 interface and minor conformational adjustments occurring upon the inhibitor binding.     

Convergent signaling pathways—interaction between methionine oxidation and serine/threonine/tyrosine O-phosphorylation

Oxidation of methionine (Met) to Met sulfoxide (MetSO) is a frequently found reversible posttranslational modification. It has been presumed that the major functional role for oxidation-labile Met residues is to protect proteins/cells from oxidative stress. However, Met oxidation has been established as a key mechanism for direct regulation of a wide range of protein functions and cellular processes. Furthermore, recent reports suggest an interaction between Met oxidation and O-phosphorylation. Such interactions are a potentially direct interface between redox sensing and signaling, and cellular protein kinase/phosphatase-based signaling. Herein, we describe the current state of Met oxidation research, provide some mechanistic insight into crosstalk between these two major posttranslational modifications, and consider the evolutionary significance and regulatory potential of this crosstalk.     

Backbone dynamics in dihydrofolate reductase complexes: role of loop flexibility in the catalytic mechanism

To elucidate the influence of local motion of the polypeptide chain on the catalytic mechanism of an enzyme, we have measured (15)N relaxation data for Escherichia coli dihydrofolate reductase in three different complexes, representing different stages in the catalytic cycle of the enzyme. NMR relaxation data were analyzed by the model-free approach, corrected for rotational anisotropy, to provide insights into the backbone dynamics. There are significant differences in the backbone dynamics in the different complexes. Complexes in which the cofactor binding site is occluded by the Met20 loop display large amplitude motions on the picosecond/nanosecond time scale for residues in the Met20 loop, the adjacent betaF-betaG loop and for residues 67-69 in the adenosine binding loop. Formation of the closed Met20 loop conformation in the ternary complex with folate and NADP(+), results in attenuation of the motions in the Met20 loop and the betaF-betaG loop but leads to increased flexibility in the adenosine binding loop. New fluctuations on a microsecond/millisecond time scale are observed in the closed E:folate:NADP(+) complex in regions that form hydrogen bonds between the Met20 and the betaF-betaG loops. The data provide insights into the changes in backbone dynamics during the catalytic cycle and point to an important role of the Met20 and betaF-betaG loops in controlling access to the active site. The high flexibility of these loops in the occluded conformation is expected to promote tetrahydrofolate-assisted product release and facilitate binding of the nicotinamide ring to form the Michaelis complex. The backbone fluctuations in the Met20 loop become attenuated once it closes over the active site, thereby stabilizing the nicotinamide ring in a geometry conducive to hydride transfer. Finally, the relaxation data provide evidence for long-range motional coupling between the adenosine binding loop and distant regions of the protein.     

Prediction of interface residues in protein-protein complexes by a consensus neural network method: test against NMR data

The number of structures of protein-protein complexes deposited to the Protein Data Bank is growing rapidly. These structures embed important information for predicting structures of new protein complexes. This motivated us to develop the PPISP method for predicting interface residues in protein-protein complexes. In PPISP, sequence profiles and solvent accessibility of spatially neighboring surface residues were used as input to a neural network. The network was trained on native interface residues collected from the Protein Data Bank. The prediction accuracy at the time was 70% with 47% coverage of native interface residues. Now we have extensively improved PPISP. The training set now consisted of 1156 nonhomologous protein chains. Test on a set of 100 nonhomologous protein chains showed that the prediction accuracy is now increased to 80% with 51% coverage. To solve the problem of over-prediction and under-prediction associated with individual neural network models, we developed a consensus method that combines predictions from multiple models with different levels of accuracy and coverage. Applied on a benchmark set of 68 proteins for protein-protein docking, the consensus approach outperformed the best individual models by 3-8 percentage points in accuracy. To demonstrate the predictive power of cons-PPISP, eight complex-forming proteins with interfaces characterized by NMR were tested. These proteins are nonhomologous to the training set and have a total of 144 interface residues identified by chemical shift perturbation. cons-PPISP predicted 174 interface residues with 69% accuracy and 47% coverage and promises to complement experimental techniques in characterizing protein-protein interfaces. .     

Beauty Is in the Eye of the Beholder: Proteins Can Recognize Binding Sites of Homologous Proteins in More than One Way

Understanding the mechanisms of protein–protein interaction is a fundamental problem with many practical applications. The fact that different proteins can bind similar partners suggests that convergently evolved binding interfaces are reused in different complexes. A set of protein complexes composed of non-homologous domains interacting with homologous partners at equivalent binding sites was collected in 2006, offering an opportunity to investigate this point. We considered 433 pairs of protein–protein complexes from the ABAC database (AB and AC binary protein complexes sharing a homologous partner A) and analyzed the extent of physico-chemical similarity at the atomic and residue level at the protein–protein interface. Homologous partners of the complexes were superimposed using Multiprot, and similar atoms at the interface were quantified using a five class grouping scheme and a distance cut-off. We found that the number of interfacial atoms with similar properties is systematically lower in the non-homologous proteins than in the homologous ones. We assessed the significance of the similarity by bootstrapping the atomic properties at the interfaces. We found that the similarity of binding sites is very significant between homologous proteins, as expected, but generally insignificant between the non-homologous proteins that bind to homologous partners. Furthermore, evolutionarily conserved residues are not colocalized within the binding sites of non-homologous proteins. We could only identify a limited number of cases of structural mimicry at the interface, suggesting that this property is less generic than previously thought. Our results support the hypothesis that different proteins can interact with similar partners using alternate strategies, but do not support convergent evolution.     

Important role of methionine 145 in dimerization of bovine β-lactoglobulin

β-Lactoglobulin (LG) contains nine β-strands (strands A-I) and one α-helix. Strands A-H form a β-barrel. At neutral pH, bovine LG (BLG) forms a dimer and the dimer interface consists of AB-loops and the I-strands of two subunits. On the other hand, equine LG (ELG) is monomeric. The residues 145-153 of BLG, which compose a dimer interface, are entirely different from those of ELG. The difference in the association states between BLG and ELG can be attributed to the residues 145-153. To confirm this, we constructed a chimeric LG, ImBLG (I-strand mutated BLG), in which the residues 145-153 were replaced with those of ELG. Gel-filtration chromatography and analytical ultracentrifugation revealed that ImBLG existed as a monomer. To identify the residues important for dimerization, we constructed several revertants and investigated their association. This experiment revealed that, in addition to the interface residues (Ile147, Leu149 and Phe151), Met145 is critical for dimerization. Although Met145 does not contact with the other protomer, it seems to be important in determining the backbone conformation of the I-strand. This was supported by the fact that all Met145-containing mutants showed circular dichroism spectra similar to BLG but different from ImBLG.     

Crystal structure of the HGF beta-chain in complex with the Sema domain of the Met receptor

The Met tyrosine kinase receptor and its ligand, hepatocyte growth factor (HGF), play important roles in normal development and in tumor growth and metastasis. HGF-dependent signaling requires proteolysis from an inactive single-chain precursor into an active alpha/beta-heterodimer. We show that the serine protease-like HGF beta-chain alone binds Met, and report its crystal structure in complex with the Sema and PSI domain of the Met receptor. The Met Sema domain folds into a seven-bladed beta-propeller, where the bottom face of blades 2 and 3 binds to the HGF beta-chain 'active site region'. Mutation of HGF residues in the area that constitutes the active site region in related serine proteases significantly impairs HGF beta binding to Met. Key binding loops in this interface undergo conformational rearrangements upon maturation and explain the necessity of proteolytic cleavage for proper HGF signaling. A crystallographic dimer interface between two HGF beta-chains brings two HGF beta:Met complexes together, suggesting a possible mechanism of Met receptor dimerization and activation by HGF.     

Charged residues at protein interaction interfaces: unexpected conservation and orchestrated divergence

Protein-protein interactions play an essential role in the functioning of cell. The importance of charged residues and their diverse role in protein-protein interactions have been well studied using experimental and computational methods. Often, charged residues located in protein interaction interfaces are conserved across the families of homologous proteins and protein complexes. However, on a large scale, it has been recently shown that charged residues are significantly less conserved than other residue types in protein interaction interfaces. The goal of this work is to understand the role of charged residues in the protein interaction interfaces through their conservation patterns. Here, we propose a simple approach where the structural conservation of the charged residue pairs is analyzed among the pairs of homologous binary complexes. Specifically, we determine a large set of homologous interactions using an interaction interface similarity measure and catalog the basic types of conservation patterns among the charged residue pairs. We find an unexpected conservation pattern, which we call the correlated reappearance, occurring among the pairs of homologous interfaces more frequently than the fully conserved pairs of charged residues. Furthermore, the analysis of the conservation patterns across different superkingdoms as well as structural classes of proteins has revealed that the correlated reappearance of charged residues is by far the most prevalent conservation pattern, often occurring more frequently than the unconserved charged residues. We discuss a possible role that the new conservation pattern may play in the long-range electrostatic steering effect.     

Complementarity of stability patches at the interfaces of protein complexes: Implication for the structural organization of energetic hot spots

A rational design of protein complexes with defined functionalities and of drugs aimed at disrupting protein–protein interactions requires fundamental understanding of the mechanisms underlying the formation of specific protein complexes. Efforts to develop efficient small‐molecule or protein‐based binders often exploit energetic hot spots on protein surfaces, namely, the interfacial residues that provide most of the binding free energy in the complex. The molecular basis underlying the unusually high energy contribution of the hot spots remains obscure, and its elucidation would facilitate the design of interface‐targeted drugs. To study the nature of the energetic hot spots, we analyzed the backbone dynamic properties of contact surfaces in several protein complexes. We demonstrate that, in most complexes, the backbone dynamic landscapes of interacting surfaces form complementary “stability patches,” in which static areas from the opposing surfaces superimpose, and that these areas are predominantly located near the geometric center of the interface. We propose that a diminished enthalpy–entropy compensation effect augments the degree to which residues positioned within the complementary stability patches contribute to complex affinity, thereby giving rise to the energetic hot spots. These findings offer new insights into the nature of energetic hot spots and the role that backbone dynamics play in facilitating intermolecular recognition. Mapping the interfacial stability patches may provide guidance for protein engineering approaches aimed at improving the stability of protein complexes and could facilitate the design of ligands that target complex interfaces.     

PPIcons: identification of protein-protein interaction sites in selected organisms

The physico-chemical properties of interaction interfaces have a crucial role in characterization of protein–protein interactions (PPI). In silico prediction of participating amino acids helps to identify interface residues for further experimental verification using mutational analysis, or inhibition studies by screening library of ligands against given protein. Given the unbound structure of a protein and the fact that it forms a complex with another known protein, the objective of this work is to identify the residues that are involved in the interaction. We attempt to predict interaction sites in protein complexes using local composition of amino acids together with their physico-chemical characteristics. The local sequence segments (LSS) are dissected from the protein sequences using a sliding window of 21 amino acids. The list of LSSs is passed to the support vector machine (SVM) predictor, which identifies interacting residue pairs considering their inter-atom distances. We have analyzed three different model organisms of Escherichia coli, Saccharomyces Cerevisiae and Homo sapiens, where the numbers of considered hetero-complexes are equal to 40, 123 and 33 respectively. Moreover, the unified multi-organism PPI meta-predictor is also developed under the current work by combining the training databases of above organisms. The PPIcons interface residues prediction method is measured by the area under ROC curve (AUC) equal to 0.82, 0.75, 0.72 and 0.76 for the aforementioned organisms and the meta-predictor respectively.     

Computational structural analysis: multiple proteins bound to DNA

Background

With increasing numbers of crystal structures of protein∶DNA and protein∶protein∶DNA complexes publically available, it is now possible to extract sufficient structural, physical-chemical and thermodynamic parameters to make general observations and predictions about their interactions. In particular, the properties of macromolecular assemblies of multiple proteins bound to DNA have not previously been investigated in detail.

Methodology/Principal Findings

We have performed computational structural analyses on macromolecular assemblies of multiple proteins bound to DNA using a variety of different computational tools: PISA; PROMOTIF; X3DNA; ReadOut; DDNA and DCOMPLEX. Additionally, we have developed and employed an algorithm for approximate collision detection and overlapping volume estimation of two macromolecules. An implementation of this algorithm is available at http://promoterplot.fmi.ch/Collision1/. The results obtained are compared with structural, physical-chemical and thermodynamic parameters from protein∶protein and single protein∶DNA complexes. Many of interface properties of multiple protein∶DNA complexes were found to be very similar to those observed in binary protein∶DNA and protein∶protein complexes. However, the conformational change of the DNA upon protein binding is significantly higher when multiple proteins bind to it than is observed when single proteins bind. The water mediated contacts are less important (found in less quantity) between the interfaces of components in ternary (protein∶protein∶DNA) complexes than in those of binary complexes (protein∶protein and protein∶DNA).The thermodynamic stability of ternary complexes is also higher than in the binary interactions. Greater specificity and affinity of multiple proteins binding to DNA in comparison with binary protein-DNA interactions were observed. However, protein-protein binding affinities are stronger in complexes without the presence of DNA.

Conclusions/Significance

Our results indicate that the interface properties: interface area; number of interface residues/atoms and hydrogen bonds; and the distribution of interface residues, hydrogen bonds, van der Walls contacts and secondary structure motifs are independent of whether or not a protein is in a binary or ternary complex with DNA. However, changes in the shape of the DNA reduce the off-rate of the proteins which greatly enhances the stability and specificity of ternary complexes compared to binary ones.     

Transient complexes of redox proteins: structural and dynamic details from NMR studies

Redox proteins participate in many metabolic routes, in particular those related to energy conversion. Protein-protein complexes of redox proteins are characterized by a weak affinity and a short lifetime. Two-dimensional NMR spectroscopy has been applied to many redox protein complexes, providing a wealth of information about the process of complex formation, the nature of the interface and the dynamic properties of the complex. These studies have shown that some complexes are non-specific and exist as a dynamic ensemble of orientations while in other complexes the proteins assume a single orientation. The binding interface in these complexes consists of a small hydrophobic patch for specificity, surrounded by polar, uncharged residues that may enhance dissociation, and, in most complexes, a ring or patch of charged residues that enhances the association by electrostatic interactions. The entry and exit port of the electrons is located within the hydrophobic interaction site, ensuring rapid electron transfer from one redox centre to the next.     

Insights into protein-DNA interactions through structure network analysis

Protein-DNA interactions are crucial for many cellular processes. Now with the increased availability of structures of protein-DNA complexes, gaining deeper insights into the nature of protein-DNA interactions has become possible. Earlier, investigations have characterized the interface properties by considering pairwise interactions. However, the information communicated along the interfaces is rarely a pairwise phenomenon, and we feel that a global picture can be obtained by considering a protein-DNA complex as a network of noncovalently interacting systems. Furthermore, most of the earlier investigations have been carried out from the protein point of view (protein-centric), and the present network approach aims to combine both the protein-centric and the DNA-centric points of view. Part of the study involves the development of methodology to investigate protein-DNA graphs/networks with the development of key parameters. A network representation provides a holistic view of the interacting surface and has been reported here for the first time. The second part of the study involves the analyses of these graphs in terms of clusters of interacting residues and the identification of highly connected residues (hubs) along the protein-DNA interface. A predominance of deoxyribose-amino acid clusters in beta-sheet proteins, distinction of the interface clusters in helix-turn-helix, and the zipper-type proteins would not have been possible by conventional pairwise interaction analysis. Additionally, we propose a potential classification scheme for a set of protein-DNA complexes on the basis of the protein-DNA interface clusters. This provides a general idea of how the proteins interact with the different components of DNA in different complexes. Thus, we believe that the present graph-based method provides a deeper insight into the analysis of the protein-DNA recognition mechanisms by throwing more light on the nature and the specificity of these interactions.