Wei Donghui, Tang Mingsheng, Zhan Chang-Guo
Department of Chemistry, Zhengzhou University, Daxue Road, Zhengzhou, Henan 450052, China.
Org Biomol Chem. 2015 Jun 28;13(24):6857-65. doi: 10.1039/c5ob00737b.
In this study, molecular dynamics (MD) simulations and first-principles quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations have been performed to uncover the fundamental reaction pathway of proteasome with a representative inhibitor syringolin A (SylA). The calculated results reveal that the reaction process consists of three steps. The first step is a proton transfer process, activating Thr1-O(γ) directly by Thr1-N(z) to form a zwitterionic intermediate. The next step is a nucleophilic attack on the olefin carbon of SylA by the negatively charged Thr1-O(γ) atom. The last step is a proton transfer from Thr1-N(z) to another olefin carbon of SylA to complete the inhibition reaction process. The calculated free energy profile demonstrates that the second step should be the rate-determining step and has the highest free energy barrier of 24.6 kcal mol(-1), which is reasonably close to the activation free energy (∼22.4-23.0 kcal mol(-1)) derived from the available experimental kinetic data. In addition, our computational results indicate that no water molecule can assist the rate-determining step, since the second step is not involved in a proton transfer process. The obtained mechanistic insights should be valuable for understanding the inhibition process of proteasome by SylA and structurally related inhibitors at a molecular level, and thus provide a solid mechanistic base and valuable clues for future rational design of novel, more potent inhibitors of proteasome.
在本研究中,已进行分子动力学(MD)模拟和第一性原理量子力学/分子力学自由能(QM/MM-FE)计算,以揭示蛋白酶体与代表性抑制剂丁香脂素A(SylA)的基本反应途径。计算结果表明,反应过程包括三个步骤。第一步是质子转移过程,由Thr1-N(z)直接激活Thr1-O(γ)形成两性离子中间体。下一步是带负电荷的Thr1-O(γ)原子对SylA的烯烃碳进行亲核攻击。最后一步是质子从Thr1-N(z)转移到SylA的另一个烯烃碳上,以完成抑制反应过程。计算得到的自由能剖面图表明,第二步应为速率决定步骤,其自由能垒最高,为24.6 kcal mol(-1),这与从现有实验动力学数据得出的活化自由能(约22.4 - 23.0 kcal mol(-1))相当接近。此外,我们的计算结果表明,由于第二步不涉及质子转移过程,因此没有水分子可以协助速率决定步骤。所获得的机理见解对于在分子水平上理解SylA和结构相关抑制剂对蛋白酶体的抑制过程应具有重要价值,从而为未来合理设计新型、更有效的蛋白酶体抑制剂提供坚实的机理基础和有价值的线索。