Cui Qizhi, Sulea Traian, Schrag Joseph D, Munger Christine, Hung Ming-Ni, Naïm Marwen, Cygler Miroslaw, Purisima Enrico O
Biotechnology Research Institute, National Research Council Canada, 6100 Royalmount Avenue, Montreal, Quebec, Canada.
J Mol Biol. 2008 Jun 13;379(4):787-802. doi: 10.1016/j.jmb.2008.04.035. Epub 2008 Apr 20.
Using the MP1-p14 scaffolding complex from the mitogen-activated protein kinase signaling pathway as model system, we explored a structure-based computational protocol to probe and characterize binding affinity hot spots at protein-protein interfaces. Hot spots are located by virtual alanine-scanning consensus predictions over three different energy functions and two different single-structure representations of the complex. Refined binding affinity predictions for select hot-spot mutations are carried out by applying first-principle methods such as the molecular mechanics generalized Born surface area (MM-GBSA) and solvated interaction energy (SIE) to the molecular dynamics (MD) trajectories for mutated and wild-type complexes. Here, predicted hot-spot residues were actually mutated to alanine, and crystal structures of the mutated complexes were determined. Two mutated MP1-p14 complexes were investigated, the p14(Y56A)-mutated complex and the MP1(L63A,L65A)-mutated complex. Alternative ways to generate MD ensembles for mutant complexes, not relying on crystal structures for mutated complexes, were also investigated. The SIE function, fitted on protein-ligand binding affinities, gave absolute binding affinity predictions in excellent agreement with experiment and outperformed standard MM-GBSA predictions when tested on the MD ensembles of Ras-Raf and Ras-RalGDS protein-protein complexes. For wild-type and mutant MP1-p14 complexes, SIE predictions of relative binding affinities were supported by a yeast two-hybrid assay that provided semiquantitative relative interaction strengths. Results on the MP1-mutated complex suggested that SIE predictions deteriorate if mutant MD ensembles are approximated by just mutating the wild-type MD trajectory. The SIE data on the p14-mutated complex indicated feasibility for generating mutant MD ensembles from mutated wild-type crystal structure, despite local structural differences observed upon mutation. For energetic considerations, this would circumvent costly needs to produce and crystallize mutated complexes. The sensitized protein-protein interface afforded by the p14(Y56A) mutation identified here has practical applications in screening-based discovery of first-generation small-molecule hits for further development into specific modulators of the mitogen-activated protein kinase signaling pathway.
以丝裂原活化蛋白激酶信号通路中的MP1-p14支架复合物为模型系统,我们探索了一种基于结构的计算方法,以探测和表征蛋白质-蛋白质界面处的结合亲和力热点。通过对复合物的三种不同能量函数和两种不同单结构表示进行虚拟丙氨酸扫描一致性预测来定位热点。通过将诸如分子力学广义玻恩表面积(MM-GBSA)和溶剂化相互作用能(SIE)等第一性原理方法应用于突变型和野生型复合物的分子动力学(MD)轨迹,对选定热点突变的结合亲和力进行精细预测。在此,将预测的热点残基实际突变为丙氨酸,并确定突变复合物的晶体结构。研究了两种突变的MP1-p14复合物,即p14(Y56A)突变复合物和MP1(L63A,L65A)突变复合物。还研究了不依赖突变复合物晶体结构生成突变复合物MD集合的替代方法。根据蛋白质-配体结合亲和力拟合的SIE函数给出的绝对结合亲和力预测与实验结果高度一致,并且在对Ras-Raf和Ras-RalGDS蛋白质-蛋白质复合物的MD集合进行测试时,其性能优于标准MM-GBSA预测。对于野生型和突变型MP1-p14复合物,酵母双杂交试验支持了SIE对相对结合亲和力的预测,该试验提供了半定量的相对相互作用强度。MP1突变复合物的结果表明,如果仅通过突变野生型MD轨迹来近似突变MD集合,SIE预测会变差。p14突变复合物的SIE数据表明,尽管在突变时观察到局部结构差异,但从突变的野生型晶体结构生成突变MD集合是可行的。出于能量考虑,这将避免生产和结晶突变复合物的高昂成本。此处鉴定的p14(Y56A)突变提供的敏感蛋白质-蛋白质界面在基于筛选发现第一代小分子命中物以进一步开发为丝裂原活化蛋白激酶信号通路的特异性调节剂方面具有实际应用。