Prediction of protein-protein binding free energies

We present an energy function for predicting binding free energies of protein-protein complexes, using the three-dimensional structures of the complex and unbound proteins as input. Our function is a linear combination of nine terms and achieves a correlation coefficient of 0.63 with experimental measurements when tested on a benchmark of 144 complexes using leave-one-out cross validation. Although we systematically tested both atomic and residue-based scoring functions, the selected function is dominated by residue-based terms. Our function is stable for subsets of the benchmark stratified by experimental pH and extent of conformational change upon complex formation, with correlation coefficients ranging from 0.61 to 0.66.     

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.     

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?     

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.     

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.     

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.     

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.     

In silico modeling of pH-optimum of protein-protein binding

Protein-protein association is a pH-dependent process and thus the binding affinity depends on the local pH. In vivo the association occurs in a particular cellular compartment, where the individual monomers are supposed to meet and form a complex. Since the monomers and the complex exist in the same micro environment, it is plausible that they coevolved toward its properties, in particular, toward the characteristic subcellular pH. Here we show that the pH at which the monomers are most stable (pH-optimum) or the pH at which stability is almost pH-independent (pH-flat) of monomers are correlated with the pH-optimum of maximal affinity (pH-optimum of binding) or pH interval at which affinity is almost pH-independent (pH-flat of binding) of the complexes made of the corresponding monomers. The analysis of interfacial properties of protein complexes demonstrates that pH-dependent properties can be roughly estimated using the interface charge alone. In addition, we introduce a parameter beta, proportional to the square root of the absolute product of the net charges of monomers, and show that protein complexes characterized with small or very large beta tend to have neutral pH-optimum. Further more, protein complexes made of monomers carrying the same polarity net charge at neutral pH have either very low or very high pH-optimum of binding. These findings are used to propose empirical rule for predicting pH-optimum of binding provided that the amino acid compositions of the corresponding monomers 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.     

Electrostatic contribution to the binding stability of protein-protein complexes

To investigate roles of electrostatic interactions in protein binding stability, electrostatic calculations were carried out on a set of 64 mutations over six protein-protein complexes. These mutations alter polar interactions across the interface and were selected for putative dominance of electrostatic contributions to the binding stability. Three protocols of implementing the Poisson-Boltzmann model were tested. In vdW4 the dielectric boundary between the protein low dielectric and the solvent high dielectric is defined as the protein van der Waals surface and the protein dielectric constant is set to 4. In SE4 and SE20, the dielectric boundary is defined as the surface of the protein interior inaccessible to a 1.4-A solvent probe, and the protein dielectric constant is set to 4 and 20, respectively. In line with earlier studies on the barnase-barstar complex, the vdW4 results on the large set of mutations showed the closest agreement with experimental data. The agreement between vdW4 and experiment supports the contention of dominant electrostatic contributions for the mutations, but their differences also suggest van der Waals and hydrophobic contributions. The results presented here will serve as a guide for future refinement in electrostatic calculation and inclusion of nonelectrostatic effects.     

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.     

Discovery of protein-protein binding disruptors using multi-component condensations small molecules

A series of small molecule compounds interfering with the binding process of VEGF and NRP1 has been identified and further optimized. Full synthetic details as well as SAR are reported which demonstrate that expeditious MCC-based syntheses may lead to valuable molecules addressing challenging targets such as protein-protein interactions. Preliminary functional assay data confirm that these compounds may be further developed toward drug candidates.     

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.     

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.     

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.     

Exploiting sequence and structure homologs to identify protein-protein binding sites

A rapid increase in the number of experimentally derived three-dimensional structures provides an opportunity to better understand and subsequently predict protein-protein interactions. In this study, structurally conserved residues were derived from multiple structure alignments of the individual components of known complexes and the assigned conservation score was weighted based on the crystallographic B factor to account for the structural flexibility that will result in a poor alignment. Sequence profile and accessible surface area information was then combined with the conservation score to predict protein-protein binding sites using a Support Vector Machine (SVM). The incorporation of the conservation score significantly improved the performance of the SVM. About 52% of the binding sites were precisely predicted (greater than 70% of the residues in the site were identified); 77% of the binding sites were correctly predicted (greater than 50% of the residues in the site were identified), and 21% of the binding sites were partially covered by the predicted residues (some residues were identified). The results support the hypothesis that in many cases protein interfaces require some residues to provide rigidity to minimize the entropic cost upon complex formation.     

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.     

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.