Alexander Balaeff

PNA versus DNA : Effects of Structural Fluctuations on Electronic Structure and Hole-Transport Mechanisms

The effects of structural fluctuations on charge transfer in double-stranded DNA and peptide nucleic acid (PNA) are investigated. A palindromic sequence with two guanine bases that play the roles of hole donor and acceptor, separated by a bridge of two adenine bases, was analyzed using combined molecular dynamics (MD) and quantum-chemical methods. Surprisingly, electronic structure calculations on individual MD snapshots show significant frontier orbital electronic population on the bridge in approximately 10% of the structures. Electron-density delocalization to the bridge is found to be gated by fluctuations of the covalent conjugated bond structure of the aromatic rings of the nucleic bases. It is concluded, therefore, that both thermal hopping and superexchange should contribute significantly to charge transfer even in short DNA/PNA fragments. PNA is found to be more flexible than DNA, and this flexibility is predicted to produce larger rates of charge transfer.

MULTISCALE METHOD FOR SIMULATING PROTEIN-DNA COMPLEXES ∗

We present a multiresolution approach to modeling complexes between protein and DNA that contain looped or coiled DNA. The approach combines a coarse-grained model of the DNA loop, based on the classical theory of elasticity, with an atom level model of proteins and protein- DNA interfaces based on molecular dynamics. The coarse-grained DNA description is controlled through the atom level protein description and vice versa. The feasibility of the resulting multiscale modeling approach is demonstrated for a protein-DNA complex in which a protein called the E. coli lac repressor forces DNA into a 76 base pair loop. The required simulation involves 230,000 atoms, a number that would triple if both protein and DNA loops were described at the atomic level.

Biological charge transfer via flickering resonance.

Significance Electron transport through DNA plays a central role in nucleic acid damage and repair, and it is usually modeled using a carrier tunneling mechanism (at short distances) and a hopping mechanism (at longer distances). We find that fluctuations into transient geometries that bring multiple bases into electronic degeneracy may support band-like transport during the resonance lifetimes over a distance of ≲15 Å, obviating the need to invoke electron tunneling at short distances. This line of research may help to reveal mechanisms of charge transport in multiheme proteins, bacterial nanowires, and synthetic nanowires, and may also assist in framing the mechanisms of coherent multipigment excitonic transport in light-harvesting proteins. 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.

Evidence for a Near-Resonant Charge Transfer Mechanism for Double-Stranded Peptide Nucleic Acid

We present evidence for a near-resonant mechanism of charge transfer in short peptide nucleic acid (PNA) duplexes obtained through electrochemical, STM break junction (STM-BJ), and computational studies. A seven base pair (7-bp) PNA duplex with the sequence (TA)(3)-(XY)-(TA)(3) was studied, in which XY is a complementary nucleobase pair. The experiments showed that the heterogeneous charge transfer rate constant (k(0)) and the single-molecule conductance (σ) correlate with the oxidation potential of the purine base in the XY base pair. The electrochemical measurements showed that the enhancement of k(0) is independent, within experimental error, of which of the two PNA strands contains the purine base of the XY base pair. 7-bp PNA duplexes with one or two GC base pairs had similar measured k(0) and conductance values. While a simple superexchange model, previously used to rationalize charge transfer in single stranded PNA (Paul et al. J. Am. Chem. Soc. 2009, 131, 6498-6507), describes some of the experimental observations, the model does not explain the absence of an enhancement in the experimental k(0) and σ upon increasing the G content in the duplexes from one to two. Moreover, the superexchange model is not consistent with other studies (Paul et al. J. Phys. Chem. B 2010, 114, 14140), that showed a hopping charge transport mechanism is likely important for PNA duplexes longer than seven base pairs. A quantitative computational analysis shows that a near-resonant charge transfer regime, wherein a mix of superexchange and hopping mechanisms are expected to coexist, can rationalize all of the experimental results.

Role of Nucleobase Energetics and Nucleobase Interactions in Single-Stranded Peptide Nucleic Acid Charge Transfer

Self-assembled monolayers of single-stranded (ss) peptide nucleic acids (PNAs) containing seven nucleotides (TTTXTTT), a C-terminus cysteine, and an N-terminus ferrocene redox group were formed on gold electrodes. The PNA monomer group (X) was selected to be either cytosine (C), thymine (T), adenine (A), guanine (G), or a methyl group (Bk). The charge transfer rate through the oligonucleotides was found to correlate with the oxidation potential of X. Kinetic measurements and computational studies of the ss-PNA fragments show that a nucleobase mediated charge transport mechanism in the deep tunneling superexchange regime can explain the observed dependence of the kinetics of charge transfer on the PNA sequence. Theoretical analysis suggests that the charge transport is dominantly hole-mediated and takes place through the filled bridge orbitals. The strongest contribution to conductance comes from the highest filled orbitals (HOMO, HOMO-1, and HOMO-2) of individual bases, with a rapid drop off in contributions from lower lying filled orbitals. Our studies further suggest that the linear correlation observed between the experimental charge transfer rates and the oxidation potential of base X arises from weak average interbase couplings and similar stacking geometries for the four TTTXTTT systems.

Solution structure of a peptide nucleic acid duplex from NMR data: features and limitations.

This paper describes the results of a 1D and 2D NMR spectroscopy study of a palindromic 8-base pair PNA duplex GGCATGCC in H2O and H2O-D2O solutions. The (1)H NMR peaks have been assigned for most of the protons of the six central base pairs, as well as for several amide protons of the backbone. The resulting 36 interbase and base-backbone distance restraints were used together with Watson-Crick restraints to generate the PNA duplex structure in the course of 10 independent simulated annealing runs followed by restrained molecular dynamics (MD) simulations in explicit water. The resulting PNA structures correspond to a P-type helix with helical parameters close to those observed in the crystal structures of PNA. Based on the current limited number of restraints obtained from NMR spectra, alternative structures obtained by MD from starting PNA models based on DNA cannot be ruled out and are also discussed.

STRUCTURE PREDICTION OF A COMPLEX BETWEEN THE CHROMOSOMAL PROTEIN HMG-D AND DNA

Non‐histone chromosomal proteins are an important part of nuclear structure and function due to their ability to interact with DNA to form and modulate chromatin structure and regulate gene expression. However, the understanding of the function of chromosomal proteins at the molecular level has been hampered by the lack of structures of chromosomal protein–DNA complexes. We have carried out a molecular dynamics modeling study to provide insight into the mode of DNA binding to the chromosomal HMG‐domain protein, HMG‐D. Three models of a complex of HMG‐D bound to DNA were derived through docking the protein to two different DNA fragments of known structure. Molecular dynamics simulations of the complexes provided data indicating the most favorable model. This model was further refined by molecular dynamics simulation and extensively analyzed. The structure of the corresponding HMG‐D‐DNA complex exhibits many features seen in the NMR structures of the sequence‐specific HMG‐domain‐DNA complexes, lymphoid enhancer factor 1 (LEF‐1) and testis determining factor (SRY). The model reveals differences from these known structures that suggest how chromosomal proteins bind to many different DNA sequences with comparable affinity. Proteins 30:113–135, 1998.

Effect of Backbone Flexibility on Charge Transfer Rates in Peptide Nucleic Acid Duplexes

Charge transfer (CT) properties are compared between peptide nucleic acid structures with an aminoethylglycine backbone (aeg-PNA) and those with a γ-methylated backbone (γ-PNA). The common aeg-PNA is an achiral molecule with a flexible structure, whereas γ-PNA is a chiral molecule with a significantly more rigid structure than aeg-PNA. Electrochemical measurements show that the CT rate constant through an aeg-PNA bridging unit is twice the CT rate constant through a γ-PNA bridging unit. Theoretical calculations of PNA electronic properties, which are based on a molecular dynamics structural ensemble, reveal that the difference in the CT rate constant results from the difference in the extent of backbone fluctuations of aeg- and γ-PNA. In particular, fluctuations of the backbone affect the local electric field that broadens the energy levels of the PNA nucleobases. The greater flexibility of the aeg-PNA gives rise to more broadening, and a more frequent appearance of high-CT rate conformations than in γ-PNA.

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.

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.Guanine-quadruplex motifs form multi-ended DNA junctions that transport electrical charges with minimal losses.