Ab Initio Ionization Energy Thresholds of DNA and RNA Bases in Gas Phase and in Aqueous Solution

Ionization Patterns and Chemical Reactivity of Cytosine-Guanine Watson-Crick Pairs

Gas-phase and aqueous vertical ionization potentials, vIPgas and vIPaq respectively and measurements of the molecular electrostatic and local ionization maps calculated at the DFT/B3LYP-D3/ 6-311+G** level of theory and the C-PCM reaction field model for single- and double-stranded CpG and 5MeCpG pairs show that the vIPaq for single- and double-stranded pairs of C-G and 5MeC-G are practically the same, in the range of 5.79 to 5.81 eV. The aqueous adiabatic ionization potentials for single-stranded CpG and 5MeCpG are 5.52 eV and 5.51 eV respectively and they reflect the nuclear reorganization that takes place after the abstraction of the electron. The aqueous adiabatic ionization energy values that correspond to the CpG+. radical cation and the hydrated electron, e-,, being at infinite distance, adIPaq+Vo, are 3.92 eV and 3.91 eV respectively with (Vo=-1.6 eV) Analysis of data suggest that the HOMO-LUMO energy gap in the hard/soft-acid/base (HSAB) concept cannot be used a priori to determine the effect of cytosine methylation on the guanine enhanced oxidative damage in DNA.

Electron interaction with DNA constituents in aqueous phase

Electron driven chemistry of biomolecules in aqueous phase presents the realistic picture to study molecular processes. In this study we have investigated the interactions of electrons with the DNA constituents in their aqueous phase in order to obtain the quantities useful for DNA damage assessment. We have computed the inelastic mean free path (IMFP), mass stopping power (MSP) and absorbed dose (D) for the DNA constituents (Adenine, Cytosine, Guanine, Thymine and Uracil) in the aqueous medium from ionisation threshold to 5000 eV. We have modified complex optical potential formalism to include band gap of the systems to calculate inelastic cross sections which are used to estimate these entities. This is the maiden attempt to report these important quantities for the aqueous DNA constituents. We have compared our results with available data in gas and other phase and have observed explicable accord for IMFP and MSP. Since these are the first results of absorbed dose (D) for these compounds, we have explored present results vis-a-vis dose absorption in water.

Effects of Intra-Base Pair Proton Transfer on Dissociation and Singlet Oxygenation of 9-Methyl-8-Oxoguanine-1-Methyl-Cytosine Base-Pair Radical Cations

8-Oxoguanosine is the most common oxidatively generated base damage and pairs with complementary cytidine within duplex DNA. The 8-oxoguanosine-cytidine lesion, if not recognized and removed, not only leads to G-to-T transversion mutations but renders the base pair being more vulnerable to the ionizing radiation and singlet oxygen (1 O2 ) damage. Herein, reaction dynamics of a prototype Watson-Crick base pair [9MOG ⋅ 1MC]⋅+ , consisting of 9-methyl-8-oxoguanine radical cation (9MOG⋅+ ) and 1-methylcystosine (1MC), was examined using mass spectrometry coupled with electrospray ionization. We first detected base-pair dissociation in collisions with the Xe gas, which provided insight into intra-base pair proton transfer of 9MOG⋅+  ⋅ 1MC← → ${{\stackrel{ {\rightarrow} } { {\leftarrow} } } }$ [9MOG - HN1 ]⋅ ⋅ [1MC+HN3' ]+ and subsequent non-statistical base-pair separation. We then measured the reaction of [9MOG ⋅ 1MC]⋅+ with 1 O2 , revealing the two most probable pathways, C5-O2 addition and HN7 -abstraction at 9MOG. Reactions were entangled with the two forms of 9MOG radicals and base-pair structures as well as multi-configurations between open-shell radicals and 1 O2 (that has a mixed singlet/triplet character). These were disentangled by utilizing approximately spin-projected density functional theory, coupled-cluster theory and multi-referential electronic structure modeling. The work delineated base-pair structural context effects and determined relative reactivity toward 1 O2 as [9MOG - H]⋅>9MOG⋅+ >[9MOG - HN1 ]⋅ ⋅ [1MC+HN3' ]+ ≥9MOG⋅+  ⋅ 1MC.

Spin-orbit charge transfer from guanine and 9-methylguanine radical cations to nitric oxide radicals and the induced triplet-to-singlet intersystem crossing

Nitric oxide (●NO) participates in many biological activities, including enhancing DNA radiosensitivity in ionizing radiation-based radiotherapy. To help understand the radiosensitization of ●NO, we report reaction dynamics between ●NO and the radical cations of guanine (a 9HG●+ conformer) and 9-methylguanine (9MG●+). On the basis of the formation of 9HG●+ and 9MG●+ in the gas phase and the collisions of the radical cations with ●NO in a guided-ion beam mass spectrometer, the charge transfer reactions of 9HG●+ and 9MG●+ with ●NO were examined. For both reactions, the kinetic energy-dependent product ion cross sections revealed a threshold energy that is 0.24 (or 0.37) eV above the 0 K product 9HG (or 9MG) + NO+ asymptote. To interrogate this abnormal threshold behavior, the reaction potential energy surface for [9MG + NO]+ was mapped out at closed-shell singlet, open-shell singlet, and triplet states using density functional and coupled cluster theories. The results showed that the charge transfer reaction requires the interaction of a triplet-state surface originating from a reactant-like precursor complex 3[9MG●+(↑)⋅(↑)●NO] with a closed-shell singlet-state surface evolving from a charge-transferred complex 1[9MG⋅NO+]. During the reaction, an electron is transferred from π∗(NO) to perpendicular π∗(9MG), which introduces a change in orbital angular momentum. The latter offsets the change in electron spin angular momentum and facilitates intersystem crossing. The reaction threshold in excess of the 0 K thermochemistry and the low charge-transfer efficiency are rationalized by the vibrational excitation in the product ion NO+ and the kinetic shift arising from a long-lived triplet intermediate.

Intramolecular and intermolecular hole delocalization rules the reducer character of isolated nucleobases and homogeneous single-stranded DNA

The use of DNA strands as nanowires or electrochemical biosensors requires a deep understanding of charge transfer processes along the strand, as well as of the redox properties. These properties are computationally assessed in detail throughout this study. By applying molecular dynamics and hybrid QM/continuum and QM/QM/continuum schemes, the vertical ionization energies, adiabatic ionization energies, vertical attachment energies, one-electron oxidation potentials, and delocalization of the hole generated upon oxidation have been determined for nucleobases in their free form and as part of a pure single-stranded DNA. We show that the reducer ability of the isolated nucleobases is explained by the intramolecular delocalization of the positively charged hole, while the enhancement of the reducer character when going from aqueous solution to the strand correlates very well with the intermolecular hole delocalization. Our simulations suggest that the redox properties of DNA strands can be tuned by playing with the balance between intramolecular and intermolecular charge delocalization.

Singlet O2 Oxidation of the Radical Cation versus the Dehydrogenated Neutral Radical of 9-Methylguanine in a Watson-Crick Base Pair. Consequences of Structural Context

In DNA, guanine is the most susceptible to oxidative damage by exogenously and endogenously produced electronically excited singlet oxygen (¹O₂). The reaction mechanism and the product outcome strongly depend on the nucleobase ionization state and structural context. Previously, exposure of a monomeric 9-methylguanine radical cation (9MG^•+, a model guanosine compound) to ¹O₂ was found to result in the formation of an 8-peroxide as the initial product. The present work explores the ¹O₂ oxidation of 9MG^•+ and its dehydrogenated neutral form [9MG - H]^• within a Watson-Crick base pair consisting of one-electron-oxidized 9-methylguanine-1-methylcytosine [9MG·1MC]^•+. Emphasis is placed on entangling the base pair structural context and intra-base pair proton transfer with and consequences thereof on the singlet oxygenation of guanine radical species. Electrospray ionization coupled with guided-ion beam tandem mass spectrometry was used to study the formation and reaction of guanine radical species in the gas phase. The ¹O₂ oxidation of both 9MG^•+ and [9MG - H]^• is exothermic and proceeds barrierlessly either in an isolated monomer or within a base pair. Single- and multi-referential theories were tested for treating spin contaminations and multi-configurations occurring in radical-¹O₂ interactions, and reaction potential energy surfaces were mapped out to support experimental findings. The work provides a comprehensive profile for the singlet oxygenation of guanine radicals in different charge states and in the absence and the presence of base pairing. All results point to an 8-peroxide as the major oxidation product in the experiment, and the oxidation becomes slightly more favorable in a neutral radical form. On the basis of a variety of reaction pathways and product profiles observed in the present and previous studies, the interplay between guanine structure, base pairing, and singlet oxygenation and its biological implications are discussed.

Nucleic Acids under Stress: Understanding and Simulating Nucleobase Fragmentation Pathways

The effects of radiations on nucleic acids and their constituents is widely studied across several research fields using different experimental and theoretical protocols. While a large number of studies were performed in this context, many fundamental physical and chemical effects are still being investigated, particularly involving the effect of the biological environment. As an example, the interpretation of experimental nucleic acid bases mass spectra, and hence inferring their reactivity in complex environment still poses great challenge. This Minireview summarizes recent theoretical advancements aiming to predict and interpret the reactivity of nucleic acid bases. We focus not only on the understanding of the inherent fragmentation pathways of isolated nucleobases but also on the modeling of a realistic nano-environments highlighting the importance of molecular dynamics simulations and the non-innocent role of the environment and also the possibility to open novel fragmentation pathways.

Efficient treatment of molecular excitations in the liquid phase environment via stochastic many-body theory

Accurate predictions of charge excitation energies of molecules in the disordered condensed phase are central to the chemical reactivity, stability, and optoelectronic properties of molecules and critically depend on the specific environment. Herein, we develop a stochastic GW method for calculating these charge excitation energies. The approach employs maximally localized electronic states to define the electronic subspace of a molecule and the rest of the system, both of which are randomly sampled. We test the method on three solute-solvent systems: phenol, thymine, and phenylalanine in water. The results are in excellent agreement with the previous high-level calculations and available experimental data. The stochastic calculations for supercells containing up to 1000 electrons representing the solvated systems are inexpensive and require ≤1000 central processing unit hrs. We find that the coupling with the environment accounts for ∼40% of the total correlation energy. The solvent-to-solute feedback mechanism incorporated in the molecular correlation term causes up to 0.6 eV destabilization of the quasiparticle energy. Simulated photo-emission spectra exhibit red shifts, state-degeneracy lifting, and lifetime shortening. Our method provides an efficient approach for an accurate study of excitations of large molecules in realistic condensed phase environments.

Quantum Chemistry Calculations for Metabolomics

A primary goal of metabolomics studies is to fully characterize the small-molecule composition of complex biological and environmental samples. However, despite advances in analytical technologies over the past two decades, the majority of small molecules in complex samples are not readily identifiable due to the immense structural and chemical diversity present within the metabolome. Current gold-standard identification methods rely on reference libraries built using authentic chemical materials ("standards"), which are not available for most molecules. Computational quantum chemistry methods, which can be used to calculate chemical properties that are then measured by analytical platforms, offer an alternative route for building reference libraries, i.e., in silico libraries for "standards-free" identification. In this review, we cover the major roadblocks currently facing metabolomics and discuss applications where quantum chemistry calculations offer a solution. Several successful examples for nuclear magnetic resonance spectroscopy, ion mobility spectrometry, infrared spectroscopy, and mass spectrometry methods are reviewed. Finally, we consider current best practices, sources of error, and provide an outlook for quantum chemistry calculations in metabolomics studies. We expect this review will inspire researchers in the field of small-molecule identification to accelerate adoption of in silico methods for generation of reference libraries and to add quantum chemistry calculations as another tool at their disposal to characterize complex samples.

Ground and Excited States of Gas-Phase DNA Nucleobase Cation-Radicals. A UV-vis Photodisociation Action Spectroscopy and Computational Study of Adenine and 9-Methyladenine

Cation radicals of adenine (A^•+) and 9-methyladenine (MA^•+) were generated in the gas phase by collision-induced intramolecular electron transfer in copper-terpyridine-nucleobase ternary complexes and characterized by collision-induced dissociation (CID) mass spectra and UV-vis photodissociation action spectroscopy in the 210-700 nm wavelength region. The action spectra of both A^•+ and MA^•+ displayed characteristic absorption bands in the near-UV and visible regions. Another tautomer of A^•+ was generated as a minor product by multistep CID of protonated 9-(2-bromoethyl)adenine. Structure analysis by density functional theory and coupled-clusters ab initio calculations pointed to the canonical 9-H-tautomer Ad1^•+ as the global energy minimum of adenine cation radicals. The canonical tautomer MA1^•+ was also calculated to be a low-energy structure among methyladenine cation radicals. However, two new noncanonical tautomers were found to be energetically comparable to MA1^•+. Vibronic absorption spectra were calculated for several tautomers of A^•+ and MA^•+ and benchmarked on equation-of-motion coupled-clusters excited-state calculations. Analysis of the vibronic absorption spectra of A^•+ tautomers pointed to the canonical tautomer Ad1^•+ as providing the best match with the action spectrum. Likewise, the canonical tautomer MA1^•+ was the unequivocal best match for the MA^•+ ion generated in the gas phase. According to potential-energy mapping, MA1^•+ was separated from energetically favorable noncanonical cation radicals by a high-energy barrier that was calculated to be above the dissociation threshold for loss of a methyl hydrogen atom, thus preventing isomerization. Structures and energetics of all four DNA nucleobase cation radicals are compared and discussed.

One-electron oxidation of ds(5'-GGG-3') and ds(5'-G(8OG)G-3') and the nature of hole distribution: a density functional theory (DFT) study

Of particular interest in radiation-induced charge transfer processes in DNA is the extent of hole localization immediately after ionization and subsequent relaxation. To address this, we considered double stranded oligomers containing guanine (G) and 8-oxoguanine (8OG), i.e., ds(5'-GGG-3') and ds(5'-G8OGG-3') in B-DNA conformation. Using DFT, we calculated a variety of properties, viz., vertical and adiabatic ionization potentials, spin density distributions in oxidized stacks, solvent and solute reorganization energies and one-electron oxidation potential (E0) in the aqueous phase. Calculations for the vertical state of the -GGG- cation radical showed that the spin was found mainly (67%) on the middle G. However, upon relaxation to the adiabatic -GGG- cation radical, the spin localized (96%) on the 5'-G, as observed in experiments. Hole localizations on the middle G and 3'-G were higher in energy by 0.5 kcal mol-1 and 0.4 kcal mol-1, respectively, than that of 5'-G. In the -G8OGG- cation radical, the spin localized only on the 8OG in both vertical and adiabatic states. The calculated vertical ionization potentials of -GGG- and -G8OGG- stacks were found to be lower than that of the vertical ionization potential of a single G in DNA. The calculated E0 values of -GGG- and -G8OGG- stacks are 1.15 and 0.90 V, respectively, which owing to stacking effects are substantially lower than the corresponding experimental E0 values of their monomers (1.49 and 1.18 V, respectively). SOMO to HOMO level switching is observed in these oxidized stacks. Consequently, our calculations predict that local double oxidations in DNA will form triplet diradical states, which are especially significant for high LET radiations.

Electron transfer in biologically important systems: Polycyclic aromatic hydrocarbons, DNA bases and free radicals

Occurrence of electron transfer was studied for different combinations of polycyclic aromatic hydrocarbons (PAHs) and DNA bases as electron donors or acceptors and free radicals only as electron acceptors. Geometries of all the molecules and radicals were optimized in aqueous medium employing the polarizable continuum model. Single electron transfer (SET) and sequential proton loss electron transfer mechanisms were investigated employing Gibbs free energies of the appropriate neutral, anionic and cationic species. Barrier energies involved in these phenomena were calculated using the Marcus theory. The SET barrier energies were found to be linearly correlated with [Formula: see text] (Electron affinities of acceptors – Ionization potentials of donors). SET barrier energies from the DNA bases to the PAHs follow the order Cy [Formula: see text] Th [Formula: see text] Ad [Formula: see text] Gu, whereas SET barrier energies from the PAHs to the DNA bases follow the order Gu [Formula: see text] Ad [Formula: see text] Th [Formula: see text] Cy. Thus, guanine, among the DNA bases, is the best electron donor to the PAHs and worst electron acceptor from the same.

Quantitative ionization energies and work functions of aqueous solutions

Despite the ubiquitous nature of aqueous solutions across the chemical, biological and environmental sciences our experimental understanding of their electronic structure is rudimentary-qualitative at best. One of the most basic and seemingly straightforward properties of aqueous solutions-ionization energies-are (qualitatively) tabulated at the water-air interface for a mere handful of solutes, and the manner in which these results are obtained assume the aqueous solutions behave like a gas in the photoelectron experiment (where the vacuum levels of the aqueous solution and of the photoelectron analyzer are equilibrated). Here we report the experimental measure of a sizeable offset (ca. 0.6 eV) between the vacuum levels of an aqueous solution (0.05 M NaCl) and that of our photoelectron analyzer, indicating a breakdown of the gas-like vacuum level alignment assumption for the aqueous solution. By quantifying the vacuum level offset as a function of solution chemical composition our measurements enable, for the first time, quantitative determination of ionization energies in liquid solutions. These results reveal that the ionization energy of liquid water is not independent of the chemical composition of the solution as is usually inferred in the literature, a finding that has important ramifications as measured ionization energies are frequently used to validate theoretical models that posses the ability to provide microscopic insight not directly available by experiment. Finally, we derive the work function, or the electrochemical potential of the aqueous solution and show that it too varies with the chemical composition of the solution.

Controlled potential electro-oxidation of genomic DNA

Exposure of mammalian cells to oxidative stress can result in DNA damage that adversely affects many cell processes. Lack of dependable DNA damage reference materials and standardized measurement methods, despite many case-control studies hampers the wider recognition of the link between oxidatively degraded DNA and disease risk. We used bulk electrolysis in an electrochemical system and gas chromatographic mass spectrometric analysis (GC/MS/MS) to control and measure, respectively, the effect of electrochemically produced reactive oxygen species on calf thymus DNA (ct-DNA). DNA was electro-oxidized for 1 h at four fixed oxidizing potentials (E = 0.5 V, 1.0 V, 1.5 V and 2 V (vs Ag/AgCl)) using a high surface area boron-doped diamond (BDD) working electrode (WE) and the resulting DNA damage in the form of oxidatively-modified DNA lesions was measured using GC/MS/MS. We have shown that there are two distinct base lesion formation modes in the explored electrode potential range, corresponding to 0.5 V < E < 1.5 V and E > 1.5 V. Amounts of all four purine lesions were close to a negative control levels up to E = 1.5 V with evidence suggesting higher levels at the lowest potential of this range (E = 0.5 V). A rapid increase in all base lesion yields was measured when ct-DNA was exposed at E = 2 V, the potential at which hydroxyl radicals were efficiently produced by the BDD electrode. The present results demonstrate that controlled potential preparative electrooxidation of double-stranded DNA can be used to purposely increase the levels of oxidatively modified DNA lesions in discrete samples. It is envisioned that these DNA samples may potentially serve as analytical control or quality assurance reference materials for the determination of oxidatively induced DNA damage.

Quantum dynamics of a hole migration through DNA: A single strand DNA model

A model predicting the behavior of a hole acting on the DNA strand was investigated. The hole-DNA interaction on the basis of a quantum-classical, non-linear DNA single strand model was described. The fact that a DNA molecule is formed by a furanose ring as its sugar, phosphate group and bases was taken into consideration. Based on the model, results were obtained for the probability of a hole location on the DNA base sequences, such as GTTGGG, GATGTGGG, GTTGTTGGG as well as on the sugar-phosphate groups mated with them.

A Functional Group Approach for Prediction of APPI Response of Organic Synthetic Targets

Atmospheric pressure photoionization (APPI) is a technique of choice for ionization of non-polar molecules in mass spectrometry (MS). Reported APPI-based studies tend to focus on a selected compound class, which may contain a variety of functional groups. These studies demonstrate that APPI response frequently differs substantially, indicating a certain dependence on the functional group present. Although this dependence could be employed for APPI response prediction, its systematic use is currently absent. Here, we apply APPI MS to a judiciously-compiled set of 63 compounds containing a number of diverse functional groups commonly utilized in synthesis, reactive functional groups, as well as those containing boron and silicon. Based on the outcome of APPI MS of these compounds, we propose and evaluate a simple guideline to estimate the APPI response for a novel compound, the key properties of which have not been characterized in the gas phase. Briefly, we first identify key functional groups in the compound and gather knowledge on the known ionization energies from the smallest analogues containing said functional groups. We then consider local inductive and resonance effects on said ionization energies for the compounds of interest to estimate the APPI response. Finally, application of APPI MS to compounds of interest considered herein demonstrated extended upper mass ionization limit of 3.5 kDa for non-polymeric compounds.

Modeling photoionization of aqueous DNA and its components

Radiation damage to DNA is usually considered in terms of UVA and UVB radiation. These ultraviolet rays, which are part of the solar spectrum, can indeed cause chemical lesions in DNA, triggered by photoexcitation particularly in the UVB range. Damage can, however, be also caused by higher energy radiation, which can ionize directly the DNA or its immediate surroundings, leading to indirect damage. Thanks to absorption in the atmosphere, the intensity of such ionizing radiation is negligible in the solar spectrum at the surface of Earth. Nevertheless, such an ionizing scenario can become dangerously plausible for astronauts or flight personnel, as well as for persons present at nuclear power plant accidents. On the beneficial side, ionizing radiation is employed as means for destroying the DNA of cancer cells during radiation therapy. Quantitative information about ionization of DNA and its components is important not only for DNA radiation damage, but also for understanding redox properties of DNA in redox sensing or labeling, as well as charge migration along the double helix in nanoelectronics applications. Until recently, the vast majority of experimental and computational data on DNA ionization was pertinent to its components in the gas phase, which is far from its native aqueous environment. The situation has, however, changed for the better due to the advent of photoelectron spectroscopy in liquid microjets and its most recent application to photoionization of aqueous nucleosides, nucleotides, and larger DNA fragments. Here, we present a consistent and efficient computational methodology, which allows to accurately evaluate ionization energies and model photoelectron spectra of aqueous DNA and its individual components. After careful benchmarking, the method based on density functional theory and its time-dependent variant with properly chosen hybrid functionals and polarizable continuum solvent model provides ionization energies with accuracy of 0.2-0.3 eV, allowing for faithful modeling and interpretation of DNA photoionization. The key finding is that the aqueous medium is remarkably efficient in screening the interactions within DNA such that, unlike in the gas phase, ionization of a base, nucleoside, or nucleotide depends only very weakly on the particular DNA context. An exception is the electronic interaction between neighboring bases which can lead to sequence-specific effects, such as a partial delocalization of the cationic hole upon ionization enabled by presence of adjacent bases of the same type.

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.

Biorecognition by DNA oligonucleotides after exposure to photoresists and resist removers

Combining biological molecules with integrated circuit technology is of considerable interest for next generation sensors and biomedical devices. Current lithographic microfabrication methods, however, were developed for compatibility with silicon technology rather than bioorganic molecules, and consequently it cannot be assumed that biomolecules will remain attached and intact during on-chip processing. Here, we evaluate the effects of three common photoresists (Microposit S1800 series, PMGI SF6, and Megaposit SPR 3012) and two photoresist removers (acetone and 1165 remover) on the ability of surface-immobilized DNA oligonucleotides to selectively recognize their reverse-complementary sequence. Two common DNA immobilization methods were compared: adsorption of 5'-thiolated sequences directly to gold nanowires and covalent attachment of 5'-thiolated sequences to surface amines on silica coated nanowires. We found that acetone had deleterious effects on selective hybridization as compared to 1165 remover, presumably due to incomplete resist removal. Use of the PMGI photoresist, which involves a high temperature bake step, was detrimental to the later performance of nanowire-bound DNA in hybridization assays, especially for DNA attached via thiol adsorption. The other three photoresists did not substantially degrade DNA binding capacity or selectivity for complementary DNA sequences. To determine whether the lithographic steps caused more subtle damage, we also tested oligonucleotides containing a single base mismatch. Finally, a two-step photolithographic process was developed and used in combination with dielectrophoretic nanowire assembly to produce an array of doubly contacted, electrically isolated individual nanowire components on a chip. Postfabrication fluorescence imaging indicated that nanowire-bound DNA was present and able to selectively bind complementary strands.

Base-specific ionization of deprotonated nucleotides by resonance enhanced two-photon detachment

The intrinsic ionization energy of a base in DNA plays a critical role in determining the energies at which damage mechanisms may emerge. Here, a two-photon resonance-enhanced ionization scheme is presented that utilizes the (1)ππ* transition, localized on the DNA base, to elucidate the base-specific ionization in a deprotonated nucleotide. In contrast to previous reports, the scheme is insensitive to competing ionization channels arising from the sugar-phosphate backbone. Using this approach, we demonstrate that for all bases except guanine, the lowest electron detachment energy corresponds to detachment from the sugar-phosphate backbone and allows us to determine the lowest adiabatic ionization energy for the other three bases for the first time in an isolated nucleotide.

Sequence and conformation effects on ionization potential and charge distribution of homo-nucleobase stacks using M06-2X hybrid density functional theory calculations

DNA is subject to oxidative damage due to radiation or by-products of cellular metabolism, thereby creating electron holes that migrate along the DNA stacks. A systematic computational analysis of the dependence of the electronic properties of nucleobase stacks on sequence and conformation was performed here, on the basis of single- and double-stranded homo-nucleobase stacks of 1-10 bases or 1-8 base pairs in standard A-, B-, and Z-conformation. First, several levels of theory were tested for calculating the vertical ionization potentials of individual nucleobases; the M06-2X/6-31G* hybrid density functional theory method was selected by comparison with experimental data. Next, the vertical ionization potential, and the Mulliken charge and spin density distributions were calculated and considered on all nucleobase stacks. We found that (1) the ionization potential decreases with the number of bases, the lowest being reached by Gua≡Cyt tracts; (2) the association of two single strands into a double-stranded tract lowers the ionization potential significantly (3) differences in ionization potential due to sequence variation are roughly three times larger than those due to conformational modifications. The charge and spin density distributions were found (1) to be located toward the 5'-end for single-stranded Gua-stacks and toward the 3'-end for Cyt-stacks and basically delocalized over all bases for Ade- and Thy-stacks; (2) the association into double-stranded tracts empties the Cyt- and Thy-strands of most of the charge and all the spin density and concentrates them on the Gua- and Ade-strands. The possible biological implications of these results for transcription are discussed.

Interactions of amino acids with oxidized guanine in the gas phase associated with the protection of damaged DNA

Density functional theory calculations were employed to study the stabilization process of the guanine radical cation through amino acid interactions as well as to understand the protection mechanisms. On the basis of our calculations, several protection mechanisms are proposed in this work subject to the type of the amino acid. Our results indicate that a series of three-electron bonds can be formed between the amino acids and the guanine radical cation which may serve as relay stations supporting hole transport. In the three-electron-bonded, π-π-stacked, and H-bonded modes, amino acids can protect guanine from oxidation or radiation damage by sharing the hole, while amino acids with reducing properties can repair the guanine radical cation through proton-coupled electron transfer or electron transfer. Another important finding is that positively charged amino acids (ArgH(+), LysH(+), and HisH(+)) can inhibit ionization of guanine through raising its ionization potential. In this situation, a negative dissociation energy for hydrogen bonds in the hole-trapped and positively charged amino acid-Guanine dimer is observed, which explains the low hole-trapping efficiency. We hope that this work provides valuable information on how to protect DNA from oxidation- or radiation-induced damages in biological systems.

Theoretical Determination of One-Electron Oxidation Potentials for Nucleic Acid Bases

The oxidation potentials for N-methyl substituted nucleic acid bases guanine, adenine, cytosine, thymine, uracil, xanthine, and 8-oxoguanine were computed using B3LYP and CBS-QB3 with the SMD solvation model. Acid-base and tautomeric equilibria present in aqueous solution were accounted for by combining standard redox potentials with calculated pKa and tautomerization energies to produce an ensemble averaged pH dependent potential. Gas phase free energies were computed using B3LYP/aug-cc-pVTZ//B3LYP/6-31+G(d,p) and CBS-QB3. Solvation free energies were computed at the SMD/B3LYP/6-31+G(d,p) level of theory. Compared to experimental results, calculations with the CBS-QB3 level of theory have a mean absolute error (MAE) of ca. 1 kcal/mol for the gas phase proton affinity/gas phase basicity and an MAE of ca. 0.04 eV for the adiabatic/vertical ionization potentials. The B3LYP calculations have a MAE of ∼2 kcal/mol for the proton affinity/gas phase basicity data but systematically underestimated ionization potentials by 0.14-0.21 eV. Solvent cavities for charged solute species were rescaled uniformly by fitting computed pKa data to experimentally measured pKa values. After solvent cavity scaling, the MAEs for computed pKa's compared to experimental results are 0.7 for B3LYP and 0.9 for CBS-QB3. In acetonitrile, the computed E°(XH(+•)/XH) redox potentials are systematically lower than experimentally measured potentials by 0.21 V for CBS-QB3 and 0.33 V for B3LYP. However, the redox potentials relative to adenine are in very good agreement with experimental results, with MAEs of 0.10 V for CBS-QB3 and 0.07 V for B3LYP. In aqueous solution, B3LYP and CBS-QB3 have MAEs of 0.21 and 0.19 V for E7(X(•),H(+)/XH). Replacing the methyl substituent with ribose changes the calculated E7 potentials by 0.1-0.2 V. The calculated difference between the guanine and adenine oxidation potentials is too large compared to experimental results, but the calculated difference between guanine and 8-oxoguanine is in good agreement with the measured values.

One-electron oxidation of neutral sugar radicals of 2'-deoxyguanosine and 2'-deoxythymidine: a density functional theory (DFT) study

One electron oxidation of neutral sugar radicals has recently been suggested to lead to important intermediates in the DNA damage process culminating in DNA strand breaks. In this work, we investigate sugar radicals in a DNA model system to understand the energetics of sugar radical formation and oxidation. The geometries of neutral sugar radicals C(1')(•), C(2')(•), C(3')(•), C(4')(•), and C(5')(•) of 2'-deoxyguanosine (dG) and 2'-deoxythymidine (dT) were optimized in the gas phase and in solution using the B3LYP and ωB97x functionals and 6-31++G(d) basis set. Their corresponding cations (C(1')(+), C(2')(+), C(3')(+), C(4')(+), and C(5')(+)) were generated by removing an electron (one-electron oxidation) from the neutral sugar radicals, and their geometries were also optimized using the same methods and basis set. The calculation predicts the relative stabilities of the neutral sugar radicals in the order C(1')(•) > C(4')(•) > C(5')(•) > C(3')(•) > C(2')(•), respectively. Of the neutral sugar radicals, C(1')(•) has the lowest vertical ionization potential (IP(vert)), ca. 6.33 eV in the gas phase and 4.71 eV in solution. C(2')(•) has the highest IP(vert), ca. 8.02 eV, in the gas phase, and the resultant C(2') cation is predicted to undergo a barrierless hydride transfer from the C(1') site to produce the C(1') cation. One electron oxidation of C(2')(•) in dG is predicted to result in a low lying triplet state consisting of G(+•) and C(2')(•). The 5',8-cyclo-2'-deoxyguanosin-7-yl radical formed by intramolecular bonding between C(5')(•) and C(8) of guanine transfers spin density from C(5') site to guanine, and this structure has IP(vert) 6.25 and 5.48 eV in the gas phase and in solution.

THEORETICAL STUDIES ON PHOTOIONIZATION OF GUANINE TAUTOMERS AND INTERCONVERSION OF CATION RADICALS

The G3MP2B3 and P3 methods have been used to calculate the adiabatic and vertical ionization potentials (IPs) of the eight most stable tautomers of guanine. The calculated energy discrepancy between adiabatic and vertical IPs are in good agreement with the changes in geometry from neutral ground state to stable cation radicals. The geometries of imino-oxo form tautomers have no obvious change in the ionization process, which results in less energy discrepancy between vertical and adiabatic IPs. In the ionization process, the geometries of the amino-oxo and amino-hydroxy form tautomers change from nonplanar to planar structures. Hence the amino-oxo and amino-hydroxy form tautomers have larger energy discrepancy between vertical and adiabatic IPs. Further studies on the interconversion of the cation radicals shed further light on the transition process between the cation radicals and the main pathways are the hydrogen migrations and internal rotations of hydroxy (OH) and imino (NH) groups. The barriers of hydrogen rotations are lower than those of hydrogen migrations. Furthermore, the barriers for the hydrogen migrations between two rings are higher, which are about 3.0 eV.

Effect of solvation on the vertical ionization energy of thymine: from microhydration to bulk

The effect of hydration on the vertical ionization energy (VIE) of thymine was characterized using equation-of-motion ionization potential coupled-cluster (EOM-IP-CCSD) and effective fragment potential (EFP) methods. We considered several microsolvated clusters as well as thymine solvated in bulk water. The VIE in bulk water was computed by averaging over solvent-solute configurations obtained from equilibrium molecular dynamics trajectories at 300 K. The effect of microsolvation was analyzed and contrasted against the combined effect of the first solvation shell in bulk water. Microsolvation reduces the ionization energy (IE) by about 0.1 eV per water molecule, while the first solvation shell increases the IE by 0.1 eV. The subsequent solvation lowers the IE, and the bulk value of the solvent-induced shift of thymine's VIE is approximately -0.9 eV. The combined effect of the first solvation shell was explained in terms of specific solute-solvent interactions, which were investigated using model structures. The convergence of IE to the bulk value requires the hydration sphere of approximately 13.5 Å radius. The performance of the EOM-IP-CCSD/EFP scheme was benchmarked against full EOM-IP-CCSD using microhydrated structures. The errors were found to be less than 0.01-0.02 eV. The relative importance of the polarization and higher multipole moments in EFP model was also investigated.

Cooperative effect in the electronic properties of human telomere sequence

The contribution of sequence elements of human telomere DNA to the interaction of DNA with electrons has been analyzed. By applying wavelength dependent low-energy photoelectron transmission and two-photon photoemission spectroscopy, we investigated the density of states of DNA oligomers with partial sequence elements of the human telomere assembled as monolayers on gold. The findings demonstrate the role of the resonance states in the DNA in accepting electrons and the effect of the sequence on these states. When guanine (G) bases are clustered together, the resonance negative ion state is stabilized, as compared to oligomers containing the same number of G bases but distributed within the sequence. The electron-capturing probability of the human telomere-like oligomer, a sequence with an additional single adenine (A) base adjacent to the G cluster, is dramatically enhanced compared to the other oligomers studied, most likely due to the enhancement of the density of states near the highest occupied molecular orbital.