Complementarity of structure ensembles in protein-protein binding

Protein-protein association is often accompanied by changes in receptor and ligand structure. This interplay between protein flexibility and protein-protein recognition is currently the largest obstacle both to our understanding of and to the reliable prediction of protein complexes. We performed two sets of molecular dynamics simulations for the unbound receptor and ligand structures of 17 protein complexes and applied shape-driven rigid body docking to all combinations of representative snapshots. The crossdocking of structure ensembles increased the likelihood of finding near-native solutions. The free ensembles appeared to contain multiple complementary conformations. These were in general not related to the bound structure. We suggest that protein-protein binding follows a three-step mechanism of diffusion, free conformer selection, and refolding. This model combines previously conflicting ideas and is in better agreement with the current data on interaction forces, time scales, and kinetics.     

High-throughput identification of interacting protein-protein binding sites

Background  

With the advent of increasing sequence and structural data, a number of methods have been proposed to locate putative protein binding sites from protein surfaces. Therefore, methods that are able to identify whether these binding sites interact are needed.     

The molecular architecture of protein-protein binding sites

The formation of specific protein interactions plays a crucial role in most, if not all, biological processes, including signal transduction, cell regulation, the immune response and others. Recent advances in our understanding of the molecular architecture of protein-protein binding sites, which facilitates such diversity in binding affinity and specificity, are enabling us to address key questions. What is the amino acid composition of binding sites? What are interface hotspots? How are binding sites organized? What are the differences between tight and weak interacting complexes? How does water contribute to binding? Can the knowledge gained be translated into protein design? And does a universal code for binding exist, or is it the architecture and chemistry of the interface that enable diverse but specific binding solutions?     

Flexibility and conformational entropy in protein-protein binding

To better understand the interplay between protein-protein binding and protein dynamics, we analyzed molecular dynamics simulations of 17 protein-protein complexes and their unbound components. Complex formation does not restrict the conformational freedom of the partner proteins as a whole, but, rather, it leads to a redistribution of dynamics. We calculate the change in conformational entropy for seven complexes with quasiharmonic analysis. We see significant loss, but also increased or unchanged conformational entropy. Where comparison is possible, the results are consistent with experimental data. However, stringent error estimates based on multiple independent simulations reveal large uncertainties that are usually overlooked. We observe substantial gains of pseudo entropy in individual partner proteins, and we observe that all complexes retain residual stabilizing intermolecular motions. Consequently, protein flexibility has an important influence on the thermodynamics of binding and may disfavor as well as favor association. These results support a recently proposed unified model for flexible protein-protein association.     

A structure-based benchmark for protein-protein binding affinity

We have assembled a nonredundant set of 144 protein-protein complexes that have high-resolution structures available for both the complexes and their unbound components, and for which dissociation constants have been measured by biophysical methods. The set is diverse in terms of the biological functions it represents, with complexes that involve G-proteins and receptor extracellular domains, as well as antigen/antibody, enzyme/inhibitor, and enzyme/substrate complexes. It is also diverse in terms of the partners' affinity for each other, with K(d) ranging between 10(-5) and 10(-14) M. Nine pairs of entries represent closely related complexes that have a similar structure, but a very different affinity, each pair comprising a cognate and a noncognate assembly. The unbound structures of the component proteins being available, conformation changes can be assessed. They are significant in most of the complexes, and large movements or disorder-to-order transitions are frequently observed. The set may be used to benchmark biophysical models aiming to relate affinity to structure in protein-protein interactions, taking into account the reactants and the conformation changes that accompany the association reaction, instead of just the final product.     

Predicting protein-protein binding sites in membrane proteins

Background  

Many integral membrane proteins, like their non-membrane counterparts, form either transient or permanent multi-subunit complexes in order to carry out their biochemical function. Computational methods that provide structural details of these interactions are needed since, despite their importance, relatively few structures of membrane protein complexes are available.     

Specific protein-protein binding in many-component mixtures of proteins

Proteins must bind to specific other proteins in vivo in order to function. The proteins are required to bind to only one or a few other proteins of the few thousand proteins typically present in vivo. To quantify this requirement we introduce a property of proteins called the capability. The capability is the maximum number of specific-binding interactions possible in a mixture, or in other words the size of largest sustainable interactome. This calculation of the maximum number possible is closely analogous to the work of Shannon and others on the maximum rate of communication through noisy channels. Using a simple model of proteins, we find specific binding to be a demanding function in the sense that it demands that the binding sites of the proteins be encoded by long sequences of elements, and the requirement for specific binding then strongly constrains these sequences.     

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.     

Zinc fingers: DNA binding and protein-protein interactions

The zinc finger domain is a very ubiquitous structural element whose hallmark is the coordination of a zinc atom by several amino acid residues (cysteines and histidines, and occasionally aspartate and glutamate). These structural elements are associated with protein-nucleic acid recognition as well as protein-protein interactions. The purpose of this review is to examine recent data on the DNA and protein binding properties of a few zinc fingers whose three dimensional structure is known.     

Understanding mechanisms governing protein-protein interactions from synthetic binding interfaces

Recent advances in methodologies and design of combinatorial library selection have enabled comprehensive characterization of sequence space for protein-protein interaction interfaces and generation of fully synthetic binding interfaces. By exhaustively introducing and quantitatively analyzing mutations in natural interfaces, new insights into their molecular architecture and plasticity have emerged. Minimalist combinatorial libraries based on a restricted amino acid code have produced synthetic interfaces that rival natural ones using a different set of rules. A two amino acid code composed of just tyrosine and serine in the context of antibody CDR loops is sufficient to produce high affinity and specific interactions with different classes of protein targets. Structural analyses highlight the dominant role of Tyr in forming productive interactions and demonstrate the dominance of conformational diversity over chemical diversity in producing na?ve binding interfaces. Synthetic binding proteins are beginning to be used as a powerful crystallization tool to attack important structural biology problems that are recalcitrant to crystallization using traditional methods.     

Kinetics of protein-protein interactions of connexins: use of enzyme linked sorbent assays

Determination of the protein-protein interactions of connexins has become a rapidly expanding field of research. While there are multiple methods of determining the identity of binding partners, determination of the strengths of interactions is not as simple. Here we describe the use of the in vitro method Enzyme Linked Sorbent Assay (ELSA) to compare binding affinities of known protein partners for Connexin43. We used the binding of Cx43 Carboxyl Terminal domain to the PDZ-2 domain of Zonula Occludens-1 and to the SH3 domain of c-Src. In the ELSA assay we found that while the binding of the SH3 domain of c-Src is pH-dependent, the interaction of the PDZ domain of ZO-1 is not. These data confirm findings using Surface Plasmon Resonance (1) and indicate that ELSA can be a useful tool in determining the kinetics of protein-protein interactions.     

Kinetics of heme binding to semi-alpha-hemoglobin

The binding of carbonmonoxyheme to semi-alpha-hemoglobin and to an apohemoglobin control was investigated using stopped-flow techniques in 0.025 M potassium phosphate buffer, pH 7 and 10 degrees C. The resultant second order kinetic data were analyzed by the classical model which assumes the existence of an intermediate complex which either redissociates to reactants or undergoes an irreversible conversion to form hemoglobin. The rate constants for the latter unimolecular process were apparently not experimentally different for semi-alpha-hemoglobin and apohemoglobin (360 ( +/- 100) s-1 and 480 ( +/- 60) s-1, respectively). However, the equilibrium dissociation constant for the intermediate of semi-alpha-hemoglobin (Kd = 9.3 ( +/- 2.6) micromolar) was approximately two fold greater than that of apohemoglobin (Kd = 4.1 ( +/- 0.5) micromolar). The reduced stability of the semi-alpha-hemoglobin complex was postulated to be due to the lower affinity of the beta pocket for heme. The studies reported here address the possible role of semi-alpha-hemoglobin as an intermediate in the assembly of hemoglobin in vivo.     

Kinetics of calmodulin binding to calcineurin

Calcineurin (CaN) binds Ca(2+)-saturated calmodulin (CaM) with relatively high affinity; however, an accurate steady-state K(d) value has not been determined. In this report, we describe, using steady-state and stopped-flow fluorescence techniques, the rates of association and dissociation of Ca(2+)-saturated CaM from CaN heterodimer (CaNA/CaNB) and CaNA only. The rate of Ca(2+)/CaM association was determined to be 4.6 x 10(7) M(-1)s(-1). The rate of Ca(2+)/CaM dissociation from CaN was slower than previously reported and was approximately 0.0012 s(-1). In preparations of CaNA alone (no regulatory CaNB subunit), the dissociation rate was slowed further to 0.00026 s(-1). From these data we calculate a K(d) for binding of Ca(2+)-saturated CaM to CaN of 28 pM. This K(d) is significantly lower than previously reported estimates of approximately 1 nM and indicates that CaN is one of the highest affinity CaM-binding proteins identified to date.     

Kinetics of ligand binding to haemoproteins

A mass of experimental data has been accumulated in the 65 years since Hartridge and Roughton made the first measurement of the rapid reaction of haemoglobin with O2 in solution on a millisecond time scale, at first by flow-mixing methods, and, for 30 years or so, by flash photolysis. Technical advances, particularly in lasers, have allowed increasingly rapid reactions to be followed and the fastest reactions now observed have half-times conveniently measured in pico-seconds. The measurements were used at first to discuss the physiology of gas transport and to describe co-operativity in haemoglobin. More recently, the process of ligand binding has been dissected into intramolecular and intermolecular components. Relating the various rates to the abundance of structural information on crystals is so difficult that the work has barely begun, but the combination of kinetic measurements with genetic engineering and crystallography has promise, as well as problems, for the future.     

In vivo protein-protein and protein-DNA crosslinking for genomewide binding microarray

Chromatin immunoprecipitation (ChrIP or ChIP) has commonly been used to map protein-DNA interaction sites at specific genomic loci through use of formaldehyde-induced crosslinking. However, formaldehyde alone has proved inadequate for crosslinking of certain proteins such as the yeast histone deacetylase Rpd3. We report here a modified crosslinking procedure that includes a protein-protein crosslinking agent in addition to formaldehyde. Using this double crosslinking method, we have successfully mapped Rpd3 binding sites in vivo. We also describe the use of ChrIP in combination with DNA microarrays (ChrIP-array) to determine the pattern of Rpd3 binding genomewide. This approach couples the versatility of ChrIP with that of microarrays to identify binding patterns that would otherwise be hidden in a gene-by-gene survey.     

Phosphorylation in protein-protein binding: effect on stability and function

Posttranslational modifications offer a dynamic way to regulate protein activity, subcellular localization, and stability. Here we estimate the effect of phosphorylation on protein binding and function for different types of complexes from human proteome. We find that phosphorylation sites tend to be located on binding interfaces in heterooligomeric and weak transient homooligomeric complexes. Analysis of molecular mechanisms of phosphorylation shows that phosphorylation may modulate the strength of interactions directly on interfaces and that binding hotspots tend to be phosphorylated in heterooligomers. Although the majority of complexes do not show significant estimated stability differences upon phosphorylation or dephosphorylation, for about one-third of all complexes it causes relatively large changes in binding energy. We discuss the cases where phosphorylation mediates the complex formation and regulates the function. We show that phosphorylation sites are more likely to be evolutionary conserved than other interfacial residues.     

Poisson-Boltzmann calculations of nonspecific salt effects on protein-protein binding free energies

The salt dependence of the binding free energy of five protein-protein hetero-dimers and two homo-dimers/tetramers was calculated from numerical solutions to the Poisson-Boltzmann equation. Overall, the agreement with experimental values is very good. In all cases except one involving the highly charged lactoglobulin homo-dimer, increasing the salt concentration is found both experimentally and theoretically to decrease the binding affinity. To clarify the source of salt effects, the salt-dependent free energy of binding is partitioned into screening terms and to self-energy terms that involve the interaction of the charge distribution of a monomer with its own ion atmosphere. In six of the seven complexes studied, screening makes the largest contribution but self-energy effects can also be significant. The calculated salt effects are found to be insensitive to force-field parameters and to the internal dielectric constant assigned to the monomers. Nonlinearities due to high charge densities, which are extremely important in the binding of proteins to negatively charged membrane surfaces and to nucleic acids, make much smaller contributions to the protein-protein complexes studied here, with the exception of highly charged lactoglobulin dimers. Our results indicate that the Poisson-Boltzmann equation captures much of the physical basis of the nonspecific salt dependence of protein-protein complexation.     

The modular organization of domain structures: insights into protein-protein binding

Domains are the building blocks of proteins and play a crucial role in protein-protein interactions. Here, we propose a new approach for the analysis and prediction of domain-domain interfaces. Our method, which relies on the representation of domains as residue-interacting networks, finds an optimal decomposition of domain structures into modules. The resulting modules comprise highly cooperative residues, which exhibit few connections with other modules. We found that non-overlapping binding sites in a domain, involved in different domain-domain interactions, are generally contained in different modules. This observation indicates that our modular decomposition is able to separate protein domains into regions with specialized functions. Our results show that modules with high modularity values identify binding site regions, demonstrating the predictive character of modularity. Furthermore, the combination of modularity with other characteristics, such as sequence conservation or surface patches, was found to improve our predictions. In an attempt to give a physical interpretation to the modular architecture of domains, we analyzed in detail six examples of protein domains with available experimental binding data. The modular configuration of the TEM1-beta-lactamase binding site illustrates the energetic independence of hotspots located in different modules and the cooperativity of those sited within the same modules. The energetic and structural cooperativity between intramodular residues is also clearly shown in the example of the chymotrypsin inhibitor, where non-binding site residues have a synergistic effect on binding. Interestingly, the binding site of the T cell receptor beta chain variable domain 2.1 is contained in one module, which includes structurally distant hot regions displaying positive cooperativity. These findings support the idea that modules possess certain functional and energetic independence. A modular organization of binding sites confers robustness and flexibility to the performance of the functional activity, and facilitates the evolution of protein interactions.     

RPA phosphorylation in mitosis alters DNA binding and protein-protein interactions

The heterotrimeric DNA-binding protein, replication protein A (RPA), consists of 70-, 34-, and 14-kDa subunits and is involved in maintaining genomic stability by playing key roles in DNA replication, repair, and recombination. RPA participates in these processes through its interaction with other proteins and its strong affinity for single-stranded DNA (ssDNA). RPA-p34 is phosphorylated in a cell-cycle-dependent fashion primarily at Ser-29 and Ser-23, which are consensus sites for Cdc2 cyclin-dependent kinase. By systematically examining RPA-p34 phosphorylation throughout the cell cycle, we have found there are distinct phosphorylated forms of RPA-p34 in different cell-cycle stages. We have isolated and purified a unique phosphorylated form of RPA that is specifically associated with the mitotic phase of the cell cycle. The mitotic form of RPA (m-hRPA) shows no difference in ssDNA binding activity as compared with recombinant RPA (r-hRPA), yet binds less efficiently to double-stranded DNA (dsDNA). These data suggest that mitotic phosphorylation of RPA-p34 inhibits the destabilization of dsDNA by RPA complex, thereby decreasing the binding affinity for dsDNA. The m-hRPA also exhibits altered interactions with certain DNA replication and repair proteins. Using highly purified proteins, m-hRPA exhibited decreased binding to ATM, DNA pol alpha, and DNA-PK as compared to unphosphorylated recombinant RPA (r-hRPA). Dephosphorylation of m-hRPA was able to restore the interaction with each of these proteins. Interestingly, the interaction of RPA with XPA was not altered by RPA phosphorylation. These data suggest that phosphorylation of RPA-p34 plays an important role in regulating RPA functions in DNA metabolism by altering specific protein-protein interactions.