Energy contribution of the solvent to the charge migration in DNA

The authors have investigated the interactions of the reaction centers, participating in the charge transfer reaction within the DNA molecule with the phosphate backbones and the solvent molecules, and have estimated the contribution of these interactions into the charge migration in DNA. They have determined the unequal shift of the energy surfaces of the initial and final transition states of the transfer reaction along the energy axis and the dependence of the magnitude of the energy shift on the nature of the reaction centers and the surrounding environment. The nonuniform distribution of the negative charge in the DNA phosphate backbones results in an increase of the positive shift of the energy surface of the DNA base pairs in the center of the structure, where the maximum density of the negative charge is concentrated. Localization of the positive charge on the guanine and the adenine in the DNA base pairs in the oxidized state results in a dependence of the free energy of reaction in the solvent on the pair sequences and their arrangement in the DNA chain. As an example, for the G-C/A-T configuration the positive charges are localized on the same strand that results in a decrease of the free energy of reaction in the solvent for charge migration from G-C to A-T pair by 0.125 eV.     

Intercalation of trioxatriangulenium ion in DNA: binding,electron transfer,x-ray crystallography,and electronic structure

Trioxatriangulenium ion (TOTA(+)) is a flat, somewhat hydrophobic compound that has a low-energy unoccupied molecular orbital. It binds to duplex DNA by intercalation with a preference for G-C base pairs. Irradiation of intercalated TOTA(+) causes charge (radical cation) injection that results in strand cleavage (after piperidine treatment) primarily at GG steps. The X-ray crystal structure of TOTA(+) intercalated in the hexameric duplex d[CGATCG](2) described here reveals that intercalation of TOTA(+) results in an unusually large extension of the helical rise of the DNA and that the orientation of TOTA(+) is sensitive to hydrogen-bonding interactions with backbone atoms of the DNA. Electronic structure calculations reveal no meaningful charge transfer from DNA to TOTA(+) because the lowest unoccupied molecular orbital of TOTA(+), (LUMO)(T), falls in the gap between the highest occupied molecular orbital, (HOMO)(D), and the (LUMO)(D) of the DNA bases. These calculations reveal the importance of backbone, water, and counterion interactions, which shift the energy levels of the bases and the intercalated TOTA(+) orbitals significantly. The calculations also show that the inserted TOTA(+) strongly polarizes the intercalation cavity where a sheet of excess electron density surrounds the TOTA(+).     

Structural thermodynamics of lamellar cationic lipid-DNA complex: DNA compressibility modulus

We have studied theoretically the compressibility modulus B of DNA and complexation adsorption isotherms of DNA and lipids, as a function of DNA spacing d(DNA) and NaCl electrolyte concentration, respectively, in isoelectric states of lamellar DNA/cationic lipid (CL) self-assemblies. The electrostatic free energy derived from the Poisson-Boltzmann theory predicts partial agreement with measured B values for interhelical separations d(DNA)>33 A when use is made of a fit of hydration repulsion from bulk DNA hexagonal phases in solution. For lower interchain separations the prediction worsens due to the hydration interaction that overcomes the electrostatic contribution. An exact match of the system's counterion electrochemical potentials and the coions of salt in aqueous phase leads to the electrostatic part of the free energy that renders isotherms of d(DNA) versus ionic strength in qualitative consistency with general trends of available experimental data of CL-DNA complexes.     

Developing a computational model that accurately reproduces the structural features of a dinucleoside monophosphate unit within B-DNA

The ability of a dinucleoside monophosphate to mimic the conformation of B-DNA was tested using a combination of different phosphate models (anionic, neutral, counterion), environments (gas, water), and density functionals (B3LYP, MPWB1K, M06-2X) with the 6-31G(d,p) basis set. Three sequences (5'-GX(Py)-3', where X(Py) = T, U or (Br)U) were considered, which vary in the (natural or modified) 3' pyrimidine nucleobase (X(Py)). These bases were selected due to their presence in natural DNA, structural similarity to T and/or applications in radical-initiated anti-tumour therapies. The accuracy of each of the 54 model, method and sequence combinations was judged based on the ability to reproduce key structural features of natural B-DNA, including the stacked base-base orientation and important backbone torsion angles. B3LYP yields distorted or tilted relative base-base orientations, while many computational challenges were encountered for MPWB1K. Despite wide use in computational studies of DNA, the neutral (protonated) phosphate model could not consistently predict a stacked arrangement for all sequences. In contrast, stacked base-base arrangements were obtained for all sequences with M06-2X in conjunction with both the anionic and (sodium) counterion phosphate models. However, comparison of calculated and experimental backbone conformations reveals the charge-neutralized counterion phosphate model best mimics B-DNA. Structures optimized with implicit solvent (water) are comparable to gas-phase structures, which suggests similar results should be obtained in an environment of intermediate polarity. We recommend the use of either gas-phase or water M06-2X optimizations with the counterion phosphate model to study the structure and/or reactivity of natural or modified DNA with a dinucleoside monophosphate.     

Lifetime regulation of the charge-separated state in DNA by modulating the oxidation potential of guanine in DNA through hydrogen bonding

A series of naphthalimide (NI)- and 5-bromocytosine ((br)C)-modified oligodeoxynucleotides (ODNs) were prepared, and their lifetimes of the charge-separated states during the photosensitized one-electron oxidation of DNA were measured. Various lifetimes of the charge-separated states were observed depending on the sequence and the incorporation sites of (br)C, and the oxidation potential of G in the (br)C:G base-pair relative to that of G in the C:G base-pair and in the GGG sequence was determined by comparing the lifetimes of the charge-separated states. The change in the cytosine C5 hydrogen to bromine resulted in a 24 mV increase in the oxidation potential of G in the (br)C:G base-pair as compared to that of G in the C:G base-pair, the value of which is comparable to a 58 mV decrease in the oxidation potential of G in the GGG sequence. These results clearly demonstrate that hole transfer in DNA can be controlled through hydrogen bonding by introducing a substituent on the cytosine.     

B—DNA分子中互补碱基对G—C、A—T水合作用的研究

本文通过对B—DNA分子中互补碱基对A—T、G—C水合作用的量子化学计算以及最优化处理求得了它们水合作用的最优配位模式。结果表明,当水分子与互补碱基对A—T、G—C共面且水分子与A—T碱基对中胸腺嘧啶(T)杂环上的O_2原子形成氢键以及与G—C碱基对中胞嘧啶(C)杂环上的O_2原子形成氢键时,它们的水合作用具有最稳定的配位模式。     

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

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

Oxidative damage to DNA: counterion-assisted addition of water to ionized DNA

Oxidative damage to DNA, implicated in mutagenesis, aging, and cancer, follows electron loss that generates a radical cation that migrates to a guanine, where it may react with water to form 8-oxo-7,8-dihydroguanine (8-OxoG). Molecular dynamics and ab initio quantum simulations on a B-DNA tetradecamer reveal activated reaction pathways that depend on the local counterion arrangement. The lowest activation barrier, 0.73 eV, is found for a reaction that starts from a configuration where a Na(+) resides in the major groove near the N7 atoms of adjacent guanines, and evolves through a transition state where a bond between a water oxygen atom and a carbon atom forms concurrently with displacement of a proton toward a neighboring water molecule. Subsequently, a bonded complex of a hydronium ion and the nearest backbone phosphate group forms. This counterion-assisted proton shuttle mechanism is supported by experiments exploiting selective substitution of backbone phosphates by methylphosphonates.     

Water induced weakly bound electrons in DNA

We have studied the effect of humidity on the electronic properties of DNA base pairs. We found that the hydrogen links of the nucleobases with water molecules lead to a shift of the pi electron density from carbon atoms to nitrogen atoms and can change the symmetry of the wave function for some nucleobases. As a result, the orbital energies are shifted which leads to a decrease in the potential barrier for the hole transfer between the G-C and A-T pairs from 0.7 eV for the dehydrated case to 0.123 eV for the hydrated. More importantly, the pi electron density redistribution activated by hydration is enhanced by the intrastrand interactions. This leads to a modification of the nucleobase chemical structures from the covalent type to a resonance structure with separated charges, where some pi electrons are not locked up into the covalent bonds. Within the (G-C)(2) sequences, there is overlapping of the electronic clouds of such unlocked electrons belonging to the stacked guanines, that significantly increases the electron coupling between them to V(DA)=0.095 eV against the V(DA)=0.025 eV for the dehydrated case. Consequently, the charge transfer between two guanines within the (G-C)(2) sequences is increased by 250 times due to hydration. The presence of nonbonded electrons suppress the band gap up to approximately 3.0 eV, that allows us to consider DNA as a narrow band gap semiconductor.     

A solid-state 23Na NMR study of monovalent cation binding to double-stranded DNA at low relative humidity

We report a solid-state (23)Na NMR study of monovalent cation (Li(+), Na(+), K(+), Rb(+), Cs(+) and NH(4) (+)) binding to double-stranded calf thymus DNA (CT DNA) at low relative humidity, ca 0-10%. Results from (23)Na--(31)P rotational echo double resonance (REDOR) NMR experiments firmly establish that, at low relative humidity, monovalent cations are directly bound to the phosphate group of CT DNA and are partially dehydrated. On the basis of solid-state (23)Na NMR titration experiments, we obtain quantitative thermodynamic parameters concerning the cation-binding affinity for the phosphate group of CT DNA. The free energy difference (DeltaG degrees ) between M(+) and Na(+) ions is as follows: Li(+) (-1.0 kcal mol(-1)), K(+) (7.2 kcal mol(-1)), NH(4) (+) (1.0 kcal mol(-1)), Rb(+) (4.5 kcal mol(-1)) and Cs(+) (1.5 kcal mol(-1)). These results suggest that, at low relative humidity, the binding affinity of monovalent cations for the phosphate group of CT DNA follows the order: Li(+) > Na(+) > NH(4) (+) > Cs(+) > Rb(+) > K(+). This sequence is drastically different from that observed for CT DNA in solution. This discrepancy is attributed to the different modes of cation binding in dry and wet states of DNA. In the wet state of DNA, cations are fully hydrated. Our results suggest that the free energy balance between direct cation-phosphate contact and dehydration interactions is important. The reported experimental results on relative ion-binding affinity for the DNA backbone may be used for testing theoretical treatment of cation-phosphate interactions in DNA.     

Laser probes of conformational isomerization in flexible molecules and complexes

Molecules with several flexible coordinates have potential energy surfaces with a large number of minima and many transition states separating them. A general experimental protocol is described that is capable of studying conformational isomerization in such circumstances, measuring the product quantum yields following conformation-specific infrared excitation, and measuring energy thresholds for isomerization of specific X --> Y reactant-product isomer pairs following excitation via stimulated emission pumping (SEP). These methods have been applied to a series of molecules of varying size and conformational complexity, including 3-indolepropionic acid (IPA), meta-ethynylstyrene, N-acetyltryptophan methyl amide (NATMA), N-acetyltryptophan amide (NATA), and melatonin. Studies of isomerization in solute-solvent complexes are also described, including a measurement of the barrier to isomerization in the IPA-H2O complex, and a unique isomerization reaction in which a single water molecule is shuttled between H-bonding sites on the trans-formanilide (TFA) molecule.     

Free energy landscape of A-DNA to B-DNA conversion in aqueous solution

The interconversion between the well-characterized A- and B-forms of DNA is a structural transition for which the intermediate states and the free energy difference between the two endpoints are not known precisely. In the present study, the difference between the Root Mean Square Distance (RMSD) from canonical A-form and B-form DNA is used as an order parameter to characterize this free energy difference using umbrella sampling molecular dynamics (MD) simulations with explicit solvent. The constraint imposed along this order parameter allows relatively unrestricted evolution of the intermediate structures away from both canonical A- and B-forms. The free energy difference between the A- and B-forms for the hexamer DNA sequence CTCGAG in aqueous solution is conservatively estimated to be at least 2.8 kcal/mol. A continuum of intermediate structures with no well-defined local minima links the two forms. The absence of any major barriers in the free energy surface is consistent with spontaneous conversion of the A-form DNA to B-form DNA in unconstrained simulations. The extensive sampling in the MD simulations (>0.1 mus) also allowed quantitative energetic characterization of local backbone conformational variables such as sugar pseudorotation angles and BI/BII state equilibria and their dependence on base identity. The absolute minimum in the calculated free energy profile corresponds closely to the crystal structure of the hexamer sequence, indicating that the present method has the potential to identify the most stable state for an arbitrary DNA sequence in water.     

Evaluation of effect of functionalized gold nanoparticles with a partial negative charge on stability of DNA molecule: a study of molecular dynamics simulation

Today, understanding the interaction between DNA molecule with nanoparticles and functionalized nanoparticles has a significant importance in medical applications and targeted drug delivery. Molecular dynamics simulation on double-stranded molecule with the structure of the double helix and sequence of (CCTCAGGCCTCC) was performed in three states. The aim was to evaluate the effect of gold nanoparticles (GNPs) with partial negative charge on the stability of a DNA molecule. During the simulation process, the GNPs become closed to the DNA molecule and phosphate groups of the DNA molecule guided the nanoparticles toward its major groove. At the end of the DNA molecule chain, the terminal nucleotide of the chain was laid flat on the surface of the GNPs due to excessive exposure to solvent molecules and occurrence of peeling and untwisting states. According to the results, proximity of the GNPs and functionalized GNPs to the DNA molecule led to increased configuration entropy. While conformational energy and van der Waals energy of the DNA molecule increased in the presence of the GNPs and functionalized GNPs, there was a reduction and an increase in the number of hydrogen bonds between complementary bases in the presence of the GNPs and functionalized GNPs, respectively. Radial distribution function was estimated for water molecules and sodium cations, compared to oxygen atoms of the phosphate group of the DNA molecule. Results were indicative of the release of water molecules from around the DNA molecule in the presence of the GNPs. In addition, the distance between sodium cations and the GNPs decreased. Nevertheless, no such phenomenon occurred in the presence of the functionalized GNPs. Therefore, according to results, it seems that GNPs decreased the stability of the DNA molecule and the functionalized GNPs with partial negative charge caused structural changes and created compression, but did not destroy the double-strand structure of the DNA molecule.     

Effects of various halogen anions and cations of alkali metals on energetics of excess charge recombination in stilbene donor-acceptor capped DNA hairpins

DNA hairpin conjugates with a stilbenedicarboxamide (Sa) hole donor and a stilbenediether (Sd) hole acceptor are considered as model systems for studying charge recombination (CR) of excess charges in DNA. Using the method of thermodynamic integration, we estimated the relative free energies of this process in hairpins with three adenine:thymine pairs between Sa and Sd surrounded by 1 M aqueous solutions of ionic compounds M(+)Cl(-) (M = Li, Na, K) and Na(+)X(-) (X = F, Cl, Br, I). The values of this quantity were calculated with respect to the free energy for the same hairpin in the 1 M NaCl aqueous solution. Based on the results obtained, we conclude that halogen anions have no significant influence on the rate of the CR reaction. By contrast, cations of other alkali metals can considerably change the potential barrier of the process, thus affecting the reaction rate. Different results obtained for cations and anions were attributed to the fundamental distinction in the electrostatic interactions of M(+) and X(-) species with negatively charged phosphate groups of the hairpin. In addition, our results show that the relative free energy of CR is larger for cations that are able to be closer to Sd and Sa structural units. The latter correlation suggests that the replacement of Na(+) by cations of other alkali metals enables one to change the CR rate modifying it in either direction.     

Surfactant encapsulated DNA: structure characterization and interaction with dye molecules in organic media

The recognition of electrostatically-bound DNA-didodecyldimethylammonium (DNA-DDDA) complex by three dye molecules, acridine orange (AO), ethidium bromide (EB) and 5,10,15,20-tetrakis(4-N-methylpyridyl) porphyrin tetra(p-toluenesulfonate) (TMPyP) in organic media was investigated through 1H NMR, UV-vis, and circular dichroism (CD) spectroscopies. When the organic solvent in which DNA-DDDA complex dissolves is changed from ethanol to chloroform, the adsorbed AO undergoes a reversible transformation from a monomer to a highly aggregated state at the interface between DNA and DDDA. EB also adsorbs at the interface between DNA and DDDA when EB interacts with the DNA-DDDA complex in organic media, but its existing state is independent of the used solvents. The third dye, TMPyP cation can intercalate into the G-C region while its anionic p-tosylate counterion remains unbound when it mixes with DNA complex in organic media. The complexes of DDDA with previously recognized DNA by the three dye molecules (DNA-dye), respectively, are also investigated. AO seems having changed its location from the grooves of DNA to the interface between DNA and DDDA after DNA-AO complex was electrostatically encapsulated with DDDA. The aggregation behavior of AO also shows a dependence on the polarity of the organic solvent. EB molecules are believed to intercalate into the base pairs of DNA in aqueous solution. The intercalation mode is still maintained after the encapsulation for DNA-EB in organic solvents, which is different from the situation between DNA-DDDA complex and EB. But in both cases, the existing states of EB are independent of the polarity of the organic solvents. Finally, TMPyP in the complex of DNA-TMPyP and DDDA is also judged to intercalate into the G-C region of DNA while its anionic p-tosylate counterion remains separated from DNA complex, which is similar to its interaction with DNA-DDDA complex in organic media. These data also strongly suggest that the intercalation state of TMPyP is more stable than its adsorption state in grooves when it is recognized with DNA. The present results are significant for the designs of both laser dye and conductive materials.     

The electronic structure of free water clusters probed by Auger electron spectroscopy

(H2O)(N) clusters generated in a supersonic expansion source with N approximately 1000 were core ionized by synchrotron radiation, giving rise to core-level photoelectron and Auger electron spectra (AES), free from charging effects. The AES is interpreted as being intermediate between the molecular and solid water spectra showing broadened bands as well as a significant shoulder at high kinetic energy. Qualitative considerations as well as ab initio calculations explain this shoulder to be due to delocalized final states in which the two valence holes are mostly located at different water molecules. The ab initio calculations show that valence hole configurations with both valence holes at the core-ionized water molecule are admixed to these final states and give rise to their intensity in the AES. Density-functional investigations of model systems for the doubly ionized final states--the water dimer and a 20-molecule water cluster--were performed to analyze the localization of the two valence holes in the electronic ground states. Whereas these holes are preferentially located at the same water molecule in the dimer, they are delocalized in the cluster showing a preference of the holes for surface molecules. The calculated double-ionization potential of the cluster (22.1 eV) is in reasonable agreement with the low-energy limit of the delocalized hole shoulder in the AES.     

HOMO energy gap dependence of hole-transfer kinetics in DNA

DNA consists of two type of base-pairs, G-C and A-T, in which the highest occupied molecular orbital (HOMO) localizes on the purine bases G and A. While the hole transfer through consecutive Gs or As occurs faster than 10(9) s(-1), a significant drop in the hole transfer rate was observed for G-C and A-T mixed random sequences. In this study, by using various natural and artificial nucleobases having different HOMO levels, the effect of the HOMO-energy gap between bases (Δ(HOMO)) on the hole-transfer kinetics in DNA was investigated. The results demonstrated that the hole transfer rate can be increased by decreasing the Δ(HOMO) and can be finely tuned over 3 orders of magnitude by varying the Δ(HOMO).     

Emergent functionality of nucleobase radical cations in duplex DNA: prediction of reactivity using qualitative potential energy landscapes

The one-electron oxidation of a series of DNA oligonucleotides was examined. Each oligomer contains a covalently linked anthraquinone (AQ) group. Irradiation of the AQ group with near-UV light results in a one-electron oxidation of the DNA that generates a radical cation (electron "hole"). The radical cation migrates through the DNA by a hopping mechanism and is trapped by reaction with water or molecular oxygen, which results in chemical reaction at particular nucleobases. This reaction is revealed as strand cleavage when the irradiated oligonucleotide is treated with piperidine. The specific oligomers examined reveal the existence of three categories of nucleobase sequences: charge shuttles, charge traps, and barriers to charge migration. The characterization of a sequence is not independent of the identity of other sequences in the oligonucleotide, and for this reason, the function of a particular sequence emerges from an analysis of the entire structure. Qualitative potential energy landscapes are introduced as a tool to assist in the rationalization and prediction of the reactions of nucleobases in oxidized DNA.     

Influence of DNA structure on the reactivity of the guanine radical cation

Oxidative damage of DNA via radical cation formation is a common cause of mutagenesis, cancer and of the physiological changes associated with aging. By using state-of-the-art ab initio molecular dynamics simulations, we study the mechanism that guides the first steps of this process. In the mechanism proposed here, guanine, which among the bases has the lowest oxidation potential, and the phosphate backbone play a crucial role. We found that the rate limiting step is the water protolysis. We illuminate the role of the local environment in considerably lowering the barrier. Of particular relevance in this respect is the role of the phosphate backbone.     

Impact of the anion on electrochemically doped regioregular and regiorandom poly(3-hexylthiophene)

Chemical and electrochemical doping of π-conjugated polymers is an important aspect in determining the performance and enabling the operation of many organic electronic devices, from organic light emitting diodes and thermoelectrics to organic electrochemical transistors. In both chemical doping and electrochemical doping an ionized dopant or counterion is present along with the doped π-conjugated polymer. This dopant or counterion is not a benign spectator, rather, its presence can significantly impact the optical, electronic, and thermoelectric properties of the resulting material. Here, we investigate how counterion structure impacts the electrochemical doping ability, oxidation potential, ionization energy, and polaron absorbance of regioregular (rr) and regiorandom (rra) P3HT. We find that in most cases the anion has a small effect on the polymer oxidation potential, except for in the case of rr-P3HT with the large tetrakis[3,5-bis(trifluoromethyl)phenyl]borate anion. We propose that this large anion is excluded from the crystalline regions and thus the oxidation potential is similar to that of rra-P3HT. The anions also result in significant differences in polaron absorbance and ionization energies, thereby emphasizing the important role of the counterion in determining the optical and electronic properties of doped π-conjugated polymers.