Smith Graham R, Sternberg Michael J E, Bates Paul A
Biomolecular Modelling Laboratory, Cancer Research UK London Research Institute, Lincoln's Inn Fields Laboratories, 44 Lincoln's Inn Fields, London WC2A 3PX, UK.
J Mol Biol. 2005 Apr 15;347(5):1077-101. doi: 10.1016/j.jmb.2005.01.058.
We investigate the extent to which the conformational fluctuations of proteins in solution reflect the conformational changes that they undergo when they form binary protein-protein complexes. To do this, we study a set of 41 proteins that form such complexes and whose three-dimensional structures are known, both bound in the complex and unbound. We carry out molecular dynamics simulations of each protein, starting from the unbound structure, and analyze the resulting conformational fluctuations in trajectories of 5 ns in length, comparing with the structure in the complex. It is found that fluctuations take some parts of the molecules into regions of conformational space close to the bound state (or give information about it), but at no point in the simulation does each protein as whole sample the complete bound state. Subsequent use of conformations from a clustered MD ensemble in rigid-body docking is nevertheless partially successful when compared to docking the unbound conformations, as long as the unbound conformations are themselves included with the MD conformations and the whole globally rescored. For one key example where sub-domain motion is present, a ribonuclease inhibitor, principal components analysis of the MD was applied and was also able to produce conformations for docking that gave enhanced results compared to the unbound. The most significant finding is that core interface residues show a tendency to be less mobile (by size of fluctuation or entropy) than the rest of the surface even when the other binding partner is absent, and conversely the peripheral interface residues are more mobile. This surprising result, consistent across up to 40 of the 41 proteins, suggests different roles for these regions in protein recognition and binding, and suggests ways that docking algorithms could be improved by treating these regions differently in the docking process.
我们研究了溶液中蛋白质的构象波动在多大程度上反映了它们形成二元蛋白质 - 蛋白质复合物时所经历的构象变化。为此,我们研究了一组41种能形成此类复合物且其三维结构已知的蛋白质,包括复合物中的结合态和未结合态。我们从每种蛋白质的未结合结构开始进行分子动力学模拟,并分析长度为5纳秒的轨迹中产生的构象波动,与复合物中的结构进行比较。结果发现,波动会使分子的某些部分进入接近结合态的构象空间区域(或给出有关该区域的信息),但在模拟过程中的任何时刻,每种蛋白质整体都不会完全处于结合态。然而,与对接未结合构象相比,在刚体对接中使用聚类MD系综中的构象仍取得了部分成功,只要未结合构象本身与MD构象一起包含在内并进行全局重新评分。对于存在子域运动的一个关键例子,即核糖核酸酶抑制剂,对MD进行了主成分分析,并且与未结合构象相比,也能够产生用于对接的构象,从而得到更好的结果。最显著的发现是,即使没有其他结合伴侣,核心界面残基的流动性(按波动大小或熵衡量)也往往比表面的其他部分小,相反,外围界面残基的流动性更大。这一令人惊讶的结果在41种蛋白质中的多达40种中都是一致的,表明这些区域在蛋白质识别和结合中具有不同的作用,并暗示了对接算法在对接过程中通过对这些区域进行不同处理而得以改进的方法。