Reaction of hypochlorous acid with imidazole: formation of 2-chloro- and 2-oxoimidazoles

Reaction of hypochlorous acid (HOCl) with imidazole (Im) taken as a model for the 5-membered ring of guanine, leading to the products 2-chloro- and 2-oxo-imidazoles was investigated at the B3LYP/6-31+G* and B3LYP/AUG-cc-pVDZ levels of density functional theory. For all cases, single point energy calculations were performed at the MP2/AUG-cc-pVDZ level of theory using the geometries optimised at the B3LYP/AUG-cc-pVDZ level. Intrinsic reaction coordinate calculations were performed to ensure genuineness of all the calculated transition states. Effect of aqueous media was investigated by solvating all the species involved in the reactions using the polarizable continuum model. It is found that 2-chloroimidazole (2-ClIm) can be formed following three different reaction schemes while 2-oxoimidazole (2-oxoIm) can be formed following two different reaction schemes. The calculated barrier energies show that formation of 2-oxoIm would be less favored than that of 2-ClIm, which explains the experimental observations on relative yields of 8-chlorodeoxyguanosine and 8-oxodeoxyguanosine.     

Catalytic involvement of CO2 in the mutagenesis caused by reactions of ONOO(-) with guanine

The catalytic role of CO2 in reactions of ONOO- with guanine, leading to the formation of the mutagenic species 8-oxoguanine (8-oxoG) and 8-nitroguanine anion (8-nitroG-), was investigated by considering the reactions of nitrosoperoxycarbonate anion (ONOOCO2-), an adduct of ONOO- and CO2, with guanine at the B3LYP/6-31G** and B3LYP/AUG-cc-pVDZ levels of density functional theory in gas phase. In order to study bulk solvent effect, single-point energy calculations in aqueous media were carried out for all the species occurring in the reactions at the B3LYP/AUG-cc-pVDZ level of theory, by use of the polarizable continuum model (PCM). Vibrational frequency analysis was performed, and zero-point-energy (ZPE)-corrected total energies and Gibbs free energy changes at 298.15 K were obtained. The genuineness of the calculated transition states was confirmed by visually examining the vibrational modes and also by intrinsic reaction coordinate (IRC) calculations. The reaction between ONOOCO2- and guanine occurring through four different mechanisms leads to the formation of 8-oxoG or its anion, while the reaction between the same two species occurring through a different scheme leads to the formation of 8-nitroG-. It has been shown that the presence of a water molecule along with ONOOCO2- would not affect the reaction mechanisms significantly. Structures of the reactant complexes, product complexes and barrier energies involved in the reactions reveal that CO2 acts as a catalyst for the reaction between ONOO- and guanine. The cause of the catalytic action of CO2 is mainly due to intermediacy of the CO3 radical anion and NO2 radical into which ONOOCO2- is fragmented while reacting with guanine. The relative stabilities of the different product complexes suggest that the mutation caused by ONOO- in the presence of CO2 would mainly involve 8-oxoG.     

Formation of 8-nitroguanine and 8-oxoguanine due to reactions of peroxynitrite with guanine

Reactions of peroxynitrite with guanine were investigated using density functional theory (B3LYP) employing 6-31G** and AUG-cc-pVDZ basis sets. Single point energy calculations were performed at the MP2/AUG-cc-pVDZ level. Genuineness of the calculated transition states (TS) was tested by visually examining the vibrational modes corresponding to the imaginary vibrational frequencies and applying the criterion that the TS properly connected the reactant and product complexes (PC). Genuineness of all the calculated TS was further ensured by intrinsic reaction coordinate (IRC) calculations. Effects of aqueous media were investigated by solvating all the species involved in the reactions using the polarizable continuum model (PCM). The calculations reveal that the most stable nitro-product complex involving the anion of 8-nitroguanine and a water molecule i.e. 8NO(2)G(-) + H(2)O can be formed according to one reaction mechanism while there are two possible reaction mechanisms for the formation of the oxo-product complex involving 8-oxoguanine and anion of the NO(2) group i.e. 8OG + NO(2)(-). The calculated relative stabilities of the PC, barrier energies of the reactions and the corresponding enthalpy changes suggest that formation of the complex 8OG + NO(2)(-) would be somewhat preferred over that of the complex 8NO(2)G(-) + H(2)O. The possible biological implications of this result are discussed.     

Mechanisms of formation of 8-oxoguanine due to reactions of one and two OH* radicals and the H2O2 molecule with guanine: A quantum computational study

Mechanisms of formation of the mutagenic product 8-oxoguanine (8OG) due to reactions of guanine with two separate OH* radicals and with H2O2 were investigated at the B3LYP/6-31G, B3LYP/6-311++G, and B3LYP/AUG-cc-pVDZ levels of theory. Single point energy calculations were carried out with the MP2/AUG-cc-pVDZ method employing the optimized geometries at the B3LYP/AUG-cc-pVDZ level. Solvent effect was treated using the PCM and IEF-PCM models. Reactions of two separate OH* radicals and H2O2 with the C2 position of 5-methylimidazole (5MI) were investigated taking 5MI as a model to study reactions at the C8 position of guanine. The addition reaction of an OH* radical at the C8 position of guanine is found to be nearly barrierless while the corresponding adduct is quite stable. The reaction of a second OH* radical at the C8 position of guanine leading to the formation of 8OG complexed with a water molecule can take place according to two different mechanisms, involving two steps each. According to one mechanism, at the first step, 8-hydroxyguanine (8OHG) complexed with a water molecule is formed ,while at the second step, 8OHG is tautomerized to 8OG. In the other mechanism, at the first step, an intermediate complexed (IC) with a water molecule is formed, the five-membered ring of which is open, while at the second step, the five-membered ring is closed and a hydrogen bonded complex of 8OG with a water molecule is formed. The reaction of H2O2 with guanine leading to the formation of 8OG complexed with a water molecule can also take place in accordance with two different mechanisms having two steps each. At the first step of one mechanism, H2O2 is dissociated into two OH* groups that react with guanine to form the same IC as that formed in the reaction with two separate OH* radicals, and the subsequent step of this mechanism is also the same as that of the reaction of guanine with two separate OH* radicals. At the first step of the other mechanism of the reaction of guanine with H2O2, the latter molecule is dissociated into a hydrogen atom and an OOH* group which become bonded to the N7 and C8 atoms of guanine, respectively. At the second step of this mechanism, the OOH* group is dissociated into an oxygen atom and an OH* group, the former becomes bonded to the C8 atom of guanine while the latter abstracts the H8 atom bonded to C8, thus producing 8OG complexed with a water molecule. Solvent effects of the aqueous medium on certain reaction barriers and released energies are appreciable. 5MI works as a satisfactory model for a qualitative study of the reactions of two separate OH* radicals or H2O2 occurring at the C8 position of guanine.     

Reactions of NO(2)Cl with imidazole: a model study for the corresponding reactions of guanine

Reactions of nitryl chloride NO(2)Cl with imidazole, taken as a model for the guanine base of DNA, leading to the formation of 2-oxoimidazole (2-oxoIm), 2-chloroimidazole (2-chloroIm), and 2-nitroimidazole (2-nitroIm) corresponding to the 8-oxo, 8-chloro, and 8-nitro derivatives of guanine were studied at the B3LYP and MP2 levels of theory employing the 6-31+G* and AUG-cc-pVDZ basis sets in gas phase. In order to incorporate solvent effect, all the B3LYP/AUG-cc-pVDZ optimized structures were solvated in aqueous media at the B3LYP/AUG-cc-pVDZ and MP2/AUG-cc-pVDZ levels of theory using the polarizable continuum model (PCM). A single mechanism was found for the formation of 2-oxoIm, while two and three mechanisms were found for the formation of 2-nitroIm and 2-chloroIm respectively. Each of these reaction mechanisms involves two steps. The calculated barrier energies show that the formation of 2-nitroIm would occur more efficiently than those of 2-oxoIm and 2-chloroIm. It suggests that formation of 8-nitroguanine would be the main DNA lesion caused by NO(2)Cl, which is consistent with experimental observations.     

A quantum chemical study of reactions of DNA bases with sulphur mustard: a chemical warfare agent

Reactions of the sulphonium ion of sulphur mustard (SM⁺¹) at the N7, N3 and O6 sites of guanine, N7, N3 and N1 sites of adenine, O2 and N3 sites of cytosine and O2 and O4 sites of thymine were studied theoretically in gas phase and aqueous media employing density functional theory (DFT) and second order Møller–Plesset perturbation (MP2) theory. The B3LYP, B3PW91 and B1B95 functionals of DFT and the 6-31+G* and AUG-cc-pVDZ basis sets were used in the calculations. Basis set superposition error was treated using the counterpoise method by single point energy calculations at the B3LYP/6-31+G* level in gas phase. The present study explains the mechanism of alkylation of the DNA bases and shows that SM⁺¹ would form stable adducts at the endocyclic nitrogen sites of the DNA bases, and at the O6 site of guanine and the O2 site of cytosine. Formation of adducts at the N7 site of guanine and N3 site of adenine are found to be most favored and next most favored respectively, which agrees with experimental observations.     

Mechanistic study of the deamination reaction of guanine: a computational study

The mechanism for the deamination of guanine with H(2)O, OH(-), H(2)O/OH(-) and for GuaH(+) with H(2)O has been investigated using ab initio calculations. Optimized geometries of the reactants, transition states, intermediates, and products were determined at RHF/6-31G(d), MP2/6-31G(d), B3LYP/6-31G(d), and B3LYP/6-31+G(d) levels of theory. Energies were also determined at G3MP2, G3MP2B3, G4MP2, and CBS-QB3 levels of theory. Intrinsic reaction coordinate (IRC) calculations were performed to characterize the transition states on the potential energy surface. Thermodynamic properties (ΔE, ΔH, and ΔG), activation energies, enthalpies, and Gibbs free energies of activation were also calculated for each reaction investigated. All pathways yield an initial tetrahedral intermediate and an intermediate in the last step that dissociates to products via a 1,3-proton shift. At the G3MP2 level of theory, deamination with OH(-) was found to have an activation energy barrier of 155 kJ mol(-1) compared to 187 kJ mol(-1) for the reaction with H(2)O and 243 kJ mol(-1) for GuaH(+) with H(2)O. The lowest overall activation energy, 144 kJ mol(-1) at the G3MP2 level, was obtained for the deamination of guanine with H(2)O/OH(-). Due to a lack of experimental results for guanine deamination, a comparison is made with those of cytosine, whose deamination reaction parallels that of guanine.     

Potential energy surface of the reaction of imidazole with peroxynitrite: Density functional theory study

This article presents a theoretical investigation of the reaction mechanism of imidazole nitration by peroxynitrite using density functional theory calculations. Understanding this reaction mechanism will help in elucidating the mechanism of guanine nitration by peroxynitrite, which is one of the assumed chemical pathways for damaging DNA in cells. This work focuses on the analysis of the potential energy surface (PES) for this reaction in the gas phase. Calculations were carried out using Hartree–Fock (HF) and density functional theory (DFT) Hamiltonians with double‐zeta basis sets ranging from 6‐31G(d) to 6‐31++G(d,p), and the triple‐zeta basis set 6‐311G(d). The computational results reveal that the reaction of imidazole with peroxynitrite in gas phase produces the following species: (i) hydroxide ion and 2‐nitroimidazole, (ii) hydrogen superoxide ion and 2‐nitrosoimidazole, and (iii) water and 2‐nitroimidazolide. The rate‐determining step is the formation of a short‐lived intermediate in which the imidazole C₂ carbon is covalently bonded to peroxynitrite nitrogen. Three short‐lived intermediates were found in the reaction path. These intermediates are involved in a proton‐hopping transport from C₂ carbon to the terminal oxygen of the ? O? O moiety of peroxynitrite via the nitroso (ON? ) oxygen. Both HF and DFT calculations (using the Becke3–Lee–Yang–Parr functional) lead to similar reaction paths for proton transport, but the landscape details of the PES for HF and DFT calculations differ. This investigation shows that the reaction of imidazole with peroxynitrite produces essentially the same types of products (nitro‐ and nitroso‐) as observed experimentally in the reaction of guanine with peroxynitrite, which makes the former reaction a good model to study by computation the essential characteristics of the latter reaction. Nevertheless, the computationally determined activation energy for imidazole nitration by peroxynitrite in the gas phase is 84.1 kcal/mol (calculated at the B3LYP/6‐31++G(d,p) level), too large for an enzymatic reaction. Exploratory calculations on imidazole nitration in solution, and on the reaction of 9‐methylguanine with peroxynitrite in the gas phase and solution, show that solvation increases the activation energy for both imidazole and guanine, and that the modest decrease (15 kcal mol^?1) in the activation energy, due to the adjacent six member ring of guanine, is counterbalanced by solvation. These results lead to the speculation that proton tunneling may be at the origin of experimentally observed high reaction rate of guanine nitration by peroxynitrite in solution. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Reactions of guanine with methyl chloride and methyl bromide: O6‐methylation versus charge transfer complex formation

Density functional theory (DFT) at the B3LYP/6‐31+G* and B3LYP/AUG‐cc‐pVDZ levels was employed to study O6‐methylation of guanine due to its reactions with methyl chloride and methyl bromide and to obtain explanation as to why the methyl halides cause genotoxicity and possess mutagenic and carcinogenic properties. Geometries of the various isolated species involved in the reactions, reactant complexes (RCs), and product complexes (PCs) were optimized in gas phase. Transition states connecting the reactant complexes with the product complexes were also optimized in gas phase at the same levels of theory. The reactant complexes, product complexes, and transition states were solvated in aqueous media using the polarizable continuum model (PCM) of the self‐consistent reaction field theory. Zero‐point energy (ZPE) correction to total energy and the corresponding thermal energy correction to enthalpy were made in each case. The reactant complexes of the keto form of guanine with methyl chloride and methyl bromide in water are appreciably more stable than the corresponding complexes involving the enol form of guanine. The nature of binding in the product complexes was found to be of the charge transfer type (O6mG⁺ · X^?, X?Cl, Br). Binding of HCl, HBr, and H₂O molecules to the PCs obtained with the keto form of guanine did not alter the positions of the halide anions in the PCs, and the charge transfer character of the PCs was also not modified due to this binding. Further, the complexes obtained due to the binding of HCl, HBr, and H₂O molecules to the PCs had greater stability than the isolated PCs. The reaction barriers involved in the formation of PCs were found to be quite high (～50 kcal/mol). Mechanisms of genotoxicity, mutagenesis and carcinogenesis caused by the methyl halides appear to involve charge transfer‐type complex formation. Thus the mechanisms of these processes involving the methyl halides appear to be quite different from those that involve the other strongly carcinogenic methylating agents. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Activation barriers for DNA alkylation by carcinogenic methane diazonium ions

Methylation reactions of the DNA bases with the methane diazonium ion, which is the reactive intermediate formed from several carcinogenic methylating agents, were examined. The SN2 transition states of the methylation reactions at N7, N3, and O6 of guanine; N7, N3, and N1 of adenine; N3 and O2 of cytosine; and O2 and O4 of thymine were calculated using the B3LYP density functional method. Solvation effects were examined using the conductor-like polarizable continuum method and the combined discrete/SCRF method. The transition states for reactions at guanine N3, adenine N7, and adenine N1 are influenced by steric interactions between the methane diazonium ion and exocyclic amino groups. Both in the gas phase and in aqueous solution, the methylation reactions at N atoms have transition states that are looser, and generally occur earlier along the reaction pathways than reactions at O atoms. The forming bonds in the transition states in water are 0.03 to 0.13 A shorter than those observed in the gas phase, and the activation energies are 13 to 35 kcal/mol higher. The combined discrete/SCRF solvation energy calculations using base-water complexes with three water molecules yield base solvation energies that are larger than those obtained from the CPCM continuum method, especially for cytosine. Reactivities calculated using barriers obtained with the discrete/SCRF method are consistent with the experimentally observed high reactivity at N7 of guanine.     

Excited state proton transfer in guanine in the gas phase and in water solution: a theoretical study

Theoretical investigations were performed to study the phenomena of ground and electronic excited state proton transfer in the isolated and monohydrated forms of guanine. Ground and transition state geometries were optimized at both the B3LYP/6-311++G(d,p) and HF/6-311G(d,p) levels. The geometries of tautomers including those of transition states corresponding to the proton transfer from the keto to the enol form of guanine were also optimized in the lowest singlet pipi* excited state using the configuration interaction singles (CIS) method and the 6-311G(d,p) basis set. The time-dependent density function theory method augmented with the B3LYP functional (TD-B3LYP) and the 6-311++G(d,p) basis set was used to compute vertical transition energies using the B3LYP/6-311++G(d,p) geometries. The TD-B3LYP/6-311++G(d,p) calculations were also performed using the CIS/6-311G(d,p) geometries to predict the adiabatic transition energies of different tautomers and the excited state proton transfer barrier heights of guanine tautomerization. The effect of the bulk aqueous environment was considered using the polarizable continuum model (PCM). The harmonic vibrational frequency calculations were performed to ascertain the nature of potential energy surfaces. The excited state geometries including that of transition states were found to be largely nonplanar. The nonplanar fragment was mostly localized in the six-membered ring. Geometries of the hydrated transition states in the ground and lowest singlet pipi* excited states were found to be zwitterionic in which the water molecule is in the form of hydronium cation (H3O(+)) and guanine is in the anionic form, except for the N9H form in the excited state where water molecule is in the hydroxyl anionic form (OH(-)) and the guanine is in the cationic form. It was found that proton transfer is characterized by a high barrier height both in the gas phase and in the bulk water solution. The explicit inclusion of a water molecule in the proton transfer reaction path reduces the barrier height drastically. The excited state barrier height was generally found to be increased as compared to that in the ground state. On the basis of the current theoretical calculation it appears that the singlet electronic excitation of guanine may not facilitate the excited state proton transfer corresponding to the tautomerization of the keto to the enol form.     

Mechanisms for the deamination reaction of cytosine with H2O/OH(-) and 2H2O/OH(-): a computational study

Mechanisms for the deamination reaction of cytosine with H 2O/OH (-) and 2H 2O/OH (-) to produce uracil were investigated using ab initio calculations. Optimized geometries of reactants, transition states, intermediates, and products were determined at MP2 and B3LYP using the 6-31G(d) basis set and at B3LYP/6-31+G(d) levels of theory. Single point energies were also determined at MP2/G3MP2Large and G3MP2 levels of theory. Thermodynamic properties (Delta E, Delta H, and Delta G), activation energies, enthalpies, and free energies of activation were calculated for each reaction pathway investigated. Intrinsic reaction coordinate (IRC) analysis was performed to characterize the transition states on the potential energy surface. Seven pathways for the deamination reaction were found. All pathways produce an initial tetrahedral intermediate followed by several conformational changes. The final intermediate for all pathways dissociates to product via a 1-3 proton shift. The activation energy for the rate-determining step, the formation of the tetrahedral intermediate for pathway D, the only pathway that can lead to uracil, is 115.3 kJ mol (-1) at the G3MP2 level of theory, in excellent agreement with the experimental value (117 +/- 4 kJ mol (-1)).     

Interactions of 1-methylimidazole with UO2(CH3CO2)2 and UO2(NO3)2: structural, spectroscopic, and theoretical evidence for imidazole binding to the uranyl ion

The first definitive high-resolution single-crystal X-ray structure for the coordination of the 1-methylimidazole (Meimid) ligand to UO2(Ac)2 (Ac = CH3CO2) is reported. The crystal structure evidence is confirmed by IR, Raman, and UV-vis spectroscopic data. Direct participation of the nitrogen atom of the Meimid ligand in binding to the uranium center is confirmed. Structural analysis at the DFT (B3LYP) level of theory showed a conformational difference of the Meimid ligand in the free gas-phase complex versus the solid state due to small energetic differences and crystal packing effects. Energetic analysis at the MP2 level in the gas phase supported stronger Meimid binding over H2O binding to both UO2(Ac)2 and UO2(NO3)2. In addition, self-consistent reaction field COSMO calculations were used to assess the aqueous phase energetics of combination and displacement reactions involving H2O and Meimid ligands to UO2R2 (R = Ac, NO3). For both UO2(NO3)2 and UO2(Ac)2, the displacement of H2O by Meimid was predicted to be energetically favorable, consistent with experimental results that suggest Meimid may bind uranyl at physiological pH. Also, log(Knitrate/KAc) calculations supported experimental evidence that the binding stoichiometry of the Meimid ligand is dependent upon the nature of the reactant uranyl complex. These results clearly demonstrate that imidazole binds to uranyl and suggest that binding of histidine residues to uranyl could occur under normal biological conditions.     

Computational study of the deamination reaction of cytosine with H2O and OH-

The mechanism for the deamination reaction of cytosine with H(2)O and OH(-) to produce uracil was investigated using ab initio calculations. Optimized geometries of reactants, transition states, intermediates, and products were determined at RHF/6-31G(d), MP2/6-31G(d), and B3LYP/6-31G(d) levels and for anions at the B3LYP/6-31+G(d) level. Single-point energies were also determined at B3LYP/6-31+G(d), MP2/GTMP2Large, and G3MP2 levels of theory. Thermodynamic properties (DeltaE, DeltaH, and DeltaG), activation energies, enthalpies, and free energies of activation were calculated for each reaction pathway that was investigated. Intrinsic reaction coordinate analysis was performed to characterize the transition states on the potential energy surface. Two pathways for deamination with H(2)O were found, a five-step mechanism (pathway A) and a two-step mechanism (pathway B). The activation energy for the rate-determining steps, the formation of the tetrahedral intermediate for pathway A and the formation of the uracil tautomer for pathway B, are 221.3 and 260.3 kJ/mol, respectively, at the G3MP2 level of theory. The deamination reaction by either pathway is therefore unlikely because of the high barriers that are involved. Two pathways for deamination with OH(-) were also found, and both of them are five-step mechanisms. Pathways C and D produce an initial tetrahedral intermediate by adding H(2)O to deprotonated cytosine which then undergoes three conformational changes. The final intermediate dissociates to product via a 1-3 proton shift. Deamination with OH(-), through pathway C, resulted in the lowest activation energy, 148.0 kJ/mol, at the G3MP2 level of theory.     

Binding to DNA purine base and structure-activity relationship of a series of structurally related Ru(II) antitumor complexes: a theoretical study

The thermodynamics of the binding of a series of structurally related Ru(II) antitumor complexes, that is, alpha-[Ru(azpy)2Cl2] 1, beta-[Ru(azpy)2Cl2] 2, alpha-[Ru(azpy)(bpy)Cl2] 3, and cis-[Ru(bpy)2Cl2] 4 to DNA purine bases (gunine, adenine at N7 site) has been studied by using the DFT method. The binding of imine form of 9-methyladenine (9-MeAde) to the Ru(II) moiety in a didentate fashion via its N6 and N7 atoms was also considered. The geometrical structures of the DNA model base adducts were obtained at the B3LYP/(LanL2DZ + 6-31G(d)) level in vacuo. The following exact single-point energy calculations were performed at the B3LYP/(LanL2DZ(f)+6-311+G(2d, 2p)) level both in vacuo and in aqueous solution using the COSMO model. The bond dissociation enthalpies and free energies, reaction enthalpies and free energies both in the gas phase and in aqueous solution for all considered Ru(II)-DNA model base adducts were obtained from the computations. The calculated bond dissociation enthalpies and free energies allow us to build a binding affinity order for the considered Ru(II)-DNA model base adducts. The theoretical results show that the guanine N7 is a preferred site for this series of complexes and support such an experimental fact that alpha-[Ru(azpy)(bpy)(9-EtGua)H2O](2+) (3-(9-EtGua)) is isomerized to alpha'-[Ru(azpy)(bpy)(9-EtGua)H2O](2+) (3'-(9-EtGua)). On the basis of structural and thermodynamical characteristics, the possible structure-activity relationship was obtained, and the distinct difference in cytotoxicities of this series of structurally related antitumor complexes was explained theoretically.     

Theoretical study of the oxidation of histidine by singlet oxygen

Herein we present a theoretical study of the reaction of singlet oxygen with histidine performed both in the gas phase and in aqueous solution. The potential energy surface of the reactive system was explored at the B3LYP/cc-pVTZ level of theory and the electronic energies were refined by means of single-point CCSD(T)/cc-pVTZ(-f) calculations. Solvent effects were taken into account by using a solvent continuum model (COSMO) and by adding explicit water molecules. The results show that the first step in the reaction mechanism corresponds to a nearly symmetric Diels-Alder addition of the singlet oxygen molecule to the imidazole ring to yield an endoperoxide, in agreement with experimental evidence. The intermediate formed can evolve along two different reaction paths leading to two isomeric hydroperoxides and, eventually, to open-chain or internally cyclised oxidised products. Water plays a significant role in stabilising the reaction structures by solvation and by acting as a bifunctional catalyst in the elimination/addition reaction steps. Our results explain why substituents at the N1-imidazole ring can hamper the evolution of the initial endoperoxide and result in Gibbs energy barriers in solution similar to those experimentally measured and suggest a likely route to the formation of peptide aggregates during the oxidation of histidine by singlet molecular oxygen.     

6-硫代黄嘌呤互变异构体的密度泛函理论计算

在密度泛函B3LYP/6-311G**水平下,对14种气相和水相中可能存在的6-硫代黄嘌呤异构体进行了几何构型的全自由度优化,并计算出它们的总能量、焓、熵、吉布斯自由能。Onsager反应场溶剂模型用于水相的计算.计算结果表明,6-硫代黄嘌呤在气相中和水相中主要以硫酮的形式存在.在气相和水相中,硫酮-N7(H)均比硫酮-N9(H)更稳定.计算结果同已有实验结果一致.6-硫代黄嘌呤异构化的熵效应小,对互变异构平衡几乎没有显著的影响,而焓变对互变异构产生了主要的影响.较详细地讨论了水溶剂化作用对异构体的能量、几何结构、电荷分布和偶极矩的影响.     

Mechanistic Study on the Monofunctional Binding of Hydrolysis Products of Non-planar Heterocyclic Complexes Trans-[PtCl_2NH_3（piperidine）] and Trans-[PtCl_2（piperidine）_2] to DNA Purine Bases

The monofunctional substitution reactions between trans-[PtCl(H2O)(NH3)(pip)]+,trans-[Pt(H2O)2(NH3)(pip)]2+,trans-[PtCl(H2O)(pip)2]+,trans-[Pt(H2O)2(pip)2]2+ (pip = piperidine) and adenine/guanine nucleotides are explored by using B3LYP hybrid functional and IEF-PCM salvation models. For the trans-[Pt(H2O)2(NH3)(pip)]2+ and trans-[PtCl(H2O)(NH3)(pip)]+ complexes,the computed barrier heights in aqueous solution are 13.5/13.5 and 11.6/11.6 kcal/mol from trans-Pt-chloroaqua complex to trans/cis-monoadduct for adenine and guanine,and the corresponding values are 20.7/20.7 and 18.8/18.8 kcal/mol from trans-Pt-diaqua complex to trans/cis-monoadduct for adenine and guanine,respectively. For trans-[PtCl(H2O)(pip)2]+ and trans-[Pt(H2O)2(pip)2]2+,the corresponding values are 21.5/21.3 and 19.4/19.4 kcal/mol,and 26.0/26.0 and 20.7/20.8 kal/mol for adenine and guanine,respectively. Our calculations demonstrate that the barrier heights of chloroaqua are lower than the corresponding values of diaqua for adenine and guanine. In addition,the free energies of activation for guanine in aqueous solution are all smaller than that for adenine,which predicts a preference of 1.9 kcal/mol when trans-[PtCl(H2O)(NH3)(pip)]+ and trans-[Pt(H2O)2(NH3)(pip)]2+ are the active agents and ～1.9 and ～ 5.3 kcal/mol when trans-[PtCl(H2O)(pip)2]+ and trans-[Pt(H2O)2(pip)2]2+ are the active agents,respectively. For the reaction of trans-Pt-chloroaqua (or diaqua) to cis-monoadduct,we obtain the same transition-state structure as from the reaction of trans-Pt-chloroaqua (or diaqua) to trans-monoadduct,which seems that the trans-Pt-chloroaqua (or diaqua) complex can generate trans-or cis-monoadduct via the same transition-state.