Evidence for a near-resonant charge transfer mechanism for double-stranded peptide nucleic acid

Electronic Properties of DNA Origami Nanostructures Revealed by In Silico Calculations

DNA origami is a pioneering approach for producing complex 2- or 3-D shapes for use in molecular electronics due to its inherent self-assembly and programmability properties. The electronic properties of DNA origami structures are not yet fully understood, limiting the potential applications. Here, we conduct a theoretical study with a combination of molecular dynamics, first-principles, and charge transmission calculations. We use four separate single strand DNAs, each having 8 bases (4 × G4C4 and 4 × A4T4), to form two different DNA nanostructures, each having two helices bundled together with one crossover. We also generated double-stranded DNAs to compare electronic properties to decipher the effects of crossovers and bundle formations. We demonstrate that density of states and band gap of DNA origami depend on its sequence and structure. The crossover regions could reduce the conductance due to a lack of available states near the HOMO level. Furthermore, we reveal that, despite having the same sequence, the two helices in the DNA origami structure could exhibit different electronic properties, and electrode position can affect the resulting conductance values. Our study provides better understanding of the electronic properties of DNA origamis and enables us to tune these properties for electronic applications such as nanowires, switches, and logic gates.

Role of Order in the Mechanism of Charge Transport across Single-Stranded and Double-Stranded DNA Monolayers in Tunnel Junctions

Deoxyribonucleic acid (DNA) has been hypothesized to act as a molecular wire due to the presence of an extended π-stack between base pairs, but the factors that are detrimental in the mechanism of charge transport (CT) across tunnel junctions with DNA are still unclear. Here we systematically investigate CT across dense DNA monolayers in large-area biomolecular tunnel junctions to determine when intrachain or interchain CT dominates and under which conditions the mechanism of CT becomes thermally activated. In our junctions, double-stranded DNA (dsDNA) is 30-fold more conductive than single-stranded DNA (ssDNA). The main reason for this large change in conductivity is that dsDNA forms ordered monolayers where intrachain tunneling dominates, resulting in high CT rates. By varying the temperature T and the length of the DNA fragments in the junctions, which determines the tunneling distance, we reveal a complex interplay between T, the length of DNA, and structural order on the mechanism of charge transport. Both the increase in the tunneling distance and the decrease in structural order result in a change in the mechanism of CT from coherent tunneling to incoherent tunneling (hopping). Our results highlight the importance of the interplay between structural order, tunneling distance, and temperature on the CT mechanism across DNA in molecular junctions.

Delocalization-Assisted Transport through Nucleic Acids in Molecular Junctions

The flow of charge through molecules is central to the function of supramolecular machines, and charge transport in nucleic acids is implicated in molecular signaling and DNA repair. We examine the transport of electrons through nucleic acids to understand the interplay of resonant and nonresonant charge carrier transport mechanisms. This study reports STM break junction measurements of peptide nucleic acids (PNAs) with a G-block structure and contrasts the findings with previous results for DNA duplexes. The conductance of G-block PNA duplexes is much higher than that of the corresponding DNA duplexes of the same sequence; however, they do not display the strong even-odd dependence conductance oscillations found in G-block DNA. Theoretical analysis finds that the conductance oscillation magnitude in PNA is suppressed because of the increased level of electronic coupling interaction between G-blocks in PNA and the stronger PNA-electrode interaction compared to that in DNA duplexes. The strong interactions in the G-block PNA duplexes produce molecular conductances as high as 3% G₀, where G₀ is the quantum of conductance, for 5 nm duplexes.

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.

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.

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.

Towards molecular electronic devices based on 'all-carbon' wires

Nascent molecular electronic devices based on linear 'all-carbon' wires attached to gold electrodes through robust and reliable C-Au contacts are prepared via efficient in situ sequential cleavage of trimethylsilyl end groups from an oligoyne, Me₃Si-(C[triple bond, length as m-dash]C)₄-SiMe₃ (1). In the first stage of the fabrication process, removal of one trimethylsilyl (TMS) group in the presence of a gold substrate, which ultimately serves as the bottom electrode, using a stoichiometric fluoride-driven process gives a highly-ordered monolayer, Au|C[triple bond, length as m-dash]CC[triple bond, length as m-dash]CC[triple bond, length as m-dash]CC[triple bond, length as m-dash]CSiMe₃ (Au|C₈SiMe₃). In the second stage, treatment of Au|C₈SiMe₃ with excess fluoride results in removal of the remaining TMS protecting group to give a modified monolayer Au|C[triple bond, length as m-dash]CC[triple bond, length as m-dash]CC[triple bond, length as m-dash]CC[triple bond, length as m-dash]CH (Au|C₈H). The reactive terminal C[triple bond, length as m-dash]C-H moiety in Au|C₈H can be modified by 'click' reactions with (azidomethyl)ferrocene (N₃CH₂Fc) to introduce a redox probe, to give Au|C₆C₂N₃HCH₂Fc. Alternatively, incubation of the modified gold substrate supported monolayer Au|C₈H in a solution of gold nanoparticles (GNPs), results in covalent attachment of GNPs on top of the film via a second alkynyl carbon-Au σ-bond, to give structures Au|C₈|GNP in which the monolayer of linear, 'all-carbon' C₈ chains is sandwiched between two macroscopic gold contacts. The covalent carbon-surface bond as well as the covalent attachment of the metal particles to the monolayer by cleavage of the alkyne C-H bond is confirmed by surface-enhanced Raman scattering (SERS). The integrity of the carbon chain in both Au|C₆C₂N₃HCH₂Fc systems and after formation of the gold top-contact electrode in Au|C₈|GNP is demonstrated through electrochemical methods. The electrical properties of these nascent metal-monolayer-metal devices Au|C₈|GNP featuring 'all-carbon' molecular wires were characterised by sigmoidal I-V curves, indicative of well-behaved junctions free of short circuits.

Effects of the Backbone and Chemical Linker on the Molecular Conductance of Nucleic Acid Duplexes

Scanning tunneling microscope break junction measurements are used to examine how the molecular conductance of nucleic acids depends on the composition of their backbone and the linker group to the electrodes. Molecular conductances of 10 base pair long homoduplexes of DNA, aeg-PNA, γ-PNA, and a heteroduplex of DNA/aeg-PNA with identical nucleobase sequence were measured. The molecular conductance was found to vary by 12 to 13 times with the change in backbone. Computational studies show that the molecular conductance differences between nucleic acids of different backbones correlate with differences in backbone structural flexibility. The molecular conductance was also measured for duplexes connected to the electrode through two different linkers, one directly to the backbone and one directly to the nucleobase stack. While the linker causes an order-of-magnitude increase in the overall conductance for a particular duplex, the differences in the electrical conductance with backbone composition are preserved. The highest molecular conductance value, 0.06G₀, was measured for aeg-PNA duplexes with a base stack linker. These findings reveal an important new strategy for creating longer and more complex electroactive, nucleic acid assemblies.

Conductance in a bis-terpyridine based single molecular breadboard circuit

Controlling charge flow in single molecule circuits with multiple electrical contacts and conductance pathways is a much sought after goal in molecular electronics. In this joint experimental and theoretical study, we advance the possibility of creating single molecule breadboard circuits through an analysis of the conductance of a bis-terpyridine based molecule (TP1). The TP1 molecule can adopt multiple conformations through relative rotations of 7 aromatic rings and can attach to electrodes in 61 possible single and multi-terminal configurations through 6 pyridyl groups. Despite this complexity, we show that it is possible to extract well defined conductance features for the TP1 breadboard and assign them rigorously to the underlying constituent circuits. Mechanically controllable break-junction (MCBJ) experiments on the TP1 molecular breadboard show an unprecedented 4 conductance states spanning a range 10 -2G0 to 10 -7G0. Quantitative theoretical examination of the conductance of TP1 reveals that combinations of 5 types of single terminal 2-5 ring subcircuits are accessed as a function of electrode separation to produce the distinct conductance steps observed in the MCBJ experiments. We estimate the absolute conductance for each single terminal subcircuit and its percentage contribution to the 4 experimentally observed conductance states. We also provide a detailed analysis of the role of quantum interference and thermal fluctuations in modulating conductance within the subcircuits of the TP1 molecular breadboard. Finally, we discuss the possible development of molecular circuit theory and experimental advances necessary for mapping conductance through complex single molecular breadboard circuits in terms of their constituent subcircuits.

Electron Transfer Proteins as Electronic Conductors: Significance of the Metal and Its Binding Site in the Blue Cu Protein, Azurin

Electron transfer (ET) proteins are biomolecules with specific functions, selected by evolution. As such they are attractive candidates for use in potential bioelectronic devices. The blue copper protein azurin (Az) is one of the most-studied ET proteins. Traditional spectroscopic, electrochemical, and kinetic methods employed for studying ET to/from the protein's Cu ion have been complemented more recently by studies of electrical conduction through a monolayer of Az in the solid-state, sandwiched between electrodes. As the latter type of measurement does not require involvement of a redox process, it also allows monitoring electronic transport (ETp) via redox-inactive Az-derivatives. Here, results of macroscopic ETp via redox-active and -inactive Az derivatives, i.e., Cu(II) and Cu(I)-Az, apo-Az, Co(II)-Az, Ni(II)-Az, and Zn(II)-Az are reported and compared. It is found that earlier reported temperature independence of ETp via Cu(II)-Az (from 20 K until denaturation) is unique, as ETp via all other derivatives is thermally activated at temperatures >≈200 K. Conduction via Cu(I)-Az shows unexpected temperature dependence >≈200 K, with currents decreasing at positive and increasing at negative bias. Taking all the data together we find a clear compensation effect of Az conduction around the Az denaturation temperature. This compensation can be understood by viewing the Az binding site as an electron trap, unless occupied by Cu(II), as in the native protein, with conduction of the native protein setting the upper transport efficiency limit.

Luminescence quenching by photoinduced charge transfer between metal complexes in peptide nucleic acids

A new scaffold for studying photoinduced charge transfer has been constructed by connecting a [Ru(Bpy)3](2+) donor to a bis(8-hydroxyquinolinate)2 copper [CuQ2] acceptor through a peptide nucleic acid (PNA) bridge. The luminescence of the [Ru(Bpy)3](2+*) donor is quenched by electron transfer to the [CuQ2] acceptor. Photoluminescence studies of these donor-bridge-acceptor systems reveal a dependence of the charge transfer on the length and sequence of the PNA bridge and on the position of the donor and acceptor in the PNA. In cases where the [Ru(Bpy)3](2+) can access the π base stack at the terminus of the duplex, the luminescence decay is described well by a single exponential; but if the donor is sterically hindered from accessing the π base stack of the PNA duplex, a distribution of luminescence lifetimes for the donor [Ru(Bpy)3](2+*) is observed. Molecular dynamics simulations are used to explore the donor-PNA-acceptor structure and the resulting conformational distribution provides a possible explanation for the distribution of electron transfer rates.

Biological charge transfer via flickering resonance

Biological electron-transfer (ET) reactions are typically described in the framework of coherent two-state electron tunneling or multistep hopping. However, these ET reactions may involve multiple redox cofactors in van der Waals contact with each other and with vibronic broadenings on the same scale as the energy gaps among the species. In this regime, fluctuations of the molecular structures and of the medium can produce transient energy level matching among multiple electronic states. This transient degeneracy, or flickering electronic resonance among states, is found to support coherent (ballistic) charge transfer. Importantly, ET rates arising from a flickering resonance (FR) mechanism will decay exponentially with distance because the probability of energy matching multiple states is multiplicative. The distance dependence of FR transport thus mimics the exponential decay that is usually associated with electron tunneling, although FR transport involves real carrier population on the bridge and is not a tunneling phenomenon. Likely candidates for FR transport are macromolecules with ET groups in van der Waals contact: DNA, bacterial nanowires, multiheme proteins, strongly coupled porphyrin arrays, and proteins with closely packed redox-active residues. The theory developed here is used to analyze DNA charge-transfer kinetics, and we find that charge-transfer distances up to three to four bases may be accounted for with this mechanism. Thus, the observed rapid (exponential) distance dependence of DNA ET rates over distances of ≲ 15 Å does not necessarily prove a tunneling mechanism.

A three-step kinetic model for electrochemical charge transfer in the hopping regime

Single-step nonadiabatic electron tunneling models are widely used to analyze electrochemical rates through self-assembled monolayer films (SAMs). For some systems, such as nucleic acids, long-range charge transfer can occur in a "hopping" regime that involves multiple charge transfer events and intermediate states. This report describes a three-step kinetic scheme to model charge transfer in this regime. Some of the features of the three-step model are probed experimentally by changing the chemical composition of the SAM. This work uses the three-step model and a temperature dependence of the charge transfer rate to extract the charge injection barrier for a SAM composed of a 10-mer peptide nucleic acid that operates in the hopping regime.

Breaking the simple proportionality between molecular conductances and charge transfer rates

A theoretical framework is presented to describe and to understand the observed relationship between molecular conductances and charge transfer rates across molecular bridges as a function of length, structure, and charge transfer mechanism. The approach uses a reduced density matrix formulation with a phenomenological treatment of system–bath couplings to describe charge transfer kinetics and a Green's function based Landauer–Buttiker method to describe steady-state currents. Application of the framework is independent of the transport regime and includes bath-induced decoherence effects. This model shows that the relationship between molecular conductances and charge transfer rates follows a power-law. The nonlinear rate–conductance relationship is shown to arise from differences in the charge transport barrier heights and from differences in environmental decoherence rates for the two experiments. This model explains otherwise puzzling correlations between molecular conductances and electrochemical kinetics.

Hole transfer kinetics of DNA

Not long after the discovery of the double-helical structure of DNA in 1952, researchers proposed that charge transfer along a one-dimensional π-array of nucleobases might be possible. At the end of the 1990s researchers discovered that a positive charge (a hole) generated in DNA migrates more than 200 Å along the structure, a discovery that ignited interest in the charge-transfer process in DNA. As a result, DNA became an interesting potential bottom-up material for constructing nanoelectronic sensors and devices because DNA can form various complex two-dimensional and three-dimensional structures, such as smiley faces and cubes. From the fundamental aspects of the hole transfer process, DNA is one of the most well-studied organic molecules with many reports on the synthesis of artificial nucleobase analogues. Thus, DNA offers a unique system to study how factors such as the HOMO energy and molecular flexibility affect hole transfer kinetics. Understanding the hole transfer mechanism requires a discussion of the hole transfer rate constants (kHT). This Account reviews the kHT values determined by our group and by Lewis and Wasielewski's group, obtained by a combination of the synthesis of modified DNA and time-resolved spectroscopy. DNA consists of G/C and A/T base pairs; the HOMO localizes on the purine bases G and A, and G has a lower oxidation potential and a higher energy HOMO. Typically, long-range hole transfer proceeded via sequential hole transfer between G/C's. The kinetics of this process in DNA sequences, including those with mismatches, is reproducible via kinetic modeling using the determined kHT for each hole transfer step between G/C's. We also determined the distance dependence parameter (β), which describes the steepness of the exponential decrease of kHT. Because of this value, >0.6 Å(-1) for hole transfer in DNA, DNA itself does not serve as a molecular wire. Interestingly, hole transfer proceeded exceptionally fast for some sequences in which G/C's are located close to each other, an observation that we cannot explain by a simple sequential hole transfer between G/C's but rather through hole delocalization over the nucleobases. To further investigate and refine the factors that affect kHT, we examined various artificial nucleobases. We clearly demonstrated that kHT depends strongly on the HOMO energy gap between the bases (ΔHOMO), and that kHT can be increased with decreasing ΔHOMO. We reduced ΔHOMO between the two type of base pairs by replacing adenines (A's) with deazaadenines ((z)A's) or diaminopurines (D's) and showed that the hole transfer rate through the G/C and A/T mix sequence increased by more than 3 orders of magnitude. We also investigated how DNA flexibility affects kHT. Locked nucleic acid (LNA) modification, which makes DNA more rigid, lowered kHT by more than 2 orders of magnitude. On the other hand, 5-Me-2'-deoxyzebularine (B) modification, which increases DNA flexibility, increased kHT by more than 1 order of magnitude. These new insights in hole transfer kinetics obtained from modified DNAs may aid in the design of new molecular-scale conducting materials.

The single-molecule conductance and electrochemical electron-transfer rate are related by a power law

This study examines quantitative correlations between molecular conductances and standard electrochemical rate constants for alkanes and peptide nucleic acid (PNA) oligomers as a function of the length, structure, and charge transport mechanism. The experimental data show a power-law relationship between conductances and charge transfer rates within a given class of molecules with the same bridge chemistry, and a lack of correlation when a more diverse group of molecules is compared, in contrast with some theoretical predictions. Surprisingly, the PNA duplexes exhibit the lowest charge-transfer rates and the highest molecular conductances. The nonlinear rate-conductance relationships for structures with the same bridging chemistries are attributed to differences in the charge-mediation characteristics of the molecular bridge, energy barrier shifts and electronic dephasing, in the two different experimental settings.

Hole transfer in LNA and 5-Me-2'-deoxyzebularine-modified DNA

We report the measurement of hole-transfer rate constants (k(ht)) in locked nucleic acid (LNA) and 5-Me-2'-deoxyzebularine (B)-modified DNA. LNA modification, which makes DNA more rigid, caused a decrease of more than 2 orders of magnitude in k(ht), whereas B modification, which increases DNA flexibility, increased k(ht) by more than 20-fold. The present results clearly showed that hole-transfer efficiency in DNA can be increased by increasing DNA flexibility.

HOMO energy gap dependence of hole-transfer kinetics in DNA

DNA consists of two type of base-pairs, G-C and A-T, in which the highest occupied molecular orbital (HOMO) localizes on the purine bases G and A. While the hole transfer through consecutive Gs or As occurs faster than 10(9) s(-1), a significant drop in the hole transfer rate was observed for G-C and A-T mixed random sequences. In this study, by using various natural and artificial nucleobases having different HOMO levels, the effect of the HOMO-energy gap between bases (Δ(HOMO)) on the hole-transfer kinetics in DNA was investigated. The results demonstrated that the hole transfer rate can be increased by decreasing the Δ(HOMO) and can be finely tuned over 3 orders of magnitude by varying the Δ(HOMO).

Charge transfer through modified peptide nucleic acids

We studied the charge transfer properties of bipyridine-modified peptide nucleic acid (PNA) in the absence and presence of Zn(II). Characterization of the PNA in solution showed that Zn(II) interacts with the bipyridine ligands, but the stability of the duplexes was not affected significantly by the binding of Zn(II). The charge transfer properties of these molecules were examined by electrochemistry for self-assembled monolayers of ferrocene-terminated PNAs and by conductive probe atomic force microscopy for cysteine-terminated PNAs. Both electrochemical and single molecular studies showed that the bipyridine modification and Zn(II) binding do not affect significantly the charge transfer of the PNA duplexes.

Formation process and solvent-dependent structure of a polyproline self-assembled monolayer on a gold surface

The formation process and structure of a self-assembled monolayer (SAM) of lipoic-acid-terminated polyproline on a gold surface in aqueous solution were investigated by several techniques. The amount of polyproline molecules on the gold surface was determined from the area of the reductive desorption peak, and orientation and thickness of the polyproline SAM were determined in situ by attenuated total reflection infrared (ATR-IR) spectroscopy and ellipsometry. The kinetics of the polyproline SAM formation process were discussed on the basis of these results. The in situ IR study confirmed that the conformation of the polyproline SAM was changed by changing the solvent from water to methanol and methanol to water, as is the case for polyproline dissolved in solution.