Castillo R, Oliva M, Martí S, Moliner V
Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castelló, Spain.
J Phys Chem B. 2008 Aug 14;112(32):10012-22. doi: 10.1021/jp8025896. Epub 2008 Jul 23.
A theoretical study of the hydride transfer between formate anion and nicotinamide adenine dinucleotide (NAD(+)) catalyzed by the enzyme formate dehydrogenase (FDH) has been carried out by a combination of two hybrid quantum mechanics/molecular mechanics techniques: statistical simulation methods and internal energy minimizations. Free energy profiles, obtained for the reaction in the enzyme active site and in solution, allow obtaining a comparative analysis of the behavior of both condensed media. Moreover, calculations of the reaction in aqueous media can be used to probe the dramatic differences between reactants state in the enzyme active site and in solution. The results suggest that the enzyme compresses the substrate and the cofactor into a conformation close to the transition structure by means of favorable interactions with the amino acid residues of the active site, thus facilitating the relative orientation of donor and acceptor atoms to favor the hydride transfer. Moreover, a permanent field created by the protein reduces the work required to reach the transition state (TS) with a concomitant polarization of the cofactor that would favor the hydride transfer. In contrast, in water the TS is destabilized with respect to the reactant species because the polarity of the solute diminishes as the reaction proceeds, and consequently the reaction field, which is created as a response to the change in the solute polarity, is also decreased. Therefore protein structure is responsible of both effects; substrate preorganization and TS stabilization thus diminishing the activation barrier. Because of the electrostatic features of the catalyzed reaction, both media preferentially stabilize the ground-state, thus explaining the small rate constant enhancement of this enzyme, but FDH does so to a much lower extent than aqueous solution. Finally, a good agreement between experimental and theoretical kinetic isotope effects is found, thus giving some credit to our results.
通过两种混合量子力学/分子力学技术的结合,即统计模拟方法和内能最小化方法,对甲酸脱氢酶(FDH)催化的甲酸根阴离子与烟酰胺腺嘌呤二核苷酸(NAD(+))之间的氢化物转移进行了理论研究。通过在酶活性位点和溶液中进行反应得到的自由能分布,能够对两种凝聚介质的行为进行比较分析。此外,在水介质中进行的反应计算可用于探究酶活性位点和溶液中反应物状态的显著差异。结果表明,酶通过与活性位点的氨基酸残基形成有利的相互作用,将底物和辅因子压缩成接近过渡结构的构象,从而促进供体和受体原子的相对取向,有利于氢化物转移。此外,蛋白质产生的永久场降低了达到过渡态(TS)所需的功,同时辅因子发生极化,这有利于氢化物转移。相比之下,在水中,相对于反应物物种,过渡态是不稳定的,因为随着反应的进行,溶质的极性降低,因此作为对溶质极性变化的响应而产生的反应场也减小。因此,蛋白质结构对这两种效应都有影响;底物预组织和过渡态稳定化从而降低了活化能垒。由于催化反应的静电特性,两种介质都优先稳定基态,从而解释了该酶的速率常数提高幅度较小的原因,但FDH的稳定程度远低于水溶液。最后,实验和理论动力学同位素效应之间取得了良好的一致性,从而为我们的结果提供了一定的可信度。