Microsecond Simulation of Electron Transfer in DNA: Bottom-Up Parametrization of an Efficient Electron Transfer Model Based on Atomistic Details

Intra-strand phosphate-mediated pathways in microsolvated double-stranded DNA

We argue that dry DNA charge transport in molecular junctions, over distances of tens of nanometers, can take place via independent intra-strand pathways involving the phosphate groups. Such pathways explain recent single-molecule experiments that compare currents in intact and nicked 100 base-pair double-stranded DNA. We explore the conditions that favor independent intra-strand transport channels with the participation of the phosphate groups, as opposed to purely base-mediated transport involving the pi-stacked bases and inter-strand transitions. Our computations demonstrate how long-distance transport pathways in DNA are tuned by the degree of solvation, which affects the level of dynamic disorder in the pi-stacking, and the energies of phosphate-group molecular orbitals.

DNA damage and repair in the nucleosome: insights from computational methods

Cellular DNA is constantly exposed to endogenous or exogenous factors that can induce lesions. Several types of lesions have been described that can result from UV/ionizing irradiations, oxidative stress, or free radicals, among others. In order to overcome the deleterious effects of such damages, i.e., mutagenicity or cytotoxicity, cells possess a highly complex DNA repair machinery, involving repair enzymes targeting specific types of lesions through dedicated cellular pathways. In addition, DNA is highly compacted in the nucleus, the first level of compaction consisting of ~ 147 DNA base pairs wrapped around a core of histones, the so-called nucleosome core particle. In this complex environment, the DNA structure is highly constrained, and fine-tuned mechanisms involving remodeling processes are required to expose the DNA to repair enzymes and to facilitate the damage removal. However, these nucleosome-specific mechanisms remain poorly understood, and computational methods emerged only recently as powerful tools to investigate DNA damages in such complex systems as the nucleosome. In this mini-review, we summarize the latest advances brought out by computational approaches in the field, opening new exciting perspectives for the study of DNA damage and repair in the nucleosome context.

Computational study of the role of counterions and solvent dielectric in determining the conductance of B-DNA

DNA naturally exists in a solvent environment, comprising water and salt molecules such as sodium, potassium, magnesium, etc. Along with the sequence, the solvent conditions become a vital factor determining DNA structure and thus its conductance. Over the last two decades, researchers have measured DNA conductivity both in hydrated and almost dry (dehydrated) conditions. However, due to experimental limitations (the precise control of the environment), it is very difficult to analyze the conductance results in terms of individual contributions of the environment. Therefore, modeling studies can help us to gain a valuable understanding of various factors playing a role in charge transport phenomena. DNA naturally has negative charges located at the phosphate groups in the backbone, which provides both the connections between the base pairs and the structural support for the double helix. Positively charged ions such as the sodium ion (Na^{+}), one of the most commonly used counterions, balance the negative charges at the backbone. This modeling study investigates the role of counterions both with and without the solvent (water) environment in charge transport through double-stranded DNA. Our computational experiments show that in dry DNA, the presence of counterions affects electron transmission at the lowest unoccupied molecular orbital energies. However, in solution, the counterions have a negligible role in transmission. Using the polarizable continuum model calculations, we demonstrate that the transmission is significantly higher at both the highest occupied and lowest unoccupied molecular orbital energies in a water environment as opposed to in a dry one. Moreover, calculations also show that the energy levels of neighboring bases are more closely aligned to ease electron flow in the solution.

Long-range charge transport in homogeneous and alternating-rigidity chains

We study the interplay of intrinsic-electronic and environmental factors in long-range charge transport across molecular chains with up to N ∼ 80 monomers. We describe the molecular electronic structure of the chain with a tight-binding Hamiltonian. Thermal effects in the form of electron decoherence and inelastic scattering are incorporated with the Landauer–Büttiker probe method. In short chains of up to ten units, we observe the crossover between coherent (tunneling, ballistic) motion and thermally-assisted conduction, with thermal effects enhancing the current beyond the quantum coherent limit. We further show that unconventional (nonmonotonic with size) transport behavior emerges when monomer-to-monomer electronic coupling is made large. In long chains, we identify a different behavior, with thermal effects suppressing the conductance below the coherent-ballistic limit. With the goal to identify a minimal model for molecular chains displaying unconventional and effective long-range transport, we simulate a modular polymer with alternating regions of high and low rigidity. Simulations show that, surprisingly, while charge correlations are significantly affected by structuring environmental conditions, reflecting charge delocalization, the electrical resistance displays an averaging effect, and it is not sensitive to this patterning. We conclude by arguing that efficient long-range charge transport requires engineering both internal electronic parameters and environmental conditions.

Fine-tuning the DNA conductance by intercalation of drug molecules

In this work we study the structure-transport property relationships of small ligand intercalated DNA molecules using a multiscale modeling approach where extensive ab initio calculations are performed on numerous MD-simulated configurations of dsDNA and dsDNA intercalated with two different intercalators, ethidium and daunomycin. DNA conductance is found to increase by one order of magnitude upon drug intercalation due to the local unwinding of the DNA base pairs adjacent to the intercalated sites, which leads to modifications of the density of states in the near-Fermi-energy region of the ligand-DNA complex. Our study suggests that the intercalators can be used to enhance or tune the DNA conductance, which opens new possibilities for their potential applications in nanoelectronics.

Predicting the DNA Conductance Using a Deep Feedforward Neural Network Model

Double-stranded DNA (dsDNA) has been established as an efficient medium for charge migration, bringing it to the forefront of the field of molecular electronics and biological research. The charge migration rate is controlled by the electronic couplings between the two nucleobases of DNA/RNA. These electronic couplings strongly depend on the intermolecular geometry and orientation. Estimating these electronic couplings for all the possible relative geometries of molecules using the computationally demanding first-principles calculations requires a lot of time and computational resources. In this article, we present a machine learning (ML)-based model to calculate the electronic coupling between any two bases of dsDNA/dsRNA and bypass the computationally expensive first-principles calculations. Using the Coulomb matrix representation which encodes the atomic identities and coordinates of the DNA base pairs to prepare the input dataset, we train a feedforward neural network model. Our neural network (NN) model can predict the electronic couplings between dsDNA base pairs with any structural orientation with a mean absolute error (MAE) of less than 0.014 eV. We further use the NN-predicted electronic coupling values to compute the dsDNA/dsRNA conductance.

Toward a First-Principles Evaluation of Transport Mechanisms in Molecular Wires

Understanding charge transport through molecular wires is important for nanoscale electronics and biochemistry. Our goal is to establish a simple first-principles protocol for predicting the charge transport mechanism in such wires, in particular the crossover from coherent tunneling for short wires to incoherent hopping for longer wires. This protocol is based on a combination of density functional theory with a polarizable continuum model introduced by Kaupp et al. for mixed-valence molecules, which we had previously found to work well for length-dependent charge delocalization in such systems. We combine this protocol with a new charge delocalization measure tailored for molecular wires, and we show that it can predict the tunneling-to-hopping transition length with a maximum error of one subunit in five sets of molecular wires studied experimentally in molecular junctions at room temperature. This suggests that the protocol is also well suited for estimating the extent of hopping sites as relevant, for example, for the intermediate tunneling-hopping regime in DNA.

Multiscale modelling reveals higher charge transport efficiencies of DNA relative to RNA independent of mechanism

In this study, we compare the charge transport properties of multiple double-stranded (ds)RNA sequences with corresponding dsDNA sequences. Recent studies have presented a contradictory picture of relative charge transport efficiencies in A-form DNA : RNA hybrids and dsDNA. Using a multiscale modelling framework, we compute conductance of dsDNA and dsRNA using Landauer formalism in the coherent limit and Marcus-Hush theory in the incoherent limit. We find that dsDNA conducts better than dsRNA in both the charge transport regimes. Our analysis shows that the structural differences in the twist angle and slide of dsDNA and dsRNA are the main reasons behind the higher conductance of dsDNA in the incoherent hopping regime. In the coherent limit however, for the same base pair length, the conductance of dsRNA is higher than that of dsDNA for the morphologies where dsRNA has a smaller end-to-end length relative to that of dsDNA.

2'-Deoxy-2'-fluoro-arabinonucleic acid: a valid alternative to DNA for biotechnological applications using charge transport

The non-biological 2'-deoxy-2'-fluoro-arabinonucleic acid (2'F-ANA) may be used as a valid alternative to DNA in biomedical and electronic applications because of its higher resistance to hydrolysis and nuclease degradation. However, the advantage of using 2'F-ANA in such applications also depends on its charge-transfer properties compared to DNA. In this study, we compare the charge conduction properties of model 2'F-ANA and DNA double-strands, using structural snapshots from MD simulations to calculate the electronic couplings and reorganization energies associated with the hole transfer steps between adjacent nucleobase pairs. Inserting these charge-transfer parameters into a kinetic model for charge conduction, we find similar conductive properties for DNA and 2'F-ANA. Moreover, we find that 2'F-ANA's enhanced chemical stability does not correspond to a reduction in the nucleobase π-stack structural flexibility relevant to both electronic couplings and reorganization free energies. Our results promote the use of 2'F-ANA in applications that can be based on charge transport, such as biosensing and chip technology, where its chemical stability and conductivity can advantageously combine.

Machine Learning Prediction of DNA Charge Transport

First-principles calculations of charge transfer in DNA molecules are computationally expensive given that conducting charge carriers interact with intra- and intermolecular atomic motion. Screening sequences, for example, to identify excellent electrical conductors, is challenging even when adopting coarse-grained models and effective computational schemes that do not explicitly describe atomic dynamics. We present a machine learning (ML) model that allows the inexpensive prediction of the electrical conductance of millions of long double-stranded DNA (dsDNA) sequences, reducing computational costs by orders of magnitude. The algorithm is trained on short DNA nanojunctions with n = 3-7 base pairs. The electrical conductance of the training set is computed with a quantum scattering method, which captures charge-nuclei scattering processes. We demonstrate that the ML method accurately predicts the electrical conductance of varied dsDNA junctions tracing different transport mechanisms: coherent (short-range) quantum tunneling, on-resonance (ballistic) transport, and incoherent site-to-site hopping. Furthermore, the ML approach supports physical observations that clusters of nucleotides regulate DNA transport behavior. The input features tested in this work could be used in other ML studies of charge transport in complex polymers in the search for promising electronic and thermoelectric materials.

A Molecular Dynamics-Quantum Mechanics Theoretical Study of DNA-Mediated Charge Transport in Hydrated Ionic Liquids

Charge transport (CT) through biomolecules is of high significance in the research fields of biology, nanotechnology, and molecular devices. Inspired by our previous work that showed the binding of ionic liquid (IL) facilitated charge transport in duplex DNA, in silico simulation is a useful means to understand the microscopic mechanism of the facilitation phenomenon. Here molecular dynamics simulations (MD) of duplex DNA in water and hydrated ionic liquids were employed to explore the helical parameters. Principal component analysis was further applied to capture the subtle conformational changes of helical DNA upon different environmental impacts. Sequentially, CT rates were calculated by a QM/MM simulation of the flickering resonance model based upon MD trajectories. Herein, MD simulation illustrated that the binding of ionic liquids can restrain dynamic conformation and lower the on-site energy of the DNA base. Confined movement among the adjacent base pairs was highly related to the increase of electronic coupling among base pairs, which may lead DNA to a CT facilitated state. Sequentially combining MD and QM/MM analysis, the rational correlations among the binding modes, the conformational changes, and CT rates illustrated the facilitation effects from hydrated IL on DNA CT and supported a conformational-gating mechanism.

Approximate DFT-based methods for generating diabatic states and calculating electronic couplings: models of two and more states

Four types of density functional theory (DFT)-based approaches are assessed in this work for the approximate construction of diabatic states and the evaluation of electronic couplings between these states. These approaches include the constrained DFT (CDFT) method, the constrained noninteracting electron (CNE) model to post-process Kohn-Sham operators, the approximate block-diagonalization (BD) of the Kohn-Sham operators, and the generalized Mulliken-Hush method. It is shown that the first three approaches provide a good description for long-distance intermolecular electron transfer (ET) reactions. On the other hand, inconsistent results were found when applying these approaches to intramolecular ET in strongly coupled, mixed-valence systems. Model analysis shows that this discrepancy is caused by the inappropriate use of the two-state model rather than the defects of the approaches themselves. The situation is much improved when more states are included in the model electronic Hamiltonian. The CNE and BD approaches can thus serve as efficient and robust alternatives for building ET models based on DFT calculations.

Remarkable similarity of force induced dsRNA conformational changes to stretched dsDNA and their detection using electrical measurements

We show the emergence of S-RNA under stretching in analogy to S-DNA and propose a method for its detection using electrical measurement.

Functional role of an unusual tyrosine residue in the electron transfer chain of a prokaryotic (6-4) photolyase

Cryptochromes and photolyases form a flavoprotein family in which the FAD chromophore undergoes light induced changes of its redox state. During this process, termed photoreduction, electrons flow from the surface via conserved amino acid residues to FAD. The bacterial (6-4) photolyase PhrB belongs to a phylogenetically ancient group. Photoreduction of PhrB differs from the typical pattern because the amino acid of the electron cascade next to FAD is a tyrosine (Tyr391), whereas photolyases and cryptochromes of other groups have a tryptophan as direct electron donor of FAD. Mutagenesis studies have identified Trp342 and Trp390 as essential for charge transfer. Trp342 is located at the periphery of PhrB while Trp390 connects Trp342 and Tyr391. The role of Tyr391, which lies between Trp390 and FAD, is however unclear as its replacement by phenylalanine did not block photoreduction. Experiments reported here, which replace Tyr391 by Ala, show that photoreduction is blocked, underlining the relevance of Tyr/Phe at position 391 and indicating that charge transfer occurs via the triad 391-390-342. This raises the question, why PhrB positions a tyrosine at this location, having a less favourable ionisation potential than tryptophan, which occurs at this position in many proteins of the photolyase/cryptochrome family. Tunnelling matrix calculations show that tyrosine or phenylalanine can be involved in a productive bridged electron transfer between FAD and Trp390, in line with experimental findings. Since replacement of Tyr391 by Trp resulted in loss of FAD and DMRL chromophores, electron transfer cannot be studied experimentally in this mutant, but calculations on a mutant model suggest that Trp might participate in the electron transfer cascade. Charge transfer simulations reveal an unusual stabilization of the positive charge on site 391 compared to other photolyases or cryptochromes. Water molecules near Tyr391 offer a polar environment which stabilizes the positive charge on this site, thereby lowering the energetic barrier intrinsic to tyrosine. This opens a second charge transfer channel in addition to tunnelling through the tyrosine barrier, based on hopping and therefore transient oxidation of Tyr391, which enables a fast charge transfer similar to proteins utilizing a tryptophan-triad. Our results suggest that evolution of the first site of the redox chain has just been possible by tuning the protein structure and environment to manage a downhill hole transfer process from FAD to solvent.