Density functional theory calculations on the molecular structures and vibration spectra of platinum(II) antitumor drugs

Systematic assessment of DFT methods for geometry optimization of mononuclear platinum-containing complexes

Platinum is used extensively as a catalyst for a wide variety of chemical reactions, though its scarcity and price present limitations to expansions of its use. To understand the origin of platinum's versatility-with the goals of both improving the efficiency of existing catalysts and mimicking its reactivity with more abundant metals-the mechanisms of platinum-catalyzed chemical reactions must be understood via structural and spectroscopic characterization of these catalysts under operando conditions. Such data, typically consisting of complex mixtures of species, often prove challenging to interpret, inviting the aid of chemical theory. DFT calculations in particular have proven successful at predicting structural and spectroscopic parameters of transition metal species, though a thorough investigation of how these methods perform for platinum-based complexes has yet to be undertaken. Herein, we evaluated the performance of geometry optimization for five commonly used functionals (BP86, PBE, B3LYP, PBE0, and TPSSh) in combination with various ligand basis sets, relativistic approximations, and solvation and dispersion models. We applied these DFT methods to a training set of 14 platinum-containing complexes with varying sizes, oxidation states, and number and type of ligands and determined that the best-performing method was the PBE0 functional together with the def2-TZVP basis set for the ligand atoms, the ZORA relativistic approximation, and solvation and dispersion corrections. The ability of this DFT methodology to accurately predict metrical parameters was confirmed using two case studies, most notably by comparing the DFT optimized geometry of a previously uncharacterized complex to newly collected EXAFS data, which showed excellent agreement.

Vibrational frequencies and intramolecular force constants for cisplatin: assessing the role of the platinum basis set and relativistic effects

The role of platinum basis set (PTBS) and relativistic effects for predicting the vibrational frequencies and intramolecular force constants for cisplatin are discussed. Nonrelativistic and relativistic computational protocols were built at B3LYP/PTBS/jorge-DZP/C-PCM and B3LYP-DKH2/PTBS/jorge-DZP-DKH/C-PCM levels, respectively, where 19 distinct PTBS were tested. As expected, the structural parameters were not very sensitive to the PTBS, however, the inclusion of relativistic effects improves the description of the cisplatin structure. When it comes to the vibrational frequencies, the results show that the PTBS, and mainly the relativistic effects, are both important. Moreover, the PBE0 functional led to better results than B3LYP in the protocols PBE0/LANL2TZ(f)/jorge-DZP/C-PCM (P20) and PBE0-DKH2/Sapporo-DKH3-DZP-2012/jorge-DZP-DKH/C-PCM (P22), which provided a mean absolute deviation (MAD) of only 10.8 cm^-1 and 9.5 cm^-1, respectively, for vibrational frequencies, which are excellent choices to study Pt complexes. Finally, a discussion of the intramolecular force constants for cisplatin is carried out, with the calculated bond and angles force constants with P20 and P22 protocols being recommended for the parameterization of the force field of cisplatin.

Remarkably Intricate Raman Spectra of Platinum(II)-Ligand Skeletal Modes in Diamminedihalido Complexes

More than 100 calculations of vibrational frequencies for cis-[Pt(NH₃)₂X₂] (X = Cl^-, Br^-) and trans-[Pt(NH₃)₂Cl₂] have been published over the past 25 years. The high-quality Raman spectra presented here have sufficient resolution in the crucial region of metal-ligand modes to permit an assessment of frequency calculations. The peaks for symmetric and antisymmetric Pt-Cl stretching modes ν_s and ν_a are resolved in the 312-333 cm^-1 region. Frequency variations due to ³⁵Cl and ³⁷Cl isotopes are shown to be comparable in size to the differences between ν_s and ν_a. Peaks corresponding to the two Pt-Br stretching frequencies are observed at 213 and 218 cm^-1 and Pt-N stretching frequencies for all compounds between 485 and 540 cm^-1. The crystal structures of trans-[Pt(NH₃)₂Cl₂] and cis-[Pt(NH₃)₂Br₂] at 300 and 120 K are reported and complement the variable-temperature vibrational frequencies.

Correct Modeling of Cisplatin: a Paradigmatic Case

Quantum chemistry is a useful tool in modern approaches to drug and material design, but only when the adopted model reflects a correct physical picture. Paradigmatic is the case of cis‐diaminodichloroplatinum(II), cis‐[Pt(NH₃)₂Cl₂], for which the correct simulation of the structural and vibrational properties measured experimentally still remains an open question. By using this molecule as a proof of concept, it is shown that state‐of‐the‐art quantum chemical calculations and a simple model, capturing the basic physical flavors, a cis‐[Pt(NH₃)₂Cl₂] dimer, can provide the accuracy required for interpretative purposes. The present outcomes have fundamental implications for benchmark studies aiming at assessing the accuracy of a given computational protocol.

DFT study of the molecular and crystal structure and vibrational analysis of cisplatin

DFT and periodic-DFT (PAW-PBE method, code VASP) calculations have been performed to study the structural and vibrational characteristics of cis-diamminedichloroplatinum(II) (cisplatin) at molecular and outside molecular level. To estimate the effect of the intermolecular interactions in crystal on the structural and vibrational properties of cisplatin, three theoretical models are considered in the present study: monomer (isolated molecule), hydrogen bonded dimer and periodic solid state structures. The work focused on the role of the theoretical models for correct modeling and prediction of geometrical and vibrational parameters of cisplatin. It has been found that the elaborate three-dimensional intermolecular hydrogen bonding network in the crystalline cisplatin significantly influences the structural and vibrational pattern of cisplatin and therefore the isolated cisplatin molecule is not the correct computational model regardless of the theoretical level used. To account for the whole intermolecular hydrogen bonding network in direction of both a and c axis and for more reliable calculations of structural and vibrational parameters periodic DFT calculations were carried out in the full crystalline periodic environment with the known lattice parameters for each cisplatin polymorph phase. The model calculations performed both at molecular level and for the periodic structures of alpha and beta cisplatin polymorph forms revealed the decisive role of the extended theoretical model for reliable prediction of the structural and vibrational characteristics of cisplatin. The powder diffraction pattern and the calculated IR and Raman spectra predicted beta polymorph form of our cisplatin sample freshly synthesized for the purposes of the present study using the Dhara's method. The various rotamers realized in the polymorph forms of cisplatin were explained by the low population of the large number of rotamers in solution as well as with the high rotamer interconversion rate due to the low energy barrier.

Computational study of the anticancer drug cisplatin

Systems containing platinum (Pt) are more challenging for reliable computations, because Pt has 78 electrons and requires relativistic treatment. However, to reduce computational demands, most previous researchers used effective core potentials. In this investigation, we perform numerous computations on the cisplatin molecule with ab initio methods and density functional theory, some of which involve all electron and zero-order relativistic approximation. Tentative conclusions on the reliability of various methods are drawn from comparison of our results with previous calculations and available experimental data. Vibrational and electronic spectra are calculated and compared with previous studies and available experimental data.

Molecular structure, IR spectra, and chemical reactivity of cisplatin and transplatin: DFT studies, basis set effect and solvent effect

Three different density functional theory (DFT) methods were employed to study the molecular structures of cis-diamminedichloroplatinum(II) (CDDP) and trans-diamminedichloroplatinum(II) (TDDP). The basis set effect on the structure was also investigated. By comparing the optimized structures with the experimental data, a relatively more accurate method was chosen for further study of the IR spectra and other properties as well as the solvent effect. Nineteen characteristic vibrational bands of the title compounds were assigned and compared with available experimental data. The number of characteristic peaks for the asymmetric stretching and deformation vibrations of N-H can serve as a judgment for the isomer between CDDP and TDDP. Significant solvent effect was observed on the molecular structures and IR spectra. The reduced density gradient analysis was performed to study the intramolecular interactions of CDDP and TDDP, and the nature of changes in the structures caused by the solvent was illustrated. Several descriptors determined from the energies of frontier molecular orbitals (HOMO and LUMO) were applied to describe the chemical reactivity of the title compounds. The molecular electrostatic potential (MESP) surfaces showed that the amino groups were the most favorable sites that nucleophilic reagents tend to attack, and CDDP was easier to be attacked by nucleophilic reagents than TDDP.

Assessment of new DFT methods for predicting vibrational spectra and structure of cisplatin: which density functional should we choose for studying platinum(II) complexes?

Ten different DFT methods, including several recently developed functionals have been tested for their performances in prediction of infrared and Raman spectra and molecular structure of cisplatin. The assessed DFT methods cover the range from meta-GGA to hybrid, double hybrid and long-range corrected hybrid models (M06-L, M06, M06-2X, PBE0, mPW1PW, B3LYP, B2PLYP, CAM-B3LYP, ωB97XD and LC-ωPBE). The calculated structural parameters and theoretical spectra have been compared to the corresponding experimental data. It is shown that the LC-ωPBE scheme is superior to other DFT methods in predicting the geometry of cisplatin. Unfortunately, the M06-L, M06-2X and B3LYP functionals are deficient in the evaluation of the strength of two Pt←NH3 coordination bonds in cisplatin (the calculated bond lengths are too long and the Pt-N stretching frequencies are underestimated). Both the PBE0 and mPW1PW functionals, in conjunction with the LanL2TZ(f) basis set for Pt give very similar theoretical results and seem to be the best methods for predicting the IR and Raman spectra of cisplatin. The long-range corrected functionals (LC-ωPBE, ωB97XD and CAM-B3LYP) have shown good performances in predicting the frequencies of Pt-ligand vibrations and are promising new tools for theoretical study of novel platinum(II) compounds.

Molecular structure and vibrational spectra studies on antipyrine derivative, 4-(2,3,4-trihydroxybenzylideneamino) antipyrine

The molecular geometry, theoretical harmonic frequencies and infrared intensities of 4-(2,3,4-trihydroxybenzylideneamino) antipyrine (THBAP) were calculated using different density functional methods (LSDA, mPW1PW91, B3LYP and HCTH) with various basic sets, including 3-21G, 6-311G, LanL2DZ, and SDD. The purpose of this research is to compare the performance of different density functional theory (DFT) methods at different basis sets in predicting geometry and vibration spectrum of THBAP. The optimized geometric band lengths and bond angles obtained by using mPW1PW91 at both LanL2DZ and SDD basic sets show the best agreement with the experimental data. A comparison between the observed fundamental vibrational frequencies of THBAP and the calculated results has been made and the result indicates that the mPW1PW91/6-311G level is superior to all the remaining levels for predicting all the vibration spectra on average for THBAP.

Polymorphism in cisplatin anticancer drug

This study reports a combined experimental and theoretical study of the solid-state polymorphism of the anticancer agent cisplatin. A complete assignment was performed for the inelastic neutron scattering (INS) and Raman spectra collected simultaneously for cisplatin, at different temperatures, with a view to obtain reliable and definitive evidence of the relative thermal stability of its α and β polymorphic species. A marked temperature-dependent hysteresis was observed, as previously reported in the literature. Theoretical calculations were carried out at the density functional theory level, using a plane-wave basis set approach and pseudopotentials. A detailed comparison with the experimental Raman and INS data showed that the α polymorph is present at the lowest temperatures, whereas the β form occurs near room temperature. Furthermore, regions of coexistence of both forms are identified, which depend on the working mode (heating or cooling). These findings imply that Raman spectroscopy allows clear identification of the α and β polymorphs at a given temperature and can unambiguously discriminate between them. Elucidation of the polymorphic equilibrium of this widely used anticancer drug is paramount for its pharmaceutical preparation and storage conditions.

Platinum and Palladium Polyamine Complexes as Anticancer Agents: The Structural Factor

Since the introduction of cisplatin to oncology in 1978, Pt(II) and Pd(II) compounds have been intensively studied with a view to develop the improved anticancer agents. Polynuclear polyamine complexes, in particular, have attracted special attention, since they were found to yield DNA adducts not available to conventional drugs (through long-distance intra- and interstrand cross-links) and to often circumvent acquired cisplatin resistance. Moreover, the cytotoxic potency of these polyamine-bridged chelates is strictly regulated by their structural characteristics, which renders this series of compounds worth investigating and their synthesis being carefully tailored in order to develop third-generation drugs coupling an increased spectrum of activity to a lower toxicity. The present paper addresses the latest developments in the design of novel antitumor agents based on platinum and palladium, particularly polynuclear chelates with variable length aliphatic polyamines as bridging ligands, highlighting the close relationship between their structural preferences and cytotoxic ability. In particular, studies by vibrational spectroscopy techniques are emphasised, allowing to elucidate the structure-activity relationships (SARs) ruling anticancer activity.