非血红素铁酶及分子间相互作用力场的研究

摘要

本论文的研究工作分为两个方面，一是单核非血红素酶的结构和催化性能的研究。单核非血红素酶是能够活化分子氧的铁氧化酶的一种，在生物、环境以及医学领域具有广泛的应用和重要的价值。非血红素酶中的一个常见的保守结构单元是2-组氨酸-1-羰基（2-His-1-carboxyl ate）三元组。我们通过分子动力学模拟和QM/MM计算从分子水平和电子水平上研究了羟乙基膦酸双加氧酶（HEPD）酶的结构及催化反应机理，分析了酶催化反应的影响因素。二是分子间弱相互作用的研究。分子间弱相互作用（比如卤键，氢键，硫键，磷键，cation-π等）在化学，生命科学和材料科学中占有重要地位。用分子动力学方法研究生物大分子和功能材料大分子的结构，形貌和性能时需要能恰当描述分子间弱相互作用的力场。但目前已有的分子力场不能合理描述分子间弱相互作用。为此我们设计了可极化椭球力场(PEff)，用来描述具有方向性的分子间弱相互作用，并用我们所建造的力场研究了生物大分子中的卤键相互作用。我们的研究结果也为设计具有方向性的配位相互作用（比如Zn酶，Fe酶中配位相互作用）的力场提供了重要思路。
　　 本论文的主要研究内容和结论如下:
　　 1.HEPD酶是一种典型的非血红素铁酶，它可以催化草铵膦(PTT)和磷霉素(Fosfomycin)等生物合成中的关键步骤。我们采用QM/MM方法，研究了HEPD酶的催化反应机理。QM/MM的计算模拟结果表明，HEPD酶中存在着不同的反应路径，整个反应过程分为四个步骤，包括氢转移，O-O键的断裂，C-C键的断裂，和产物HMP的生成。我们的研究结果表明，谷氨酸残基Glu176可以在一定程度上通过质子转移过程参与催化反应。对于整个反应过程来说，不同的反应路径具有一定的可调节性，并且它们与蛋白质环境有明显的相关性。我们提出的反应机理在很大程度上解释一些实验上观察到的比较矛盾的现象，修正了已有的羟基过氧化机理(hydroperoxylation)和羟基化机理(hydroxyl ation)，这些结果为进一步理解和调节非血红素铁酶中的催化反应机理提供了一定的指导。
　　 2.目前的一些实验研究表明，水分子在HEPD的催化反应过程中扮演着非常重要的角色。分子动力学模拟结果也表明，活性口袋附近存在一些运动受限的水分子。我们通过QM/MM方法研究了水分子参与的催化反应机理，水分子可以把羟基自由基(·OH)束缚在活性口袋中，并且可以把羟基自由基通过特定的反应转换成铁的氢氧化物（Fe-OH）。这个反应路径对于防止羟基自由基对酶的活性损伤有一定的意义。我们还发现，水分子可以作为非血红素铁酶活性中心的氧原子的一个来源，我们的计算结果表明，最终产物HMP中的氧元素大约有40％来自水分子，因此，我们提出的反应机理可以定量地解释O18同位素标记的实验结果。
　　 3.我们设计了一个新的力场函数，称之为可极化椭球力场(PolarizableEllipsoidal force field，PEff)用来描述方向性较强的分子间弱相互作用。卤键是一种非常重要的弱键相互作用，而目前大部分经典力场都不能正确地模拟卤键作用，鉴于这种情况，我们首先将PEff力场用于卤键研究。静电势分布的各向异性和短程作用是卤键的两个基本特征，PEff力场较好地描述了这两个特征。在PEff力场中，共价卤原子的各向异性电荷分布采用一个球形的负电荷分布和一个椭球形的正电荷电荷分布的组合来描述，而极化能的贡献则通过一个诱导偶极模型来描述。通过模拟卤键势能面获得了PEff力场参数，最后得到的力场中包含一些有确定物理意义的组成部分，包括静电作用，排斥/色散作用以及极化相互作用等。我们的研究结果表明，PEff模型基本上可以正确地再现MP2水平计算得到的卤键势能面特点。
　　 4.可极化椭球力场PEff与现在广泛应用的分子动力学模拟软件包在很大程度上兼容，因此很容易整合到现有的程序中。我们拟合得到的参数具有一定的可迁移性，同时这些参数与AMBER通用力场(GAFF)的参数是基本上兼容的。进一步地，我们采用PEff模型预测了人类组织蛋白酶L(hcatL)的卤键结构，所得结果与该蛋白酶的晶体结构非常吻合。这些结果表明，我们提出的PEff力场可以作为设计配位相互作用（比如Zn酶，Fe酶的活性中心）力场的重要基础。

目录

声明

摘要

第一章 绪论

1．1 非血红素铁酶

1．2 分子间弱相互作用

1．3 卤键

1．4 本论文的研究内容及意义

参考文献

第二章 理论基础和计算方法

2．1 量子化学概述

2．1．1 Born-Oppenhemer近似

2．1．2 分子轨道理论

2．1．3 电子相关与相关能

2．1．4 密度泛函理论

2．1．5 能量分解方法

2．1．6 分子静电势

2．1．7 复合物相互作用能计算

2．2 分子动力学模拟

2．2．1 分子动力学

2．2．2 分子力场

2．2．3 分子力场的参数化

2．3 QM／MM方法

2．3．1 QM／MM方法的计算模型

2．3．2 QM／MM方法的基本理论框架

2．3．3 QM／MM计算程序

参考文献

第三章 HEPD酶的分子动力学模拟

3．1 前言

3．2 计算方法

3．2．1 模型构建

3．2．2 分子动力学模拟

3．2．3 QM／MM分子动力学模拟

3．3 结果与讨论

3．4 结论

参考文献

第四章 HEPD酶的催化反应机理

4．1 前言

4．2 计算方法

4．2．1 团簇模型的计算

4．2．2 QM／MM计算模型

4．3 结果与讨论

4．3．1 HEPD酶的催化反应机理

4．3．2 铁超氧化物中间体的性质

4．3．3 末端和前端氢转移反应

4．3．4 羟基过氧化物中间体

4．3．5 C-C键断裂的反应机理

4．3．6 酶环境的重要作用

4．3．7 碳碳键断裂反应的可调节性

4．4 结论

参考文献

第五章 HEPD酶中水分子的催化作用

5．1 前言

5．2 计算细节

5．3 结果与讨论

5．3．1 QM／MM反应机理总结

5．3．2 酶活性中心附近受限水分子的性质

5．3．3 配体与水分子交换取代反应的检验

5．3．4 水分子参与的HEPD酶催化反应机理

5．3．5 与实验结果的比较

5．4 结论

参考文献

第六章 可极化椭球力场在卤键中的应用

6．1 卤键简介

6．2 卤键势能面的特点

6．3 卤键的可极化椭球模型

6．3．1 理论和方法

6．3．2 PEff力场参数拟合方法

6．4 结果和讨论

6．4．1 PEff模型中的参数

6．4．2 PEff模型描述卤键的效果

6．5 结论

参考文献

展望

致谢

攻读博士期间发表的论文