Bjelic Sinisa, Brandsdal Bjørn O, Aqvist Johan
Department of Cell and Molecular Biology, Uppsala University Biomedical Center, Box 596, SE-751 24 Uppsala, Sweden.
Biochemistry. 2008 Sep 23;47(38):10049-57. doi: 10.1021/bi801177k. Epub 2008 Aug 30.
A major issue for organisms living at extreme temperatures is to preserve both stability and activity of their enzymes. Cold-adapted enzymes generally have a reduced thermal stability, to counteract freezing, and show a lower enthalpy and a more negative entropy of activation compared to mesophilic and thermophilic homologues. Such a balance of thermodynamic activation parameters can make the reaction rate decrease more linearly, rather than exponentially, as the temperature is lowered, but the structural basis for rate optimization toward low working temperatures remains unclear. In order to computationally address this problem, it is clear that reaction simulations rather than standard molecular dynamics calculations are needed. We have thus carried out extensive computer simulations of the keto-enol(ate) isomerization steps in differently adapted citrate synthases to explore the structure-function relationships behind catalytic rate adaptation to different temperatures. The calculations reproduce the absolute rates of the psychrophilic and mesophilic enzymes at 300 K, as well as the lower enthalpy and more negative entropy of activation of the cold-adapted enzyme, where the latter simulation result is obtained from high-precision Arrhenius plots. The overall catalytic effect originates from electrostatic stabilization of the transition state and enolate and the reduction of reorganization free energy. The simulations, however, show psychrophilic, mesophilic, and hyperthermophilic citrate synthases to have increasingly stronger electrostatic stabilization of the transition state, while the energetic penalty in terms of internal protein interactions follows the reverse order with the cold-adapted enzyme having the most favorable energy term. The lower activation enthalpy and more negative activation entropy observed for cold-adapted enzymes are found to be associated with a decreased protein stiffness. The origin of this effect is, however, not localized to the active site but to other regions of the protein structure.
对于生活在极端温度下的生物体来说,一个主要问题是如何保持其酶的稳定性和活性。与嗜温和嗜热的同源酶相比,冷适应酶通常热稳定性降低,以应对结冰情况,并且活化焓较低,活化熵更负。这种热力学活化参数的平衡会使反应速率在温度降低时更呈线性而非指数下降,但针对低温工作温度进行速率优化的结构基础仍不清楚。为了通过计算解决这个问题,显然需要进行反应模拟而非标准的分子动力学计算。因此,我们对不同适应性的柠檬酸合酶中的酮 - 烯醇(酸根)异构化步骤进行了广泛的计算机模拟,以探索催化速率适应不同温度背后的结构 - 功能关系。这些计算重现了嗜冷和嗜温酶在300 K时的绝对速率,以及冷适应酶较低的活化焓和更负的活化熵,后者的模拟结果是通过高精度阿伦尼乌斯图获得的。整体催化效应源于过渡态和烯醇酸根的静电稳定以及重组自由能的降低。然而,模拟结果表明,嗜冷、嗜温和超嗜热柠檬酸合酶对过渡态的静电稳定作用越来越强,而从内部蛋白质相互作用角度来看的能量代价则相反,冷适应酶具有最有利的能量项。冷适应酶观察到的较低活化焓和更负的活化熵与蛋白质刚性降低有关。然而,这种效应的起源并非局限于活性位点,而是蛋白质结构的其他区域。