Charge transport in DNA via thermally induced hopping.

In this contribution we advance and explore the thermally induced hopping (TIH) mechanism for long-range charge transport (CT) in DNA and in large-scale chemical systems. TIH occurs in donor-bridge-acceptor systems, which are characterized by off-resonance donor-bridge interactions (energy gap DeltaE > 0), involving thermally activated donor-bridge charge injection followed by intrabridge charge hopping. We observe a "transition" from superexchange to TIH with increasing the bridge length (i.e., the number N of the bridge constituents), which is manifested by crossing from the exponential N-dependent donor-acceptor CT rate at low N (< N(X)) to a weakly (algebraic) N-dependent CT rate at high N (>N(X)). The "critical" bridge size N(X) is determined by the energy gap, the nearest-neighbor electronic couplings, and the temperature. Experimental evidence for the TIH mechanism was inferred from our analysis of the chemical yields for the distal/proximal guanine (G) triplets in the (GGG)(+)TTXTT(GGG) duplex (X = G, azadine (zA), and adenine (A)) studied by Nakatani, Dohno and Saito [J. Am. Chem. Soc. 2000, 122, 5893]. The TIH sequential model, which involves hole hopping between (GGG) and X, is analyzed in terms of a sequential process in conjunction with parallel reactions of (GGG)(+) with water, and provides a scale of (free) energy gaps (relative to (GGG)(+)) of Delta = 0.21-0.24 eV for X = A, Delta = 0.10-0.14 eV for X = zA, and Delta = 0.05-0.10 eV for X = G. We further investigated the chemical yields for long-range TIH in (G)l(+)Xn(G)l (l = 1-3) duplexes, establishing the energetic constraints (i.e., the donor - bridge base (X) energy gap Delta), the bridge structural constraints (i.e., the intrabridge X-X hopping rates k(m)), and the kinetic constraints (i.e., the rate k(d) for the reaction of with water). Effective TIH is expected to prevail for Delta less than or approximately equal to 0.20 eV with a "fast" water reaction (k(d)/k(m) approximately 10(-3)) and for Delta < 0.30 eV with a "slow" water reaction (k(d)/k(m) approximately 10(-5)). We conclude that (T)n bridges (for which Delta approximately equals 0.6 eV) cannot act in TIH of holes. From an analysis based on the energetics of the electronic coupling matrix elements in G(+)(T-A)n(GGG) duplexes we conclude that the superexchange mechanism is expected to dominate for n = 1-4. For long (A)n bridges (n > or approximately equal to 4) the TIH prevails, provided that the water side reaction is slow, raising the issue of chemical control of TIH through long (A)n bridges in DNA attained by changing the solution composition.     

Dynamics and energetics of single-step hole transport in DNA hairpins

The dynamics of single-step hole transport processes have been investigated in a number of DNA conjugates possessing a stilbenedicarboxamide electron acceptor, a guanine primary donor, and several secondary donors. Rate constants for both forward and return hole transport between the primary and secondary donor are obtained from kinetic modeling of the nanosecond transient absorption decay profiles of the stilbene anion radical. The kinetic model requires that the hole be localized on either the primary or the secondary donor and not delocalized over both the primary and the secondary donor. Rate constants for hole transport are found to be dependent upon the identity of the secondary donor, the intervening bases, and the location of the secondary donor in the same strand as the primary donor or in the complementary strand. Rate constants for hole transport are much slower than those for the superexchange process used to inject the hole on the primary donor. This difference is attributed to the larger solvent reorganization energy for charge transport versus charge separation. The hole transport rate constants obtained in these experiments are consistent with experimental data for single-step hole transport from other transient absorption studies. Their relevance to long-distance hole migration over tens of base pairs remains to be determined. The forward and return hole transport rate constants provide equilibrium constants and free energies for hole transport equilibria. Secondary GG and GGG donors are found to form very shallow hole traps, whereas the nucleobase deazaguanine forms a relatively deep hole trap. This conclusion is in accord with selected strand cleavage data and thus appears to be representative of the behavior of holes in duplex DNA. Our results are discussed in the context of current theoretical models of hole transport in DNA.     

Absolute rates of hole transfer in DNA

Absolute rates of hole transfer between guanine nucleobases separated by one or two A:T base pairs in stilbenedicarboxamide-linked DNA hairpins were obtained by improved kinetic analysis of experimental data. The charge-transfer rates in four different DNA sequences were calculated using a density-functional-based tight-binding model and a semiclassical superexchange model. Site energies and charge-transfer integrals were calculated directly as the diagonal and off-diagonal matrix elements of the Kohn-Sham Hamiltonian, respectively, for all possible combinations of nucleobases. Taking into account the Coulomb interaction between the negative charge on the stilbenedicarboxamide linker and the hole on the DNA strand as well as effects of base pair twisting, the relative order of the experimental rates for hole transfer in different hairpins could be reproduced by tight-binding calculations. To reproduce quantitatively the absolute values of the measured rate constants, the effect of the reorganization energy was taken into account within the semiclassical superexchange model for charge transfer. The experimental rates could be reproduced with reorganization energies near 1 eV. The quantum chemical data obtained were used to discuss charge carrier mobility and hole-transport equilibria in DNA.     

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.     

Electronic couplings and on-site energies for hole transfer in DNA: systematic quantum mechanical/molecular dynamic study

The electron hole transfer (HT) properties of DNA are substantially affected by thermal fluctuations of the pi stack structure. Depending on the mutual position of neighboring nucleobases, electronic coupling V may change by several orders of magnitude. In the present paper, we report the results of systematic QM/molecular dynamic (MD) calculations of the electronic couplings and on-site energies for the hole transfer. Based on 15 ns MD trajectories for several DNA oligomers, we calculate the average coupling squares V(2) and the energies of basepair triplets XG(+)Y and XA(+)Y, where X, Y=G, A, T, and C. For each of the 32 systems, 15,000 conformations separated by 1 ps are considered. The three-state generalized Mulliken-Hush method is used to derive electronic couplings for HT between neighboring basepairs. The adiabatic energies and dipole moment matrix elements are computed within the INDO/S method. We compare the rms values of V with the couplings estimated for the idealized B-DNA structure and show that in several important cases the couplings calculated for the idealized B-DNA structure are considerably underestimated. The rms values for intrastrand couplings G-G, A-A, G-A, and A-G are found to be similar, approximately 0.07 eV, while the interstrand couplings are quite different. The energies of hole states G(+) and A(+) in the stack depend on the nature of the neighboring pairs. The XG(+)Y are by 0.5 eV more stable than XA(+)Y. The thermal fluctuations of the DNA structure facilitate the HT process from guanine to adenine. The tabulated couplings and on-site energies can be used as reference parameters in theoretical and computational studies of HT processes in DNA.     

Charge on the move: how electron-transfer dynamics depend on molecular conformation

This tutorial review illustrates the many facets whereby the molecular conformation helps to control the rates of through-bond electron transfer. A brief introduction to Marcus theory is given, highlighting the importance of the coupling element and the super-exchange mechanism, before considering the reasons why the coupling element might depend on the molecular geometry. The methods currently available for determination of both the coupling element and the geometry are reviewed and various examples are given for systems where the structure controls the degree of electronic coupling along the molecular axis. The role of the "bridge" connecting the donor and acceptor is emphasized.     

Systematic study of the influence of base-step parameters on the electronic coupling between base-pair dimers: comparison of A-DNA and B-DNA forms

The electronic coupling is one of the key parameters governing electron hole transfer along DNA helices. In this study, we established the first comprehensive data base of electronic coupling elements, calculated at the ab initio level. The data set comprises all possible Watson-Crick base pair dimers, both in standard A-DNA and B-DNA geometries. We also quantified the sensitivity of the coupling elements with respect to geometry changes by varying each of the six standard base step parameters, which specify the relative orientation of neighboring base pairs. We compare the couplings in a systematic way by discussing variations in the coupling magnitude due to geometry or nucleotide sequence in the dimer, and we analyze how the structure affects the electronic coupling in terms of general and dimer-specific trends. Furthermore, we studied how the coupling changes when one introduces the chemically modified base 7-deazaguanine in the corresponding base-pair dimers. Finally, on the basis of the calculated coupling elements, we suggest a model duplex with an enhanced capacity for hole transfer.     

B-DNA to zip-DNA: simulating a DNA transition to a novel structure with enhanced charge-transport characteristics

The forced extension of a DNA segment is studied in a series of steered molecular dynamics simulations, employing a broad range of pulling forces. Throughout the entire force range, the formation of a zipper-like (zip-) DNA structure is observed. In that structure, first predicted by Lohikoski et al., the bases of the DNA strands interdigitate with each other and form a single-base aromatic stack. Similar motifs, albeit only a few base pairs in extent, have been observed in experimental crystal structures. Analysis of the dynamics of structural changes in pulled DNA shows that S-form DNA, thought to be adopted by DNA under applied force, serves as an intermediate between B-DNA and zip-DNA. Therefore, the phase transition plateau observed in force-extension curves of DNA is suggested to reflect the B-DNA to zip-DNA structural transition. Electronic structure analysis of purine bases in zip-DNA indicates a several-fold to order of magnitude increase in the π-π electronic coupling among nearest-neighbor nucleobases, compared to B-DNA. We further observe that zip-DNA does not require base pair complementarity between DNA strands, and we predict that the increased electronic coupling in zip-DNA will result in a much higher rate of charge transfer through an all-purine zip-DNA compared to B-DNA of equal length.     

Conformational study of 2-phenylethylamine by molecular-beam Fourier transform microwave spectroscopy

The conformational preferences of the simplest amine neurotransmitter 2-phenylethylamine have been investigated using molecular beam Fourier transform microwave (MB-FTMW) spectroscopy. Two new conformers have been observed together with the two previously reported by Godfrey et al. [J. Am. Chem. Soc., 1995, 117, 8204]. The (14)N nuclear quadrupole hyperfine structure has been resolved for all four conformers. Comparison of the experimental rotational and quadrupole coupling constants with those calculated theoretically provides a conclusive test for the identification of all conformers. The two most stable conformers present a gauche (folded) disposition of the alkyl-amine chain and are stabilised by a weak NH...pi interaction between the amino group and the aromatic ring. The other two conformers show an anti (extended) arrangement of the alkyl-amine chain. Tunnelling splittings have been observed in the spectrum of one of the anti conformers. The post expansion relative abundances in the supersonic jet have been also investigated and related to the conformer energies.     

Femtosecond electron-transfer reactions in mono- and polynucleotides and in DNA

Quenching of redox active, intercalating dyes by guanine bases in DNA can occur on a femtosecond time scale both in DNA and in nucleotide complexes. Notwithstanding the ultrafast rate coefficients, we find that a classical, nonadiabatic Marcus model for electron transfer explains the experimental observations, which allows us to estimate the electronic coupling (330 cm(-1)) and reorganization (8070 cm(-1)) energies involved for thionine-[poly(dG-dC)](2) complexes. Making the simplifying assumption that other charged, pi-stacked DNA intercalators also have approximately these same values, the electron-transfer rate coefficients as a function of the driving force, DeltaG, are derived for similar molecules. The rate of electron transfer is found to be independent of the speed of molecular reorientation. Electron transfer to the thionine singlet excited state from DNA obtained from calf thymus, salmon testes, and the bacterium, micrococcus luteus (lysodeikticus) containing different fractions of G-C pairs, has also been studied. Using a Monte Carlo model for electron transfer in DNA and allowing for reaction of the dye with the nearest 10 bases in the chain, the distance dependence scaling parameter, beta, is found to be 0.8 +/- 0.1 A(-1). The model also predicts the redox potential for guanine dimers, and we find this to be close to the value for isolated guanine bases. Additionally, we find that the pyrimidine bases are barriers to efficient electron transfer within the superexchange limit, and we also infer from this model that the electrons do not cross between strands on the picosecond time scale; that is, the electronic coupling occurs predominantly through the pi-stack and is not increased substantially by the presence of hydrogen bonding within the duplex. We conclude that long-range electron transfer in DNA is not exceptionally fast as would be expected if DNA behaved as a "molecular wire" but nor is it as slow as is seen in proteins, which do not benefit from pi-stacking.     

Pathway analysis of super-exchange electronic couplings in electron transfer reactions using a multi-configuration self-consistent field method

We present a novel pathway analysis of super-exchange electronic couplings in electron transfer reactions using localized molecular orbitals from multi-configuration self-consistent field (MCSCF) calculations. In our analysis, the electronic coupling and the tunneling pathways can be calculated in terms of the configuration interaction (CI) Hamiltonian matrix obtained from the localized MCSCF wave function. Making use of the occupation restricted multiple active spaces (ORMAS) method can effectively produce the donor, acceptor, and intermediate configuration state functions (CSFs) and CIs among these CSFs. In order to express the electronic coupling as a sum of individual tunneling pathways contributions, we employed two perturbative methods: L?wdin projection-iteration method and higher-order super-exchange method. We applied them to anion couplings of butane-1,4-diyl and pentane-1,5-diyl. The results were (1) the electronic couplings calculated from the two perturbative methods were in reasonable agreement with those from a non-perturbative method (one-half value of the energy difference between the ground and first excited states), (2) the main tunneling pathways consisted of a small number of lower-order super-exchange pathways where bonding, anti-bonding, or extra-valence-shell orbitals were used once or twice, and (3) the interference among a huge number of higher-order super-exchange pathways significantly contributed to the overall electronic coupling, whereas each of them contributed only fractionally. Our method can adequately take into account both effects of non-dynamical electron correlation and orbital relaxation. Comparing with the analyses based on the Koopmans' theorem (ignoring both effects) and the ORMAS-CIs from frozen localized reference orbitals (ignoring the effect of orbital relaxation), we discuss these effects.     

Rational design of hetero-ring-expanded guanine analogs with enhanced properties for modified DNA building blocks

The properties and modes of recognition of physiological DNAs associated with the four natural nucleobases might be extended, in principle, by the design of non-natural nucleobase derivatives. The goal is an expansion of the genetic alphabet, with the possible outcome of producing new DNAs with improved physical or biological properties. In this work, a new series of hetero-ring-expanded guanine analogs are proposed, and their relevant structural characteristics and electronic properties are determined by density functional theory. The stabilities of the decamer DNA duplexes (dn.dC)10 (where n represents the corresponding expanded guanine analog designed here) are also examined, using molecular dynamics. The simulations show that the designed motifs can form stable DNA-like structures. We determined the pairing energies for the Watson-Crick (WC) hydrogen-bonded dimers between the expanded G-analogs and the natural C, and found that the pairing energies are close to those of the natural GC pair. The calculated adiabatic ionization potentials (IPs) of the size-expanded guanine analogs and their base pairs, and the corresponding vertical ionization potentials, show that some are distinctly smaller than the corresponding natural versions. The HOMO-LUMO energy gaps for most of the size-expanded guanine analogs and their WC base pairs are considerably lower than those of the corresponding natural base and base pairs. Thus, the expanded G bases may be considered as DNA genetic motifs, and they may serve as building blocks for potential biological applications and the development of molecular electronic devices.     

DFT and MP2 investigations on the hydrogen bonding interaction between 5,6-dihydrothymine and dna bases: A, C, G, T

5,6-Dihydrothymine (DHT) is a nucleobase lesion induced by the action of ionizing radiation on thymine residue in DNA. In this work, we present the hydrogen bonding base pairs involving 5,6-dihydrothymine bound to four bases in DNA: adenine (A), cytosine (C), guanine (G), and thymine (T). Full geometry optimizations are performed for the studied complexes by the B3LYP method. Interaction energies are corrected for the basis set superposition error, using the full Boys-Bernardi counterpoise correction scheme. Hydrogen bonding patterns of these base pairs are characterized using NBO and AIM analysis. According to the calculated binding energies and structural parameters, the stability of the base pairs decrease in the following order: DHT:G ～ DHT:A > DHT:C > DHT:T.     

Balanced carrier transports of electrons and holes in silole-based compounds--a theoretical study

By employing a diabatic model and a first-principle direct method, we have investigated the carrier transport properties of the highly efficient organic light-emitting materials 1,1,2,3,4,5-hexaphenylsilole (HPS) and 1-methyl-1,2,3,4,5-pentaphenylsilole (MPPS). The electronic coupling constants and reorganization energies are calculated for a wide variety of nearest-neighbor charge transfer pathways. The theoretical calculations show that (i) the electron mobility is very close to that of the hole, which indicates a balanced carrier transport in these materials, and (ii) the carrier mobilities for MPPS are larger than those for HPS. These results help explain the underlying microscopic mechanism for the high electroluminescence efficiency.     

核酸碱基互变异构体的结构、稳定性及其物理化学性质

核酸碱基是DNA及RNA分子的重要组成部分, 在基因遗传信息的传递方面起着主导作用. 核酸碱基存在多种互变异构体, 它们在DNA及RNA分子中主要以最稳定的异构体形式存在, 但是在气相或凝聚相中也有少量的其他异构体形式存在. 核酸碱基的稀有互变异构体往往能够引起碱基对的错配对, 这可能会导致DNA及RNA分子形成不规则的结构, 并进一步导致DNA或RNA双螺旋的自发突变. 因此, 对核酸碱基的互变异构体进行系统的研究, 有助于人们深入认识DNA和RNA分子的结构和性质. 国际上有很多研究小组已经通过实验和理论对核酸碱基互变异构体的结构、相对能量及其性质进行了研究. 本文对文献中有关核酸碱基互变异构体的实验和理论研究进行了综述. 在对前人研究进行归纳总结的基础上, 我们利用密度泛函计算对核酸碱基的互变异构体进行了排序, 得到的最优异构体结构参数和相对能量与实验值相比较为一致. 此外, 因为核酸碱基的物理化学性质可以为生物、化学、材料等方面的研究提供重要的基础性信息, 因此我们还对它们的电子亲和能、电离能、质子亲和能等研究进行了总结.     

Theoretical investigation of the mechanism and dynamics of intramolecular coherent resonance energy transfer in soft molecules: a case study of dithia-anthracenophane

A computational study is conducted on dithia-anthracenophane (DTA), for which there is experimental evidence for coherent resonance energy transfer dynamics, and on dimethylanthracene (DMA), a molecule representing the energy donor and the acceptor in DTA. Electronic excitation energies are calculated by configuration interaction singles (CIS) and time-dependent density functional theory (TD-DFT) methods and are compared to experimental ones. Electronic coupling constants are calculated between two DMAs embedded into the ground-state structure of DTA employing methods based on transition densities. The resulting values of electronic coupling provide a more consistent interpretation of experiments than those based on one-half the level spacing of DTA excitation energies. Solvation effects are studied based on the polarizable continuum model (PCM). Solvent-induced polarization and screening effects are shown to make opposite contributions, and the net electronic coupling is little different from the value in a vacuum. The likelihood of coherent population transfer is assessed on the basis of a recently developed theory of coherent resonance energy transfer. The time scale of bath is shown to have an important role in sustaining the quantum coherence. The combination of quantum chemical and dynamical data suggests that the electronic coupling in DTA is in the range of 50-100 cm(-1). The presence of oscillatory excitation population dynamics can be understood from the picture of polaronic excitation moderately dressed with dispersive vibrational modes. The effect of torsional modulation on the excitation energies of DTA and electronic coupling is examined on the basis of optimized structures with the torsional angle constrained. The result suggests that inelastic effect due to torsional motion cannot be disregarded in DTA.