Hendsch Z S, Tidor B
Department of Chemistry, Massachusetts Institute of Technology, Cambridge 02139-4307, USA.
Protein Sci. 1999 Jul;8(7):1381-92. doi: 10.1110/ps.8.7.1381.
The GCN4 leucine zipper is a peptide homodimer that has been the subject of a number of experimental and theoretical investigations into the determinants of affinity and specificity. Here, we utilize this model system to investigate electrostatic effects in protein binding using continuum calculations. A particularly novel feature of the computations made here is that they provide an interaction-by-interaction breakdown of the electrostatic contributions to the free energy of docking that includes changes in the interaction of each functional group with solvent and changes in interactions between all pairs of functional groups on binding. The results show that (1) electrostatic effects disfavor binding by roughly 15 kcal/mol due to desolvation effects that are incompletely compensated in the bound state, (2) while no groups strongly stabilize binding, the groups that are most destabilizing are charged and polar side chains at the interface that have been implicated in determining binding specificity, and (3) attractive intramolecular interactions (e.g., backbone hydrogen bonds) that are enhanced on binding due to reduced solvent screening in the bound state contribute significantly to affinity and are likely to be a general effect in other complexes. A comparison is made between the results obtained in an electrostatic analysis carried out calculationally and simulated results corresponding to idealized data from a scanning mutagenesis experiment. It is shown that scanning experiments provide incomplete information on interactions and, if overinterpreted, tend to overestimate the energetic effect of individual side chains that make attractive interactions. Finally, a comparison is made between the results available from a continuum electrostatic model and from a simpler surface-area dependent solvation model. In this case, although the simpler model neglects certain interactions, on average it performs rather well.
GCN4亮氨酸拉链是一种肽同二聚体,已成为许多关于亲和力和特异性决定因素的实验和理论研究的对象。在此,我们利用这个模型系统,通过连续介质计算来研究蛋白质结合中的静电效应。这里所做计算的一个特别新颖之处在于,它们提供了对接自由能静电贡献的逐个相互作用分解,其中包括每个官能团与溶剂相互作用的变化以及结合时所有官能团对之间相互作用的变化。结果表明:(1)由于在结合状态下未完全补偿的去溶剂化效应,静电效应使结合不利,大约为15千卡/摩尔;(2)虽然没有基团强烈稳定结合,但最不稳定的基团是界面处的带电和极性侧链,这些侧链与确定结合特异性有关;(3)由于结合状态下溶剂屏蔽减少,结合时增强的有吸引力的分子内相互作用(例如主链氢键)对亲和力有显著贡献,并且可能是其他复合物中的普遍效应。将计算进行的静电分析结果与对应于扫描诱变实验理想化数据的模拟结果进行了比较。结果表明,扫描实验提供的相互作用信息不完整,如果过度解读,往往会高估产生有吸引力相互作用的单个侧链的能量效应。最后,对连续介质静电模型和更简单的依赖表面积的溶剂化模型的结果进行了比较。在这种情况下,虽然更简单的模型忽略了某些相互作用,但平均而言它表现相当不错。