Ultrafast deactivation of an excited cytosine-guanine base pair in DNA

Multiconfigurational ab initio calculations and QM/MM molecular dynamics simulations of a photoexcited cytosine-guanine base pair in both gas phase and embedded in the DNA provide detailed structural and dynamical insights into the ultrafast radiationless deactivation mechanism. Photon absorption promotes transfer of a proton from the guanine to the cytosine. This proton transfer is followed by an efficient radiationless decay of the excited state via an extended conical intersection seam. The optimization of the conical intersection revealed that it has an unusual topology, in that there is only one degeneracy-lifting coordinate. This is the central mechanistic feature for the decay both in vacuo and in the DNA. Radiationless decay occurs along an extended hyperline nearly parallel to the proton-transfer coordinate, indicating the proton transfer itself is not directly responsible for the deactivation. The seam is displaced from the minimum energy proton-transfer path along a skeletal deformation of the bases. Decay can thus occur anywhere along the single proton-transfer coordinate, accounting for the remarkably short excited-state lifetime of the Watson-Crick base pair. In vacuo, decay occurs after a complete proton transfer, whereas in DNA, decay can also occur much earlier. The origin of this effect lies in the temporal electrostatic stabilization of dipole in the charge-transfer state in DNA.     

Pt(Ⅱ), Pd(Ⅱ) and Ni(Ⅱ) Complexes Binding to the N(7) Position of Guanine: Influence on the Guanine and Watson-crick GC Base Pair

Comprehensive ab initio calculations were performed on the coordination of Pt(Ⅱ),Pd(Ⅱ) and Ni(Ⅱ) adducts to the N(7) of guanine and guanine-cytosine (GC) base pair at the DFT level. The fully optimized geometries of the metal complexes were obtained and the stabilization energies of the interaction between metal adducts and nucleobase were calculated with B3LYP method by using 6-31^* basis set for the light atom. While the effective core potential (ECP) is used for metal cation. The results show that both cispalladium and cisnickel cause similar geometric changes of the base pair as cisplatin. For the coordination of metal adducts to guanine, platinum adduct possesses the highest stabilization energy; but the interaction between metal-guanine and cytosine for nickel is larger than that for platinum and palladium. It is worthy to note that hydrolysis effect can also cause significant changes in H-bonds.     

The models of proton assisted and the unassisted formation of CGC base triplets

The triple helix is formed by combining a double and a single strand DNAs in low pH and dissociates in high pH. Under such conditions, protonation of cytosine in the single strand is necessary for triplex formation where cytosine-guanine-cytosine (CGC+) base triplet stabilizes the triple helix. The mechanism of CGC+ triplet formation from guanine-cytosine (GC) and a protonated cytosine (C+) shows the importance of N3 proton. Similarly in the case of CGC (unprotonated) triplet, the donor acceptor H-bond at N3 hydrogen of the cytosine analog (C) initiates the interaction with GC. The correspondence between the two models of triplets, CGC+ and CGC, unambiguously assigned that protonation at N3 cytosine in low pH to be the first step in triplet formation, but a donor acceptor triplet (CGC) can be designed without involving a proton in the Hoogsteen H-bond. Further, the bases of cytosine analogue also show the capability of forming Watson Crick (WC) H-bonds with guanine.     

The Protonated Guanine–Cytosine Base Pair

Protonated base pairs were recently implicated in the context of DNA proton transfer and charge migration. The effects of protonating different sites of the guanine–cytosine (GC) base pair are studied here by using the DZP++ B3LYP density functional method. Optimized structures for the protonated GC base pair are compared with those of parent GC and the neutral hydrogenated GC radical (GCH). Proton and hydrogen‐atom additions significantly disturb the structure of the GC base pair. However, the structural perturbations arising from protonation are often less than those arising from hydrogenation of GC. Protonation of the GC base pair causes significant strengthening of the interstrand hydrogen bonds and a concomitant increase in the base dissociation energies. The adiabatic ionization potentials (AIPs), vertical ionization potentials (VIPs), and proton affinities (PAs) for the different protonation sites of the GC base pair are predicted. The N7 site of guanine is the preferred site for protonation of the GC base pair.     

Base-pair interactions in the gas-phase proton-bonded complexes of C+G and C+GC

Interactions involved in the formation of gas-phase proton-bonded molecular complexes of cytosine (C) and guanine (G) were theoretically investigated for the case of C(+)G and C(+)GC using B3LYP density functional theory. In this study, particular focus was on the dimeric interaction of proton-bonded C(+)G, where a proton bond and a hydrogen bond are cooperatively involved. The dimer interaction energy in terms of dissociation energy (D(e)) was predicted to be 41.8 kcal/mol. The lowest (frozen) energy structure for the C(+)G dimeric complex was found to be CH(+)...G rather than C...H(+)G in spite of the lower proton affinity of the cytosine moiety, which was more stable by 3.3 kcal/mol. The predicted harmonic vibrational frequencies and bond lengths suggest that the combined contributions of proton and hydrogen bonding may determine the resultant stability of each complex structure. In contrast to the dimer case, in the case of the isolated C(+)GC triplet, the two minimum energy structures of CH(+)...GC and C...H(+)GC were predicted to be almost equivalent in total energy. The dissociation energy (D(e)) for the C(+)G pairing in the C(+)GC triplet was 43.7 kcal/mol. Other energetics are also reported. As for the proton-transfer reaction in the proton-bond axis, the forward proton-transfer barriers for the dimer and trimer complexes were also predicted to be very low, 3.6 and 1.5 kcal/mol (DeltaE(e)(PT)), respectively.     

Effects of OH radical addition on proton transfer in the guanine-cytosine base pair

Double proton transfer (PT) reactions in guanine-cytosine OH radical adducts are studied by the hybrid density functional B3LYP approach. Concerted and stepwise proton-transfer processes are explored between N1(H) on guanine (G) and N3 on cytosine (C), and between N4(H) on C and O6 on G. All systems except GC6OH display a concerted mechanism. 8OHGC has the highest dissociation energy and is 1.2 kcal/mol more stable than the nonradical GC base pair. The origin of the interactions are investigated through the estimation of intrinsic acid-basic properties of the *OH-X monomer (X = G or C). Solvent effects play a significant role in reducing the dissociation energy. The reactions including *OH-C adducts have significantly lower PT barriers than both the nonradical GC pair and the *OH-G adducts. All reactions are endothermic, with the GC6OH --> GC6OHPT reaction has the lowest reaction energy (4.6 kcal/mol). In accordance with earlier results, the estimated NBO charges show that the G moiety carries a slight negative charge (and C a corresponding positive one) in each adduct. The formation of a partial ion pair may be a potential factor leading to the PT reactions being thermodynamically unfavored.     

Interactions of Solvated Electrons with Nucleobases: The Effect of Base Pairing

We have investigated the effect of base pairing on the electron attachment to nucleobases in bulk water, taking the guanine−cytosine (GC) base pair as a test case. The presence of the complementary base reinforces the stabilization effect provided by water and preferentially stabilizes the anion by hydrogen bonding. The electron attachment in bulk-solvated GC happens through a doorway mechanism, where the initial electron attached state is water bound, and it subsequently gets converted to a GC bound state. The additional electron in the final GC bound state is localized on the cytosine, similar to that in the gas phase. The transfer of the electron from the initial water-bound state to the final GC bound state happens due to the mixing of electronic and nuclear degrees of freedom and takes place at a picosecond time scale.     

Cu2+/+ cation coordination to adenine--thymine base pair. Effects on intermolecular proton-transfer processes

Intermolecular proton-transfer processes in the Watson & Crick adenine-thymine Cu+ and Cu2+ cationized base pairs have been studied using the density functional theory (DFT) methods. Cationized systems subject to study are those resulting from cation coordination to the main basic sites of the base pair, N7 and N3 of adenine and O2 of thymine. For Cu+ coordinated to N7 or N3 of adenine, only the double proton-transferred product is found to be stable, similarly to the neutral system. However, when Cu+ interacts with thymine, through the O2 carbonyl atom, the single proton transfer from thymine to adenine becomes thermodynamically spontaneous, and thus rare forms of the DNA bases may spontaneously appear. For Cu2+ cation, important effects on proton-transfer processes appear due to oxidation of the base pair, which stabilizes the different single proton-transfer products. Results for hydrated systems show that the presence of the water molecules interacting with the metal cation (and their mode of coordination) can strongly influence the ability of Cu2+ to induce oxidation on the base pair.     

B型DNA中鸟嘌呤-胞嘧啶碱基对内双质子转移反应

采用ONIOM(M06-2X/6-31G^*:PM3)方法研究了单个鸟嘌呤-胞嘧啶(GC)碱基对和含GC碱基对的四种排序的DNA三聚体(dATGCAT, dGCGCGC, dTAGCTA, dCGGCCG)的双质子转移反应. 通过分析其双质子转移方式、质子转移过程中各结构的能量和氢键变化, 总结出环境因素对GC碱基对双质子转移机理的影响. 气相中, dCGGCCG三聚体中发生分步双质子转移, 其它四种模型中均发生协同双质子转移. 分析发现质子转移方式受上下相邻碱基对的静电相互作用和质子接受位的质子亲和势影响, dATGCAT和dGCGCGC排序有助于质子H4a转移, 而dTAGCTA和dCGGCCG排序有助于质子H1转移, 胞嘧啶的N3位较高的质子亲和势有助于质子H1转移. 水溶剂中, 上下相邻碱基对的静电相互作用被减弱, 水溶剂稳定了分步转移过程中的单质子转移产物, 因此分步转移机理占据优势, 五种模型中均出现分步双质子转移, 在此过程中能量变化趋势相似. 溶剂效应有利于单质子转移, 却增加了双质子转移反应的反应能.     

Effect of nucleobase sequence on the proton-transfer reaction and stability of the guanine-cytosine base pair radical anion

The formation of base pair radical anions is closely related to many fascinating research fields in biology and chemistry such as radiation damage to DNA and electron transport in DNA. However, the relevant knowledge so far mainly comes from studies on isolated base pair radical anions, and their behavior in the DNA environment is less understood. In this study, we focus on how the nucleobase sequence affects the properties of the guanine-cytosine (GC) base pair radical anion. The energetic barrier and reaction energy for the proton transfer along the N(1)(G)-H···N(3)(C) hydrogen bond and the stability of GC˙(-) (i.e., electron affinity of GC) embedded in different sequences of base-pair trimer were evaluated using density functional theory. The computational results demonstrated that the presence of neighboring base pairs has an important influence on the behavior of GC˙(-) in the gas phase. The excess electron was found to be localized on the embedded GC and the charge leakage to neighboring base pairs was very minor in all of the investigated sequences. Accordingly, the sequence behavior of the proton-transfer reaction and the stability of GC˙(-) is chiefly governed by electrostatic interactions with adjacent base pairs. However, the effect of base stacking, due to its electrostatic nature, is severely screened upon hydration, and thus, the sequence dependence of the properties of GC˙(-) in aqueous environment becomes relatively weak and less than that observed in the gas phase. The effect of geometry relaxation associated with neighboring base pairs as well as the possibility of proton transfer along the N(2)(G)-H···O(2)(C) channel have also been investigated. The implications of the present findings to the electron transport and radiation damage of DNA are discussed.     

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.     

Cu2+对腺嘧啶-胸腺嘧啶碱基对阴离子质子转移的影响

采用B3LYP/DZP++方法研究了腺嘌呤-胸腺嘧啶(A-T)碱基对阴离子(AT)-的单质子转移机理以及金属离子Cu2+对(AT)-碱基对质子转移的影响.(AT)-碱基对的单质子转移路径是由胸腺嘧啶N25位上的质子H26沿分子间的氢键N25-H26…N10转移到腺嘌呤的N10位.金属Cu2+可通过络合作用分别吸附在(AT)-碱基对O24、O28、N4、N13上,从而影响(AT)-碱基对中质子转移过程.Cu2+络合作用在胸腺嘧啶(T)的O24、O28上时,发生了从胸腺嘧啶到腺嘌呤方向上的单质子转移反应;而作用在腺嘌呤(A)的N4、N13上时,得到了双质子转移的稳定产物.     

Effects of protonation on proton-transfer processes in guanine–cytosine Watson–Crick base pairs

Intermolecular proton-transfer processes in guanine–cytosine Watson–Crick base pairs have been studied using the B3LYP density functional method. Protonation of the base pair was carried out both at the N7 and at the O6 atoms of guanine. It is found that protonation induces a strengthening of the base pair and facilitates the N1–N3 single-proton-transfer reaction. The double-proton-transfer reaction, however, turns out to be unfeasible when the system is protonated at these sites. Mutagenic implications of these proton-transfer processes are discussed.Proceedings of the 11th International Congress of Quantum Chemistry satellite meeting in honor of Jean-Louis Rivail     

Role of the guanine N1 imino proton in the migration and reaction of radical cations in DNA oligomers

Oxidation of a guanine nucleobase to its radical cation in DNA oligomers causes an increase in the acidity of the N1 imino proton that may lead to its spontaneous transfer to N3 of the paired cytosine. This proton transfer is suspected of playing an important role in long-distance radical cation hopping in DNA and the decisive product-determining role in the reaction of the radical cation with H2O or O2. We prepared and investigated DNA oligomers in which certain deoxycytidines are replaced by 5-fluoro-2'-deoxycytidines (F5dC). The pKa of F5C was determined to be 1.7 units below that of dC, which causes proton transfer from the guanine radical cation to be thermodynamically unfavorable. Photoinitiated one-electron oxidation of the DNA by UV irradiation of a covalently attached anthraquinone derivative introduces a radical cation that hops throughout the oligomer and is trapped selectively at GG steps. The introduction of F5dC does not affect the efficiency of charge hopping, but it significantly reduces the amount of reaction at the GG sites, as revealed by subsequent reaction with formamidopyrimidine glycosylase. These findings suggest that transfer of the guanine radical cation N1 proton to cytosine does not play a significant role in long-range charge transfer, but this process does influence the reactions with H2O and/or O2.     

Efficient and accurate calculations on the electronic structure of B-type poly(dG).poly(dC) DNA by elongation method: first step toward the understanding of the biological properties of aperiodic DNA

Elongation method was applied to determine the electronic structures of B-type poly(dG).poly(dC) DNA at the ab initio molecular orbital level as a first step toward the calculation of aperiodic DNA. The discrepancy in total energy between the elongation method and a conventional calculation was negligibly small in the order of 10(-8) hartreeat. for 14 G-C base pair model. The local density of states for 10 G-C base pair model estimated by the elongation method well reproduced the results by the conventional calculation. It was found that the band gap of the whole system is mainly due to the energy difference between the valence band of guanine and the conduction band of cytosine. Moreover, the electron transfer path through stacking G-C base pairs rather than sugar-phosphate backbones has been confirmed by the authors' calculations.     

Electron attachment to the hydrogenated Watson-Crick guanine cytosine base pair (GC+H): conventional and proton-transferred structures

The anionic species resulting from hydride addition to the Watson-Crick guanine-cytosine (GC) DNA base pair are investigated theoretically. Proton-transferred structures of GC hydride, in which proton H1 of guanine or proton H4 of cytosine migrates to the complementary base-pair side, have been studied also. All optimized geometrical structures are confirmed to be minima via vibrational frequency analyses. The lowest energy structure places the additional hydride on the C6 position of cytosine coupled with proton transfer, resulting in the closed-shell anion designated 1T (G(-)C(C6)). Energetically, the major groove side of the GC pair has a greater propensity toward hydride/hydrogen addition than does the minor grove side. The pairing (dissociation) energy and electron-attracting ability of each anionic structure are predicted and compared with those of the neutral GC and the hydrogenated GC base pairs. Anion 8T (G(O6)C(-)) is a water-extracting complex and has the largest dissociation energy. Anion 2 (GC(C4)(-)) and the corresponding open-shell radical GC(C4) have the largest vertical electron detachment energy and adiabatic electron affinity, respectively. From the difference between the dissociation energy and electron-removal ability of the normal GC anion and the most favorable structure of GC hydride, it is clear that one may dissociate the GC anion and maintain the integrity of the GC hydride.     

Probing irradiation induced DNA damage mechanisms using excited state Car-Parrinello molecular dynamics

Photoinduced proton transfer in the Watson-Crick guanine (G)-cytosine (C) base pair has been studied using Car-Parrinello molecular dynamics (CP-MD). A flexible mechanical constraint acting on all three hydrogen bonds in an unbiased fashion has been devised to explore the free energy profile along the proton transfer coordinate. The lowest barrier has been found for proton transfer from G to C along the central hydrogen bond. The resulting charge transfer excited state lies energetically close to the electronic ground state suggesting the possibility of efficient radiationless decay. It is found that dynamic, finite temperature fluctuations significantly reduce the energy gap between the ground and excited states for this charge transfer product, promoting the internal conversion process. A detailed analysis of the internal degrees of freedom reveals that the energy gap is considerably reduced by out-of-plane molecular vibrations, in particular. Consequently, it appears that considering only the minimum energy path provides an upper-bound estimate of the associated energy gap compared to the full-dimension dynamical reaction coordinate. Furthermore, the first CP-MD simulations of the G-C base pair in liquid water are presented, and the effects of solvation on its electronic structure are analyzed.     

Regulation of one-electron oxidation rate of guanine by base pairing with cytosine derivatives

Effects of base pairing on the one-electron oxidation rate of guanine derivatives, guanine, 8-bromoguanine, and 8-oxo-7,8-dihydroguanine have been studied. The one-electron oxidation rate of guanine derivatives was determined by triplet-quenching experiments, using N,N'-dibutylnaphthaldiimide (NDI) in the triplet excited state (3NDI*) and fullerene (C(60)) in the triplet excited state ((3)C(60*)) as oxidants. In all three guanine derivatives studied here, acceleration of the one-electron oxidation was observed upon hydrogen bonding with cytosine, which demonstrates lowering of the oxidation potential of guanine derivatives by base pairing with cytosine. When a methyl or bromo group was introduced to the C5 position of cytosine, acceleration or suppression of the one-electron oxidation relative to the guanine:cytosine base pair was observed, respectively. The results demonstrate that the one-electron oxidation rate of guanine in DNA can be regulated by introducing a substituent on base pairing cytosine.     

Low‐Energy‐Barrier Proton Transfer Induced by Electron Attachment to the Guanine⋅⋅⋅Cytosine Base Pair

The photoelectron spectrum of the anion of the guanine ??? cytosine base apair (GC)^.? is recorded for the first time. The observed variation in the spectral peak‐height ratios with the source conditions suggests the presence of two or more anionic isomers. Two maxima of the broad bands in the photoelectron spectrum were measured at about 1.9 and about 2.6 eV. These values are very well reproduced by the vertical detachment energies of the B3LYP/6‐31++G(d,p) calculated low‐energy anionic structures, which are 1) the Watson–Crick base‐pair anion with proton transferred from N1 of guanine to N3 of cytosine, 2) its analogue in which the proton is transferred from N9 of guanine to N7 of guanine, and 3) the global minimum geometry, which is formed from the latter anion by rotation of guanine about the axis approximately defined by C2 of guanine and C4 of cytosine. Furthermore, a minor difference in the stabilities of the two lowest energy anions explains the experimentally observed source (temperature) dependence of the PES spectrum. A rational procedure, based on the chemistry involved in the formation of anionic dimers, which enables the low‐energy anions populated in the photoelectron spectrum to be identified is proposed. In contrast to the alternative combinatorial approach, which in the studied case would lead to carrying out quantum chemical calculations for 2000–2500 structures, the procedure described here reduces the computational problem to only 15 geometries.