Proton Transfer in Guanine−Cytosine Radical Anion Embedded in B-Form DNA

Tunneling Effect in Proton Transfer: Transfer Matrix Approach

The transfer matrix (TM) method was applied to calculate the transmission probability (TP) for proton transfer reactions. The tunneling factors in the reaction rate constants were also evaluated using the TPs. To test this method, TPs for the Eckart potentials modeled as a guanine-cytosine base pair were calculated by the TM method and compared to TPs by the analytical solution. As a result, the errors in the TPs by the TM method were quite small. The tunneling factors for the guanine-thymine (G-T) and adenine-cytosine (A-C) mispair reactions were then evaluated by the TM method. A shoulder appears on each potential for these reactions [Odai, K.; Umesaki,K. J. Phys. Chem. A. 2021, 125, 8196-8204]. As a result, the shoulder for the G-T mispair reaction contributes significantly to the tunneling, while the shoulder for the A-C mispair reaction contributes little to the tunneling. These results are difficult to obtain with methods such as Wigner's tunneling factor.

The DNA Radical Code. Resolution of Identity in Dissociations of Trinucleotide Codon Cation Radicals in the Gas Phase

Sixty DNA trinucleotide cation radicals covering a large part of the genetic code alphabet were generated by electron transfer in the gas phase, and their chemistry was studied by collision-induced dissociation tandem mass spectrometry and theoretical calculations. The major dissociations involved loss of nucleobase molecules and radicals, backbone cleavage, and cross-ring fragmentations that depended on the nature and position of the nucleobases. Mass identity in dissociations of symmetrical trinucleotide cation radicals of the (XXX+2H)^+• and (XYX+2H)^+• type was resolved by specific ¹⁵N labeling. The specific features of trinucleotide cation radical dissociations involved the dominant formation of d₂⁺ ions, hydrogen atom migrations accompanying the formation of (w₂+H)^+•, (w₂+2H)⁺, and (d₂+2H)⁺ sequence ions, and cross-ring cleavages in the 3'- and 5'-deoxyribose moieties that depended on the nucleobase type and its position in the ion. Born-Oppenheimer molecular dynamics (BOMD) and density functional theory calculations were used to obtain structures and energies of several cation-radical protomers and conformers for (AAA+2H)^+•, (CCC+2H)^+•, (GGG+2H)^+•, (ACA+2H)^+•, and (CAA+2H)^+• that were representative of the different types of backbone dissociations. The ion electronic structure, protonation and radical sites, and hydrogen bonding were used to propose reaction mechanisms for the dissociations.

Collision-induced dissociation of homodimeric and heterodimeric radical cations of 9-methylguanine and 9-methyl-8-oxoguanine: correlation between intra-base pair proton transfer originating from the N1-H at a Watson-Crick edge and non-statistical dissociation

It has been shown previously in protonated, deprotonated and ionized guanine-cytosine base pairs that intra-base pair proton transfer from the N1-H at the Watson-Crick edge of guanine to the complementary nucleobase prompts non-statistical dissociation of the base-pair system, and the dissociation of a proton-transferred base-pair structure is kinetically more favored than that of the starting, conventional base-pair structure. However, the fundamental chemistry underlying this anomalous and intriguing kinetics has not been completely revealed, which warrants the examination of more base-pair systems in different structural contexts in order to derive a generalized base-pair structure-kinetics correlation. The purpose of the present work is to expand the investigation to the non-canonical homodimeric and heterodimeric radical cations of 9-methylguanine (9MG) and 9-methyl-8-oxoguanine (9MOG), i.e., [9MG·9MG]˙⁺, [9MOG·9MG]˙⁺ and [9MOG·9MOG]˙⁺. Experimentally, collision-induced dissociation tandem mass spectrometry coupled with an electrospray ionization (ESI) source was used for the formation of base-pair radical cations, followed by detection of dissociation product ions and cross sections in the collisions with Xe gas under single ion-molecule collision conditions and as a function of the center-of-mass collision energy. Computationally, density functional theory and coupled cluster theory were used to calculate and identify probable base-pair structures and intra-base pair proton transfer and hydrogen transfer reactions, followed by kinetics modeling to explore the properties of dissociation transition states and kinetic factors. The significance of this work is twofold: it provides insight into base-pair opening kinetics in three biologically-important, non-canonical systems upon oxidative and ionization damage; and it links non-statistical dissociation to intra-base pair proton-transfer originating from the N1-H at the Watson-Crick edge of 8-oxoguanine, enhancing understanding towards the base-pair fragmentation assisted by proton transfer.

Unveiling the complex pattern of intermolecular interactions responsible for the stability of the DNA duplex

Herein, we provide new insights into the intermolecular interactions responsible for the intrinsic stability of the duplex structure of a large portion of human B-DNA by using advanced quantum mechanical methods. Our results indicate that (i) the effect of non-neighboring bases on the inter-strand interaction is negligibly small, (ii) London dispersion effects are essential for the stability of the duplex structure, (iii) the largest contribution to the stability of the duplex structure is the Watson-Crick base pairing - consistent with previous computational investigations, (iv) the effect of stacking between adjacent bases is relatively small but still essential for the duplex structure stability and (v) there are no cooperativity effects between intra-strand stacking and inter-strand base pairing interactions. These results are consistent with atomic force microscope measurements and provide the first theoretical validation of nearest neighbor approaches for predicting thermodynamic data of arbitrary DNA sequences.

Experimental and theoretical assessment of protonated Hoogsteen 9-methylguanine-1-methylcytosine base-pair dissociation: kinetics within a statistical reaction framework

We investigated the collision-induced dissociation (CID) reactions of a protonated Hoogsteen 9-methylguanine-1-methylcytosine base pair (HG-[9MG·1MC + H]+), which aims to address the mystery of the literature reported "anomaly" in product ion distributions and compare the kinetics of a Hoogsteen base pair with its Watson-Crick isomer WC-[9MG·1MC + H]+ (reported recently by Sun et al.; Phys. Chem. Chem. Phys., 2020, 22, 24986). Product ion cross sections and branching ratios were measured as a function of center-of-mass collision energy using guided-ion beam tandem mass spectrometry, from which base-pair dissociation energies were determined. Product structures and energetics were assessed using various theories, of which the composite DLPNO-CCSD(T)/aug-cc-pVTZ//ωB97XD/6-311++G(d,p) was adopted as the best-performing method for constructing a reaction potential energy surface. The statistical Rice-Ramsperger-Kassel-Marcus theory was found to provide a useful framework for rationalizing the dominating abundance of [1MC + H]+ over [9MG + H]+ in the fragment ions of HG-[9MG·1MC + H]+. The kinetics analysis proved the necessity for incorporating into kinetics modeling not only the static properties of reaction minima and transition states but more importantly, the kinetics of individual base-pair conformers that have formed in collisional activation. The analysis also pinpointed the origin of the statistical kinetics of HG-[9MG·1MC + H]+vs. the non-statistical behavior of WC-[9MG·1MC + H]+ in terms of their distinctively different intra-base-pair hydrogen-bonds and consequently the absence of proton transfer between the N1 position of 9MG and the N3' of 1MC in the Hoogsteen base pair. Finally, the Hoogsteen base pair was examined in the presence of a water ligand, i.e., HG-[9MG·1MC + H]+·H2O. Besides the same type of base-pair dissociation as detected in dry HG-[9MG·1MC + H]+, secondary methanol elimination was observed via the SN2 reaction of water with nucleobase methyl groups.

Tackling a Curious Case: Generation of Charge-Tagged Guanosine Radicals by Gas-Phase Electron Transfer and Their Characterization by UV-vis Photodissociation Action Spectroscopy and Theory

We report the generation of gas-phase riboguanosine radicals that were tagged at ribose with a fixed-charge 6-(trimethylammonium)hexane-1-aminocarbonyl group. The radical generation relied on electron transfer from fluoranthene anion to noncovalent dibenzocrown-ether dication complexes which formed nucleoside cation radicals upon one-electron reduction and crown-ether ligand loss. The cation radicals were characterized by collision-induced dissociation (CID), photodissociation (UVPD), and UV-vis action spectroscopy. Identification of charge-tagged guanosine radicals was challenging because of spontaneous dissociations by loss of a hydrogen atom and guanine that occurred upon storing the ions in the ion trap without further excitation. The loss of H proceeded from an exchangeable position on N-7 in guanine that was established by deuterium labeling and was the lowest energy dissociation of the guanosine radicals according to transition-state energy calculations. Rate constant measurements revealed an inverse isotope effect on the loss of either hydrogen or deuterium with rate constants k_H = 0.25-0.26 s^-1 and k_D = 0.39-0.54 s^-1. We used time-dependent density functional theory calculations, including thermal vibronic effects, to predict the absorption spectra of several protomeric radical isomers. The calculated spectra of low-energy N-7-H guanine-radical tautomers closely matched the action spectra. Transition-state-theory calculations of the rate constants for the loss of H-7 and guanine agreed with the experimental rate constants for a narrow range of ion effective temperatures. Our calculations suggest that the observed inverse isotope effect does not arise from the isotope-dependent differences in the transition-state energies. Instead, it may be caused by the dynamics of post-transition-state complexes preceding the product separation.

Ion transfer channel network formed by flower and rod shape crystals of hair hydrolysate in poly(vinyl alcohol) matrix and its application as anion exchange membrane in fuel cells

Ion transfer phenomena occurring in nature are known to be most efficient. Many efforts have been made to mimic such phenomena, especially in the area of energy transfer. Proteins consisting of various amino acids are known to be the fundamental materials behind these phenomena. In the current study, an effort was made to extract proteinaceous material from human hair bio-waste by a green chemical-free thermal hydrolysis process. A simple heat treatment of the human hairs in presence of water led to the formation of a water soluble material, which was called hair hydrolysate (HH), contains 70 wt% proteinaceous material. It was utilized for the fabrication of poly(vinyl alcohol) (PVA) matrix-based anion exchange membrane (AEM). Presence of 27 wt% charged amino acids and 19 wt% polarizable amino acids in the HH provided effective charge transfer sites. 7 wt% arginine present in the HH, having continuous delocalized net positive charge helped the membrane to be stable in highly alkaline conditions, which was confirmed by an indirect analysis of alkaline stability. Formation of rod and flower shaped crystal morphology by the HH in glutaraldehyde crosslinked PVA matrix, created a continuous channel network at higher loadings, which provided a simple path for ion transfer, achieving OH^- conductivity of 7.46 mS/cm at 70 °C. Swelling of the PVA matrix was minimized by annealing of the HH loaded sample, which resulted in reduction of ionic conductivity to 6.16 mS/cm (at 70 °C). At the same time, improvement in the properties like increase in thermal, mechanical and thermo-mechanical stability, reduction in water uptake, %swelling and methanol permeability was observed. The selectivity of the membrane was increased to almost a decimal place. Thus, the HH obtained from simple green thermal hydrolysis of human hair bio-waste is a cheap material, which is found to be suitable as ion conductive material for alkaline fuel cells.

Density functional theory studies on cytosine analogues for inducing double-proton transfer with guanine

To induce double-proton transfer (DPT) with guanine in a biological environment, 12 cytosine analogues (Ca) were formed by atomic substitution. The DPT reactions in the Watson-Crick cytosine-guanine model complex (Ca₀G) and 12 modified cytosine-guanine complexes (Ca_1-12G) were investigated using density functional theory methods at the M06-2X/def2svp level. The intramolecular proton transfers within the analogues are not facile due to high energy barriers. The hydrogen bond lengths of the Ca_1-12G complexes are shorter than those in the Ca₀G complex, which are conducive to DPT reactions. The DPT energy barriers of Ca_1-12G complexes are also lower than that of the Ca₀G complex, in particular, the barriers in the Ca₇G and Ca₁₁G complexes were reduced to -1.33 and -2.02 kcal/mol, respectively, indicating they are significantly more prone to DPT reactions. The DPT equilibrium constants of Ca_1-12G complexes range from 1.60 × 10⁰ to 1.28 × 10⁷, among which the equilibrium constants of Ca₇G and Ca₁₁G are over 1.0 × 10⁵, so their DPT reactions may be adequate. The results demonstrate that those cytosine analogues, especially Ca₇ and Ca₁₁, are capable of inducing DPT with guanine, and then the guanine tautomer will form mismatches with thymine during DNA replication, which may provide new strategies for gene therapy.

Is non-statistical dissociation a general feature of guanine–cytosine base-pair ions? Collision-induced dissociation of a protonated 9-methylguanine–1-methylcytosine Watson–Crick base pair, and comparison with its deprotonated and radical cation analogues

Non-statistical dissociation of a protonated guanine–cytosine Watson–Crick base pair.

Electron and hole interactions with P, Z, and P:Z and the formation of mutagenic products by proton transfer reactions

Z would act as an electron acceptor and P would capture a hole in the unnatural DNA. The latter process would produce mutagenic products via a proton transfer reaction.

The Role of Proton Transfer on Mutations

Hydrogen bonds play a critical role in nucleobase studies as they encode genes, map protein structures, provide stability to the base pairs, and are involved in spontaneous and induced mutations. Proton transfer mechanism is a critical phenomenon that is related to the acid-base characteristics of the nucleobases in Watson-Crick base pairs. The energetic and dynamical behavior of the proton can be depicted from these characteristics and their adjustment to the water molecules or the surrounding ions. Further, new pathways open up in which protonated nucleobases are generated by proton transfer from the ionized water molecules and elimination of a hydroxyl radical in this review, the analysis will be focused on understanding the mechanism of untargeted mutations in canonical, wobble, Hoogsteen pairs, and mutagenic tautomers through the non-covalent interactions. Further, rare tautomer formation through the single proton transfer (SPT) and the double proton transfer (DPT), quantum tunneling in nucleobases, radiation-induced bystander effects, role of water in proton transfer (PT) reactions, PT in anticancer drugs-DNA interaction, displacement and oriental polarization, possible models for mutations in DNA, genome instability, and role of proton transfer using kinetic parameters for RNA will be discussed.

Accurate Base Pair Energies of Artificially Expanded Genetic Information Systems (AEGIS): Clues for Their Mutagenic Characteristics

Recently, several artificial nucleobases, such as B, S, J, V, X, K, P, and Z, have been proposed to help in the expansion of the genetic information system and diagnosis of diseases. Among these bases, P and Z were identified to form stable DNA and to participate in the replication. However, the stabilities of P:Z and other artificial base pairs are not fully understood. The abilities of these unnatural nucleobases in mispairing with themselves and with natural bases are also not known. Here, the ωB97X-D dispersion-corrected density functional theoretical and complete basis set (CBS-QB3) methods are used to obtain accurate structural and energetic data related to base pair interactions involving these unnatural nucleobases. The roles of protonation and deprotonation of certain artificial bases in inducing mutations are also studied. It is found that each artificial purine has a complementary artificial pyrimidine, the base pair interactions between which are similar to those of the natural Watson-Crick base pairs. Hence, these base pairs will function naturally and would not impart mutagenicity. Among these base pairs, the J:V complex is found to be the most stable and promising artificial base pair. Remarkably, the noncomplementary artificial nucleobases are found to form stable mispairs, which may generate mutagenic products in DNA. Similarly, the misinsertions of natural bases opposite artificial bases are also found to be mutagenic. The mechanisms of these mutations are explained in detail. These results are in agreement with earlier biochemical studies. It is thus expected that this study would aid in the advancement of the synthetic biology to design more robust artificial nucleotides.

The prediction of intermolecular proton-transfer of guanine-cytosine base pair under the influence of fragments from decomposed MOFs

Metal-organic frameworks (MOFs) can be decomposed into various fragments, including negative/positive charges, Zn⁺ or Cu²⁺ when used as drug delivery materials. To evaluate the safety of MOFs, different mechanisms of intermolecular proton-transfer in guanine-cytosine (GC) base pair under the influence of such fragments were investigated by density functional theory methods. In a vacuum, calculation results show that an excess electron assists proton transfer in the anionic GC radical, and a hole assists proton transfer in the cationic GC radical with small energy barriers. The mechanism for Zn⁺-GC transfer is that the located hole assists proton transfer from G to C. All proton-transfers of Cu²⁺-GC become spontaneous with stable proton-transferred structures, and the driving force is the Cu²⁺ due to its electrostatic and oxidative effects. However, in a micro-water environment, the average energy barrier of all proton-transfer processes increases by 2.8 kcal mol^-1 because of the redistribution of charges. Water molecules play a very important role in buffering, and the influence of fragments on intermolecular proton-transfer processes of GC is reduced.

Deprotonated guanine·cytosine and 9-methylguanine·cytosine base pairs and their "non-statistical" kinetics: a combined guided-ion beam and computational study

We report a guided-ion beam mass spectrometric study on collision-induced dissociation (CID) of deprotonated guanine(G)·cytosine(C) base pairs and their 9-methylguanine (9MG) analogue with Xe, including measurements of product cross sections as a function of collision energy and determination of dissociation thresholds. DFT, RI-MP2 and DLPNO-CCSD(T) calculations and Rice-Ramsperger-Kassel-Marcus (RRKM) modeling were performed to elucidate structures and kinetics. The experiment and theoretical study have provided considerable insight into tautomerization, intra-base-pair proton transfer and dissociation of deprotonated G·C and 9MG·C. In contrast to the previously reported lowest-energy deprotonated base pair structure G·[C-H1]^- that consists of H-bonded neutral guanine and N1-deprotonated cytosine, we found that proton transfer from guanine N1 to cytosine N3 within G·[C-H1]^- (or 9MG·[C-H1]^-) leads to another slightly more stable conformer denoted as G·[C-H1]^-_PT1 (or 9MG·[C-H1]^-_PT1). The conventional (non-proton-transferred) and the proton-transferred conformers are close in energy and interconvert quickly, but they can be distinguished by dissociation products. The conventional structure dissociates into deprotonated cytosine and neutral guanine, while the other dissociates into deprotonated guanine and neutral cytosine. The two dissociation asymptotes have similar threshold energies, but surprisingly the CID product mass spectra of deprotonated G·C and 9MG·C are both overwhelmingly dominated by deprotonated G or 9MG, with their branching ratios greater than RRKM predictions by one to two orders of magnitude. The proton-transferred structures of deprotonated base pairs and the "unexpected" non-statistical kinetics provide new leads for understanding purine-pyrimidine interactions, forming rare nucleobase tautomers, and base pair opening.

Theoretical Investigations of the Chiral Transition of α-Amino Acid Confined in Various Sized Armchair Boron-Nitride Nanotubes

We computationally study the chiral transition process of the α-Ala molecule under confined different sizes of armchair SWBNNTs to explore the confinement effect. We find that the influence of a confinement environment (in armchair SWBNNTs) on the α-Ala molecule would lead to different reaction pathways. Meanwhile, the preferred reaction pathway is also different in various sizes of armchair SWBNNTs, and their energy barriers for the rate-limiting step decrease rapidly with the decreasing of the diameters of the nanotubes. It is obvious that significant decrease of the chiral transition energy barrier occurs compared with the isolated α-Ala molecule chirality conversion mechanism, by ∼15.6 kcal mol^-1, highlighting the improvement in the activity the enantiomers of α-Ala molecule. We concluded that the confinement environment has a significant impact at the nanoscale on the enantiomer transformation process of the chiral molecule.

Charge transfer dynamics in DNA revealed by time-resolved spectroscopy

In the past few decades, charge transfer in DNA has attracted considerable attention from researchers in a wide variety of fields, including bioscience, physical chemistry, and nanotechnology. Charge transfer in DNA has been investigated using various techniques. Among them, time-resolved spectroscopic methods have yielded valuable information on charge transfer dynamics in DNA, providing an important basis for numerical practical applications such as development of new therapy applications and nanomaterials. In DNA, holes and excess electrons act as positive and negative charge carriers, respectively. Although hole transfer dynamics have been investigated in detail, the dynamics of excess electron transfer have only become clearer relatively recently. In the present paper, we summarize studies on the dynamics of hole and excess electron transfer conducted by several groups including our own.

Influence of Coupling and Embedding Schemes on QM Size Convergence in QM/MM Approaches for the Example of a Proton Transfer in DNA

The influence of embedding and coupling schemes on the convergence of the QM size in the QM/MM approach is investigated for the transfer of a proton in a DNA base pair. We find that the embedding scheme (mechanical or electrostatic) has a much greater impact on the convergence behavior than the coupling scheme (additive QM/MM or subtractive ONIOM). To achieve size convergence, QM regions with up to 6000 atoms are necessary for pure QM or mechanical embedding. In contrast, electrostatic embedding converges faster: for the example of the transfer of a proton between DNA base pairs, we recommend including at least five base pairs and 5 Å of solvent (including counterions) into the QM region, i.e., a total of 1150 atoms.

Exposing the G-quadruplex to electric fields: the role played by telomeres in the propagation of DNA errors

We use quantum calculations to assess the impact of external electric fields on the stability of G-quadruplex, a key structure in telomere functionality.

Electron interaction with a DNA duplex: dCpdC:dGpdG

Electron attachment to double-stranded cytosine-rich DNA, dCpdC:dGpdG, has been studied by density functional theory. This system represents a minimal descriptive unit of a cytosine-rich double-stranded DNA helix. A significant electron affinity for the formation of a cytosine-centered radical anion is revealed to be about 2.2 eV. The excess electron may reside on the nucleobase at the 5' position (dC˙(-)pdC:dGpdG) or at the 3' position (dCpdC˙(-):dGpdG). The inter-strand proton transfer between the radical anion centered cytosine (N3) and the paired guanine (HN1) results in the formation of radical anion center separated complexes dC1H˙pdC:dG2-H(-)pdG and dCpdC2H˙:dGpdG1-H(-). These distonic radical anions are found to be approximately 1 to 4 kcal mol(-1) more stable than the normal radical anions. Intra-strand cytosine π→π transition energies are below the electron detachment energy. Inter-strand π→π transitions of the excess electron from C to G are predicted to be less than 2.79 eV. Electron transfer might also be possible through the inter-strand base-jumping mode. An analysis of absorption visible spectra reveals the absorption bands ranging from 500 nm to 700 nm for the cytosine-rich radical anions of the DNA duplex. Electron attachment to cytidine oligomers might add color to the DNA duplex.

How Does Guanine-Cytosine Base Pair Affect Excess-Electron Transfer in DNA?

Charge transfer and proton transfer in DNA have attracted wide attention due to their relevance in biological processes and so on. Especially, excess-electron transfer (EET) in DNA has strong relation to DNA repair. However, our understanding on EET in DNA still remains limited. Herein, by using a strongly electron-donating photosensitizer, trimer of 3,4-ethylenedioxythiophene (3E), and an electron acceptor, diphenylacetylene (DPA), two series of functionalized DNA oligomers were synthesized for investigation of EET dynamics in DNA. The transient absorption measurements during femtosecond laser flash photolysis showed that guanine:cytosine (G:C) base pair affects EET dynamics in DNA by two possible mechanisms: the excess-electron quenching by proton transfer with the complementary G after formation of C(•-) and the EET hindrance by inserting a G:C base pair as a potential barrier in consecutive thymines (T's). In the present paper, we provided useful information based on the direct kinetic measurements, which allowed us to discuss EET through oligonucleotides for the investigation of DNA damage/repair.

Mutagenic effects induced by the attack of NO2 radical to the guanine-cytosine base pair

We investigate the attack of the nitrogen dioxide radical (NO^•₂) to the guanine—cytosine (GC) base pair and the subsequent tautomeric reactions able to induce mutations, by means of density functional theory (DFT) calculations. The conducted simulations allow us to identify the most reactive sites of the GC base pair. Indeed, the computed relative energies demonstrate that the addition of the NO^•₂ radical to the C8 position of the guanine base forms to the most stable adduct. Although the initial adducts might evolve to non-canonical structures via inter-base hydrogen bonds rearrangements, the probability for the proton exchange to occur lies in the same range as that observed for undamaged DNA. As a result, tautomeric errors in NO₂-attacked DNA arises at the same rate as in canonical DNA, with no macroscopic impact on the overall stability of DNA. The potential mutagenic effects of the GC–NO^•₂ radical adducts likely involve side reactions, e.g., the GC deprotonation to the solvent, rather than proton exchange between guanine and cytosine basis.

DNA spontaneous mutation and its role in the evolution of GC-content: assessing the impact of the genetic sequence

We use theoretical tools to investigate the possible role played by a DNA sequence in the base pair tautomerization phenomena.

DFT reinvestigation of DNA strand breaks induced by electron attachment

The benchmark study of DFT methods on the activation energies of phosphodiester C3'-O and C5'-O bond ruptures and glycosidic C1'-N bond ruptures induced by electron attachment was performed. While conventional pure and hybrid functionals provide a relatively reasonable description for the C1'-N bond rupture, they significantly underestimate the energy barriers of the C-O bond ruptures. This is because the transition states of the later reactions, which are characterized by an electron distribution delocalized from the nucleobase to sugar-phosphate backbone, suffer from a severe self-interaction error in common DFT methods. CAM-B3LYP, M06-2X, and ωB97XD are the top three methods that emerged from the benchmark study; the mean absolute errors relative to the CCSD(T) values are 1.7, 1.9, and 2.2 kcal/mol, respectively. The C-O bond cleavages of 3'- and 5'-dXMP(•-), where X represents four nucleobases, were then recalculated at the M06-2X/6-31++G**//M06-2X/6-31+G* level, and it turned out that the C-O bond cleavages do not proceed as easily as previously predicted by the B3LYP calculations. Our calculations revealed that the C-O bonds of purine nucleotides are more susceptible than pyrimidine nucleotides to the electron attachment. The energies of electron attachment to nucleotides were calculated and discussed as well.

Environmental-confinement-induced stability enhancement of chiral molecules

We computationally study the transition process of a chiral difluorobenzo[c]phenanthrene (DFBcPh) molecule within non-polar fullerene C(260) to explore the confinement effect. We find blue-shifts in the infrared and Raman spectra of the molecule inside the fullerene relative to those of isolated systems. Six types of spectrum features of the molecule appear in the 0-60 cm(-1) band. Interestingly, the energy barrier of the chiral transformation of the molecule is elevated by 15.88 kcal mol(-1) upon the confinement by the fullerene, indicating improvement in the stability of the enantiomers. The protection by C(260) lowers the highest occupied molecular orbital energy level and lifts the lowest unoccupied molecular orbital energy level of the chiral molecule such that the chiral molecule is further chemically stabilized. We concluded that the confinement environment has an impact at the nanoscale on the enantiomer transformation process of the chiral molecule.

Assessing the importance of proton transfer reactions in DNA

Although engineered by millions of years of evolution, the cellular machinery is not flawless, and errors regularly appear during DNA replication. The subsequent alteration of the stored genetic message results in a mutation and might be the starting point of important health disorders. The question therefore is what causes DNA mutations? All living organisms are constantly exposed to a number of external agents such as free radicals and to radiation, which may lead to induced mutations. There are also mutations happening without invoking the action of any exogenous element, the so-called spontaneous mutations. The former can be partially controlled by avoiding exposure to high-risk environments, while the latter are more intriguing because their origin is unclear and difficult to determine. As noted by Watson and Crick when they first discovered the DNA structure, the correct replication of DNA rests on the assumption that the base pairs remain in their most stable, canonical form. However, protons along the interbase hydrogen-bond network are not static entities. They can in fact interchange their positions in DNA bases through proton transfer (PT) reactions before strands unwind, giving rise to noncanonical structures defined as rare tautomers. The importance of these rare tautomers was also cleverly anticipated by Watson and Crick and some years later claimed by Löwdin to be a source of spontaneous mutations. In Watson and Crick's words: "It would be of interest to know the precise difference in free energy between the various tautomeric forms under physiological conditions." Unfortunately, rare tautomeric forms are very difficult to detect, so no direct and accurate free energy measure has been discerned. In contrast, theoretical chemistry is making good progress toward the quantification of PT reactions in DNA and their biological consequences. This Account touches upon the theoretical studies devoted to appraising the importance of rare tautomers as promoters of spontaneous mutations. We focus in particular on the crucial role played by the biological environment on DNA stability. It has now been demonstrated that valuable macroscopic predictions require not only highly accurate theories but also refined chemical models. Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations performed on short but complete DNA sequence fragments emerge in this context as the most adequate tools. In addition, these methods can be used to quantify the effect of different external agents on the PT tautomeric equilibria and, eventually, to conveniently handle them. This is the case for the possible alteration of the naturally observed mutation rate by exposure to intense electric fields. Theoretical predictions envision in this respect promising applications of ultrashort electric pulses in medicine to selectively modify the mutated/canonical ratio in DNA.

Coupling between hydrogen atoms transfer and stacking interaction in adenine-thymine/guanine-cytosine complexes: a theoretical study

Four different complexes of two base pairs, an adenine-thymine and a guanine-cytosine one, have been studied in order to understand the modifications induced by the staking interaction between the two base pairs on the hydrogen atoms transfers between the bases in either base pair. The inclusion of these two kinds of interactions allows us to clarify if some properties, as the mechanism of hydrogen transfer, is exclusively a local effect of a base pair or can be modified by a more long-range interaction between the base pairs. The results on these four complexes are compared with those of the monomeric systems, the A-T and G-C base pair, and with those of the A-T and G-C dimers. The specificity of each complex and of each hydrogen bond has been analyzed.

Hydrogen-bonded proton transfer in the hydrated adenine–thymine anion

The proton transfer processes of microhydrated adenine–thymine anions are studied using density functional theory with the B3LYP method and DZP++ basis set. The microhydration effects on the geometrical structures, adsorption site, and the proton transfer reaction of the adenine–thymine anion are investigated. The site N10 atom of the adenine moiety has a larger proton affinity than the site O24 atom of thymine, which facilitates the proton H26 transfers from the N25 site of thymine to the N10 site of adenine. Therefore, the first single-proton transfer pathway (SPT1) is found for the all of the monohydrated adenine–thymine anions (AN4T)−·H2O, (AN13T)−·H2O, (ATO24)−·H2O, and (ATO28)−·H2O and tetrahydrated adenine–thymine anions (AT)−·4H2O. The proton H9 at the N7 site of adenine is also found to transfer to the O24 site of thymine for (AN4T)−·H2O and (AN13T)−·H2O in the gas phase. The double-proton transferred pathway is found when one water molecule interacts with the O28 atom of thymine. The reactant structures before the proton transfer are more stable than the product structures, and the structural changes mainly occur in thymine. The reaction energies of the microhydrated adenine–thymine anion in the gas phase and in the aqueous environment predict that the proton transfer process of the microhydrated adenine–thymine anion are more favorable in the gas phase than in aqueous solution.

Theoretical investigation of the coupling between hydrogen-atom transfer and stacking interaction in adenine-thymine dimers

Three different dimers of the adenine-thymine (A-T) base pair are studied to point out the changes of important properties (structure, atomic charge, energy and so on) induced by coupling between the movement of the atoms in the hydrogen bonds and the stacking interaction. The comparison of these results with those for the A-T monomer system explains the role of the stacking interaction in the hydrogen-atom transfer in this biologically important base pair. The results support the idea that this coupling depends on the exact dimer considered and is different for the N-N and N-O hydrogen bonds. In particular, the correlation between the hydrogen transfer and the stacking interaction is more relevant for the N-N bridge than for the N-O one. Also, the two different mechanisms of two-hydrogen transfer (step by step and concerted) can be modified by the stacking interaction between the base pairs.