The role of the basis set and the level of quantum mechanical theory in the prediction of the structure and reactivity of cisplatin

Abstract

In this article, we conducted an extensive ab initio study on the importance of the level of theory and the basis set for theoretical predictions of the structure and reactivity of cisplatin [cis-diamminedichloroplatinum(II) (cDDP)]. Initially, the role of the basis set for the Pt atom was assessed using 24 different basis sets, including three all-electron basis sets (ABS). In addition, a modified all-electron double zeta polarized basis set (mDZP) was proposed by adding a set of diffuse d functions onto the existing DZP basis set. The energy barrier and the rate constant for the first chloride/water exchange ligand process, namely, the aquation reaction, were taken as benchmarks for which reliable experimental data are available. At the B3LYP/mDZP/6-31+G(d) level (the first basis set is for Pt and the last set is for all of the light atoms), the energy barrier was 22.8 kcal mol(-1), which is in agreement with the average experimental value, 22.9 ± 0.4 kcal mol(-1). For the other accessible ABS (DZP and ADZP), the corresponding values were 15.4 and 24.5 kcal mol(-1), respectively. The ADZP and mDZP are notably similar, raising the importance of diffuse d functions for the prediction of the kinetic properties of cDDP. In this article, we also analyze the ligand basis set and the level of theory effects by considering 36 basis sets at distinct levels of theory, namely, Hartree-Fock, MP2, and several DFT functionals. From a survey of the data, we recommend the mPW1PW91/mDZP/6-31+G(d) or B3PW91/mDZP/6-31+G(d) levels to describe the structure and reactivity of cDDP and its small derivatives. Conversely, for large molecules containing a cisplatin motif (for example, the cDDP-DNA complex), the lower levels B3LYP/LANL2DZ/6-31+G(d) and B3LYP/SBKJC-VDZ/6-31+G(d) are suggested. At these levels of theory, the predicted energy barrier was 26.0 and 25.9 kcal mol(-1), respectively, which is only 13% higher than the actual value.