Calculated pKa’s of the DNA Base Radical Ions

pH-dependent complex formation with TAR RNA and DNA: application towards logic gates

Nucleic acid-based logic gates have shown great potential in biotechnology, medicine as well as diagnostics. Herein, we have constructed pH-responsive logic devices by utilizing HIV-1 TAR hairpins in combination with a thiazole peptide that exhibits turn-on fluorescence upon interacting with TAR RNA or DNA. Based on this, INHIBIT-AND and YES-INHIBIT-AND logic gates were constructed in parallel. The pH alteration leads to conformational changes of the hairpin structure, enabling the construction of a multi-reset reusable logic system which could be developed for in vitro sensing of the HIV-1 viral RNA.

Shape Resonances in DNA: Nucleobase Release, Reduction, and Dideoxynucleoside Products Induced by 1.3 to 2.3 eV Electrons

Understanding the details of DNA damage caused by high-energy particles or photons is complicated by the multitude of reactive species, arising from the ionization and dissociation of H₂O, DNA, and protein. In this work, oligonucleotides (ODNs) are irradiated with a beam of low-energy electrons of 1.3 to 2.3 eV, which can only induce damage via the decay of shape resonances into various dissociative electron attachment channels. Using LC-MS/MS analysis, the major products are the release of nonmodified nucleobases (NB; Cyt ≫ Thy ∼ Ade > Gua). Additional damage includes 5,6-dihydropyrimidines (dHT > dHU) and eight nucleosides with modified sugar moieties consisting of 2',3'- and 2',5'-dideoxynucleosides (ddG > ddA ∼ ddC > ddT). The distribution of products is remarkably different in a 16-mer ODN compared to that observed previously with thymidylyl-(3'-5')-thymidine. This difference is explained by electron delocalization occurring within a sufficiently long strand, the DEA theory of O'Malley, and recent time-dependent density functional theory calculations.

Proton-Transfer Reactions in One-Electron-Oxidized G-Quadruplexes: A Density Functional Theory Study

Recently, G-quadruplexes (Gq) formed in B-DNA as secondary structures are found to be important therapeutic targets and material for developing nanodevices. Gq are guanine-rich and thus susceptible to oxidative damage by producing short-lived intermediate radicals via proton-transfer reactions. Understanding the mechanisms of radical formation in Gq is of fundamental interest to understand the early stages of DNA damage. Herein, we used density functional theory including aqueous phase (ωB97XD-PCM/6-31++G**) and considered single layer of Gq [G-quartets (G4): association of four guanines in a cyclic Hoogsteen hydrogen-bonded arrangement (Scheme 1)] to unravel the mechanisms of formation of intermediates by calculating the relative Gibbs free energies and spin density distributions of one-electron-oxidized G4 and its various proton-transfer states: G^•+, G(N₁-H)^•, G(N₂-H')^•, G(N₂-H″)^•, G(N₁-H)^•-(H⁺O₆)G, and G(N₂-H)^•-(H⁺N₇)G. The present calculation predicts the formation of G(N₂-H)^•-(H⁺N₇)G, which is only ca. 0.8 kcal/mol higher in energy than the initially formed G^•+. The formation of G(N₂-H)^•-(H⁺N₇)G plays a key role in explaining the formation of 8-OG along with G(N₁-H)^• formation via tautomerization from G(N₂-H)^•, as proposed recently.

Deadlocks of adenine ribonucleotide synthesis: evaluation of adsorption and condensation reactions in a zeolite micropore space

Herein, we report on adenine, d-ribose, and monophosphate adsorption/co-adsorption into the synthetic analog of the zeolite mineral mordenite followed by drying at 50 °C and thermal activation at 150 °C under an argon atmosphere.

Reaction mechanism and dynamics for C8-hydroxylation of 9-methylguanine radical cation by water molecules

In contrast to their spontaneous deprotonation in aqueous solution, reactions of guanine and guanosine radical cations with water in the gas phase are exclusively initiated by hydration of the radical cations as reported in recent work (Y. Sun et al., Phys. Chem. Chem. Phys., 2018, 20, 27510). As gas-phase hydration reactions closely mimic the actual scenario for guanine radical cations in double-stranded DNA, exploration of subsequent reactions within their water complexes can provide an insight into the resulting oxidative damage to nucleosides. Herein guided-ion beam mass spectrometry experiment and direct dynamics trajectory simulations were carried out to examine prototype complexes of the 9-methylguanine radical cation with one and two water ligands (i.e., 9MG˙⁺·(H₂O)_1-2) in the gas phase, wherein the complexes were activated by collisional activation in the experiment and by thermal excitation at high temperatures in the simulations. Guided by mass spectroscopic measurements, trajectory results and reaction potential energy surface, three reaction pathways were identified. The first two reaction pathways start with H-atom abstraction from water by the O6 and N7 atoms in 9MG˙⁺ and are referred to as HA_O6 and HA_N7, respectively. The primary products of HA_O6 and HA_N7 reactions, including [9MG + H_O6]⁺/[9MG + H_N7]⁺ and ˙OH, react further to either form [8OH-9MG + H_O6]˙⁺ and [8OH-9MG + H_N7]˙⁺via C8-hydroxylation or form radical cations of 6-enol-guanine (6-enol-G˙⁺) and 7H-guanine (7HG˙⁺) via S_N2-type methanol elimination. The third reaction pathway corresponds to the formation of 8OH-9MG⁺ by H elimination from the complex, referred to as HE. Among these product channels, [8OH-9MG + H_N7]˙⁺ has the most favorable formation probability, especially in the presence of additional water molecules. This product may serve as a preceding structure to the 8-oxo-7,8-dihydroguanine lesion in DNA and has implications for health effects of radiation exposure and radiation therapy.

The Significance of Nonlinear Screening and the pH Interference Mechanism in Field-Effect Transistor Molecular Sensors

Electrolyte screening is well known for its detrimental impact on the sensitivity of liquid-gated field-effect transistor (FET) molecular sensors and is mostly described by the linearized Debye-Hückel model. However, charged and pH-sensitive FET sensing surfaces can limit the FET molecular sensitivity beyond the Debye-Hückel screening formalism. Pre-existing surface charges can lead to the breakdown of Debye-Hückel screening and induce enhanced nonlinear Poisson-Boltzmann screening. Moreover, the charging of the pH-sensitive surface groups interferes with biomolecule sensing resulting in a pH interference mechanism. With analytical equations and TCAD simulations, we highlight that the Debye-Hückel approximation can underestimate screening and overestimate FET molecular sensitivity by more than an order of magnitude. Screening strengthens significantly beyond Debye-Hückel in the proximity of even moderately charged surfaces and biomolecule charge densities (≥1 × 10¹² q/cm²). We experimentally show the strong impact of both nonlinear screening and the pH interference effect on charge-based biomolecular sensing using a model system based on the covalent binding of single-stranded DNA on silicon FET sensors. The DNA signal increases from 24 mV at pH 7 to 96 mV at pH 3 in 1.5 mM PBS for a DNA density of 7 × 10¹² DNA/cm². Our model quantitatively explains the signal's pH dependence with roughly equal nonlinear screening and pH interference contributions. This work shows the importance of reducing the net charge and the pH sensitivity of the sensing surface to improve molecular sensing. Therefore, tailoring the gate dielectric and functional layer of FET sensors is a promising route to strong silicon FET molecular sensitivity boosts.

Thymine-modified chitosan with broad-spectrum antimicrobial activities for wound healing

Avoiding wound infections is a major challenge in wound care management, and new materials are urgently needed to address these problems. Herein, four water-soluble thymine-modified chitosan (TC) derivatives with the degree of substitution (DS) ranging from 0.23 to 0.62 were synthesized and freeze-dried to obtain porous sponge dressings. TC derivatives had broad-spectrum antibacterial activities against gram-negative bacteria, gram-positive bacteria, fungi, drug-resistance bacteria, Pseudomonas aeruginosa and Acinetobacter baumannii. Especially, the antioxidant and antibacterial properties of TC derivatives increased with increasing DS. Furthermore, TC derivatives showed excellent biocompatibility and blood compatibility. TC sponges could significantly accelerate the wound healing process than gauze and chitosan sponge. The histological analysis and immunohistochemical staining showed that the wounds treated with TC sponges displayed fewer inflammatory cells, and faster regeneration of epithelial tissue, collagen deposition and new blood vessel formation speed. Therefore, TC sponges can shed new light for wound dressing design.

Theoretical insight into 7,8-dihydrogen-8-oxoguanine radical cation deprotonation

The pK_a values of reactive protons in 8-oxoG˙⁺ and potential energy profiles for 8-oxoG radical cation deprotonation reaction (N1–H and N7–H) were firstly calculated.

Gold Nanoparticles as Colorimetric Sensors for the Detection of DNA Bases and Related Compounds

Results regarding interaction of colloidal gold solutions with nucleobases, including uracil (U), as well as its sulfur derivatives, 2-thiouracil (2TU) and 4-thiouracil (4TU), cytosine (C), adenine (A), and guanine (G), as well as urea and thiourea (TU), are reported. Anionic stabilized citrate gold nanoparticles (AuNPs) were synthesized by reducing the tetrachloroaurate (III) trihydrate with trisodium citrate. The surface plasmon resonance (SPR) band was used in the characterization of synthesized AuNPs, as well as transmission electron microscope (TEM) imaging, which was used in the characterization of dispersed and aggregated gold nanoparticles. Interactions of nucleobases with the gold surface was analyzed by following the plasmon absorbance band red shift of the AuNPs. The sulfur-containing compounds adsorbed to the nanoparticle surfaces by chemisorption-type interactions; with TU and 4TU, the process is accompanied by a sudden change in color; in contrast, 2TU forms stable functionalized gold nanoparticles. Urea and U do not adsorb to nanoparticle surfaces, but the other heterocyclic bases containing nitrogen interact effectively with the gold surface, causing the assembly of nanoparticles, even though the interparticle self-aggregation process was slower than that mediated by either TU or 4TU. The method is efficient in the colorimetric detection of nucleobases and derivatives at concentration levels on the order of 1 µM.

One-electron oxidation of TAT-motif triplex DNA and the ensuing Hoogsteen hydrogen-bonding dissociation

One-electron oxidation of adenine (A) leads initially to the formation of adenine radical cation (A^•+). Subsequent deprotonation of A^•+ can provoke deoxyribonucleic acid (DNA) damage, which further causes senescence, cancer formation, and even cell death. However, compared with considerable reports on A^•+ reactions in free deoxyadenosine (dA) and duplex DNA, studies in non-B-form DNA that play critical biological roles are rare at present. It is thus of vital importance to explore non-B-form DNA, among which the triplex is an emerging topic. Herein, we investigate the deprotonation behavior of A^•+ in the TAT triplex with continuous A bases by time-resolved laser flash photolysis. The rate constants for the one-oxidation of triplex 8.4 × 10⁸ M^-1 s^-1 and A^•+ deprotonation 1.3 × 10⁷ s^-1 are obtained. The kinetic isotope effect of A^•+ deprotonation in the TAT triplex is 1.8, which is characteristic of a direct release of the proton into the solvent similar to free base dA. It is thus elucidated that the A^•+ proton bound with the third strand is most likely to be released into the solvent because of the weaker Hoogsteen H-bonding interaction and the presence of the highly mobile hydration waters within the third strand. Additionally, it is confirmed through Fourier transform infrared spectroscopy that the deprotonation of A^•+ results in the dissociation of the third strand and disruption of the secondary structure of the triplex. These results provide valuable kinetic data and in-depth mechanistic insights for understanding the adenine oxidative DNA damage in the triplex.

Proton Transfer Accompanied by the Oxidation of Adenosine

Despite numerous experimental and theoretical studies, the proton transfer accompanying the oxidation of 2'-deoxyadenosine 5'-monophosphate 2'-deoxyadenosine 5'-monophosphate (5'-dAMP, A) is still under debate. To address this issue, we have investigated the oxidation of A in acidic and neutral solutions by using transient absorption (TA) and time-resolved resonance Raman (TR³ ) spectroscopic methods in combination with pulse radiolysis. The steady-state Raman signal of A was significantly affected by the solution pH, but not by the concentration of adenosine (2-50 mm). More specifically, the A in acidic and neutral solutions exists in its protonated (AH⁺ (N1+H⁺ )) and neutral (A) forms, respectively. On the one hand, the TA spectral changes observed at neutral pH revealed that the radical cation (A^.+ ) generated by pulse radiolysis is rapidly converted into A^. (N6-H) through the loss of an imino proton from N6. In contrast, at acidic pH (<4), AH^.2+ (N1+H⁺ ) generated by pulse radiolysis of AH⁺ (N1+H⁺ ) does not undergo the deprotonation process owing to the pK_a value of AH^.2+ (N1+H⁺ ), which is higher than the solution pH. Furthermore, the results presented in this study have demonstrated that A, AH⁺ (N1+H⁺ ), and their radical species exist as monomers in the concentration range of 2-50 mm. Compared with the Raman bands of AH⁺ (N1+H⁺ ), the TR³ bands of AH^.2+ (N1+H⁺ ) are significantly down-shifted, indicating a decrease in the bond order of the pyrimidine and imidazole rings due to the resonance structure of AH^.2+ (N1+H⁺ ). Meanwhile, A^. (N6-H) does not show a Raman band corresponding to the pyrimidine+NH₂ scissoring vibration due to diprotonation at the N6 position. These results support the final products generated by the oxidation of adenosine in acidic and neutral solutions being AH^.2+ (N1+H⁺ ) and A^. (N6-H), respectively.

Direct observation of the oxidation of DNA bases by phosphate radicals formed under radiation: a model of the backbone-to-base hole transfer

In irradiated DNA, by the base-to-base and backbone-to-base hole transfer processes, the hole (i.e., the unpaired spin) localizes on the most electropositive base, guanine. Phosphate radicals formed via ionization events in the DNA-backbone must play an important role in the backbone-to-base hole transfer process. However, earlier studies on irradiated hydrated DNA, on irradiated DNA-models in frozen aqueous solution and in neat dimethyl phosphate showed the formation of carbon-centered radicals and not phosphate radicals. Therefore, to model the backbone-to-base hole transfer process, we report picosecond pulse radiolysis studies of the reactions between H2PO4˙ with the DNA bases - G, A, T, and C in 6 M H3PO4 at 22 °C. The time-resolved observations show that in 6 M H3PO4, H2PO4˙ causes the one-electron oxidation of adenine, guanine and thymine, by forming the cation radicals via a single electron transfer (SET) process; however, the rate constant of the reaction of H2PO4˙ with cytosine is too low (<107 L mol-1 s-1) to be measured. The rates of these reactions are influenced by the protonation states and the reorganization energies of the base radicals and of the phosphate radical in 6 M H3PO4.

Modeling p K Shift in DNA Triplexes Containing Locked Nucleic Acids

The protonation states for nucleic acid bases are difficult to assess experimentally. In the context of DNA triplex, the protonation state of cytidine in the third strand is particularly important, because it needs to be protonated in order to form Hoogsteen hydrogen bonds. A sugar modification, locked nucleic acid (LNA), is widely used in triplex forming oligonucleotides to target sites in the human genome. In this study, the parameters for LNA are developed in line with the CHARMM nucleic acid force field and validated toward the available structural experimental data. In conjunction, two computational methods were used to calculate the protonation state of the third strand cytidine in various DNA triplex environments: λ-dynamics and multiple pH regime. Both approaches predict p K of this cytidine shifted above physiological pH when cytidine is in the third strand in a triplex environment. Both methods show an upshift due to cytidine methylation, and a small downshift when the sugar configuration is locked. The predicted p K values for cytidine in DNA triplex environment can inform the design of better-binding oligonucleotides.

Photoactivated proton coupled electron transfer in DNA: insights from quantum mechanical calculations

The energetics of the two main proton coupled electron transfer processes that could occur in DNA are determined by means of time dependent-DFT calculations, using the M052X functional and the polarizable continuum model to include solvent effect.

Aminyl Radical Generation via Tandem Norrish Type I Photocleavage, β-Fragmentation: Independent Generation and Reactivity of the 2'-Deoxyadenosin- N6-yl Radical

Formal hydrogen atom abstraction from the nitrogen-hydrogen bonds in purine nucleosides produces reactive intermediates that are important in nucleic acid oxidation. Herein we describe an approach for the independent generation of the purine radical resulting from hydrogen atom abstraction from the N6-amine of 2'-deoxyadenosine (dA•). The method involves sequential Norrish Type I photocleavage of a ketone (7b) and β-fragmentation of the initially formed alkyl radical (8b) to form dA• and acetone. The formation of dA• was followed by laser flash photolysis, which yields a transient with λ_max ≈ 340 nm and a broader weaker absorption centered at ∼560 nm. This transient grows in at ≥2 × 10⁵ s^-1; however, computations and reactivity data suggest that β-fragmentation occurs much faster, implying the consumption of dA• as it is formed. Continuous photolysis of 7b in the presence of ferrous ion or thiophenol produces good yields of dA, whereas less reactive thiols afford lower yields presumably due to a polarity mismatch. This tandem photochemical, β-fragmentation method promises to be useful for site-specific production of dA• in nucleic acid oligomers and/or polymers and also for the production of aminyl radicals, in general.

DNA-stabilized Ag-Au bimetallic clusters: the effects of alloying and embedding on optical properties

Global geometry optimization and time-dependent density functional theory calculations have been used to study the structural evolution and optical properties of AgnAun (n = 2-6) nanoalloys both as individual clusters and as clusters stabilized with the fragments of DNA of different size. We show that alloying can be used to control and tune the level of interaction between the metal atoms of the cluster and the organic fragments of the DNA ligands. For instance, gold and silver atoms are shown to exhibit synergistic effects in the process of charge transfer from the nucleobase to the cluster, with the silver atoms directly connected to the nitrogen atoms of cytosine increasing their positive partial charge, while their more electronegative neighbouring gold atoms host the excess negative charge. This allows the geometrical structures and optical absorption spectra of small bimetallic clusters to retain many of their main features upon aggregation with relatively large DNA fragments, such as a cytosine-based 9-nucleotide hairpin loop, which suggests a potential synthetic route to such hybrid metal-organic compounds, and opens up the possibility of bringing the unique tunable properties of bimetallic nanoalloys to biological applications.

Calculations of the Energetics of Oxidation of Aqueous Nucleosides and the Effects of Prototropic Equilibria

Recently the calculated standard reduction potentials of the radical-cations of N-methyl substituted DNA bases have been reported that agree fairly well with the experimental results. However, there are issues reflecting the fact that the experimental results usually relate to the couple E(o)(Nuc(•),H(+)/NucH(+)), whereas the calculated results are for the E(o)(Nuc(•+)/Nuc) couple. To calculate the midpoint reduction potential at pH 7 (Em7), it is important to have accurate acid dissociation constants (pKs) for both ground-state bases and their radicals, and the effects of uncertainty in some of these values (e.g., that of the adenosine radical) must be considered. Calculations of the pKs of the radicals of the nucleic acid bases (as nucleosides) have been performed to explore the effects the various pKs have on calculating the values of Em7 and to see what improvements can be made with the accuracy of the calculations.

Gamma and Ion-Beam Irradiation of DNA: Free Radical Mechanisms, Electron Effects, and Radiation Chemical Track Structure

The focus of our laboratory's investigation is to study the direct-type DNA damage mechanisms resulting from γ-ray and ion-beam radiation-induced free radical processes in DNA which lead to molecular damage important to cellular survival. This work compares the results of low LET (γ-) and high LET (ion-beam) radiation to develop a chemical track structure model for ion-beam radiation damage to DNA. Recent studies on protonation states of cytosine cation radicals in the N1-substituted cytosine derivatives in their ground state and 5-methylcytosine cation radicals in ground as well as in excited state are described. Our results exhibit a radical signature of excitations in 5-methylcytosine cation radical. Moreover, our recent theoretical studies elucidate the role of electron-induced reactions (low energy electrons (LEE), presolvated electrons (e_pre^-), and aqueous (or, solvated) electrons (e_aq^-)). Finally DFT calculations of the ionization potentials of various sugar radicals show the relative reactivity of these species.

Efficient and Substantial DNA Lesions From Near 0 eV Electron-Induced Decay of the O4-Hydrogenated Thymine Nucleotides: A DFT Study

Possible electron-induced ruptures of C3'-O3', C5'-O5', and N1-C1' bonds in O4-hydrogenated 2'-deoxythymidine-3'-monophosphate (3'-dT(O4H)MPH) and 2'-deoxythymidine-5'-monophosphate (5'-dT(O4H)MPH) are investigated using density functional theory calculations, and efficient pathways are proposed. Electron attachment causes remarkable structural relaxation in the thymine C6 site. A concerted process of intramolecular proton transfer (IPT) from the C2' site of 2'-deoxyribose to the C6 site and the C3'-O3' bond rupture is observed in [3'-dT(O4H)MPH](-). A low activation barrier (9.32 kcal/mol) indicates that this pathway is the most efficient one as compared to other known pathways leading to backbone breaks of a single strand DNA at the non-3'-end thymine, which prevents the N1-C1' bond cleavage in [3'-dT(O4H)MPH](-). However, essentially spontaneous N1-C1' bond cleavage following similar IPT is predicted in [5'-dT(O4H)MPH](-). A moderate activation barrier (13.02 kcal/mol) for the rate-controlling IPT step suggests that base release from the N1-C1' cleavage arises readily at the 3'-end of single strand DNA with the strand ended by a thymine. The C5'-O5' bond has only an insignificant change in the IPT process. Solvent effects are found to increase slightly the energy requirements for either bond ruptures (11.23 kcal/mol (C3'-O3') vs 16.18 kcal/mol (N1-C1')), but not change their relative efficiencies.

π-Radical to σ-Radical Tautomerization in One-Electron-Oxidized 1-Methylcytosine and Its Analogs

In this work, iminyl σ-radical formation in several one-electron-oxidized cytosine analogs, including 1-MeC, cidofovir, 2'-deoxycytidine (dCyd), and 2'-deoxycytidine 5'-monophosphate (5'-dCMP), were investigated in homogeneous, aqueous (D2O or H2O) glassy solutions at low temperatures by employing electron spin resonance (ESR) spectroscopy. Upon employing density functional theory (DFT) (DFT/B3LYP/6-31G* method), the calculated hyperfine coupling constant (HFCC) values of iminyl σ-radical agree quite well with the experimentally observed ones, thus confirming its assignment. ESR and DFT studies show that the cytosine iminyl σ-radical is a tautomer of the deprotonated cytosine π-cation radical [cytosine π-aminyl radical, C(N4-H)(•)]. Employing 1-MeC samples at various pHs ranging from ca. 8 to 11, ESR studies show that the tautomeric equilibrium between C(N4-H)(•) and the iminyl σ-radical at low temperature is too slow to be established without added base. ESR and DFT studies agree that, in the iminyl σ-radical, the unpaired spin is localized on the exocyclic nitrogen (N4) in an in-plane pure p-orbital. This gives rise to an anisotropic nitrogen hyperfine coupling (Azz = 40 G) from N4 and a near isotropic β-nitrogen coupling of 9.7 G from the cytosine ring nitrogen at N3. Iminyl σ-radical should exist in its N3-protonated form, as the N3-protonated iminyl σ-radical is stabilized in solution by over 30 kcal/mol (ΔG = -32 kcal/mol) over its conjugate base, the N3-deprotonated form. This is the first observation of an isotropic β-hyperfine ring nitrogen coupling in an N-centered DNA radical. Our theoretical calculations predict that the cytosine iminyl σ-radical can be formed in double-stranded DNA by a radiation-induced ionization-deprotonation process that is only 10 kcal/mol above the lowest energy path.

Feasibility of occurrence of different types of protonated base pairs in RNA: a quantum chemical study

Protonated nucleobases have significant roles in facilitating catalytic functions of RNA, and in stabilizing different structural motifs. Reported pKa values of nucleobase protonation suggest that the population of neutral nucleobases is 10(3)-10(4) times higher than that of protonated nucleobases under physiological conditions (pH ∼ 7.4). Therefore, a molecular level understanding of various putative roles of protonated nucleobases cannot be achieved without addressing the question of how their occurrence propensities and stabilities are related to the free energy costs associated with the process of protonation under physiological conditions. With water as the proton donor, we use advanced QM methods to evaluate the site specific protonation propensities of nucleobases in terms of their associated free energy changes (ΔGprot). Quantitative follow up on the energetics of base pair formation and database search for evaluating their occurrence frequencies, reveal a lack of correlation between base pair stability and occurrence propensities on the one hand, and ease of protonation on the other. For example, although N7 protonated adenine (ΔGprot = 40.0 kcal mol(-1)) is found to participate in stable base pairing, base pairs involving N7 protonated guanine (ΔGprot = 36.8 kcal mol(-1)), on geometry optimization, converge to a minima where guanine transfers its extra proton to its partner base. Such observations, along with examples of weak base pairs involving N3 protonation of cytosine (ΔGprot = 37.0 kcal mol(-1)) are rationalized by analysing the protonation induced charge redistributions which are found to significantly influence, both positively and negatively, the hydrogen bonding potentials of different functional sites of individual nucleobases. Protonation induced charge redistribution is also found to strongly influence (i) the aromatic character of the rings of the participating bases and (ii) hydrogen bonding potential of the free edges of the protonated base pair. Comprehensive analysis of a non-redundant RNA crystal structure dataset further reveals that, while availability of stabilization possibilities determine the feasibility of occurrence of protonated bases, their occurrence context and specific functional roles are important factors determining their occurrence propensities.

Calculations of pKa's and redox potentials of nucleobases with explicit waters and polarizable continuum solvation

The SMD implicit solvation model augmented with one and four explicit water molecules was used to calculate pKa's and redox potentials of N-methyl-substituted nucleic acid bases guanine, adenine, cytosine, thymine, and uracil. Calculations were carried out with the B3LYP/6-31+G(d,p) level of theory. The same numbers of water molecules were hydrogen bonded to the neutral, protonated, and deprotonated nucleobases in their unoxidized and oxidized forms. The improvement in pKa1 involving neutrals and cations was modest. By contrast, the improvement in pKa2 involving neutrals and anions was quite significant, reducing the mean absolute error from 4.6 pKa units with no waters, to 2.6 with one water and 1.7 with four waters. For the oxidation of nucleobases, adding explicit waters did little to improve E(X(•),H(+)/XH), possibly because both species in the redox couple are neutral molecules at pH 7.

Proton transfer induced SOMO-to-HOMO level switching in one-electron oxidized A-T and G-C base pairs: a density functional theory study

In the present study, we show that for one-electron oxidized A-T or G-C base pairs the singly occupied molecular orbital (SOMO) is located on A or G and is lower in energy than the doubly occupied highest-occupied molecular orbital (HOMO) localized to the pyrimidines, T or C. This directs second ionizations to the pyrimidine bases resulting in triplet state diradical dications, (A^•+-T^•+) and (G^•+-C^•+). On interbase proton transfer, the SOMO and HOMO levels switch and the second oxidation is redirected to G and A. For G-C, the doubly oxidized singlet G(-H)⁺-C(H⁺) is more stable than its triplet (G^•+-C^•+); however, for A-T, the triplet (A^•+-T^•+) lies lowest in energy. The study demonstrates that double ionization of the A-T base pair results in a triplet dication diradical, which is more stable than the proton-transferred triplet or singlet species; whereas, double ionization of the G-C base pair, the proton transferred doubly oxidized singlet, G(-H)⁺-C(H⁺), is more stable and has both oxidations on guanine. In DNA, with both A-T and G-C, multiple oxidations would transfer to the guanine base alone.

π- vs σ-radical states of one-electron-oxidized DNA/RNA bases: a density functional theory study

As a result of their inherent planarity, DNA base radicals generated by one-electron oxidation/reduction or bond cleavage form π- or σ-radicals. While most DNA base systems form π-radicals, there are a number of nucleobase analogues such as one-electron-oxidized 6-azauraci1, 6-azacytosine, and 2-thiothymine or one-electron reduced 5-bromouracil that form more reactive σ-radicals. Elucidating the availability of these states within DNA, base radical electronic structure is important to the understanding of the reactivity of DNA base radicals in different environments. In this work, we address this question by the calculation of the relative energies of π- and σ-radical states in DNA/RNA bases and their analogues. We used density functional theory B3LYP/6-31++G** method to optimize the geometries of π- and σ-radicals in Cs symmetry (i.e., planar) in the gas phase and in solution using the polarized continuum model (PCM). The calculations predict that σ- and π-radical states in one-electron-oxidized bases of thymine, T(N3-H)(•), and uracil, U(N3-H)(•), are very close in energy; i.e., the π-radical is only ca. 4 kcal/mol more stable than the σ-radical. For the one-electron-oxidized radicals of cytosine, C(•+), C(N4-H)(•), adenine, A(•+), A(N6-H)(•), and guanine, G(•+), G(N2-H)(•), G(N1-H)(•), the π-radicals are ca. 16-41 kcal/mol more stable than their corresponding σ-radicals. Inclusion of solvent (PCM) is found to stabilize the π- over σ-radical of each of the systems. U(N3-H)(•) with three discrete water molecules in the gas phase is found to form a three-electron σ bond between the N3 atom of uracil and the O atom of a water molecule, but on inclusion of full solvation and discrete hydration, the π-radical remains most stable.