查尔酮合成酶催化机理中的分子力场参数拟合

【摘要】：査尔酮合成酶广泛存在于植物体内，对体系内其它酶的形成起关键作用，其催化机理的研究已是理论化学领域的热点。查尔酮合成酶的催化性能与查尔酮分子有关,对查尔酮进行分子力场优化能够为后继的分子动力学计算提供基础。 本文利用GAUSSIAN03程序、CHARMM27力场,研究了查尔酮分子及模型化合物F的力场,从几何构型、原子电荷、键长、键角、二面角方面进行优化,对与水相互作用、振动频率、势能曲线进行了分析。主要围绕以下几方面展开: (1)拆分目标分子确定模型化合物F,研究了MP2/6-31G(d)水平下的几何构型,调整了新型原子的电荷并对F与水的相互作用进行了计算,结果显示:O2与水分子相互作用最强,H6与水的相互作用最弱,空间位阻等因素导致H6的作用距离比H7的作用距离小,O2的作用距离最小。 (2)对模型化合物F的振动频率进行了计算,在力场中优化平衡状态下的键长、键角,得到了与QM数据一致的键长、键角;调整了键长、键角、二面角的力常数,得到了2组与QM数据比较符合的振动频率。从0度到360度,间隔15度进行势能扫描,含双键的二面角扭转得到了势能面1,另一个得到势能面2。调整二面角参数,结果显示连接羰基碳原子与烯烃碳原子的键的扭转对势能面2的影响较大,其它二面角参数对势能面1的影响较大,得到的势能面1与QM结果趋势一致,而势能面2与QM结果完全一致。优化的过程表明:电荷的变化对相互作用影响大,键长和键角项主要影响振动频率,二面角项对势能变化起关键作用。 (3)建立了查尔酮分子力场,模型化合物两端的甲基被苯酚替代,在整体电荷为零的情况下对电荷进行微小调整,重新优化了分子各二面角项,发现连接羰基碳原子和苯酚的键、连接烯烃碳原子和苯酚的键的扭转对整体势能的影响较大,研究了其它12组二面角的各项参数,得到了与QM趋势一致的势能优化曲线。