Computation of electron transfer properties for heme-containing oxidoreductases

含血红素氧化还原酶的电子转移特性的计算

• 批准号: EP/F004699/1
• 负责人: Jochen Blumberger
• 金额: $ 28.05万
• 依托单位: University of Cambridge
• 依托单位国家: 英国
• 项目类别: Research Grant
• 财政年份: 2008
• 资助国家: 英国
• 起止时间: 2008 至 无数据
• 项目状态: 已结题
• 来源: https://gtr.ukri.org/projects?ref=EP%2FF004699%2F1
• 关键词: Computation; electron; transfer; properties; heme

Atomistic computer simulations are proposed to advance aquantitative understanding of long-range electron transfer (ET) reactionsin electron transport proteins. This work is carried out to aid and complementexperimental work in this field. Existing computer simulation techniques areextended to compute key parameters that govern the rate of biologicalET reactions. These parameters are often unknown and usually ratherdifficult to measure. Using numerical methods we aim at providingquantitative estimates that can be used in electron tunneling simulations ofbiological electron carriers. Moreover, exploiting the microscopic information ofmolecular simulations, the contributions of single amino acid residues andthe surrounding solvent to the ET rate are analyzed and used to interpret effects ofprotein mutations that could guide the design of efficient biomimetics. In the short term the simulation methods are validated on simple and experimentally well characterized electron transport proteins (specific aim I). In the medium to long term redox and ET reactions in the more complex heme catalases are investigated (specific aim II). Research programme (specific aim I): Reliable experimental ET parameters areavailable for only very few biological ET reactions where donor and acceptor havea well defined structure. Among these systems are ruthenium modified cytochrome (cyt) proteins such as cyt c, myoglobin, cyt b5 and cyt b562. A wealth of structural, thermodynamic and kinetic data available for these proteins, which makes them possibly the best benchmark systems for assessment of the accuracy ofcomputer simulations. First objective is the computation of ET parametersfor ET from the heme group located inside the proteinto the ruthenium complex located at the surface of the protein. The simulationsof the proteins are carried out for models with increasing degree of complexity. Depending on the deviation with experimental data, new ideas for improvement of the simulation methodology are explored. The finite temperature motion of the different cytochromes are analyzed and used to investigate the influence of the different protein folds and heme groups on the ET parameters. Thereafter the focus of research will shift to the more complicated heme catalases described below.Research programme (specific aim II):Heme catalases prevent cells from oxidative damage by catalyzing the decomposition of hydrogen peroxide in oxygen and water. The catalytic activity is drasticallyreduced or even lost under certain conditions, if reaction intermediatecompound I undergoes one-electron reduction to compound II.Recent crystallographic and computational studies have given evidencethat the oxidized form of heme b containing catalase of H. pylori (HPC) formsthe catalytically less active compound II whereas the heme d containing catalase of P. vitale remains in theactive compound I form. Naturally, the question arises whether this difference is related to the different heme groups or to the different protein structure of the two catalases. Using the methods validated on simple cytochrome proteins (see above), we propose to calculate ET parameters for HPC and PVC and to determine the corresponding contributions of cofactor, protein and solvent.In the next step ionizable groups in the vicinity of the heme groups are identified and the ET parameters computed for ET from the ionizable group to the oxidized heme center. We hope that on the basis of these calculations one can explain the different tendencies of HPC and PVC to form the catalytically less active compound II.

原子计算机模拟被提出来促进对电子传递蛋白质中的远程电子传递（ET）反应的定量理解。这项工作是为了辅助和补充这一领域的实验工作。现有的计算机模拟技术被扩展到计算控制生物ET反应速率的关键参数。这些参数通常是未知的，通常是难以测量的。利用数值方法，我们的目标是提供可用于生物电子载流子的电子隧穿模拟的定量估计.此外，利用分子模拟的微观信息，分析单个氨基酸残基和周围溶剂对ET速率的贡献，并用于解释蛋白质突变的影响，从而指导高效仿生的设计。在短期内，模拟方法在简单的和实验上充分表征的电子传递蛋白（具体目标I）上进行了验证。在中长期的氧化还原和ET反应在更复杂的血红素过氧化氢酶进行了研究（具体目标II）。研究计划（具体目标I）：只有极少数供体和受体具有明确结构的生物ET反应可获得可靠的实验ET参数。这些系统中有钌修饰的细胞色素（cyt）蛋白，如cyt c、肌红蛋白、cyt b5和cyt b562。这些蛋白质有丰富的结构、热力学和动力学数据，这使它们可能成为评估计算机模拟准确性的最佳基准系统。第一个目标是计算ET参数，从位于蛋白质内部的血红素基团到位于蛋白质表面的钌络合物。蛋白质的模拟是针对复杂程度不断增加的模型进行的。根据与实验数据的偏差，探索了改进模拟方法的新思路。分析了不同细胞色素在有限温度下的运动规律，并研究了不同蛋白质折叠和血红素基团对ET参数的影响。研究计划（具体目标II）：血红素过氧化氢酶通过催化过氧化氢在氧气和水中的分解，防止细胞受到氧化损伤。如果反应中间体化合物I发生单电子还原，在一定条件下催化活性会急剧降低甚至丧失。最近的晶体学和计算研究表明，含血红素B的H. pylori（HPC）形成催化活性较低的化合物II，而P.vitale的含血红素的过氧化氢酶保持活性化合物I形式。当然，问题是这种差异是否与不同的血红素基团或两种过氧化氢酶的不同蛋白质结构有关。使用在简单细胞色素蛋白上验证的方法（见上文），我们建议计算HPC和PVC的ET参数，并确定辅因子、蛋白质和溶剂的相应贡献。下一步，确定血红素基团附近的可电离基团，并计算从可电离基团到氧化血红素中心的ET参数。我们希望，在这些计算的基础上，可以解释HPC和PVC形成催化活性较低的化合物II的不同倾向。

项目成果

Electronic coupling matrix elements from charge constrained density functional theory calculations using a plane wave basis set.

• DOI: 10.1063/1.3507878
• 发表时间: 2010-12
• 期刊: The Journal of chemical physics
• 作者: H. Oberhofer;J. Blumberger

Insight into the mechanism of the Ru2+-Ru3+ electron self-exchange reaction from quantitative rate calculations.

从定量速率计算洞察Ru2-Ru3电子自交换反应的机理。

• DOI: 10.1002/anie.200906455
• 发表时间: 2010
• 期刊: Angewandte Chemie (International ed. in English)
• 作者: Oberhofer H

Prediction of reorganization free energies for biological electron transfer: a comparative study of Ru-modified cytochromes and a 4-helix bundle protein.

• DOI: 10.1021/ja107876p
• 发表时间: 2010-11
• 期刊: Journal of the American Chemical Society
• 影响因子: 15
• 作者: V. Tipmanee;H. Oberhofer;Mina Park;Kwang Soo Kim;J. Blumberger