学术文献库

Understanding the molecular mechanisms driving the binding between bio-molecules is a crucial challenge in molecular biology. In this respect, characteristics like the preferentially hydrophobic composition of the binding interfaces, the role of van der Waals interactions (short range forces), and the consequent shape complementarity between the interacting molecular surfaces are well established. However, no consensus has yet been reached on how and how much electrostatic participates in the various stages of protein-protein interactions. Here, we perform extensive analyses on a large dataset of protein complexes for which both experimental binding affinity and pH data were available. We found that (i) although different classes of dimers do not present marked differences in the amino acid composition and charges disposition in the binding region, (ii) homodimers with identical binding region show higher electrostatic compatibility with respect to both homodimers with non-identical binding region and heterodimers. The level of electrostatic compatibility also varies with the pH of the complex, reaching the lowest values for low pH. Interestingly, (iii) shape and electrostatic complementarity behave oppositely when one stratifies the complexes by their binding affinity. Conversely, complexes with low values of binding affinity exploit Coulombic complementarity to acquire specificity, suggesting that electrostatic complementarity may play a greater role in transient (or less stable) complexes. In light of these results, (iv) we provide a fast and efficient method to measure electrostatic complementarity without the need of knowing the complex structure. Expanding the electrostatic potential on a basis of 2D orthogonal polynomials, we can discriminate between transient and permanent protein complexes with an AUC of the ROC of 0.8.