Futera Zdenek, Klenko Julia, Sponer Judit E, Sponer Jirí, Burda Jaroslav V
Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 121 16 Prague 2, Czech Republic.
J Comput Chem. 2009 Sep;30(12):1758-70. doi: 10.1002/jcc.21179.
Piano stool ruthenium complexes of the composition Ru(II)(eta6-arene)(en)Cl (en = ethylenediamine) represent an emerging class of cisplatin-analogue anticancer drug candidates. In this study, we use computational quantum chemistry to characterize the structure, stability and reactivity of these compounds. All these structures were optimized at DFT(B3LYP)/6-31G(d) level and their single point properties were determined by the MP2/6-31++G(2df,2pd) method. Thermodynamic parameters and rate constants were determined for the aquation process, as a replacement of the initial chloro ligand by water and subsequent exchange reaction of aqua ligand by nucleobases. The computations were carried out at several levels of DFT and ab initio theories (B3LYP, MP2 and CCSD) utilizing a range of bases sets (from 6-31G(d) to aug-cc-pVQZ). Excellent agreement with experimental results for aquation process was obtained at the CCSD level and reasonable match was achieved also with the B3LYP/6-31++G(2df,2pd) method. This level was used also for nucleobase-water exchange reaction where a smaller rate constant for guanine exchange was found in comparison with adenine. Although adenine follows a simple replacement mechanism, guanine complex passes by a two-step mechanism. At first, Ru-O6(G) adduct is formed, which is transformed through a chelate TS2 to the Ru-N7(G) final complex. In case of guanine, the exchange reaction is more favorable thermodynamically (releasing in total by about 8 kcal/mol) but according to our results, the rate constant for guanine substitution is slightly smaller than the analogous constant in adenine case when reaction course from local minimum is considered.
组成式为Ru(II)(η6 - 芳烃)(乙二胺)Cl(乙二胺 = en)的钢琴凳钌配合物代表了一类新兴的顺铂类似物抗癌药物候选物。在本研究中,我们使用计算量子化学来表征这些化合物的结构、稳定性和反应活性。所有这些结构均在DFT(B3LYP)/6 - 31G(d)水平上进行了优化,并通过MP2/6 - 31++G(2df,2pd)方法确定了它们的单点性质。测定了水合过程的热力学参数和速率常数,即初始氯配体被水取代以及随后水合配体与核碱基的交换反应。计算是在多个DFT和从头算理论水平(B3LYP、MP2和CCSD)上进行的,使用了一系列基组(从6 - 31G(d)到aug - cc - pVQZ)。在CCSD水平上获得了与水合过程实验结果的极佳一致性,并且使用B3LYP/6 - 31++G(2df,2pd)方法也实现了合理匹配。该水平也用于核碱基 - 水交换反应,其中发现鸟嘌呤交换的速率常数比腺嘌呤的小。尽管腺嘌呤遵循简单的取代机制,但鸟嘌呤配合物通过两步机制进行。首先,形成Ru - O6(G)加合物,其通过螯合过渡态TS2转化为Ru - N7(G)最终配合物。就鸟嘌呤而言,交换反应在热力学上更有利(总共释放约8千卡/摩尔),但根据我们的结果,当考虑从局部最小值开始的反应过程时,鸟嘌呤取代的速率常数略小于腺嘌呤情况中的类似常数。