The solvent effect on the aquation processes of the cis-dichloro(ethylenediammine)platinum(II) using continuum solvation models

Interaction of 5-fluorouracil anticancer drug with nucleobases: insight from DFT, TD-DFT, and AIM calculations

In this study, the interaction of 5-fluorouracil (5FU) drug with adenine (A), guanine(G), cytosine(C), uracil (U), and thymine (T) nucleobases of DNA and RNA are surveyed at the ωB97XD/LANL2DZ, M06-2X/6-31G (d, p), MPW1PWQ1/6-31G(d, p), PBEPBE/6-31(d, p) and ωB97XD/6-31G(d, p) levels of density functional theory (DFT). The considered complexes of 5FU drug with nucleobases are optimized at the above level of theories. Max Force and RMS of optimization criteria are 0.00035 (Ha), and 0.0003 respectively. From optimized structures, the adsorption energy, thermodynamic parameters in gas and solvent media, quantum theory atom in molecule (QTAIM), electron localized function (ELF), and reduced density gradient (RDG) are calculated at ωB97XD/LANL2DZ and M06-2X/6-31G (d, p) level of DFT theory. The QTAIM, ELF, and RDG results confirm that the nature of bonding between 5FU drug with A, C, G, U, and T nucleobases is electrostatic or hydrogen bond type. The adsorption and thermodynamic energy results demonstrate that the interaction of the 5FU drug with C and G nucleobases is stronger than other nucleobases. The results of this study can be suggested the mechanism of interaction of the 5FU drug with nucleobases of DNA and RNA.Communicated by Ramaswamy H. Sarma.

The conformation effect of the diamine bridge on the stability of dinuclear platinum(II) complexes and their hydrolysis

In this paper, the hydrolysis process of a bisplatinum complex containing the flexible chain 1,6-hexanediamine between the two metal centers was investigated through the use of density functional theory (DFT) with the analysis of the role of the spacing group arrangement on the values of free energy activation barrier. All structures were fully optimized in aqueous solution using implicit model for solvent at DFT level. The energy profiles for the hydrolysis reaction were determined by using the supermolecule approach. Five transition states were proposed differing by the conformation of the bridge group, and the activation free energy calculated as a weighted average within the selected forms. The Gibbs population for reactant was used as a statistical weight leading to the predicted value of 23.1kcalmol(-1), in good accordance with experiment, 23.8kcalmol(-1). Our results suggests that for 1,6-hexanediamine bridge ligand, the extend forms with average torsional angle over the carbon chain larger than 130° have the greatest contribution to the hydrolysis kinetics. The results presented here point out that the hydrolysis mechanism might follow different paths for each conformation and each of these contributes to the observed energy barrier.

Palladium(II) Complexes Containing Mixed Nitrogen-Sulphur Donor Ligands: Interaction of [Pd(Methionine Methyl Ester)(H2O)2](2+) with Biorelevant Ligands

Pd(MME)Cl₂ complex (MME = methionine methyl ester) was synthesised and characterized by physicochemical measurements. The reaction of [Pd(MME)(H₂O)₂]²⁺ with amino acids, peptides, or dicarboxylic acids was investigated at 25°C and 0.1 M ionic strength. Amino acids and dicarboxylic acids form 1 : 1 complexes. Peptides form both 1 : 1 complexes and the corresponding deprotonated amide species. The stability of the complexes formed was determined and the binding centres of the ligands were assigned. Effect of solvent on the stability constant of Pd(MME)-CBDCA complex, taken as a representative example, shows that the complex is more favoured in a medium of low dielectric constant. The concentration distribution diagrams of the complexes were evaluated.

Theoretical analysis of trans-[PtCl2(NH3)(thiazole)] and trans-[PtCl2(thiazole)2] binding to biological targets — Factors influence binding kinetics and adduct stability

Full reaction energy profiles for trans-[PtCl2(NH3)(thiazole)] and trans-[PtCl2(thiazole)2] binding to sulfur- and nitrogen-containing biorelevant ligands were constructed by the density functional theory (DFT) method. Calculated results demonstrate that trans-platinum complexes can interact with biological targets, affording cis and trans products via very similar transition states. For different substituents, sulfur-containing ligands constitute kinetically preferred targets for platination, whereas the platination of nitrogen-containing ligands is more favorable thermodynamically. This is consistent with previous experimental studies. Calculated results also suggest that the trans effect, the influence of the ligand, the size of the ligand, and hydrogen bonding play important roles in binding kinetics and stabilizing adducts.

Cisplatin interaction with amino acids cysteine and methionine from gas phase to solutions with constant pH

This work is focused on the computational studies of reactions of hydrated forms of cisplatin with sulphur-containing amino acids cysteine and methionine. First, the appropriate model for solvation of the examined complexes was searched for. The suggested procedure employs the B3LYP density functional, 6-311++G(2df,2pd) basis set with Stuttgart-Dresden pseudopotentials on heavy atoms, the D-PCM solvation model and the UAKS cavity which uses more realistic NPA partial charges instead of formal partial charges for platinum ligands. In the second part this model is applied to the evaluation of the Legendre transformed reaction Gibbs free energy of cisplatin with cysteine and methionine in solution at constant pH.

An experimental and theoretical investigation into the hydrolysis of dichloro(ethylenediamine)platinum(II) via electrospray mass spectrometry and density functional theory

Dichloro(ethylenediamine)platinum(II), Pt(en)Cl(2), was dissolved in H(2)O and D(2)O, and the resulting aqueous solutions were electrosprayed into a quadrupole ion-trap mass spectrometer. A series of major and minor ionic hydrolysis products were detected. These ions were then subjected to collision-induced dissociation. As an aid in interpreting the experimental results, density functional theory calculations were carried out. These computations permitted the structures and energetics associated with the hydrolysis products to be determined. An understanding of the hydrolysis of PtenCl(2) and related coordination complexes is essential in the rational design of metal-based drugs.

Synthesis, characterization, DFT studies and DNA binding of mixed platinum (II) complexes containing quinoxaline and 1,2-dithiolate ligands

The complexes Pt(pq)Cl2(1) and Pt(pq)(bdt) (2) (where pq = 2-(2'pyridyl)quinoxaline and bdt=benzene-1,2-dithiolate) have been synthesized and fully characterized by UV-visible (UV-Vis), Fourier Transformer Infrared Spectra (FTIR), 1 and 2D NMR and cyclic voltammetry (CV). Interactions of the tested systems (the aforementioned complexes 1 and 2) and the free ligands pq and bdt with double stranded calf thymus DNA (CT-DNA) were studied by UV-spectrophotometric (melting curves) and circular dichroism (CD) measurements. The results suggest that both complexes 1 and 2, are able to form adducts with DNA and to distort the double helix by changing the base stacking. Complex 2 forms stronger adducts to CT-DNA than complex 1 and this is probably due to the substitution of the chlorine atoms of 1 by the 1,2-dithiolate ligand (bdt) in 2. The latter induces an extensive distortion in the planarity of 2 as density functional theory (DFT) calculations reveal. Besides, the light absorbing complex 2 possess intense mixed metal ligand to ligand charge transfer (MM'LLCT) transition in the visible region of the spectrum and could act as photoluminescent metal-based probe for the study of DNA binding. Thus, the photocleavage of DNA by 2 has been studied by UV-Vis and CD spectra and monitored by agarose gel electrophoresis. Under our experimental conditions, it is unclear that complex 2 can photocleave DNA. Furthermore, the ability of 2 to inhibit proliferation of human tumor cell lines was tested and the results indicate some cytoxytic effect on the SF-286 cells.

Multiple Proton-Transfer Reactions in DNA Base Pairs by Coordination of Pt Complex

The possibility of multiple proton-transfer reactions in DNA base pairs because of coordination of cisplatin is theoretically elucidated by density functional theory (DFT) and by quantum mechanics/molecular mechanics (QM/MM) methods with an ONIOM method. From the energetics of two base pairs with the cisplatin, it is theoretically confirmed that the Pt complex is likely to bind in the order cis-(CG)-Pt-(GC), cis-(CG)-Pt-(AT), cis-(TA)-Pt-(AT), where G, C, A, and T are guanine, cytosine, adenine, and thymine, respectively, and the Pt atom bonds to the N7 site of G and A. This result supports the experimental evidence, where the structure cis-A-Pt-A is seldom observed at room temperature. The single proton-transfer reaction occurs in one of the two GC pairs. No simultaneous single proton-transfer reaction can occur in both base pairs. Two different single proton-transferred structures (cis-(CG*)(d)-Pt-(GC)(p) and cis-(CG)(d)-Pt-(G*C)(p), where the asterisk means a proton donor of G) are as stable as the original structures (CG)(d)-Pt-(GC)(p). The same tendency was observed with cis-(CG*)-Pt-(AT). In contrast to cisplatin, multiple single proton-transfer reactions occur in the system consisting of two base pairs with transplatin. The optimized structure agrees with the experimental data for Pt-G coordination except for the hydrogen-bonding length.

Monte Carlo Simulation of Cisplatin Molecule in Aqueous Solution

The Lennard-Jones (12-6) parameters were obtained for all atoms of cisplatin molecule using the ab initio quantum mechanical potential energy surface for the water-cisplatin interaction as reference data. The parameters found were (epsilon/kcal.mol(-1) and sigma/angstroms) 1.0550, 3.6590 (Pt); 0.0381, 4.6272 (Cl); 0.0455, 3.3783 (N); and 0.0185, 0.0936 (H) and provided very good results for the description of the aqueous solution of cisplatin through Monte Carlo simulation. From statistical analysis of solute-solvent interactions, we observed that the NH3 groups are involved in 53% of the calculated hydrogen bonds with a significant contribution from chlorides (41%) and only 6% involving the Pt center. This is in agreement with the expected behavior for such molecules. Two hydration shells with 22 and 81 water molecules, respectively, centered around 4.6 and 7.3 angstroms were found from the center of mass pair correlation function analysis. The cisplatin atomic Lennard-Jones parameters are reported for the first time, and they might be useful for studying the structure, properties, and processes of cisplatin-like molecules in aqueous solution, including explicitly the solvent effect.

Modeling anticancer drug–DNA interactions via mixed QM/MM molecular dynamics simulations

The development of anticancer drugs started over four decades ago, with the serendipitous discovery of the antitumor activity of cisplatin and its successful use in the treatment of various cancer types. Despite the efforts made in unraveling the mechanism of the action of cisplatin, as well as in the rational design of new anticancer compounds, in many cases detailed structural and mechanistic information is still lacking. Many of these drugs exert their anticancer activity by covalently binding to DNA inducing a distortion or simply impeding replication, thus triggering a cellular response, which eventually leads to cell death. A detailed understanding of the structural and electronic properties of drug-DNA complexes and their mechanism of binding is the key step in elucidating the principles of their anticancer activity. At the theoretical level, the description of covalent drug-DNA complexes requires the use of state-of-the-art computer simulation techniques such as hybrid quantum/classical molecular dynamics simulations. In this review we provide a general overview on: drugs which covalently bind to DNA duplexes, the basic concepts of quantum mechanics/molecular mechanics (QM/MM), molecular dynamics methods and a list of selected applications of these simulations to the study of drug-DNA adducts. Finally, the potential and the limitations of this approach to the study of such systems are critically evaluated.

Cisplatin interaction with cysteine and methionine, a theoretical DFT study

Interactions of hydrated cisplatin complexes with sulphur-containing amino acids cysteine and methionine were explored. The square-planar cis-[Pt(NH3)2(H2O)X]+ complexes (where X=Cl- and OH-) were chosen as mono- and dihydrated reactants. Calculations using density functional theory (DFT) techniques with B3LYP functional were performed. The isolated molecules and the supermolecular approaches were employed for the determination of the reaction energies. Bond dissociation energies (BDE) were estimated in the model of isolated molecules and supermolecules were used for the determination of the association energies between the two interacting parts. Formation of monodentate complexes by replacing the aqua-ligand with the S, N, and O-sites of both amino acids represents an exothermic process. The highest BDE was found in cysteine structures for the Pt-S coordination. The bonding energy is about 114 kcal/mol, which is comparable with cisplatin-guanine adducts. Analogous BDE for methionine complexes is smaller by about 40 kcal/mol. This correlates well with the known fact that cysteine forms irreversible cisplatin adducts while similar adducts in the methionine case are reversible. The formation of chelate structures is an exothermic reaction only for the hydroxo-form of reactants in the supermolecular approach where additional association interactions between the released water and chelate molecules sufficiently stabilize the final product.

Study on Conformation Interconversion of 3-Alkyl-4-acetyl-3,4-dihydro-2H-1,4-benzoxazines from Dynamic NMR Experiments and ab Initio Density Functional Calculations

Variable-temperature NMR experiments and ab initio density functional calculations were carried out to investigate the conformation interconversion of novel chiral 3-alkyl-3,4-dihydro-2H-benzo[1,4]oxazine derivatives. With CDCl3 as the solvent, the coalescence temperatures of H2, H3, H11, and H19 of product 1 are about 289, 304, 292, and 316 K, with the corresponding activation free energies at 58.0 +/- 6.7, 60.9 +/- 7.1, 58.3 +/- 6.8, and 59.6 +/- 6.9 kJ.mol(-1), respectively. When dimethyl sulfoxide (DMSO-d6) was used as the solvent, 1H and 13C NMR signals were completely assigned at 375 K. The effects of solvent and temperature were investigated through a polarizable continuum model. At each theoretical level (MP2 or B3LYP), the changing tendencies of the calculated activation free energies and interconversion rates agree well with those of the NMR results. In addition, the interconversion rate at each specified temperature was calculated to be about 1.5 times faster in DMSO-d6 than in CDCl3. Accordingly, we failed to observe the coalescence phenomena of H3 and H19 in DMSO-d6 by NMR measurements from 296 to 375 K. The substitution effect at the R1-R5 positions was considered using density functional calculations, with the activation barriers decreasing as follows: product 6 > 3 > 1 > 7 > 2. This sequence is consistent with that of the reaction heats, except for product 7, implying that the interconversion processes may be thermodynamically controlled. Surprisingly, the substituted groups near the acetyl group in product 2 and 7 do not elevate the activation barrier but, instead, lower it somewhat, with the possible reasons for this provided in the paper.

Linear free energy relationship for 4-substituted (o-phenylenediamine)platinum(II) dichloride derivatives using quantum mechanical descriptors

In the present work quantum mechanical methods were used to calculate the rate constants for the first step of the aquation of a set of 4-substituted (o-phenylenediamine)platinum(II) dichloride derivatives containing electron-donating and withdrawing substituents at the 4-position of the aromatic ring. A linear free energy relationship was obtained for log(k(X)/k(H)), k being the rate constant for the first step of hydrolysis, and the electronic Hammett constants sigma(m) and sigma(p). The results showed that electron-donating groups promote the hydrolysis reaction. The quantitative models described here may be useful for the rational design of new, less mutagenic drugs based on platinum complexes.

Hydration process as an activation of trans- and cisplatin complexes in anticancer treatment. DFT andab initio computational study of thermodynamic and kinetic parameters

The thermodynamic and kinetic aspects of hydration reactions of cis-/transplatin were explored. The polarizable continuum model was used for estimation of solvent effects. Using the B3LYP/6-31+G(d) method, the structures were optimized and vibrational frequencies estimated. Interaction energies and activation barriers were determined at the CCSD(T)/6-31++G(d,p) level within the COSMO approach. An associative mechanism was assumed with a trigonal-bipyramidal structure of the transition state. Within the applied model, all the hydration reactions are slightly endothermic. The Gibbs energies of cisplatin hydration amount to 7.0 and 14.2 kcal/mol for the chloride and ammonium replacement, respectively. Analogous values for the transplatin reactions are 6.8 and 11.9 kcal/mol. The determined rate constants are by several (three to four) orders of magnitude larger for the dechlorination process than for deammination. The cisplatin dechlorination rate constant was established as 1.3 x 10(-4) s(-1) in excellent accord with the experiment.