On the electrostatic component of protein-protein binding free energy
© Talley et al; 2008
Received: 25 July 2008
Accepted: 05 November 2008
Published: 05 November 2008
Calculations of electrostatic properties of protein-protein complexes are usually done within framework of a model with a certain set of parameters. In this paper we present a comprehensive statistical analysis of the sensitivity of the electrostatic component of binding free energy (ΔΔGel) with respect with different force fields (Charmm, Amber, and OPLS), different values of the internal dielectric constant, and different presentations of molecular surface (different values of the probe radius). The study was done using the largest so far set of entries comprising 260 hetero and 2148 homo protein-protein complexes extracted from a previously developed database of protein complexes (ProtCom). To test the sensitivity of the energy calculations with respect to the structural details, all structures were energy minimized with corresponding force field, and the energies were recalculated. The results indicate that the absolute value of the electrostatic component of the binding free energy (ΔΔGel) is very sensitive to the force field parameters, the minimization procedure, the values of the internal dielectric constant, and the probe radius. Nevertheless our results indicate that certain trends in ΔΔGel behavior are much less sensitive to the calculation parameters. For instance, the fraction of the homo-complexes, for which the electrostatics was found to oppose binding, is 80% regardless of the force fields and parameters used. For the hetero-complexes, however, the percentage of the cases for which electrostatics opposed binding varied from 43% to 85%, depending on the protocol and parameters employed. A significant correlation was found between the effects caused by raising the internal dielectric constant and decreasing the probe radius. Correlations were also found among the results obtained with different force fields. However, despite of the correlations found, the absolute ΔΔGel calculated with different force field parameters could differ more than tens of kcal/mol in some cases. Set of rules of obtaining confident predictions of absolute ΔΔGel and ΔΔGel sign are provided in the conclusion section.
PACS codes: 87.15.A-, 87.15. km
Proteins are essential components of the living cell. They function by interacting with other biological macromolecules, including other proteins. Therefore, understanding protein-protein interactions is a very important step in learning about cell function, which in turn requires understanding the forces and effects that influence these interactions [2–10]. Protein-protein binding is a complex process mostly driven by the hydrophobic effect and van der Waals interactions, with significant contribution from entropy and electrostatics. While these energies can not be easy experimentally separated, the electrostatic energy is the energy mostly affected by pH and salt concentration. The interest of modeling pH and salt dependent phenomena inspired many researchers to model the electrostatic component of the binding free energy alone . Nowadays, the 3D structures of a significant number of protein-protein complexes are available, which allows researchers to model the contribution of electrostatic energy to the binding [12–19]. Experimental [20–23] and computational mutagenesis studies have found that most of the ionizable residues at protein-protein interfaces contribute significantly to the binding free energy (i.e., replacement of those residues with the alanine residues critically affects the protein binding affinity). In many cases, the formation of a complex creates favorable pair-wise electrostatic interactions across the interface, as was demonstrated in a series of work on the barnase-barstar complex[17, 25–28] and on other complexes[19, 29–33]. Authors of several studies[7, 34, 35] have concluded that electrostatic interactions have a dominant contribution into total binding free energy for the complexes with small interfaces. The contribution of electrostatic energy to the binding affinity of the particular, Rap/Raf, complex was the subject of a series of investigations[36, 37]. Despite the fact that all of the abovementioned studies agreed that there are many specific pair-wise electrostatic interactions across the interface, the conclusions about the role of electrostatics on binding affinity remain controversial. It was found that, in some cases, electrostatic energy favors binding, while in other cases, it opposes it. Since the electrostatic component of the binding free energy is the difference between two large and similar numbers (energy of pair-wise interactions and the desolvation penalty), the calculations outcome inevitably depends on parameters of simulations such as the force field parameters, the choice of the internal dielectric constant of proteins, and the molecular surface representation. In addition, structural refinement also affects the outcome of electrostatic calculations.
The internal dielectric constant is an important characteristic of protein's polarizability [38–41]. However, it was repeatedly shown in the literature that the dielectric constant is neither a constant nor is there an optimal value that works for all cases. Molecular dynamic simulations with implicit solvent models use a dielectric constant of either 1.0 or 2.0, the molecular mechanics Poisson-Boltzmann (MMPB) method usually employs a dielectric constant of 2.0[43, 44], and pKa's calculations are usually carried out with a dielectric constant larger or equal to 4.0[45, 46]. In our study the dielectric constant will be considered a parameter, and we will vary its value to reveal the sensitivity of the calculations of the electrostatic component of binding. The other parameter that will be varied is the probe radius. As pointed out by Zhou and co-workers[17, 47, 48], the electrostatic component of the binding free energy is very sensitive to the method used to build the molecular surface. Recently, the effect of different presentations of the molecular surface was tested on a set of 64 mutants, which demonstrated that calculations done with van der Waals molecular surface give better agreement with experimental data delivered from a double mutant cycle. The charge distribution and atomic radii greatly affect numerical calculations as well. Currently, there are many force fields, each having different partial charge distributions and sets of atomic radii. Three of the most popular force fields, Charmm, Amber, and OPLS will be used in this study to reveal the sensitivity of the calculations with respect to the force field parameters. Previous work has shown that optimal hydrogen bonding is crucial for correct predictions of pKa's of ionizable groups[50, 51]. However, different force fields generate different hydrogen networks. In addition, the calculated energies strongly depend on structural details, and many studies have been done on refined structures, i.e., on structures that were subjected to energy minimization. Therefore, it is desirable to know what will be the effect of minimization on the calculations of the electrostatic component of binding, if minimization is done with different force fields.
In this study we test the sensitivity of the electrostatic component of the binding free energy with respect to different force field parameters, structural relaxation, values of the internal dielectric constant, and probe radii. The study will be done on a large set of protein-protein complexes comprised of 260 hetero- and 2148 homo-complexes. Such a large dataset has never been investigated before, especially with energy minimized structures.
Protein-protein complexes used in the study and their parameters
Protein-protein hetero-complexes subjected to the study were extracted from the ProtCom database (as of November 2006) http://www.ces.clemson.edu/compbio/protcom, which contains 1771 entries at 95% sequence identity level. To avoid the bias toward overrepresented complexes, the entries were purged with CD-hit at 40% sequence identity level for all components of the hetero-complexes, including monomers that belong to the same protein-protein complex. This resulted in 299 structures out of which 39 structures were excluded from further consideration due to large defects in the PDB files (like large segments of missing polypeptide chains, for example), which constitute untreatable obstacles for the correct protonation of the structures. The structures of homo-complexes were taken from 40% sequence identity of the ProtCom database. Some structures were excluded because of large structural defects, which reduced our homo-complex set to 2617 structures at 40% sequence identity level.
This resulted in an initial set of 260 hetero- and 2617 homo-protein complexes, which after removing structures for which some of the computational procedures could not produce the desired accuracy was reduced to 260 hetero- and 2148 homo-protein complexes. The interfacial area was calculated with the surfv program by subtracting the accessible surface area of the complex from the accessible surface area of the free monomers. The net charge of the complex and the free monomers were calculated assuming that all titratable amino acids were in their charged state.
Hydrogen placement and energy minimization
Some of the structures in our initial data set had structural defects, and thus all of the structures were subjected to the profix program from the Jackal package developed in Honig's lab (http://wiki.c2b2.columbia.edu/honiglab_public/index.php/Software: Jackal) in order to add missing atoms and/or sequence fragments. All of the complexes that had a missing segment chain longer than fifteen amino acids were removed from the set because building such long sequence stretches could introduce significant structural errors. The rest of the structures were subjected to the TINKER software to add missing hydrogens using the package's pdbxyz.x and the xyzpdb.z modules with three different force field parameters: Charmm27, AMBER, and OPLS. We will refer to this set of structures as the non-minimized set.
A common approach in computing biophysical quantities using the 3D structures of biological macromolecules is to refine these structures by performing energy minimization with a particular force field. Following this strategy, we created another set of protonated structures (hereafter referred to as the minimized set) for all of the PDB files described above by running the minimize.x module of the TINKER software between running the pdbxyz.x and the xyzpdb.z modules. The minimize.x module performs energy minimization using the Limited Memory BFGS Quasi-Newton Optimization algorithm. The implicit solvent was modeled using the Still Generalized Born model, and the internal dielectric constant was set to 1.0 to be consistent with the corresponding force field parameters used in the calculations. To test the sensitivity of the results, the minimizations were done with Charmm27, AMBER, and OPLS parameters. A weak convergence criteria (RMS gradient per atom = 0.1) was applied to make computation tractable. Even in such a case, minimizing 2408 protein-protein complexes, some of which were larger than 50,000 atoms, was quite a challenge from a computational point of view (monomers were not minimized separately and their structures were taken from minimized complexes). Therefore, for these calculations we utilized a high throughput distributed computing resource, CONDOR, originally developed at the University of Wisconsin-Madison http://www.cs.wisc.edu/condor, which is available at Clemson University with more than 2,000 single CPUs of computational power.
Calculation of the electrostatic component of the binding free energy
where ΔG el (X)is the electrostatic component of the folding energy of the complex, monomer A and monomer B, respectively (X = AB, A or B). The 3D structures of the monomers were taken from the corresponding complex, i.e. no conformational changes associated with the binding were modeled. Since the 3D structures of the bound and unbound monomers are the same, their internal mechanical energies are the same as well and do not affect the calculations. These energies were calculated as a sum of Coulombic interactions and reaction field energies, as discussed in detail in Ref.. The Coulombic energy was calculated in homogeneous media with a dielectric constant equal to the internal dielectric constant of the protein. The reaction field energy was calculated as the energy of the interaction of permanent charges with induced charges at the surface of the corresponding entity – protein complex or a monomer. These energy terms were calculated with the Delphi[64, 65] program using a grid size of 129. The calculations were performed as zero ionic strength. The external dielectric constant was kept at 80. The convergence criteria was rmsc = 0.0001 kT/e (see the Delphi manual for details at http://wiki.c2b2.columbia.edu/honiglab_public/index.php/Software: Delphi).
The calculations were performed assuming that all of the Arg, Asp, Glu, and Lys residues were ionized in both free and bound states, while the His residues were considered neutral. In order to reduce the complexity of the problem, all ionizable residues were kept in their default charge state. Thus, ionization changes induced by complex formation (suggested either by experimental data[66, 67] or by theoretical simulations [68–71]) were not considered.
Parameters that will be varied
Force field parameters
Internal dielectric constant
What the value of the internal dielectric constant (ε(in)) should be for a given protein is the subject of much debate in the scientific community. In this study we adopt a pragmatic approach and test the sensitivity of the results at ε(in) = 1, 2, 4, 8, 20.
The molecular surface in FDPB algorithms is determined by rolling a water probe with a particular radius to find the molecular surface according to Richard's algorithm. A typical probe radius is 1.4A. However, Zhou and co-workers have introduced the concept of a zero probe radius, and here we will explore the effect of different probe radii on the electrostatic component of the binding free energy. The probe radius will vary taking values of R = 0.0, 0.5, 1.0, and 1.4A.
Macroscopic parameters of hetero- and homo-complexes
Max interfacial area (A**2)
Min interfacial area (A**2)
Mean interfacial area (A**2)
Max interfacial residues
Min interfacial residues
Mean interfacial residues
Max/min net charge
Max/min interfacial charge
Percentage having opposite net charge of the monomers
Percentage having opposite interfacial change of the monomers (in parentheses are the cases for which at least one of the interfaces had zero net charge)
58% (23% zero)
8% (33% zero)
Effects of the internal dielectric constant
Why were the calculated distributions of ΔΔGel using the low internal dielectric constant quite similar for the hetero- and homo-complexes but significantly different using high values of the internal dielectric constant? A possible explanation is offered by the numbers reported in the last two rows of Table 1. From the results in Table 1, it can be concluded that in the vast majority of the cases, the homo-complexes form interfaces that carry the same polarity charge, while in the hetero-complex cases, more than 58% of the complexes have oppositely charged interfaces. Thus, in the vast majority of the cases, the Coulombic energy has a large positive magnitude for homo-complexes (interactions between the same polarity charges across interfaces), while this percentage is much smaller for hetero-complexes. This results in a larger amplitude of the calculated ΔΔGel for homo-complexes compared with the ΔΔGel calculated for the hetero-complexes. Increasing the internal dielectric constant combined with minimizing the energy of the structures could turn some of the small positive ΔΔGel calculated for the hetero-complexes into small negative ΔΔGel. However, ΔΔGel are large positive numbers in the homo-complexes set, and an increase of the internal dielectric constant combined with the refinement of the structures cannot dramatically change the percentage of the cases in which electrostatics was calculated to oppose binding.
The observation made for the hetero-complexes, that their structural refinement through the energy minimization of corresponding complexes results in a change of the sign of the calculated ΔΔGel, deserves attention. The electrostatics is one of the components of the total energy and the total energy is minimized in the minimization protocol. There is no reason to expect that minimizing the total energy should result in the minimization of the energy components. However, apparently the electrostatic component was actually optimized in the energy minimization protocol. Perhaps this indicates that if a given energy component favors a particular energy state, the minimization of the total energy will most likely result in the enhancement of this energy component as well.
Effects of the probe radius
Slopes of the fitting lines and the corresponding correlation coefficients
Type of complex
Probe radius R = 1.4A and different ε(in)
ε(in) = 2 and different probe radius R
Slope of the fitting line
The effects of the force field
The parameters of the distributions of ΔΔΔGel calculated for hetero-complexes.
The parameters of the distributions of ΔΔΔGel calculated for homo-complexes.
where X, Y stand for either the Charmm27, Amber98 or OPLS force fields. These differences were calculated for all types of complexes, both non-minimized and minimized. The parameters of the resulting distributions are reported in Tables 3, 4. The corresponding graphs are shown in Additional file 1.
It can be seen that both the correlation coefficients and the slope of the fitting lines are close to 1.0, indicating that ΔΔGels calculated with different force fields are quite similar. However, a closer look at Tables 3, 4 shows significant differences as well. The mean of ΔΔGel calculated with Charmm27 is 5.8 and is 10.3 kcal/mol less positive than that of the corresponding means calculated with Amber98 in the case of non-minimized and minimized hereto-complexes, respectively. The same offset is even more pronounced when comparing the results obtained with the Amber and OPLS force fields. In the hetero-complex cases, the means differ by 13.0 and 19.0 kcal/mol for non-minimized and minimized complexes, respectively. Lastly, comparing the results obtained with the Charmm27 and OPLS force fields, we observed that the mean of OPLS is less positive as compared to the mean of the energy distribution calculated with Charmm27. Thus, these results suggest that the OPLS force field resulted in less positive ΔΔGel (less unfavorable energies), and the Charmm27 and Amber98 force fields resulted in the most unfavorable ΔΔGel in the hetero-complex cases.
In addition to the above observations, in some exceptional cases the ΔΔGel calculated with different force fields differed by more than 100 kcal/mol. Particular example is tail-associated lysozyme bound to baseplate structural protein GP27 (PDB ID 1k28) for which ΔΔGel calculated with Charmm27, e(in) = 1.0 and probe radius 1.4A was 261.7 kcal/mol while it was calculated to be 391.1 kcal/mol with Amber98 force field. We have not investigated these exceptions in detail; however, the differences reported in Tables 3, 4 as "Min" and "Max" indicate such cases. In many cases, even a difference of several kcal/mol could be crucial for getting the correct biophysics of the binding.
We have done extensive testing of the sensitivity of the electrostatic component of the binding free energy with respect to internal dielectric constant values, probe radius, and several commonly used force fields. The goal was also to reveal the general trends and to elaborate on the role of electrostatics on the binding for hetero- and homo-complexes. The study was done on a very large set of protein-protein complexes: 260 hetero-complexes and 2148 homo-complexes. Energy minimization and full-scale electrostatic calculations have never been done on such a large dataset before. This ensures that the obtained results and conclusions are statistically significant.
A common practice is to validate the results of numerical simulations against experimental data. However, it is impossible to experimentally determine the electrostatic component of the binding free energy. The binding free energy is made of a variety of energy terms whose interplay results in the observed affinity. Even the pH-dependence of the binding is not a pure electrostatic event. The changes of the protonation states of ionizable groups induced by either pH changes or the binding in turn cause conformational changes and thus alter the contribution of non-electrostatic energy terms to the binding affinity as demonstrated in a recent study. The most successful attempt so far to compare electrostatic calculations to experimental data was done by Dong and Zhou. They carried out Poisson-Boltzmann calculations on a set of 64 mutations over six protein-protein complexes and compared their predictions with experimental data. The mutants were selected to involve charge-charge interactions across the interface of the corresponding complex and thus suggested the putative dominance of electrostatic interactions over the rest of the energy terms contributing to the binding. It was shown that the calculations done with a zero probe radius and an internal dielectric constant of 4.0 provided the closest agreement with experimental data. However, as pointed out by the authors, other contributions, such as van der Waals interactions and hydrophobic effect, may play important role. Modeling the non-electrostatic contributions on such a large set of data as ours could introduce significant noise and obscure the results. That is our reasoning for not attempting to compare our calculations to cases in which the total binding free energy was experimentally measured. Proper modeling of the absolute binding free energy requires that the entropy be also accounted for, a task that is not easily achieved, especially on a large set of data.
Despite the fact that electrostatics is only one of the many forces contributing to the binding, in many cases, its involvement could be related to biologically important effects. For this reason, many researchers employ electrostatic calculations to find out the electrostatic component of the binding free energy. However, the value of the dielectric constant, the method of determining the molecular surface (probe radius), and the force field parameters are still a personal preference. Our study demonstrated that the absolute value of the calculated ΔΔGel is very sensitive to all of the above-mentioned parameters. While this is expected to be the case for different internal dielectric constants and to some extent for different probe radii, the observation that the results are quite sensitive to the force field parameters deserves attention. The average difference of the ΔΔGel calculated with different force fields can be as large as 20 kcal/mol and more. It seems to us that the energy minimization with the corresponding force field does not make much difference. Since the minimization protocol minimizes the total energy and not just the electrostatic component, the outcome for electrostatic energy will vary case by case. All of these observations indicate that the absolute values of the ΔΔGel should be considered with precaution.
This study revealed a major difference between hetero- and homo-complexes with respect to the calculated ΔΔGel. The calculated electrostatic component of the binding free energy opposes binding for approximately 80% of the homo-complex cases, despite the internal dielectric constant, probe radius, and force field. In contrast, the role of electrostatics on the binding of hetero-complexes depends on all of the factors above. Using a low internal dielectric constant and a probe radius of 1.4A, despite the force field used, most of the calculated ΔΔGel opposes binding. However, increasing the internal dielectric constant to 20 or decreasing the probe radius to R = 0.0A results in 60% of the cases having a ΔΔGel favoring the binding. These findings offer a pragmatic prescription for assessing the confidence of the predicted role of ΔΔGel on the binding affinity. It seems that despite the differences between the homo- and hetero-complexes, if the binding free energy is calculated to be a large positive quantity, then the conclusions could be considered independent from the choice of the parameters. However, cases in which the ΔΔGel are calculated with a particular set of parameters and corresponding magnitude is a small number close to zero are more complicated. Our study indicates that the calculated ΔΔGel could easily change sign upon changing the parameters of the protocol or upon switching the force field. Therefore, the conclusions made for the role of electrostatics on the binding in case of a small ΔΔGel should be carefully investigated by varying the force fields and the parameters of the computational protocol.
Our study, being applied on 2408 protein-protein complexes, offers the possibility to statistically address the question of what role electrostatics play in protein-protein binding. A similar question of what role salt bridges play on protein stability was extensively studied by Tidor and co-workers [78–80]. It was concluded that in most cases, salt bridges destabilize proteins. They also investigated the effect of electrostatic interactions on protein binding since salt bridges are also formed across protein-protein interfaces. It was found that charge interactions in barnase-barstar complexes are optimized to favor the binding[25, 26]. Schreiber and co-workers[22, 82, 83] employed experimental and computational methods to investigate the role of electrostatic interactions on barnase-barstar association. It was shown that charge-charge interactions can be re-designed to result in tighter binding. The association rate was also studied on a set 68 transient hetero-complexes, demonstrating that electrostatics plays a marginal role. Recently the association rates of four protein complexes and twenty three mutants were modeled by Zhou and co-workers . It was shown that the calculations with vdW surface presentation (zero probe radius) and non-linear PB equation can successfully reproduce experimental data, while the models with molecular surface presentation (probe radius 1.4A) gave nonphysical results. The role of electrostatic interactions on the stability of monomeric protein was also systematically investigated. Makhatadze and co-workers have shown that so called "complex" salt bridges, in which the anchor residue forms hydrogen bonds with two or more opposite charged residues simultaneously, contribute significantly to protein stability . This observation is valid not only for buried but for surface exposed groups as well [87–90]. It terms of different types of hydrophilic group, Pace and co-workers demonstrated that Asp and Glu amino acids contribute the most to protein stability , however the individual contributions are strongly context-dependent . These observations and the successes in re-engineering of more stable proteins and tighter complexes than the wild type, demonstrate that the current computational methods can adequately model electrostatic interactions in biological macromolecules. However, it should be pointed out that the effect of re-designing is calculated and measured in respect to a reference state, the wild type protein or protein complex. Thus, an increased stability or tighter binding due to electrostatic optimization does not necessary mean that electrostatics plays favorable role. The only conclusion that can be made is that the re-engineering makes electrostatic contribution either more favorable or less unfavorable.
The results of our study indicate that the role electrostatics have on protein-protein binding in hetero-complexes cases is the same as the role salt bridges have on protein stability; i.e., in some cases, it will favor the binding, but in other cases it will oppose the binding. In our previous study  we have shown that electrostatic component of the binding free energy tends to be optimized in respect to random shuffling of the amino acid sequence of the corresponding partners. It was also demonstrated that the salt dependence of the binding is not correlated with macroscopic parameters of the monomers [63, 94]. Taking together these observations, the findings of above mentioned studies and the results presented in this work, it could be stated that in majority of the cases electrostatics favor hetero-complexes formation. However, the role of electrostatics for each individual case depends on the fine details at the interface as formation of "complex" salt bridges, for example. The situation is quite different for homo-complexes. As was shown, all of the computational protocols predicted that, in majority of the cases, the electrostatics opposes binding. Since it was stated that macroscopic parameters do not matter much for protein-protein binding, this is a surprising observation. However, Table 1 provides an explanation. The homo-complexes differ from hetero-complexes not only by their macroscopic parameters as being made of two monomer carrying the same charge, but their interfaces are also different. Only 8% of homo-complexes have opposite charged interfaces, while this percentage is 58% for hetero-complexes. In remaining 92% of the cases, homo-complexes are formed between interfaces carrying either the same net charge or not having net charge at all. Perhaps, formation of "complex" salt bridges  and favorable polar interaction is much more difficult in this case as compared with hetero-complexes.
In this work we have not considered possible protonation changes induced by the binding. It was done in order to reduce computational demand of the numerical protocol. However, protonation changes may be common phenomena in protein binding as demonstrated recently by Jensen and co-workers . Our recent investigation  also indicates that from statistical point of view there is 70% chance that binding will cause proton uptake/release of more than 0.5 proton units at pH = 7.0. Adjusting the ionization states according to the predicted pKa's of ionizable groups will definitely make the calculated ΔΔGel more favorable (or less unfavorable) compared with calculations done with default charge states. However, most of the predicted ionization changes are fractional numbers and modeling of non-integer charge changes requires ensemble approach in computing ΔΔGel, a task which is quite difficult to accomplish on a set of 2408 protein complexes.
The study revealed that for a significant fraction of hetero-complexes and in vast majority of homo-complexes the electrostatics opposes binding. Since most of the protein-protein complexes in the cell are homo-complexes, the overall role of electrostatics is predicted to oppose binding. Instead, electrostatics provides the necessary specificity and steering, processes equally important for protein-protein association.
The analysis of the sensitivity of the ΔΔGel, calculated with different force fields, internal dielectric constants, probe radius value and minimization protocols, gives us the opportunity to suggest a set of rules of calculating (a) the absolute value of ΔΔGel and (b) the sign of ΔΔGel: (a1) if there is no prior knowledge what the effective value of both internal dielectric constant and probe radius are for a given protein complex and a particular protocol of calculating the ΔΔGel, then the absolute value of calculated ΔΔGel is meaningless; (a2) provided that the effective value of both internal dielectric constant and probe radius are known for a given protein complex and a particular protocol of calculating the ΔΔGel, then the absolute value of ΔΔGel should be obtained as an average of ΔΔGel calculated with different force fields; (a3) energy minimization does not help in obtaining consistent ΔΔGel's, since it minimizes the total energy, not just the electrostatic component. Therefore, the absolute value of ΔΔGel should be obtained as an average of ΔΔGel calculated with different force fields, provided that other parameters are known; (b1) if there is no prior knowledge what the effective value of both internal dielectric constant and probe radius are for a given protein complex and a particular protocol of calculating the ΔΔGel, then the sign of ΔΔGel (determining the role of electrostatics on binding) is meaningful only if the absolute ΔΔGel calculated with either high internal dielectric constant or zero probe radius is larger than 1 kcal/mol. In case of hetero-complexes, if ΔΔGel is calculated to be positive and larger than 1 kcal/mol with ε(in) = 20 and probe radius 1.4A with given force field, then the probability of remaining positive in the calculations performed with different force fields is 0.99 (or 0.97 with ε(in) = 2). In case of homo-complexes, this probability is 0.95 (or 0.94 with ε(in) = 2); (b2) if there is a prior knowledge of the effective value of internal dielectric constant and probe radius, then calculations with a particular force field with and without minimization provide a good estimate of ΔΔGel sign in about 90% of the cases. In case of hetero-complexes, if ΔΔGel is calculated to be positive and larger than 1 kcal/mol with ε(in) = 20 and probe radius 1.4A with given force field, then the probability of remaining positive in the calculations performed with different dielectric constant is 0.98. Probability of change of the sign ΔΔGel using different probe radii is about 0.1. In case of homo-complexes, these probabilities are 0.96 and 0.87.
The minimization of the structures used in this work was made possible by a Condor pool deployed and maintained by Clemson Computing and Information Technology. The authors would like to acknowledge the support of the staff from the Cyber Infrastructure Technology Integration group, especially Barr Von Oehsen. This research was supported by an award to Clemson University from the Howard Hughes Medical Institute Undergraduate Science Education Program.
- Alberts B, Bray D, Lewis J, Raff M, Roberts K, Watson J: Molecular Biology of the Cell. 1994, Garland PublishingGoogle Scholar
- Alexov E: Protein-protein interactions. Curr Pharm Biotechnol. 2008, 9 (2): 55-56.View ArticleGoogle Scholar
- Sham Y, Chu Z, Tao H, Warshel A: Examining Methods for Calculations of Binding Free Energies: LRA, LIE, PDLD-LRA, and PDLD/S-LRA Calculations of Ligands Binding to an HIV Protease. Proteins. 2000, 39: 393-407.View ArticleGoogle Scholar
- Keskin O, Ma BY, Rogale K, Gunasekaran K, Nussinov R: Protein-protein interactions: organization, cooperativity and mapping in a bottom-up Systems Biology approach. Phys Biol. 2005, 2 (2): S24-S35.View ArticleADSGoogle Scholar
- McDonald IK, Thornton JM: Satisfying hydrogen bonding potential in proteins. J Mol Biol. 1994, 238: 777-793.View ArticleGoogle Scholar
- Jones S, Thornton J: Principles of protein-protein interactions derived from structural studies. Proceedings of the National Academy of Sciences. 1996, 93: 13-20.View ArticleADSGoogle Scholar
- Sheinerman F, Norel R, Honig B: Electrostatics Aspects of Protein-Protein Interactions. Current Opinion in Structural Biology. 2000, 10: 153-159.View ArticleGoogle Scholar
- Kundrotas PJ, Alexov E: Electrostatic properties of protein-protein complexes. Biophys J. 2006, 91 (5): 1724-1736.View ArticleGoogle Scholar
- Lensink MF, Mendez R: Recognition-induced conformational changes in protein-protein docking. Curr Pharm Biotechnol. 2008, 9 (2): 77-86.View ArticleGoogle Scholar
- Keskin O, Gursoy A, Ma B, Nussinov R: Principles of Protein-Protein Interactions: What are the Preferred Ways For Proteins To Interact?. Chem Rev. 2008, 108 (4): 1225-1244.View ArticleGoogle Scholar
- Jensen JH: Calculating pH and salt dependence of protein-protein binding. Curr Pharm Biotechnol. 2008, 9 (2): 96-102.View ArticleGoogle Scholar
- Bordner AJ, Abagyan R: Statistical analysis and prediction of protein-protein interfaces. Proteins. 2005, 60 (3): 353-366.View ArticleGoogle Scholar
- Jones S, Thornton J: Prediction of Protein-Protein Interaction Sites using Patch Analysis. J Mol Biol. 1997, 272: 113-143.Google Scholar
- Jones S, Thornton J: Principles of protein-protein interactions. PNAS (USA). 1996, 93: 13-20.View ArticleADSGoogle Scholar
- Nooren IMA, Thornton JM: Structural characterisation and functional significance of transient protein-protein interactions. Journal of Molecular Biology. 2003, 325 (5): 991-1018.View ArticleGoogle Scholar
- Albeck S, Schreiber G: Biophysical characterization of the interaction of the beta-lactamase TEM-1 with its protein inhibitor BLIP. Biochemistry. 1999, 38 (1): 11-21.View ArticleGoogle Scholar
- Dong F, Vijayakumar M, Zhou HX: Comparison of calculation and experiment implicates significant electrostatic contributions to the binding stability of barnase and barstar. Biophys J. 2003, 85 (1): 49-60.View ArticleGoogle Scholar
- Elcock A, Sept D, McCammon J: Computer simulation of protein-protein interactions. J Phys Chem. 2000, 105: 1504-1518.View ArticleGoogle Scholar
- Elcock AH, Sept D, McCammon JA: Computer simulation of protein-protein interactions. J Phys Chem B. 2001, 105 (8): 1504-1518.View ArticleGoogle Scholar
- Clackson T, Wells JA: A Hot-Spot of Binding-Energy in a Hormone-Receptor Interface. Science. 1995, 267 (5196): 383-386.View ArticleADSGoogle Scholar
- Wells JA: Systematic Mutational Analyses of Protein Protein Interfaces. Method Enzymol. 1991, 202: 390-411.View ArticleGoogle Scholar
- Schreiber G, Fersht AR: Interaction of barnase with its popypeptide inhibitor barstar studied by protein engineering. Biochemistry. 1993, 32: 5145-5150.View ArticleGoogle Scholar
- Schreiber G, Frisch C, Fersht AR: The role of Glu73 of barnase in catalysis and the binding of barstar. J Mol Biol. 1997, 270: 111-122.View ArticleGoogle Scholar
- Massova I, Kollman PA: Computational alanine scanning to probe protein-protein interactions: A novel approach to evaluate binding free energies. J Am Chem Soc. 1999, 121 (36): 8133-8143.View ArticleGoogle Scholar
- Lee LP, Tidor B: Optimization of binding electrostatics: charge complementarity in the barnase-barstar protein complex. Protein Sci. 2001, 10 (2): 362-377.View ArticleGoogle Scholar
- Lee LP, Tidor B: Barstar is electrostatically optimized for tight binding to barnase. Nat Struct Biol. 2001, 8 (1): 73-76.View ArticleGoogle Scholar
- Schreiber G, Ferst A: Interaction of Barnase with Its Polypeptide Inhibitor Barstar Studied by Protein Engineering. Biochemistry. 1993, 32: 5145-5150.View ArticleGoogle Scholar
- Frisch C, Schreiber G, Johnson CM, Fersht AR: Thermodynamics of the interactions of barnase and barstar: changes in free energy versus changes in enthalpy on mutation. J Mol Biol. 1997, 267: 696-706.View ArticleGoogle Scholar
- Lee L, Tidor B: Optimization of binding electrostatics: Charge complementarity in the barnase-barstar protein complex. Protein Sci. 2001, 10: 362-377.View ArticleGoogle Scholar
- Seizer T, Albeck S, Schreiber G: Rational design of faster associating and tigher binding protein complexes. Nature Struc Biol. 2000, 7: 537-541.View ArticleGoogle Scholar
- Ma CS, Baker NA, Joseph S, McCammon JA: Binding of aminoglycoside antibiotics to the small ribosomal subunit: A continuum electrostatics investigation. J Am Chem Soc. 2002, 124 (7): 1438-1442.View ArticleGoogle Scholar
- Sims PA, Wong CF, McCammon JA: Charge optimization of the interface between protein kinases and their ligands. J Comput Chem. 2004, 25 (11): 1416-1429.View ArticleGoogle Scholar
- Sims PA, Wong CF, Vuga D, McCammon JA, Sefton BM: Relative contributions of desolvation, inter- and intramolecular interactions to binding affinity in protein kinase systems. J Comput Chem. 2005, 26 (7): 668-681.View ArticleGoogle Scholar
- Norel R, Sheinerman F, Petrey D, Honig B: Electrostatic contribution to protein-protein interactions: Fast energetic filters for docking and their physical basis. Prot Sci. 2001, 10: 2147-2161.View ArticleGoogle Scholar
- Sheinerman FB, Honig B: On the role of electrostatic interactions in the design of protein-protein interfaces. Journal of Molecular Biology. 2002, 318 (1): 161-177.View ArticleGoogle Scholar
- Muegge I, Schweins T, Warshel A: Electrostatic Contribution to Protein-Protein Binding Affinities: Application to Rap/Raf Interaction. Proteins. 1998, 30: 407-423.View ArticleGoogle Scholar
- Gohlke H, Kiel C, Case D: Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes. J Mol Biol. 2003, 330: 891-913.View ArticleGoogle Scholar
- Sharp KA: Electrostatic Interactions in Macromolecules. Curr Opin Struct Biol. 1994, 4: 234-239.View ArticleMathSciNetGoogle Scholar
- Gilson M, Honig B: The Dielectric Constant of a Folded Protein. Biopolymers. 1986, 25: 2097-2119.View ArticleGoogle Scholar
- Simonson T, Brooks C: Charge Screening and the Dielectric Constant of Proteins: Insights from Molecular Dynamics. J Am Chem Soc. 1996, 118: 8452-8458.View ArticleGoogle Scholar
- Schulz C, Warshel A: What Are the Dielectric "Constants" of Proteins and How To Validate Electrostatic Models. Proteins. 2001, 44: 400-417.View ArticleGoogle Scholar
- Brooks BR, Bruccoleri RE, Olafson BD, States DJ, Swaminathan S, Karplus M: CHARMM: A program for macromolecular energy, minimization and dynamic calculations. J Comp Chem. 1983, 4: 187-217.View ArticleGoogle Scholar
- Gilson MK, Honig B: The inclusion of electrostatic hydration energies in molecular mechanics calculations. J Comput Aided Mol Des. 1991, 5 (1): 5-20.View ArticleGoogle Scholar
- Schiffer CA, Caldwell J, Stroud R, Kollman P: Inclusion of Solvation Free Energy with Molecular Mechanics Energy: Alanine Dipeptide as a Test Case. Prot Sci. 1992, 1: 396-400.View ArticleGoogle Scholar
- Georgescu R, Alexov E, Gunner M: Combining Conformational Flexibility and Continuum Electrostatics for Calculating Residue pKa's in Proteins. Biophys J. 2002, 83: 1731-1748.View ArticleGoogle Scholar
- Antosiewicz J, McCammon J, Gilson M: Prediction of pH dependent properties of proteins. JMol Bio. 1994, 238: 415-436.View ArticleGoogle Scholar
- Dong F, Zhou H-X: Electrostatic contribution to the binding stability of protein-protein complexes. Proteins. 2006, 65: 87-102.View ArticleGoogle Scholar
- Qin S, Zhou HX: Do electrostatic interactions destabilize protein-nucleic acid binding?. Biopolymers. 2007Google Scholar
- Antosiewicz J, Briggs J, Elcock A, Gilson M, McCammon J: Computing the ionization states of proteins with a detail charge model. J Comp Chem. 1996, 17: 1633-1644.View ArticleGoogle Scholar
- Nielsen J, Vriend G: Optimizing the Hydrogen-Bond Network in Poisson-Boltzmann Equation-Based pKa Calculations. Proteins. 2001, 43: 403-412.View ArticleGoogle Scholar
- Nielsen J, Andersen K, Honig B, Hooft R, Klebe G, Vriend G, Wade R: Improving macromolecular electrostatic calculations. Protein Eng. 1999, 12: 657-662.View ArticleGoogle Scholar
- Forrest L, Honig B: An assessment of the accuracy of methods for predicting hydrogen positions in protein structures. Proteins. 2005, 61: 296-309.View ArticleGoogle Scholar
- Nielsen J, McCammon A: On the evaluation and optimization of protein X-ray structures for pKa calculations. Protein Science. 2003, 12: 313-326.View ArticleGoogle Scholar
- Kundrotas PJ, Alexov E: PROTCOM: searchable database of protein complexes enhanced with domain-domain structures. Nucleic Acids Res. 2007, 35: D575-D579.View ArticleGoogle Scholar
- Li W, Jaroszewski L, Godzik A: Clustering of highly homologous sequences to reduce the size of large protein databases. Bioinformatics. 2001, 17 (3): 282-283.View ArticleGoogle Scholar
- Rocchia W, Sridharan S, Nicholls A, Alexov E, Chiabrera A, Honig B: Rapid grid-based construction of the molecular surface and the use of induced surface charge to calculate reaction field energies: applications to the molecular systems and geometric objects. J Comput Chem. 2002, 23 (1): 128-137.View ArticleGoogle Scholar
- Ponder JW: TINKER-software tools for molecular design. 1999, St. Luis: Washington University, 3.7Google Scholar
- Brooks BR, Bruccoleri RE, Olafson BD, States DJ, Swaminathan S, Karplus M: CHARMM: a program for macromolecular energy, minimization, and dynamics calculations. J Comput Chem. 1983, 4: 187-217.View ArticleGoogle Scholar
- Wang J, Wolf RM, Caldwell JW, Kollman PA, Case DA: Development and testing of a general amber force field. J Comput Chem. 2004, 25 (9): 1157-1174.View ArticleGoogle Scholar
- Jorgensen WL, Tirado-Rives J: The OPLS potential function for proteins. J Am Chem Soc. 1988, 110: 1657-1666.View ArticleGoogle Scholar
- Still WC, Tempczyk A, Hawley RC, Hendrickson T: Semianalytical Treatment of Solvation for Molecular Mechanics and Dynamics. J Am Chem Soc. 1990, 112: 6127-6129.View ArticleGoogle Scholar
- MacKerell AD, Bashford D, Bellot M, Dunbrack RL, Evanseck JD, Field MJ, Fischer S, Gao J, Guo H, Ha S, et al: All-atom empirical potential for molecular modeling and dynamics studies of proteins. J Phys Chem. 1998, 102: 3586-3616.View ArticleGoogle Scholar
- Bertonati C, Honig B, Alexov E: Poisson-Boltzmann calculations of nonspecific salt effects on protein-protein binding free energies. Biophys J. 2007, 92 (6): 1891-1899.View ArticleGoogle Scholar
- Rocchia W, Alexov E, Honig B: Extending the applicability of the nonlinear Poisson-Boltzmann equation: Multiple dielectric constants and multivalent ions. J Phys Chem. 2001, 105 (85): 6507-6514.View ArticleGoogle Scholar
- Rocchia W, Sridharan S, Nicholls A, Alexov E, Chiabrera A, Honig B: Rapid Grid-based Construction of the Molecular Surface and the Use of Induced Surface Charges to Calculate Reaction Field Energies: Applications to the Molecular Systems and Geometrical Objects. J Comp Chem. 2002, 23: 128-137.View ArticleGoogle Scholar
- Nielsen JE, McCammon JA: Calculating pKa values in enzyme active sites. Protein Science. 2003, 12 (9): 1894-1901.View ArticleGoogle Scholar
- Gomez J, Freire E: Thermodynamic Mapping of the Inhibitor Site of the Aspartic Protease Endothiapepsin. J Mol Biol. 1995, 252: 337-350.View ArticleGoogle Scholar
- Elcock A, McCammon A: Electrostatic Contributions to the Stability of Halophilic Proteins. J Mol Biol. 1998, 280: 731-748.View ArticleGoogle Scholar
- Sharp KA: Electrostatic interactions in hirudin-thrombin binding. Biophys Chem. 1996, 61: 37-49.View ArticleGoogle Scholar
- Trylska J, Antosiewich J, Geller M, Hodge C, Klabe R, Head M, Gilson M: Thermodynamic linkage between the binding of protons and inhibitors to HIV-1 protease. Protein Science. 1999, 8: 180-195.View ArticleGoogle Scholar
- Alexov E: Calculating Proton Uptake/Release and the Binding Free Energy Taking into Account Ionization and Conformation Changes Induced by Protein-Inhibitor Association. Application to Plasmepsin, Cathepsin D and Endothiapepsin-Pepstatin Complexes. Proteins. 2004, 56: 572-584.View ArticleGoogle Scholar
- Richards FM: Areas, volumes, packing and protein structure. Ann Rev Biophys Bioeng. 1977, 6: 151-176.View ArticleGoogle Scholar
- Vijayakumar M, Zhou H: Salt Bridges Stabilize the Folded State of Barnase. J Phys Chem. 2001, 105: 7334-7340.View ArticleGoogle Scholar
- Sham Y, Chu Z, Warshel A: Consistent Calculations of pKa's of Ionizable Residues in Proteins: Semi-microscopic amd Microscopic Approaches. J Phys Chem. 1997, 101: 4458-4472.View ArticleGoogle Scholar
- Roca M, Messer B, Warshel A: Electrostatic contributions to protein stability and folding energy. FEBS Lett. 2007, 581 (10): 2065-2071.View ArticleGoogle Scholar
- Warshel A, Sharma PK, Kato M, Parson WW: Modeling electrostatic effects in proteins. Biochim Biophys Acta. 2006, 1764 (11): 1647-1676.View ArticleGoogle Scholar
- Alexov E: The role of the protein side chain fluctuations on the strength of pair wise electrostatic interactions. Comparing experimental with computed pKa's. Proteins. 2003, 50: 94-103.View ArticleGoogle Scholar
- Hendsch Z, Tidor B: Do salt bridges stabilize proteins? A continuum electrostatics analysis. Prot Sci. 1994, 3: 211-226.View ArticleGoogle Scholar
- Luisi DL, Snow CD, Lin JJ, Hendsch ZS, Tidor B, Raleigh DP: Surface salt bridges, double-mutant cycles, and protein stability: an experimental and computational analysis of the interaction of the Asp 23 side chain with the N-terminus of the N-terminal domain of the ribosomal protein l9. Biochemistry. 2003, 42 (23): 7050-7060.View ArticleGoogle Scholar
- Hendsch ZS, Sindelar CV, Tidor B: Parameter dependence in continuum electrostatic calculations: a study using protein salt bridges. J Phys Chem. 1998, 102: 4404-4410.View ArticleGoogle Scholar
- Kangas E, Tidor B: Charge optimization leads to favorable electrostatic binding free energy. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics. 1999, 59 (5 Pt B): 5958-5961.Google Scholar
- Martinez JC, Filimonov VV, Mateo PL, Schreiber G, Fersht AR: A calorimetric study of the thermal-stability of barstar and its interaction with barnase. Biochemistry. 1995, 34: 5224-5233.View ArticleGoogle Scholar
- Schreiber G, Fersht AR: Energetics of protein-protein interactions: analysis of the barnase-barstar interface by single mutations and double mutant cycles. J Mol Biol. 1995, 248: 478-486.Google Scholar
- Shaul Y, Schreiber G: Exploring the charge space of protein-protein association: A proteomic study. Proteins-Structure Function and Bioinformatics. 2005, 60 (3): 341-352.View ArticleGoogle Scholar
- Alsallaq R, Zhou HX: Electrostatic rate enhancement and transient complex of protein-protein association. Proteins. 2008, 71 (1): 320-335.View ArticleGoogle Scholar
- Gvritishvili AG, Gribenko AV, Makhatadze GI: Cooperativity of complex salt bridges. Protein Sci. 2008, 17 (7): 1285-1290.View ArticleGoogle Scholar
- Schweiker KL, Zarrine-Afsar A, Davidson AR, Makhatadze GI: Computational design of the Fyn SH3 domain with increased stability through optimization of surface charge charge interactions. Protein Sci. 2007, 16 (12): 2694-2702.View ArticleGoogle Scholar
- Gribenko AV, Makhatadze GI: Role of the charge-charge interactions in defining stability and halophilicity of the CspB proteins. J Mol Biol. 2007, 366 (3): 842-856.View ArticleGoogle Scholar
- Strickler SS, Gribenko AV, Gribenko AV, Keiffer TR, Tomlinson J, Reihle T, Loladze VV, Makhatadze GI: Protein stability and surface electrostatics: a charged relationship. Biochemistry. 2006, 45 (9): 2761-2766.View ArticleGoogle Scholar
- Makhatadze GI, Loladze VV, Ermolenko DN, Chen X, Thomas ST: Contribution of surface salt bridges to protein stability: guidelines for protein engineering. J Mol Biol. 2003, 327 (5): 1135-1148.View ArticleGoogle Scholar
- Trevino SR, Scholtz JM, Pace CN: Amino acid contribution to protein solubility: Asp, Glu, and Ser contribute more favorably than the other hydrophilic amino acids in RNase Sa. J Mol Biol. 2007, 366 (2): 449-460.View ArticleGoogle Scholar
- Takano K, Scholtz JM, Sacchettini JC, Pace CN: The contribution of polar group burial to protein stability is strongly context-dependent. J Biol Chem. 2003, 278 (34): 31790-31795.View ArticleGoogle Scholar
- Brock K, Talley K, Coley K, Kundrotas P, Alexov E: Optimization of electrostatic interactions in protein-protein complexes. Biophys J. 2007, 93 (10): 3340-3352.View ArticleGoogle Scholar
- Talley K, Alexov E: Modelling Salt Dependence of Protein-Protein Association:Linear vs Non-Linear Poisson-Bolzmann Equation. Comm in Computational Physics. 2008, 4: 169-179.Google Scholar
- Mason AC, Jensen JH: Protein-protein binding is often associated with changes in protonation state. Proteins. 2007Google Scholar
- Mitra R, Shyam R, Mitra I, Miteva MA, Alexov E: Calculation of the proponation states of proteins and small molecules: Inmplications to ligand-receptor interactions. Current computer-aided drug design. 2008, 4: 169-179.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.