Protonation states of molecules significantly influence the thermodynamics and kinetics of chemical reactions. This is especially important in biochemical processes, where appropriate protonation states of amino acids control the exo/endoergicity of practically all biochemical cycles. This paper is focused on appraisal of the impact of DFT functionals and PCM solvation models on the accuracy of pKa evaluations for all proteinogenic amino acids. Eight functionals (B3LYP, PBE0, revPBE0, M06-2X, M11, M11-L, TPSSh, and ωB97X-D) and four basis sets are considered, together with four kinds of implicit solvation models when additional attention is paid to a cavity construction. An influence of nonelectrostatic contributions and Wertz's corrections on Gibbs free energy is investigated together with accuracy of provided proton solvation energy. The best model is based on the M06-2X/6-311++G**/D-PCM/UAKS computational level. The fitting procedure is utilized to improve the accuracy of the evaluated models. All of these results are also compared with values obtained from the COSMOtherm program and CCSD(T) calculations. Results for cysteine and histidine are discussed individually, as they can be found in different protonation states at neutral pH.

• Publikační typ
• časopisecké články MeSH

The kinetics of the hydration reaction on trans-[Pt(NH3)2(pyrX)Cl]+ (pyr = pyridine) complexes (X = OH-, Cl-, F-, Br-, NO2 -, NH2, SH-, CH3, C≡CH, and DMA) was studied by density functional theory calculations in the gas phase and in water solution described by the implicit polarizable continuum model method. All possible positions ortho, meta, and para of the substituent X in the pyridine ring were considered. The substitution of the pyr ligand by electron-donating X's led to the strengthening of the Pt-N1(pyrX) (Pt-NpyrX) bond and the weakening of the trans Pt-Cl or Pt-Ow bonds. The electron-withdrawing X's have exactly the opposite effect. The strengths of these bonds can be predicted from the basicity of sigma electrons on the NpyrX atom determined on the isolated pyrX ligand. As the pyrX ring was oriented perpendicularly with respect to the plane of the complex, the nature of the X···Cl electrostatic interaction was the decisive factor for the transition-state (TS) stabilization which resulted in the highest selectivity of ortho-substituted systems with respect to the reaction rate. Because of a smaller size of X's, the steric effects influenced less importantly the values of activation Gibbs energies ΔG ⧧ but caused geometry changes such as the elongation of the Pt-NpyrX bonds. Substitution in the meta position led to the highest ΔG ⧧ values for most of the X's. The changes of ΔG ⧧ because of electronic effects were the same in the gas phase and the water solvent. However, as the water solvent dampened electrostatic interactions, 2200 and 150 times differences in the reaction rate were observed between the most and the least reactive mono-substituted complexes in the gas phase and the water solvent, respectively. An additional NO2 substitution of the pyrNO2 ligand further decelerated the rate of the hydration reaction, but on the other hand, the poly-NH2 complexes were no more reactive than the fastest o-NH2 system. In the gas phase, the poly-X complexes showed the additivity of the substituent effects with respect to the Pt-ligand bond strengths and the ligand charges.

• Publikační typ
• časopisecké články MeSH

Interaction of cisplatin in activated diaqua-form with His-Met dipeptide is explored using DFT approach with PCM model. First the conformation space of the dipeptide is explored to find the most stable structure (labeled 0683). Several functionals with double-zeta basis set are used for optimization and obtained order of conformers is confirmed by the CCSD(T) single-point calculations. Supermolecular model is used to determine reaction coordinate for the replacement of aqua ligands consequently by N-site of histidine and S-site of methionine and reversely. Despite the monoadduct of Pt-S(Met) is thermodynamically less stable this reaction passes substantially faster (by several orders of magnitude) than coordination of cisplatin to histidine. The consequent chelate formation occurs relatively fast with energy release up to 12 kcal mol-1.

• Klíčová slova
• Anticancer drug, Computational chemistry, Density functional theory, Heavy metal, Thermodynamics,
• MeSH
• chelátory chemie   MeSH
• cisplatina chemie   MeSH
• dipeptidy chemie   MeSH
• histidin chemie   MeSH
• kinetika MeSH
• methionin chemie   MeSH
• protinádorové látky chemie   MeSH
• teorie funkcionálu hustoty * MeSH
• termodynamika MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• chelátory MeSH
• cisplatina MeSH
• dipeptidy MeSH
• histidin MeSH
• methionin MeSH
• protinádorové látky MeSH