DNA Charge Transfer: An Atomistic Model

Ultrafast excitonic dynamics in DNA: Bridging correlated quantum dynamics and sequence dependence

After photoexcitation of DNA, the excited electron (in the LUMO) and the remaining hole (in the HOMO) localized on the same DNA base form a bound pair, called the Frenkel exciton, due to their mutual Coulomb interaction. In this study, we demonstrate that a tight-binding (TB) approach, using TB parameters for electrons and holes available in the literature, allows us to correlate relaxation properties, average charge separation, and dipole moments to a large ensemble of double-stranded DNA sequences (all 16384 possible sequences with 14 nucleobases). This way, we are able to identify a relatively small subensemble of sequences responsible for long-lived excited states, high average charge separation, and high dipole moment. Further analysis shows that these sequences are particularly T rich. By systematically screening the impact of electron-hole interaction (Coulomb forces), we verify that these correlations are relatively robust against finite-size variations of the interaction parameter, not directly accessible experimentally. This methodology combines simulation methods from quantum physics and physical chemistry with statistical analysis known from genetics and epigenetics, thus representing a powerful bridge to combine information from both fields.

A quantum physics layer of epigenetics: a hypothesis deduced from charge transfer and chirality-induced spin selectivity of DNA

BACKGROUND

Epigenetic mechanisms are informational cellular processes instructing normal and diseased phenotypes. They are associated with DNA but without altering the DNA sequence. Whereas chemical processes like DNA methylation or histone modifications are well-accepted epigenetic mechanisms, we herein propose the existence of an additional quantum physics layer of epigenetics.

RESULTS

We base our hypothesis on theoretical and experimental studies showing quantum phenomena to be active in double-stranded DNA, even under ambient conditions. These phenomena include coherent charge transfer along overlapping pi-orbitals of DNA bases and chirality-induced spin selectivity. Charge transfer via quantum tunneling mediated by overlapping orbitals results in charge delocalization along several neighboring bases, which can even be extended by classical (non-quantum) electron hopping. Such charge transfer is interrupted by flipping base(s) out of the double-strand e.g., by DNA modifying enzymes. Charge delocalization can directly alter DNA recognition by proteins or indirectly by DNA structural changes e.g., kinking. Regarding sequence dependency, charge localization, shown to favor guanines, could influence or even direct epigenetic changes, e.g., modification of cytosines in CpG dinucleotides. Chirality-induced spin selectivity filters electrons for their spin along DNA and, thus, is not only an indicator for quantum coherence but can potentially affect DNA binding properties.

CONCLUSIONS

Quantum effects in DNA are prone to triggering and manipulation by external means. By the hypothesis put forward here, we would like to foster research on "Quantum Epigenetics" at the interface of medicine, biology, biochemistry, and physics to investigate the potential epigenetic impact of quantum physical principles on (human) life.

The Dynamics of Hole Transfer in DNA

High-energy radiation and oxidizing agents can ionize DNA. One electron oxidation gives rise to a radical cation whose charge (hole) can migrate through DNA covering several hundreds of Å, eventually leading to irreversible oxidative damage and consequent disease. Understanding the thermodynamic, kinetic and chemical aspects of the hole transport in DNA is important not only for its biological consequences, but also for assessing the properties of DNA in redox sensing or labeling. Furthermore, due to hole migration, DNA could potentially play an important role in nanoelectronics, by acting as both a template and active component. Herein, we review our work on the dynamics of hole transfer in DNA carried out in the last decade. After retrieving the thermodynamic parameters needed to address the dynamics of hole transfer by voltammetric and spectroscopic experiments and quantum chemical computations, we develop a theoretical methodology which allows for a faithful interpretation of the kinetics of the hole transport in DNA and is also capable of taking into account sequence-specific effects.

Modeling DNA oxidation in water

A novel set of hole-site energies and electronic coupling parameters to be used, in the framework of the simplest tight-binding approximation, for predicting DNA hole trapping efficiencies and rates of hole transport in oxidized DNA is proposed. The novel parameters, significantly different from those previously reported in the literature, have been inferred from reliable density functional calculations, including both the sugar-phosphate ionic backbone and the effects of the aqueous environment. It is shown that most of the experimental oxidation free energies of DNA tracts and of oligonucleotides available from photoelectron spectroscopy and voltammetric measurements are reproduced with great accuracy, without the need for introducing sequence dependent parameters.

Thermodynamic integration network study of electron transfer: from proteins to aggregates

We describe electron transfer through the NrfHA nitrite reductase using a thermodynamic integration scheme. Driving forces are hardly affected by dimerization, but the transport mechanism only emerges simulating the dimer.

Sensing organic molecules by charge transfer through aptamer-target complexes: theory and simulation

Aptamers, i.e., short sequences of RNA and single-stranded DNA, are capable of specificilly binding objects ranging from small molecules over proteins to entire cells. Here, we focus on the structure, stability, dynamics, and electronic properties of oligonucleotides that interact with aromatic or heterocyclic targets. Large-scale molecular dynamics simulations indicate that aromatic rings such as dyes, metabolites, or alkaloides form stable adducts with their oligonucleotide host molecules at least on the simulation time scale. From molecular dynamics snapshots, the energy parameters relevant to Marcus' theory of charge transfer are computed using a modified Su-Schrieffer-Heeger Hamiltonian, permitting an estimate of the charge transfer rates. In many cases, aptamer binding seriously influences the charge transfer kinetics and the charge carrier mobility within the complex, with conductivities up to the nanoampere range for a single complex. We discuss the conductivity properties with reference to potential applications as biosensors.

Quantum transport anomalies in DNA containing mispairs

The effect of mispairs on charge transport in DNA of sequence (GC)(TA)(N)(GC)(3) connected to platinum electrodes is studied using the tight-binding model. With parameters derived from an ab initio density functional result, we calculate the current versus bias voltage for DNA with and without a mispair and for different numbers of (TA) basepairs N between the single and triple (GC) basepairs. The current decays exponentially with N under low bias but reaches a minimum under high bias when a multichannel transport mechanism is established. A (GA) mispair substituting a (TA) basepair near the middle of the (TA)(N) sequence usually enhances the current by one order due to its low ionization energy but may decrease the current significantly when an established multichannel mechanism is broken.

Electronic parameters for charge transfer along DNA

We systematically examine all the tight-binding parameters pertinent to charge transfer along DNA. The pi molecular structure of the four DNA bases (adenine, thymine, cytosine, and guanine) is investigated by using the linear combination of atomic orbitals method with a recently introduced parametrization. The HOMO and LUMO wave functions and energies of DNA bases are discussed and then used for calculating the corresponding wave functions of the two B-DNA base-pairs (adenine-thymine and guanine-cytosine). The obtained HOMO and LUMO energies of the bases are in good agreement with available experimental values. Our results are then used for estimating the complete set of charge transfer parameters between neighboring bases and also between successive base-pairs, considering all possible combinations between them, for both electrons and holes. The calculated microscopic quantities can be used in mesoscopic theoretical models of electron or hole transfer along the DNA double helix, as they provide the necessary parameters for a tight-binding phenomenological description based on the pi molecular overlap. We find that usually the hopping parameters for holes are higher in magnitude compared to the ones for electrons. Our findings are also compared with existing calculations from first principles.

A Direct Simulation of Adiabatic Charge Transfer Through Bridged Organic Molecules

We approach the dynamics of adiabatic charge transfer through bridged triarylamine cations by a direct molecular dynamics simulation involving classical and quantum mechanical degrees of freedom. Within a simple yet chemically specific model, the quantum mechanical subsystem is described by a tight-binding Hamiltonian, which is coupled to a classical force field. From a population analysis of the quantum part, the charge transfer rate can be readily extracted, including the influence of memory effects. The direct computation of the associated thermodynamic potential establishes a close link to analytical rate concepts. The theoretical data are compared to experiments, and the limits and possible extensions of our approach are discussed.

Solvent reorganization energy of hole transfer in DNA

We present an application of a molecular-dynamics-based scheme to evaluate the solvent reorganization energy of hole transfer in DNA. The obtained parameters can be used for simulations of hole transfer in DNA by means of Marcus' theory. Also, we perform an analysis of the reorganization energies, including the case of transfer of a delocalized hole.

Direct simulation of electron transfer reactions in DNA radical cations

The electron transfer properties of DNA radical cations are important in DNA damage and repair processes. Fast long-range charge transfer has been demonstrated experimentally, but the subtle influences that experimental conditions as well as DNA sequences and geometries have on the details of electron transfer parameters are still poorly understood. In this work, we employ an atomistic QM/MM approach, based on a one-electron tight binding Hamiltonian and a classical molecular mechanics forcefield, to conduct nanosecond length MD simulations of electron holes in DNA oligomers. Multiple spontaneous electron transfer events were observed in 100 ns simulations with neighboring adenine or guanine bases. Marcus parameters of charge transfer could be extracted directly from the simulations. The reorganization energy lambda for hopping between neighboring bases was found to be ca. 25 kcal/mol and charge transfer rates of 4.1 x 10(9) s(-1) for AA hopping and 1.3 x 10(9) s(-1) for GG hopping were obtained.

Local Aromaticity in Natural Nucleobases and Their Size-Expanded Benzo-Fused Derivatives

The influence of the insertion/addition of a benzene ring to the natural nucleic acid bases on the local aromaticity of the so-called size-expanded (xN, with N being adenine, guanine, cytosine, and thymine) bases is examined. To this end, the local aromaticity of the six- and five-membered rings in both the natural bases and their benzoderivatives is determined using HOMA, NICS, aromatic fluctuation index (FLU), and para-delocalization index (PDI) descriptors. In general, there is a good correspondence between the different indices, so that ring moieties with more negative NICS values also have larger HOMA and PDI measures and lower FLU indices. The results also point out notable differences in the aromatic character of the natural and size-expanded bases, which generally are hardly affected upon hydrogen bonding. The differences in the highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) gap determined for the size-expanded nucleobases show an inverse correlation with the aromaticity of the fused benzene ring, so that the larger the HOMO-LUMO gap is, the lower the destabilization experienced by the benzene upon insertion/addition to the natural bases. This finding suggests that the introduction of suitable chemical modifications in the benzene ring might be useful to modulate the HOMO-LUMO gap while enabling the design of modified DNA duplexes that are able to act as molecular wires.

Electronic conductivity in polyaromatic hydrocarbon glasses: A theoretical perspective

Based upon Monte Carlo simulations of amorphous molecular glasses, we have computed the electronic structure of five prototypical polyaromatic hydrocarbons using an extended Su-Schrieffer-Heeger model [J. R. Schrieffer, W. P. Su, and A. J. Heeger, Phys. Rev. Lett. 42, 1698 (1979)]. In the presence of excess charges, the resulting potential energy surfaces have been analyzed using Marcus' [Annu. Rev. Phys. Chem. 15, 155 (1964)] theory of charge transfer to yield reaction coefficients and--via the application of linear response theory--local conductivities. Applying Kirchhoff's rules, the emerging random resistor network problem leads to global conductivities of the order of 10(-1)-1 Scm, which correlate with the structural characteristics of the underlying geometry.

Environmental Effect on the Relative Contribution of the Charge-Transfer Mechanisms within a Short DNA Sequence

Time evolution of the charge-transfer site population is studied in a short DNA sequence to determine the type of governing charge-transfer mechanism. The system consists of a 5'-GAGGG-3' nucleobase sequence coupled with a dissipative bath that represents the DNA phosphate backbone and solvents. Relative contribution of transfer mechanisms to the whole charge-transfer process has been obtained using the on-the-fly filtered propagator functional path integral method with the density matrix decomposition. Partial density matrixes of the incoherent hopping and coherent superexchange pathways as well as the full reduced density matrix have been evaluated and discussed for both debye and ohmic baths. It was found that the relative contribution of the transfer mechanisms is rather sensitive to the frequency-dependent environmental description.

Static and dynamic aspects of DNA charge transfer: a theoretical perspective

In this work, we approach the impact of dynamic and static disorder on DNA charge transfer from a theoretical and numerical perspective. Disordered or defect geometries are either realized via molecular dynamics simulations using a classical force field or by experimentally determined DNA bulge structures. We apply a chemically specific, atomically resolved extended Su-Schrieffer-Heeger model to compute the energy parameters relevant to DNA charge transfer. For both models studied here, the effective donor-acceptor couplings--and hence the charge transfer rates--significantly depend upon the geometry. Dynamic disorder leads to a correlation time in this quantity of the order of 30 fs, and the transfer rates universally exhibit a broad, yet well-defined, exponential distribution. For DNA bulges, the angle characterizing the defect controls the charge transfer efficiency. The results are discussed and extensively compared to experimental findings and other calculations.