物理化学学报 >> 2009, Vol. 25 >> Issue (04): 668-676.doi: 10.3866/PKU.WHXB200904271

研究论文 上一篇    下一篇

新型二氟甲基磷酸类酪氨酸蛋白磷酸酯酶1B抑制剂的分子动力学模拟和结合自由能计算

崔巍 张怀 计明娟   

  1. 中国科学院研究生院化学与化学工程学院, 北京 100049; 中国科学院研究生院计算地球动力学重点实验室, 北京 100049
  • 收稿日期:2008-10-10 修回日期:2008-12-04 发布日期:2009-03-31
  • 通讯作者: 计明娟 E-mail:jmj@gucas.ac.cn

Molecular Dynamics Simulations and Free Energy Calculations of a Novel Series of Protein Tyrosine Phosphatase 1B Difluoromethylenephosphonic Acid Inhibitors

 CUI Wei, ZHANG Huai, JI Ming-Juan   

  1. College of Chemistry and Chemical Engineering, Graduate University of the Chinese Academy of Sciences, Beijing 100049, P. R. China; Laboratory of ComputationalGeodynamics, Graduate University of the Chinese Academy of Sciences, Beijing 100049, P. R. China
  • Received:2008-10-10 Revised:2008-12-04 Published:2009-03-31
  • Contact: JI Ming-Juan E-mail:jmj@gucas.ac.cn

摘要: 通过分子对接建立了一系列含二氟甲基磷酸基团(DFMP)或二氟甲基硫酸基团(DFMS)的抑制剂与酪氨酸蛋白磷酸酯酶1B(PTP1B)的相互作用模式, 并通过1 ns的分子动力学模拟和molecular mechanics/generalized Born surface area (MM/GBSA)方法计算了其结合自由能. 计算获得的结合自由能排序和抑制剂与靶酶间结合能力排序一致; 通过基于主方程的自由能计算方法, 获得了抑制剂与靶酶残基间相互作用的信息, 这些信息显示DFMP/DFMS基团的负电荷中心与PTP1B的221位精氨酸正电荷中心之间的静电相互作用强弱决定了此类抑制剂的活性, 进一步的分析还显示位于DFMP/DFMS基团中的氟原子或其他具有适当原子半径的氢键供体原子会增进此类抑制剂与PTP1B活性位点的结合能力.

关键词: 酪氨酸蛋白磷酸酯酶1B, 分子动力学模拟, 自由能计算, 自由能分解, MM/GBSA

Abstract: Binding models for a series of difluoromethylenephosphonic (DFMP) and difluoromethylenesulfonic (DFMS) acids to protein tyrosine phosphatase 1B (PTP1B) were studied by molecular docking, molecular dynamics (MD) simulations, and free energy calculations. Binding free energies were computed using the molecular mechanics/generalized Born surface area (MM/GBSA) methodology based on 1 ns MD simulations. The order of affinities for the studied inhibitors can be accurately predicted using previously predicted binding free energies. Inhibitor/residue interaction profiles for all inhibitors were systematically generated using MM/GBSA free energy decomposition analysis. Inhibitor/residue interaction profiles demonstrated that electrostatic interactions between the negative charge center of DFMP/DFMS groups and Arg221 of PTP1B are a crucial part of the studied molecule affinities. Furthermore, the fluorine atom or other hydrogen bonding donor atoms with appropriate radii will improve inhibitor binding to the primary binding site of PTP1B.

Key words: Protein tyrosine phosphatase 1B, Molecular dynamics simulation, Free energy calculation, Free energy decomposition, MM/GBSA

MSC2000: 

  • O641