This Article Full Text Full Text (PDF) Supplement All Versions of this Article: biophysj.107.112367v1 93/10/3340    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Brock, K. Articles by Alexov, E. Search for Related Content PubMed PubMed Citation Articles by Brock, K. Articles by Alexov, E.

Optimization of Electrostatic Interactions in Protein-Protein Complexes

Kelly Brock ^*, Kemper Talley ^*, Kacey Coley ^*, Petras Kundrotas  and Emil Alexov 

^* South Carolina Governor School for Science and Mathematics, Hartsville, South Carolina; and Computational Biophysics and Bioinformatics, Department of Physics, Clemson University, Clemson, South Carolina

Correspondence: Address reprint requests to E. Alexov, E-mail: ealexov{at}clemson.edu.

In this article, we present a statistical analysis of the electrostatic properties of 298 protein-protein complexes and 356 domain-domain structures extracted from the previously developed database of protein complexes (ProtCom, http://www.ces.clemson.edu/compbio/protcom). For each structure in the dataset we calculated the total electrostatic energy of the binding and its two components, Coulombic and reaction field energy. It was found that in a vast majority of the cases (>90%), the total electrostatic component of the binding energy was unfavorable. At the same time, the Coulombic component of the binding energy was found to favor the complex formation while the reaction field component of the binding energy opposed the binding. It was also demonstrated that the components in a wild-type (WT) structure are optimized/anti-optimized with respect to the corresponding distributions, arising from random shuffling of the charged side chains. The degree of this optimization was assessed through the Z-score of WT energy in respect to the random distribution. It was found that the Z-scores of Coulombic interactions peak at a considerably negative value for all 654 cases considered while the Z-score of the reaction field energy varied among different types of complexes. All these findings indicate that the Coulombic interactions within WT protein-protein complexes are optimized to favor the complex formation while the total electrostatic energy predominantly opposes the binding. This observation was used to discriminate WT structures among sets of structural decoys and showed that the electrostatic component of the binding energy is not a good discriminator of the WT; while, Coulombic or reaction field energies perform better depending upon the decoy set used.