Kinetics of desolvation-mediated protein-protein binding

The role of desolvation in protein binding kinetics is investigated using Brownian dynamics simulations in complexes in which the electrostatic interactions are relatively weak. We find that partial desolvation, modeled by a short-range atomic contact potential, is not only a major contributor to the binding free energy but also substantially increases the diffusion-limited rate for complexes in which long-range electrostatics is weak. This rate enhancement is mostly due to weakly specific pathways leading to a low free-energy attractor, i.e., a precursor state before docking. For alpha-chymotrypsin and human leukocyte elastase, both interacting with turkey ovomucoid third domain, we find that the forward rate constant associated with a collision within a solid angle phi around their corresponding attractor approaches 10(7) and 10(6) M(-1)s(-1), respectively, in the limit phi approximately 2 degrees. Because these estimates agree well with experiments, we conclude that the final bound conformation must be preceded by a small set of well-defined diffusion-accessible precursor states. The inclusion of the otherwise repulsive desolvation interaction also explains the lack of aggregation in proteins by restricting nonspecific association times to approximately 4 ns. Under the same reaction conditions but without short range forces, the association rate would be only approximately 10(3) M(-1)s(-1). Although desolvation increases these rates by three orders of magnitude, desolvation-mediated association is still at least 100-fold slower than the electrostatically assisted binding in complexes such as barnase and barstar.     

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.     

The role of weak protein-protein interactions in multivalent lectin-carbohydrate binding: crystal structure of cross-linked FRIL

Binding of multivalent glycoconjugates by lectins often leads to the formation of cross-linked complexes. Type I cross-links, which are one-dimensional, are formed by a divalent lectin and a divalent glycoconjugate. Type II cross-links, which are two or three-dimensional, occur when a lectin or glycoconjugate has a valence greater than two. Type II complexes are a source of additional specificity, since homogeneous type II complexes are formed in the presence of mixtures of lectins and glycoconjugates. This additional specificity is thought to become important when a lectin interacts with clusters of glycoconjugates, e.g. as is present on the cell surface. The cryst1al structure of the Glc/Man binding legume lectin FRIL in complex with a trisaccharide provides a molecular snapshot of how weak protein-protein interactions, which are not observed in solution, can become important when a cross-linked complex is formed. In solution, FRIL is a divalent dimer, but in the crystal FRIL forms a tetramer, which allows for the formation of an intricate type II cross-linked complex with the divalent trisaccharide. The dependence on weak protein-protein interactions can ensure that a specific type II cross-linked complex and its associated specificity can occur only under stringent conditions, which explains why lectins are often found forming higher-order oligomers.     

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.     

Secondary structure based analysis and classification of biological interfaces: identification of binding motifs in protein-protein interactions

MOTIVATION: The increasing amount of data on protein-protein interaction needs to be rationalized for deriving guidelines for the alteration or design of an interface between two proteins. RESULTS: We present a detaild structural analysis and comparison of homo- versus heterodimeric protein-protein interfaces. Regular secondary structures (helices and strands) are the main components of the former, whereas non-regular structures (turns, loops, etc.) frequently mediate interactions in the latter. Interface helices get longer with increasing interface area, but only in heterocomplexes. On average, the homodimers have longer helical segments and prominent helix-helix pairs. There is a surprising distinction in the relative orientation of interface helices, with a tendency for aligned packing in homodimers and a clear preference for packing at 90 degrees in heterodimers. Arg and the aromatic residues have a higher preference to occur in all secondary structural elements (SSEs) in the interface. Based on the dominant SSE, the interfaces have been grouped into four classes: alpha, beta, alphabeta and non-regular. Identity between protein and interface classes is the maximum for alpha proteins, but rather mediocre for the other protein classes. The interface classes of the two chains forming a heterodimer are often dissimilar. Eleven binding motifs can capture the prominent architectural features of most of the interfaces.     

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?     

ProMate: a structure based prediction program to identify the location of protein-protein binding sites

Is the whole protein surface available for interaction with other proteins, or are specific sites pre-assigned according to their biophysical and structural character? And if so, is it possible to predict the location of the binding site from the surface properties? These questions are answered quantitatively by probing the surfaces of proteins using spheres of radius of 10 A on a database (DB) of 57 unique, non-homologous proteins involved in heteromeric, transient protein-protein interactions for which the structures of both the unbound and bound states were determined. In structural terms, we found the binding site to have a preference for beta-sheets and for relatively long non-structured chains, but not for alpha-helices. Chemically, aromatic side-chains show a clear preference for binding sites. While the hydrophobic and polar content of the interface is similar to the rest of the surface, hydrophobic and polar residues tend to cluster in interfaces. In the crystal, the binding site has more bound water molecules surrounding it, and a lower B-factor already in the unbound protein. The same biophysical properties were found to hold for the unbound and bound DBs. All the significant interface properties were combined into ProMate, an interface prediction program. This was followed by an optimization step to choose the best combination of properties, as many of them are correlated. During optimization and prediction, the tested proteins were not used for data collection, to avoid over-fitting. The prediction algorithm is fully automated, and is used to predict the location of potential binding sites on unbound proteins with known structures. The algorithm is able to successfully predict the location of the interface for about 70% of the proteins. The success rate of the predictor was equal whether applied on the unbound DB or on the disjoint bound DB. A prediction is assumed correct if over half of the predicted continuous interface patch is indeed interface. The ability to predict the location of protein-protein interfaces has far reaching implications both towards our understanding of specificity and kinetics of binding, as well as in assisting in the analysis of the proteome.     

Conformer selection and induced fit in flexible backbone protein-protein docking using computational and NMR ensembles

Accommodating backbone flexibility continues to be the most difficult challenge in computational docking of protein-protein complexes. Towards that end, we simulate four distinct biophysical models of protein binding in RosettaDock, a multiscale Monte-Carlo-based algorithm that uses a quasi-kinetic search process to emulate the diffusional encounter of two proteins and to identify low-energy complexes. The four binding models are as follows: (1) key-lock (KL) model, using rigid-backbone docking; (2) conformer selection (CS) model, using a novel ensemble docking algorithm; (3) induced fit (IF) model, using energy-gradient-based backbone minimization; and (4) combined conformer selection/induced fit (CS/IF) model. Backbone flexibility was limited to the smaller partner of the complex, structural ensembles were generated using Rosetta refinement methods, and docking consisted of local perturbations around the complexed conformation using unbound component crystal structures for a set of 21 target complexes. The lowest-energy structure contained > 30% of the native residue-residue contacts for 9, 13, 13, and 14 targets for KL, CS, IF, and CS/IF docking, respectively. When applied to 15 targets using nuclear magnetic resonance ensembles of the smaller protein, the lowest-energy structure recovered at least 30% native residue contacts in 3, 8, 4, and 8 targets for KL, CS, IF, and CS/IF docking, respectively. CS/IF docking of the nuclear magnetic resonance ensemble performed equally well or better than KL docking with the unbound crystal structure in 10 of 15 cases. The marked success of CS and CS/IF docking shows that ensemble docking can be a versatile and effective method for accommodating conformational plasticity in docking and serves as a demonstration for the CS theory—that binding-competent conformers exist in the unbound ensemble and can be selected based on their favorable binding energies.     

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.     

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.     

Selection and flexible optimization of binding modes from conformation ensembles

This paper describes a new semi-flexible docking approach named Fado (flexible alignment and docking), which incorporates flexibility by using an ensemble of precomputed ligand conformers. A primary ligand is defined as a linear combination over all input conformers. An optimization with regard to the linear coefficients makes the ligand flexible. Initially, a point matching problem utilizing the Merck Molecular Force Field (MMFF) is modeled in order to compute the correct orientation of the ligand with respect to the target. The problem is then solved through a local optimization approach (RPROP). This is done for 20 randomized ligand orientations, yielding 20 binding modes per complex. Evaluating these modes illustrates that our method is able to reproduce the binding modes of molecules within a few minutes of CPU time. A representative dataset of diverse protein-ligand complexes could be reproduced with 78% accuracy below 2A RMSD distance to the reference crystal structure. Fado is available upon request to the authors (see also http://www.zib.eu/Numerik/projects/docking/projectlong.en.html).     

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.     

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.     

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.     

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.