Conformationally Gated Charge Transfer in DNA Three-Way Junctions

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.

Role of intercalation in the electrical properties of nucleic acids for use in molecular electronics

Intercalating ds-DNA/RNA with small molecules can play an essential role in controlling the electron transmission probability for molecular electronics applications such as biosensors, single-molecule transistors, and data storage. However, its applications are limited due to a lack of understanding of the nature of intercalation and electron transport mechanisms. We addressed this long-standing problem by studying the effect of intercalation on both the molecular structure and charge transport along the nucleic acids using molecular dynamics simulations and first-principles calculations coupled with the Green's function method, respectively. The study on anthraquinone and anthraquinone-neomycin conjugate intercalation into short nucleic acids reveals some universal features: (1) the intercalation affects the transmission by two mechanisms: (a) inducing energy levels within the bandgap and (b) shifting the location of the Fermi energy with respect to the molecular orbitals of the nucleic acid, (2) the effect of intercalation was found to be dependent on the redox state of the intercalator: while oxidized anthraquinone decreases, reduced anthraquinone increases the conductance, and (3) the sequence of the intercalated nucleic acid further affects the transmission: lowering the AT-region length was found to enhance the electronic coupling of the intercalator with GC bases, hence yielding an increase of more than four times in conductance. We anticipate our study to inspire designing intercalator-nucleic acid complexes for potential use in molecular electronics via creating a multi-level gating effect.

Revisiting the Hole Size in Double Helical DNA with Localized Orbital Scaling Corrections

The extent of electronic wave function delocalization for the charge carrier (electron or hole) in double helical DNA plays an important role in determining the DNA charge transfer mechanism and kinetics. The size of the charge carrier's wave function delocalization is regulated by the solvation induced localization and the quantum delocalization among the π stacked base pairs at any instant of time. Using a newly developed localized orbital scaling correction (LOSC) density functional theory method, we accurately characterized the quantum delocalization of the hole wave function in double helical B-DNA. This approach can be used to diagnose the extent of delocalization in fluctuating DNA structures. Our studies indicate that the hole state tends to delocalize among 4 guanine-cytosine (GC) base pairs and among 3 adenine-thymine (AT) base pairs when these adjacent bases fluctuate into degeneracy. The relatively small delocalization in AT base pairs is caused by the weaker π-π interaction. This extent of delocalization has significant implications for assessing the role of coherent, incoherent, or flickering coherent carrier transport in DNA.

Tuning molecular fluctuation to boost the conductance in DNA based molecular wires

Inherent molecular fluctuations are known to have a significant influence on the charge transport properties of biomolecules like DNA, PNA and proteins. In this work, we show ways to control these fluctuations and further demonstrate their use to enhance the conductance of two widely studied molecular wires, namely dsDNA (DNA) and G4 Quadruplex (G4-Quad). We quantify the molecular fluctuation in terms of the root mean square deviation (RMSD) of the molecule. In the case of DNA, we use temperature to control the fluctuations, while in the case of G4-Quad the fluctuations are tuned by the ions inside the pore. The electronic coupling between the bases of dsDNA and G4-Quad, which measures the conductance of these molecular wires, shows a non-monotonic behaviour with the increase in fluctuation. We find values of fluctuation which give rise to maximum electronic coupling and hence high conductivity for both the cases. In the case of DNA, these optimal fluctuations (∼2.5 Å) are achieved at a temperature of 210 K, which gives rise to an electronic coupling of 0.135 eV between the DNA bases. The optimal fluctuations in G4-Quad are achieved (∼7 Å) in a 4 base pair long system with 2 Na⁺ ions inside the pore, giving rise to an electronic coupling of 0.09 eV.

The Role of Proton Transfer on Mutations

Hydrogen bonds play a critical role in nucleobase studies as they encode genes, map protein structures, provide stability to the base pairs, and are involved in spontaneous and induced mutations. Proton transfer mechanism is a critical phenomenon that is related to the acid-base characteristics of the nucleobases in Watson-Crick base pairs. The energetic and dynamical behavior of the proton can be depicted from these characteristics and their adjustment to the water molecules or the surrounding ions. Further, new pathways open up in which protonated nucleobases are generated by proton transfer from the ionized water molecules and elimination of a hydroxyl radical in this review, the analysis will be focused on understanding the mechanism of untargeted mutations in canonical, wobble, Hoogsteen pairs, and mutagenic tautomers through the non-covalent interactions. Further, rare tautomer formation through the single proton transfer (SPT) and the double proton transfer (DPT), quantum tunneling in nucleobases, radiation-induced bystander effects, role of water in proton transfer (PT) reactions, PT in anticancer drugs-DNA interaction, displacement and oriental polarization, possible models for mutations in DNA, genome instability, and role of proton transfer using kinetic parameters for RNA will be discussed.

TacoxDNA: A user-friendly web server for simulations of complex DNA structures, from single strands to origami

Simulations of nucleic acids at different levels of structural details are increasingly used to complement and interpret experiments in different fields, from biophysics to medicine and materials science. However, the various structural models currently available for DNA and RNA and their accompanying suites of computational tools can be very rarely used in a synergistic fashion. The tacoxDNA webserver and standalone software package presented here are a step toward a long-sought interoperability of nucleic acids models. The webserver offers a simple interface for converting various common input formats of DNA structures and setting up molecular dynamics (MD) simulations. Users can, for instance, design DNA rings with different topologies, such as knots, with and without supercoiling, by simply providing an XYZ coordinate file of the DNA centre-line. More complex DNA geometries, as designable in the cadnano, CanDo and Tiamat tools, can also be converted to all-atom or oxDNA representations, which can then be used to run MD simulations. Though the latter are currently geared toward the native and LAMMPS oxDNA representations, the open-source package is designed to be further expandable. TacoxDNA is available at http://tacoxdna.sissa.it. © 2019 Wiley Periodicals, Inc.

Why Are DNA and Protein Electron Transfer So Different?

The corpus of electron transfer (ET) theory provides considerable power to describe the kinetics and dynamics of electron flow at the nanoscale. How is it, then, that nucleic acid (NA) ET continues to surprise, while protein-mediated ET is relatively free of mechanistic bombshells? I suggest that this difference originates in the distinct electronic energy landscapes for the two classes of reactions. In proteins, the donor/acceptor-to-bridge energy gap is typically several-fold larger than in NAs. NA ET can access tunneling, hopping, and resonant transport among the bases, and fluctuations can enable switching among mechanisms; protein ET is restricted to tunneling among redox active cofactors and, under strongly oxidizing conditions, a few privileged amino acid side chains. This review aims to provide conceptual unity to DNA and protein ET reaction mechanisms. The establishment of a unified mechanistic framework enabled the successful design of NA experiments that switch electronic coherence effects on and off for ET processes on a length scale of multiple nanometers and promises to provide inroads to directing and detecting charge flow in soft-wet matter.

Detection of Mg2+-dependent, coaxial stacking rearrangements in a bulged three-way DNA junction by single-molecule FRET

Three-way helical junctions (3WJs) arise in genetic processing, and they have architectural and functional roles in structured nucleic acids. An internal bulge at the junction core allows the helical domains to become oriented into two possible, coaxially stacked conformers. Here, the helical stacking arrangements for a series of bulged, DNA 3WJs were examined using ensemble fluorescence resonance energy transfer (FRET) and single-molecule FRET (smFRET) approaches. The 3WJs varied according to the GC content and sequence of the junction core as well as the pyrimidine content of the internal bulge. Mg²⁺ titration experiments by ensemble FRET show that both stacking conformations have similar Mg²⁺ requirements for folding. Strikingly, smFRET experiments reveal that a specific junction sequence can populate both conformers and that this junction undergoes continual interconversion between the two stacked conformers. These findings will support the development of folding principles for the rational design of functional DNA nanostructures.

Charge splitters and charge transport junctions based on guanine quadruplexes

Self-assembling circuit elements, such as current splitters or combiners at the molecular scale, require the design of building blocks with three or more terminals. A promising material for such building blocks is DNA, wherein multiple strands can self-assemble into multi-ended junctions, and nucleobase stacks can transport charge over long distances. However, nucleobase stacking is often disrupted at junction points, hindering electric charge transport between the two terminals of the junction. Here, we show that a guanine-quadruplex (G4) motif can be used as a connector element for a multi-ended DNA junction. By attaching specific terminal groups to the motif, we demonstrate that charges can enter the structure from one terminal at one end of a three-way G4 motif, and can exit from one of two terminals at the other end with minimal carrier transport attenuation. Moreover, we study four-way G4 junction structures by performing theoretical calculations to assist in the design and optimization of these connectors.

Conformational Heterogeneity in a Fully Complementary DNA Three-Way Junction with a GC-Rich Branchpoint

DNA three-way junctions (3WJs) are branched structures that serve as important biological intermediates and as components in DNA nanostructures. We recently derived the global structure of a fully complementary 3WJ and found that it contained unpaired bases at the branchpoint, which is consistent with previous observations of branch flexibility and branchpoint reactivity. By combining high-resolution single-molecule Förster resonance energy transfer, molecular modeling, time-resolved ensemble fluorescence spectroscopy, and the first ¹⁹F nuclear magnetic resonance observations of fully complementary 3WJs, we now show that the 3WJ structure can adopt multiple distinct conformations depending upon the sequence at the branchpoint. A 3WJ with a GC-rich branchpoint adopts an open conformation with unpaired bases at the branch and at least one additional conformation with an increased number of base interactions at the branchpoint. This structural diversity has implications for branch interactions and processing in vivo and for technological applications.

Reaction Path Averaging: Characterizing the Structural Response of the DNA Double Helix to Electron Transfer

A polarizable environment, prominently the solvent, responds to electronic changes in biomolecules rapidly. The knowledge of conformational relaxation of the biomolecule itself, however, may be scarce or missing. In this work, we describe in detail the structural changes in DNA undergoing electron transfer between two adjacent nucleobases. We employ an approach based on averaging of tens to hundreds of thousands of nonequilibrium trajectories generated with molecular dynamics simulation, and a reduction of dimensionality suitable for DNA. We show that the conformational response of the DNA proceeds along a single collective coordinate that represents the relative orientation of two consecutive base pairs, namely, a combination of helical parameters shift and tilt. The structure of DNA relaxes on time scales reaching nanoseconds, contributing marginally to the relaxation of energies, which is dominated by the modes of motion of the aqueous solvent. The concept of reaction path averaging (RPA), conveniently exploited in this context, makes it possible to filter out any undesirable noise from the nonequilibrium data, and is applicable to any chemical process in general.