Search

We address the issue of whether chemical alterations of nucleobases are an effective tool to modulate charge transfer through DNA molecules. Our investigation uses a multilevel computational approach based on classical molecular dynamics and quantum chemistry. We find yet another piece of evidence that structural fluctuations are a key factor to determine the electronic structure of double-stranded DNA. We argue that the electronic structure and charge transfer ability of flexible polymers is the result of a complex intertwining of various structural, dynamical and chemical factors. Chemical intuition may be used to design molecular wires, but this is not the sole component in the complex charge transfer mechanism through DNA.

Among the distinct strategies proposed to expand the genetic alphabet, size-expanded nucleobases are promising for the development of modified DNA duplexes with improved biotechnological properties. In particular, duplexes built up by replacing canonical bases with the corresponding benzo-fused counterparts could be valuable as molecular nanowires. In this context, this study reports the results of classical molecular dynamics simulations carried out to examine the structural and dynamical features of size-expanded DNAs, including both hybrid duplexes containing mixed pairs of natural and benzo-fused bases (xDNA) and pure size-expanded (xxDNA) duplexes. Furthermore, the electronic structure of both natural and size-expanded duplexes is examined by means of density functional computations. The results confirm that the structural and flexibility properties of the canonical DNA are globally little affected by the presence of benzo-fused bases. The most relevant differences are found in the enhanced size of the grooves, and the reduction in the twist. However, the analysis also reveals subtle structural effects related to the nature and sequence of benzo-fused bases in the duplex. On the other hand, electronic structure calculations performed for xxDNAs confirm the reduction in the HOMO-LUMO gap predicted from the analysis of the natural bases and their size-expanded counterparts, which suggests that pure size-expanded DNAs can be good conductors. A more complex situation is found for xDNAs, where fluctuations in the electrostatic interaction between base pairs exerts a decisive influence on the modulation of the energy gap.

Quantum chemical calculations reveal that neither steric hindrance nor less polar solvent medium is able to reduce the otherwise high hydrolysis rate of arsenate-esters. These results question the stability of As-DNA not only in aqueous but also in non-aqueous environments.

Naphtho-homologated DNA bases have been recently used to build a new type of size-expanded DNA known as yyDNA. We have used theoretical techniques to investigate the structure, tautomeric preferences, base-pairing ability, stacking interactions, and HOMO-LUMO gaps of the naphtho-bases. The structure of these bases is found to be similar to that of the benzo-fused predecessors (y-bases) with respect to the planarity of the aromatic rings and amino groups. Tautomeric studies reveal that the canonical-like forms of naphtho-thymine (yyT) and naphtho-adenine (yyA) are the most stable tautomers, leading to hydrogen-bonded dimers with the corresponding natural nucleobases that mimic the Watson-Crick pairing. However, the canonical-like species of naphtho-guanine (yyG) and naphtho-cytosine (yyC) are not the most stable tautomers, and the most favorable hydrogen-bonded dimers involve wobble-like pairings. The expanded size of the naphtho-bases leads to stacking interactions notably larger than those found for the natural bases, and they should presumably play a dominant contribution in modulating the structure of yyDNA duplexes. Finally, the HOMO-LUMO gap of the naphtho-bases is smaller than that of their benzo-base counterparts, indicating that size-expansion of DNA bases is an efficient way of reducing their HOMO-LUMO gap. These results are examined in light of the available experimental evidence reported for yyT and yyC.

xDNA is an artificial duplex made of natural and benzo-homologated bases. The latter can be seen as a fusion between benzene and a natural base. We have used two different ab initio techniques, one based on B3LYP and a Gaussian expansion of the wave functions, and the other on GGA and plane-waves, to investigate the electronic properties of an xDNA duplex and a natural one with an analogous sequence. The calculations were performed in dry conditions, i.e., H atoms were used to neutralize the charge. It is found that the HOMO-LUMO gap of xDNA is about 0.5 eV smaller than that of B-DNA, independent of the technique used. The pi-pi* gap of xDNA is 1.3 or 1.0 eV smaller than that of B-DNA, depending on whether one uses B3LYP/6-31G or GGA/plane-waves, respectively. An analysis of how saturation changes the electronic properties of the nucleotide pairs that make up these duplexes suggests that different saturation schemes significantly affect the HOMO-LUMO gap value of xDNA and B-DNA. The same is not true for the pi-pi* gap. That xDNA has a smaller pi-pi* gap than B-DNA suggests that xDNA could be a plausible candidate for molecular-wire applications.

M-DNA is a type of metalated DNA that forms at high pH and in the presence of Zn, Ni, and Co, with the metals placed in between each base pair, as in G-Zn-C. Experiments have found that M-DNA could be a promising candidate for a variety of nanotechnological applications, as it is speculated that the metal d-states enhance the conductivity, but controversy still clouds these findings. In this paper, we carry out a comprehensive ab initio study of eight G-Zn-C models in the gas phase to help discern the structure and electronic properties of Zn-DNA. Specifically, we study whether a model prefers to be planar and has electronic properties that correlate with Zn-DNA having a metallic-like conductivity. Out of all the studied models, there is only one which preserves its planarity upon full geometry optimization. Nevertheless, starting from this model, one can deduce a parallel Zn-DNA architecture only. This duplex would contain the imino proton, in contrast to what has been proposed experimentally. Among the nonplanar models, there is one that requires less than 8 kcal/mol to flatten (both in gas and solvent conditions), and we propose that it is a plausible model for building an antiparallel duplex. In this duplex, the imino proton would be replaced by Zn, in accordance with experimental models. Neither planar nor nonplanar models have electronic properties that correlate with Zn-DNA having a metallic-like conductivity due to Zn d-states. To understand whether density functional theory (DFT) can describe appropriately the electronic properties of M-DNAs, we have investigated the electronic properties of G-Co-C base pairs. We have found that when self-interaction corrections (SIC) are not included the HOMO state contains Co d-levels, whereas these levels are moved below the HOMO state when SIC are considered. This result indicates that caution should be exercised when studying the electronic properties of M-DNAs with functionals that do not account for strong electronic correlations.

Characterization of the electrical properties of the DNA bases (adenine, cytosine, guanine, and thymine), in addition to building the basic knowledge on these fundamental constituents of a DNA, is a crucial step in developing a DNA sequencing technology. We present a first-principles study of the current-voltage characteristics of nucleotidelike molecules of the DNA bases, placed in a 1.5 nm gap formed between gold nanoelectrodes. The quantum transport calculations in the tunneling regime are shown to vary strongly with the electrode-molecule geometry and the choice of the density-functional theory exchange-correlation functionals. Analysis of the results in the zero-bias limit indicates that distinguishable current-voltage characteristics of different DNA bases are dominated by the geometrical conformations of the bases and nanoelectrodes.

First-principles calculation of the transverse conductance across DNA fragments placed between gold nanoelectrodes reveals that such conductance describes electron tunneling that depends critically on geometrical rather than electronic-structure properties. By factoring the first-principles result into two simple and approximately independent tunneling factors, we show that the conductances of the A, C, G, and T fragments differ only because of their sizes: the larger is the DNA base, the smaller its distance to the electrode, and the larger its conductance. Because the geometrical factors are difficult to control in an experiment, the direct-current measurements across DNA with gold contact electrodes may not be a convenient approach to DNA sequencing.

The size-expanded DNA bases, xA, xC, xG, and xT, are benzo-homologue forms of the natural DNA bases; i.e., their structure can be seen as the fusion of a natural base and a benzene ring. Recently, a variety of DNAs, known as xDNAs, have been synthesized in which size-expanded and natural bases are paired. In this paper we use second-order Møller-Plesset perturbation theory and density functional theory to investigate the structural and electronic properties of xA, xC, xG, and xT and their natural counterparts. We find that whereas natural and size-expanded bases have both nonplanar amino groups the latter have also nonplanar aromatic rings. When density functional theory is used to investigate the electronic properties of size-expanded and natural bases, it is found that the HOMO-LUMO gap of the size-expanded bases is smaller than that of the natural bases. Also, xG should be easier to oxidize than G.

Molecular dynamics simulations of the DNA duplex d(CCAACGTTGG)2 were used to study the relationship between DNA sequence and structure in a crystal environment. Three different force fields were used: a traditional description based on atomic point charges, a polarizable force field, and an “extra-point” force field (with additional charges on extranuclear sites). It is found that all the force fields reproduce fairly well the sequence-dependent features of the experimental structure. The polarizable force field, however, provides the most accurate representation of the crystal structure and the sequence-dependent effects observed in the experiment. These results point out to the need of the inclusion of polarization for accurate descriptions of DNA. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Size-expanded DNA bases are analogues of natural bases that can be described as a synthesis between benzene and a natural base. Size-expanded bases have been combined with natural bases to form xDNA and yDNA, a new class of synthetic nucleic acids. We are interested in xDNA and yDNA because they might function as molecular wires. Recently, we also became intrigued by the possibility of altering the electronic conductivity of xDNA and yDNA by means of structural changes in the constituent bases. This possibility appeared after we noticed that the highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO–LUMO) gap of the base yG can be increased dramatically, ∼0.73 eV, by changing the aromaticity of its benzene ring. Therefore, if one is able to alter the HOMO–LUMO gap of size-expanded bases, it should be possible to change the electronic conductivity of xDNAs and yDNAs as well. In the present work, we extend our study on aromaticity-induced changes on the electronic properties of size-expanded bases by investigating the HOMO–LUMO gap of all possible tautomers of xC. We have found that, as for yG, the HOMO–LUMO gap of xC can be modified by ∼ 0.74 eV, and that this can be accomplished by changing the aromaticity of its benzene ring. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006 [ABSTRACT FROM AUTHOR]