de Hatten Xavier, Cournia Zoe, Huc Ivan, Smith Jeremy C, Metzler-Nolte Nils
Department for Chemistry and Biochemistry, University of Bochum, Universitätstrasse 150, 44809 Bochum, Germany.
Chemistry. 2007;13(29):8139-52. doi: 10.1002/chem.200700358.
The increasing importance of hydrogenase enzymes in the new energy research field has led us to examine the structure and dynamics of potential hydrogenase mimics, based on a ferrocene-peptide scaffold, using molecular dynamics (MD) simulations. To enable this MD study, a molecular mechanics force field for ferrocene-bearing peptides was developed and implemented in the CHARMM simulation package, thus extending the usefulness of the package into peptide-bioorganometallic chemistry. Using the automated frequency-matching method (AFMM), optimized intramolecular force-field parameters were generated through quantum chemical reference normal modes. The partial charges for ferrocene were derived by fitting point charges to quantum-chemically computed electrostatic potentials. The force field was tested against experimental X-ray crystal structures of dipeptide derivatives of ferrocene-1,1'-dicarboxylic acid. The calculations reproduce accurately the molecular geometries, including the characteristic C2-symmetrical intramolecular hydrogen-bonding pattern, that were stable over 0.1 micros MD simulations. The crystal packing properties of ferrocene-1-(D)alanine-(D)proline-1'-(D)alanine-(D)proline were also accurately reproduced. The lattice parameters of this crystal were conserved during a 0.1 micros MD simulation and match the experimental values almost exactly. Simulations of the peptides in dichloromethane are also in good agreement with experimental NMR and circular dichroism (CD) data in solution. The developed force field was used to perform MD simulations on novel, as yet unsynthesized peptide fragments that surround the active site of [Ni-Fe] hydrogenase. The results of this simulation lead us to propose an improved design for synthetic peptide-based hydrogenase models. The presented MD simulation results of metallocenes thereby provide a convincing validation of our proposal to use ferrocene-peptides as minimal enzyme mimics.
氢化酶在新能源研究领域的重要性日益凸显,这促使我们基于二茂铁 - 肽支架,运用分子动力学(MD)模拟来研究潜在氢化酶模拟物的结构与动力学。为开展此MD研究,我们开发了一种用于含二茂铁肽的分子力学力场,并将其应用于CHARMM模拟软件包中,从而将该软件包的应用范围拓展至肽 - 生物有机金属化学领域。采用自动频率匹配方法(AFMM),通过量子化学参考正则模式生成了优化的分子内力场参数。二茂铁的部分电荷通过将点电荷拟合到量子化学计算的静电势来推导。该力场针对二茂铁 - 1,1'- 二羧酸二肽衍生物的实验X射线晶体结构进行了测试。计算结果准确再现了分子几何结构,包括在0.1微秒MD模拟中稳定的特征性C2对称分子内氢键模式。二茂铁 - 1 - (D)丙氨酸 - (D)脯氨酸 - 1'- (D)丙氨酸 - (D)脯氨酸的晶体堆积性质也得到了准确再现。在0.1微秒MD模拟过程中,该晶体的晶格参数保持不变,且与实验值几乎完全匹配。二氯甲烷中肽的模拟结果也与溶液中的实验核磁共振(NMR)和圆二色性(CD)数据高度吻合。所开发的力场用于对围绕[Ni - Fe]氢化酶活性位点的新型、尚未合成的肽片段进行MD模拟。该模拟结果使我们提出了一种改进的基于合成肽的氢化酶模型设计。二茂金属的MD模拟结果有力地验证了我们使用二茂铁 - 肽作为最小化酶模拟物的提议。