Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, Jilin, People's Republic of China.
Department of Hepato-Biliary-Pancreatic Surgery, The Second Hospital of Jilin University, Changchun 130041, Jilin, People's Republic of China.
J Phys Chem B. 2022 Jul 7;126(26):4914-4924. doi: 10.1021/acs.jpcb.2c01492. Epub 2022 Jun 22.
Protein-protein interaction plays an important role in the development of almost all cells. Elucidating the dynamic binding and affinity of a protein-protein complex is essential for understanding the biological functions of proteins. EphA6 and Odin proteins are members of the Eph (erythropoietin-producing hepatocyte) receptor family and the Anks (ankyrin repeat and sterile α motif domain-containing) family, respectively. Odin significantly functions in regulating endocytosis, degradation, and stability of EphA receptors. In this work, the key residues of the interaction interface were determined through a hydrogen bond, root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), and dynamic correlation analysis of the conventional molecular dynamics (MD) simulations. The calculated standard binding free energy, -7.92 kcal/mol, between EphA6 and Odin is quite consistent with the experimental measurement value, -8.73 kcal/mol. By the combination of several MD simulation techniques, our investigation of the binding process reveals the detailed representative characteristics of the entire binding pathway at the molecular level. Based on the obtained potential of the mean force (PMF) curve, the analysis of the simulation trajectories shows that the residue Arg1013 in the receptor EphA6 is responsible for capturing Asp739 and Asp740 in the ligand Odin during the initial stage of binding. In the later stage of binding, the hydrogen bonds and salt bridges between a series of residues Lys973, Leu1007, Gly1009, His1010, and Arg1012 in the receptor and residues Leu735, Asn736, Asp739, Asp740, and Asp753 in the ligand mainly contribute to the stability of the protein complex. In addition, the specific change process of the receptor-ligand-binding mode is also clarified during the binding process. Our present simulation will promote a deep understanding of the protein-protein interaction, and the identified key interresidue interaction will be theoretical guidance for the design of protein drugs.
蛋白质-蛋白质相互作用在几乎所有细胞的发育中都起着重要作用。阐明蛋白质-蛋白质复合物的动态结合和亲和力对于理解蛋白质的生物学功能至关重要。EphA6 和 Odin 蛋白分别是 Eph(促红细胞生成素产生肝细胞)受体家族和 Anks(锚蛋白重复和无菌α基序结构域)家族的成员。Odin 显著调节 EphA 受体的内吞作用、降解和稳定性。在这项工作中,通过氢键、均方根偏差(RMSD)、均方根波动(RMSF)和传统分子动力学(MD)模拟的动态相关分析确定了相互作用界面的关键残基。计算得出 EphA6 和 Odin 之间的标准结合自由能为-7.92 kcal/mol,与实验测量值-8.73 kcal/mol 非常吻合。通过几种 MD 模拟技术的结合,我们对结合过程的研究揭示了整个结合途径在分子水平上的详细代表性特征。基于获得的平均力势能(PMF)曲线分析模拟轨迹表明,受体 EphA6 中的残基 Arg1013 在结合的初始阶段负责捕获配体 Odin 中的残基 Asp739 和 Asp740。在结合的后期阶段,受体中一系列残基 Lys973、Leu1007、Gly1009、His1010 和 Arg1012 与配体中残基 Leu735、Asn736、Asp739、Asp740 和 Asp753 之间的氢键和盐桥主要有助于蛋白质复合物的稳定性。此外,在结合过程中还阐明了受体-配体结合模式的具体变化过程。我们目前的模拟将促进对蛋白质-蛋白质相互作用的深入理解,并且确定的关键残基相互作用将为蛋白质药物的设计提供理论指导。