Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes

Absolute binding free energy calculations and free energy decompositions are presented for the protein-protein complexes H-Ras/C-Raf1 and H-Ras/RalGDS. Ras is a central switch in the regulation of cell proliferation and differentiation. In our study, we investigate the capability of the molecular mechanics (MM)-generalized Born surface area (GBSA) approach to estimate absolute binding free energies for the protein-protein complexes. Averaging gas-phase energies, solvation free energies, and entropic contributions over snapshots extracted from trajectories of the unbound proteins and the complexes, calculated binding free energies (Ras-Raf: -15.0(+/-6.3)kcal mol(-1); Ras-RalGDS: -19.5(+/-5.9)kcal mol(-1)) are in fair agreement with experimentally determined values (-9.6 kcal mol(-1); -8.4 kcal mol(-1)), if appropriate ionic strength is taken into account. Structural determinants of the binding affinity of Ras-Raf and Ras-RalGDS are identified by means of free energy decomposition. For the first time, computationally inexpensive generalized Born (GB) calculations are applied in this context to partition solvation free energies along with gas-phase energies between residues of both binding partners. For selected residues, in addition, entropic contributions are estimated by classical statistical mechanics. Comparison of the decomposition results with experimentally determined binding free energy differences for alanine mutants of interface residues yielded correlations with r(2)=0.55 and 0.46 for Ras-Raf and Ras-RalGDS, respectively. Extension of the decomposition reveals residues as far apart as 25A from the binding epitope that can contribute significantly to binding free energy. These "hotspots" are found to show large atomic fluctuations in the unbound proteins, indicating that they reside in structurally less stable regions. Furthermore, hotspot residues experience a significantly larger-than-average decrease in local fluctuations upon complex formation. Finally, by calculating a pair-wise decomposition of interactions, interaction pathways originating in the binding epitope of Raf are found that protrude through the protein structure towards the loop L1. This explains the finding of a conformational change in this region upon complex formation with Ras, and it may trigger a larger structural change in Raf, which is considered to be necessary for activation of the effector by Ras.     

蛋白质-蛋白质对接中打分函数的研究

通过分析蛋白质－蛋白质间的静电、疏水作用和熵效应与相对于晶体结构的蛋白质主链原子的均方根偏差（RMSD）的相关性,定量地考查了它们在蛋白质－蛋白质对接中作为打分函数评价近天然构象的能力。对7个蛋白质复合物体系的分析表明,就水化能而言,原子接触势模型（ACE）优于原子水化参数模型（ASP）,且修正的ACE模型具有更好的评价近天然构象的能力;水化能与静电能结合对评价能力有进一步的提高。最后,我们将静电和修正的ACE水化能结合作为打分函数用于36个蛋白质复合物体系的对接研究,进一步证实了这两种能量项的组合能有效地将近天然结构从分子对接模式中区分出来。     

Hydrophobic regions on protein surfaces. Derivation of the solvation energy from their area distribution in crystallographic protein structures.

For the first time, a direct approach for the derivation of an atomic solvation parameter from macromolecular structural data alone is presented. The specific free energy of solvation for hydrophobic surface regions of proteins is delineated from the area distribution of hydrophobic surface patches. The resulting value is 18 cal/(mol.A2), with a statistical uncertainty of +/-2 cal/mol.A2) at the 5% significance level. It compares favorably with the parameters for carbon obtained by other authors who use the the crystal geometry of succinic acid or energies of transfer from hydrophobic solvent to water for small organic compounds. Thus, the transferability of atomic solvation parameters for hydrophobic atoms to macromolecules has been directly demonstrated. A careful statistical analysis demonstrates that surface energy parameters derived from thermodynamic data of protein mutation experiments are clearly less confident.     

Quantifying the effect of burial of amino acid residues on protein stability

The average contribution of individual residue to folding stability and its dependence on buried accessible surface area (ASA) are obtained by two different approaches. One is based on experimental mutation data, and the other uses a new knowledge-based atom-atom potential of mean force. We show that the contribution of a residue has a significant correlation with buried ASA and the regression slopes of 20 amino acid residues (called the buriability) are all positive (pro-burial). The buriability parameter provides a quantitative measure of the driving force for the burial of a residue. The large buriability gap observed between hydrophobic and hydrophilic residues is responsible for the burial of hydrophobic residues in soluble proteins. Possible factors that contribute to the buriability gap are discussed.     

Nine hydrophobic side chains are key determinants of the thermodynamic stability and oligomerization status of tumour suppressor p53 tetramerization domain.

The contribution of almost each amino acid side chain to the thermodynamic stability of the tetramerization domain (residues 326-353) of human p53 has been quantitated using 25 mutants with single-residue truncations to alanine (or glycine). Truncation of either Leu344 or Leu348 buried at the tetramer interface, but not of any other residue, led to the formation of dimers of moderate stability (8-9 kcal/mol of dimer) instead of tetramers. One-third of the substitutions were moderately destabilizing (<3.9 kcal/mol of tetramer). Truncations of Arg333, Asn345 or Glu349 involved in intermonomer hydrogen bonds, Ala347 at the tetramer interface or Thr329 were more destabilizing (4.1-5.7 kcal/mol). Strongly destabilizing (8.8- 11.7 kcal/mol) substitutions included those of Met340 at the tetramer interface and Phe328, Arg337 and Phe338 involved peripherally in the hydrophobic core. Truncation of any of the three residues involved centrally in the hydrophobic core of each primary dimer either prevented folding (Ile332) or allowed folding only at high protein concentration or low temperature (Leu330 and Phe341). Nine hydrophobic residues per monomer constitute critical determinants for the stability and oligomerization status of this p53 domain.     

Comparison of atomic solvation parametric sets: applicability and limitations in protein folding and binding.

Atomic solvation parameters (ASP) are widely used to estimate the solvation contribution to the thermodynamic stability of proteins as well as the free energy of association for protein-ligand complexes. They are also included in several molecular mechanics computer programs. In this work, a total of eight atomic solvation parametric sets has been employed to calculate the solvation contribution to the free energy of folding delta Gs for 17 proteins. A linear correlation between delta Gs and the number of residues in each protein was found for each ASP set. The calculations also revealed a great variety in the absolute value and in the sign of delta Gs values such that certain ASP sets predicted the unfolded state to be more stable than the folded, whereas others yield precisely the opposite. Further, the solvation contribution to the free energy of association of helix pairs and to the disassociation of loops (connection between secondary structural elements in proteins) from the protein tertiary structures were computed for each of the eight ASP sets and discrepancies were evident among them.     

Magnitude of the hydrophobic effect at central versus peripheral sites in protein-protein interfaces

Hydrophobic interactions are essential for stabilizing protein-protein complexes, whose interfaces generally consist of a central cluster of hot spot residues surrounded by less important peripheral residues. According to the O-ring hypothesis, a condition for high affinity binding is solvent exclusion from interacting residues. This hypothesis predicts that the hydrophobicity at the center is significantly greater than at the periphery, which we estimated at 21 cal mol(-1) A(-2). To measure the hydrophobicity at the center, structures of an antigen-antibody complex where a buried phenylalanine was replaced by smaller hydrophobic residues were determined. By correlating structural changes with binding free energies, we estimate the hydrophobicity at this central site to be 46 cal mol(-1) A(-2), twice that at the periphery. This context dependence of the hydrophobic effect explains the clustering of hot spots at interface centers and has implications for hot spot prediction and the design of small molecule inhibitors.     

Studies of protein-protein interfaces: a statistical analysis of the hydrophobic effect.

Data sets of 362 structurally nonredundant protein-protein interfaces and of 57 symmetry-related oligomeric interfaces have been used to explore whether the hydrophobic effect that guides protein folding is also the main driving force for protein-protein associations. The buried nonpolar surface area has been used to measure the hydrophobic effect. Our analysis indicates that, although the hydrophobic effect plays a dominant role in protein-protein binding, it is not as strong as that observed in the interior of protein monomers. Comparison of interiors of the monomers with those of the interfaces reveals that, in general, the hydrophobic amino acids are more frequent in the interior of the monomers than in the interior of the protein-protein interfaces. On the other hand, a higher proportion of charged and polar residues are buried at the interfaces, suggesting that hydrogen bonds and ion pairs contribute more to the stability of protein binding than to that of protein folding. Moreover, comparison of the interior of the interfaces to protein surfaces indicates that the interfaces are poorer in polar/charged than the surfaces and are richer in hydrophobic residues. The interior of the interfaces appears to constitute a compromise between the stabilization contributed by the hydrophobic effect on the one hand and avoiding patches on the protein surfaces that are too hydrophobic on the other. Such patches would be unfavorable for the unassociated monomers in solution. We conclude that, although the types of interactions are similar between protein-protein interfaces and single-chain proteins overall, the contribution of the hydrophobic effect to protein-protein associations is not as strong as to protein folding. This implies that packing patterns and interatom, or interresidue, pairwise potential functions, derived from monomers, are not ideally suited to predicting and assessing ligand associations or design. These would perform adequately only in cases where the hydrophobic effect at the binding site is substantial.     

Attempts to delineate the relative contributions of changes in hydrophobicity and packing to changes in stability of ribonuclease S mutants

While the hydrophobic driving force is thought to be a major contributor to protein stability, it is difficult to experimentally dissect out its contribution to the overall free energy of folding. We have made large to small substitutions of buried hydrophobic residues at positions 8 and 13 in the peptide/protein complex, RNase-S, and have characterized the structures by X-ray crystallography. The thermodynamics of association of these mutant S peptides with S protein was measured in the presence of different concentrations of methanol and ethanol. The reduction in the strength of the hydrophobic driving force in the presence of these organic solvents was estimated from surface-tension data as well as from the dependence of the DeltaC(p) of protein/peptide binding on the alcohol concentration. The data indicated a decrease in the strength of the hydrophobic driving force of about 30-40% over a 0-30% range of the alcohol concentration. We observe that large to small substitutions destabilize the protein. However, the amount of destabilization, relative to the wild type, is independent of the alcohol concentration over the range of alcohol concentrations studied. The data clearly indicate that decreased stability of the mutants is primarily due to the loss of packing interactions rather than a reduced hydrophobic driving force and suggest a value of the hydrophobic driving force of less than 18 cal mol(-)(1) A(2).     

Theoretical studies on solvation contribution to the thermodynamic stability of mutants of lysozyme T4

Atomic solvation parameters (ASPs) are widely used to estimate the solvation contribution to the thermodynamic stability of proteins as well as the free energy of association for protein-ligand complexes. In view of discrepancies in the results of free energies of solvation of folding for various proteins obtained using different atomic solvation parameter sets, systematic studies have been carried out for the calculation of accessible surface area and the changes in free energy of solvation of folding (deltaG(s,f)) for mutants of lysozyme T4 where threonine 157 is replaced by amino acids: cysteine, aspartate, glutamate, phenylalanine, glycine, histidine, isoleucine, leucine, asparagine, arginine, serine and valine. The deviations of the calculated results from the experimental results are discussed to highlight the discrepancies in the atomic solvation parameter sets and possible reasons for them. The results are also discussed to throw light on the effect of chain free energy and hydrogen bonding on the stability of mutants. The octanol to water-based ASP sets 'Sch1' and 'EM' perform better than the vacuum to water-based ASP sets. The vacuum to water-based ASP sets 'Sch3' and 'WE' can be used to predict the stability of mutants if a proper method to calculate the hydrogen bond contribution to overall stability is in place.     

Empirical free energy calculation: comparison to calorimetric data.

An effective free energy potential, developed originally for binding free energy calculation, is compared to calorimetric data on protein unfolding, described by a linear combination of changes in polar and nonpolar surface areas. The potential consists of a molecular mechanics energy term calculated for a reference medium (vapor or nonpolar liquid), and empirical terms representing solvation and entropic effects. It is shown that, under suitable conditions, the free energy function agrees well with the calorimetric expression. An additional result of the comparison is an independent estimate of the side-chain entropy loss, which is shown to agree with a structure-based entropy scale. These findings confirm that simple functions can be used to estimate the free energy change in complex systems, and that a binding free energy evaluation model can describe the thermodynamics of protein unfolding correctly. Furthermore, it is shown that folding and binding leave the sum of solute-solute and solute-solvent van der Waals interactions nearly invariant and, due to this invariance, it may be advantageous to use a nonpolar liquid rather than vacuum as the reference medium.     

A procedure for detection and quantitation of cavity volumes proteins. Application to measure the strength of the hydrophobic driving force in protein folding

Accurate identification of cavities is important in the study of protein structure, stability, design, and ligand binding. Identification and quantitation of cavities is a nontrivial problem because most cavities are connected to the protein exterior. We describe a computational procedure for quantitating cavity volumes and apply this to derive an estimate of the hydrophobic driving force in protein folding. A grid-based Monte Carlo procedure is used to position water molecules on the surface of a protein. A Voronoi procedure is used to identify and quantitate empty space within the solvated protein. Additional cavities not detected by other existing procedures can be identified. Most of these are close to surface concavities. Residue volumes for both the interior and the surface residues as well as cavity volumes are in good agreement with volumes calculated from fully hydrated protein structures obtained from molecular dynamic simulations. We show that the loss of stability because of cavity-creating mutations correlates better with cavity volumes determined by this procedure than with cavity volumes determined by other methods. Available structural and thermodynamic data for a number of cavity-containing mutants were analyzed to obtain estimates of 26.1 cal x mol(-1) x A(-3) and 18.5 cal x mol(-1) x A(-2) for the relative contributions of cavity formation and the hydrophobic effect to the observed stability changes. The present estimate for the hydrophobic driving force is at the lower end of estimates derived from model compound studies and considerably lower than previous estimates of approximately 50 cal x mol(-1) x A(-2) derived from protein mutational data. In the absence of structural rearrangement, on average, deletion of a single methylene group is expected to result in losses in stability of 0.41 and 0.70 kcal x mol(-1) resulting from decrease in hydrophobicity and packing, respectively.     

Hydrophilic framework in proteins?

The spatial neighborhood composition of residues was determined in a 511-structure set by taking only side-chain atoms into account to generate a hydrophobicity scale. This scale is symmetrical and has been divided into seven functional groups. Hydrophobic (LIVFMCAWYG) and hydrophilic (PTHSQRNKED) residues obey an equipartition rule: not only are they found in equal proportions, but they play equivalent roles in many of their properties. The nearest neighbors of all residues are always hydrophilic. However, hydrophobic residues are mostly surrounded by other hydrophobic residues located at a peak at 3.9 Å, while hydrophilic residues show three peaks at 5.0, 6.5, and 8.0 Å, suggesting a hydrophilic structural framework. This leads us to question the importance of hydrophobic cores believed to be at the origin of protein folding.     

Estimation of the hydrophobic effect in an antigen-antibody protein-protein interface

Antigen-antibody complexes provide useful models for analyzing the thermodynamics of protein-protein association reactions. We have employed site-directed mutagenesis, X-ray crystallography, and isothermal titration calorimetry to investigate the role of hydrophobic interactions in stabilizing the complex between the Fv fragment of the anti-hen egg white lysozyme (HEL) antibody D1.3 and HEL. Crystal structures of six FvD1.3-HEL mutant complexes in which an interface tryptophan residue (V(L)W92) has been replaced by residues with smaller side chains (alanine, serine, valine, aspartate, histidine, and phenylalanine) were determined to resolutions between 1.75 and 2.00 A. In the wild-type complex, V(L)W92 occupies a large hydrophobic pocket on the surface of HEL and constitutes an energetic "hot spot" for antigen binding. The losses in apolar buried surface area in the mutant complexes, relative to wild-type, range from 25 (V(L)F92) to 115 A(2) (V(L)A92), with no significant shifts in the positions of protein atoms at the mutation site for any of the complexes except V(L)A92, where there is a peptide flip. The affinities of the mutant Fv fragments for HEL are 10-100-fold lower than that of the original antibody. Formation of all six mutant complexes is marked by a decrease in binding enthalpy that exceeds the decrease in binding free energy, such that the loss in enthalpy is partly offset by a compensating gain in entropy. No correlation was observed between decreases in apolar, polar, or aggregate (sum of the apolar and polar) buried surface area in the V(L)92 mutant series and changes in the enthalpy of formation. Conversely, there exist linear correlations between losses of apolar buried surface and decreases in binding free energy (R(2) = 0.937) as well as increases in the solvent portion of the entropy of binding (R(2) = 0.909). The correlation between binding free energy and apolar buried surface area corresponds to 21 cal mol(-1) A(-2) (1 cal = 4.185 J) for the effective hydrophobicity at the V(L)92 mutation site. Furthermore, the slope of the line defined by the correlation between changes in binding free energy and solvent entropy approaches unity, demonstrating that the exclusion of solvent from the binding interface is the predominant energetic factor in the formation of this protein complex. Our estimate of the hydrophobic contribution to binding at site V(L)92 in the D1.3-HEL interface is consistent with values for the hydrophobic effect derived from classical hydrocarbon solubility models. We also show how residue V(L)W92 can contribute significantly less to stabilization when buried in a more polar pocket, illustrating the dependence of the hydrophobic effect on local environment at different sites in a protein-protein interface.     

Thermodynamics of hydrogen bond and hydrophobic interactions in cyclodextrin complexes.

Values of K, delta G(o), delta H(o), delta S(o) and delta C(po) for the binding reaction of small organic ligands forming 1:1 complexes with either alpha- or beta-cyclodextrin were obtained by titration calorimetry from 15 degrees C to 45 degrees C. A hydrogen bond or hydrophobic interaction was introduced by adding a single functional group to the ligand. The thermodynamics of binding with and without the added group are compared to estimate the contribution of the hydrogen bond or hydrophobic interaction. A change in the environment of a functional group is required to influence the binding thermodynamics, but molecular size-dependent solute-solvent interactions have no effect. For phenolic O-H-O hydrogen bond formation, delta H(o) varies from -2 to -1.4 kcal mol(-1) from 15 degrees C to 45 degrees C, and delta C(p) is increased by 18 cal K(-1) mol(-1). The hydrophobic interaction has an opposite effect: in alpha-cyclodextrin, delta C(po) = -13.3 cal K(-1) mol(-1) per ligand -CH(2)-, identical to values found for the transfer of a -CH(2)-group from water to a nonpolar environment. At room temperature, the hydrogen bond and the -CH(2)-interaction each contribute about -600 cal mol(-1) to the stability (delta G(o)) of the complex. With increased temperature, the hydrogen bond stability decreases (i.e., hydrogen bonds "melt"), but the stability of the hydrophobic interaction remains essentially constant.     

Conformational stability and mechanism of folding of ribonuclease T1

Urea and thermal unfolding curves for ribonuclease T1 (RNase T1) were determined by measuring several different physical properties. In all cases, steep, single-step unfolding curves were observed. When these results were analyzed by assuming a two-state folding mechanism, the plots of fraction unfolded protein versus denaturant were coincident. The dependence of the free energy of unfolding, delta G (in kcal/mol), on urea concentration is given by delta G = 5.6 - 1.21 (urea). The parameters characterizing the thermodynamics of unfolding are: midpoint of the thermal unfolding curve, Tm = 48.1 degrees C, enthalpy change at Tm, delta Hm = 97 kcal/mol, and heat capacity change, delta Cp = 1650 cal/mol deg. A single kinetic phase was observed for both the folding and unfolding of RNase T1 in the transition and post-transition regions. However, two slow kinetic phases were observed during folding in the pre-transition region. These two slow phases account for about 90% of the observed amplitude, indicating that a faster kinetic phase is also present. The slow phases probably result from cis-trans isomerization at the 2 proline residues that have a cis configuration in folded RNase T1. These results suggest that RNase T1 folds by a highly cooperative mechanism with no structural intermediates once the proline residues have assumed their correct isomeric configuration. At 25 degrees C, the folded conformation is more stable than the unfolded conformations by 5.6 kcal/mol at pH 7 and by 8.9 kcal/mol at pH 5, which is the pH of maximum stability. At pH 7, the thermodynamic data indicate that the maximum conformational stability of 8.3 kcal/mol will occur at -6 degrees C.     

Determination of partial unfolding enthalpy for lysozyme upon interaction with dodecyltrimethylammonium bromide using an extended solvation model

The interactions of dodecyltrimethylammonium bromides (DTABs) with hen egg lysozyme have been investigated at pH = 7.0 and 27 degrees C in phosphate buffer by isothermal titration calorimetry. DTAB interacts endothermically and activate lysozyme. The endothermicity of the lysozyme-DTAB interaction is in marked contrast to the exothermic interactions between sodium dodecyl sulphate (SDS) and lysozyme which have been attributed to specific binding between the anionic sulphate head groups and cationic amino acid residues. The enthalpies of interaction between the cationic surfactant (DTAB) and lysozyme are dominated by the endothermic unfolding of the native structure followed by an exothermic solvation of the lysozyme-DTAB complex by the addition of extra DTAB. A new direct calorimetric method to follow protein denaturation, and the effect of surfactants on the stability of proteins was introduced. The extended solvation model was used to reproduce the enthalpies of lysozyme-DTAB interaction over the whole range of DTAB concentrations. The solvation parameters recovered from the new equation, attributed to the structural change of lysozyme and its biological activity. At low concentrations of DTAB, the binding is mainly electrostatic, with some simultaneous interaction of the hydrophobic tail with nearby hydrophobic patches on the lysozyme. These initial interactions presumably cause some protein unfolding and expose additional hydrophobic sites. The DTAB-induced denaturation enthalpy of lysozyme is 86.46 +/- 0.02 kJ mol(-1).     

Polar and nonpolar atomic environments in the protein core: Implications for folding and binding

Hydrophobic interactions are believed to play an important role in protein folding and stability. Semi-empirical attempts to estimate these interactions are usually based on a model of solvation, whose contribution to the stability of proteins is assumed to be proportional to the surface area buried upon folding. Here we propose an extension of this idea by defining an environment free energy that characterizes the environment of each atom of the protein, including solvent, polar or nonpolar atoms of the same protein or of another molecule that interacts with the protein. In our model, the difference of this environment free energy between the folded state and the unfolded (extended) state of a protein is shown to be proportional to the area buried by nonpolar atoms upon folding. General properties of this environment free energy are derived from statistical studies on a database of 82 well-refined protein structures. This free energy is shown to be able to discriminate misfolded from correct structural models, to provide an estimate of the stabilization due to oligomerization, and to predict the stability of mutants in which hydrophobic residues have been substituted by site-directed mutagenesis, provided that no large structural modifications occur. © 1994 Wiley-Liss, Inc.     

Frequencies of amino acid strings in globular protein sequences indicate suppression of blocks of consecutive hydrophobic residues

Patterns of hydrophobic and hydrophilic residues play a major role in protein folding and function. Long, predominantly hydrophobic strings of 20-22 amino acids each are associated with transmembrane helices and have been used to identify such sequences. Much less attention has been paid to hydrophobic sequences within globular proteins. In prior work on computer simulations of the competition between on-pathway folding and off-pathway aggregate formation, we found that long sequences of consecutive hydrophobic residues promoted aggregation within the model, even controlling for overall hydrophobic content. We report here on an analysis of the frequencies of different lengths of contiguous blocks of hydrophobic residues in a database of amino acid sequences of proteins of known structure. Sequences of three or more consecutive hydrophobic residues are found to be significantly less common in actual globular proteins than would be predicted if residues were selected independently. The result may reflect selection against long blocks of hydrophobic residues within globular proteins relative to what would be expected if residue hydrophobicities were independent of those of nearby residues in the sequence.