A DFT study of 2-aminopurine-containing dinucleotides: prediction of stacked conformations with B-DNA structure

Duplex DNA Retains the Conformational Features of Single Strands: Perspectives from MD Simulations and Quantum Chemical Computations

Molecular dynamics simulations and geometry optimizations carried out at the quantum level as well as by quantum mechanical/molecular mechanics methods predict that short, single-stranded DNA oligonucleotides adopt conformations very similar to those observed in crystallographic double-stranded B-DNA, with rise coordinates close to ≈3.3 Å. In agreement with the experimental evidence, the computational results show that DNA single strands rich in adjacent purine nucleobases assume more regular arrangements than poly-thymine. The preliminary results suggest that single-stranded poly-cytosine DNA should also retain a substantial helical order in solution. A comparison of the structures of single and double helices confirms that the B-DNA motif is a favorable arrangement also for single strands. Indeed, the optimal geometry of the complementary single helices is changed to a very small extent in the formation of the duplex.

Fundamental photophysics of isomorphic and expanded fluorescent nucleoside analogues

Fluorescent nucleoside analogues (FNAs) are structurally diverse mimics of the natural essentially non-fluorescent nucleosides which have found numerous applications in probing the structure and dynamics of nucleic acids as well as their interactions with various biomolecules. In order to minimize disturbance in the labelled nucleic acid sequences, the FNA chromophoric groups should resemble the natural nucleobases in size and hydrogen-bonding patterns. Isomorphic and expanded FNAs are the two groups that best meet the criteria of non-perturbing fluorescent labels for DNA and RNA. Significant progress has been made over the past decades in understanding the fundamental photophysics that governs the spectroscopic and environmentally sensitive properties of these FNAs. Herein, we review recent advances in the spectroscopic and computational studies of selected isomorphic and expanded FNAs. We also show how this information can be used as a rational basis to design new FNAs, select appropriate sequences for optimal spectroscopic response and interpret fluorescence data in FNA applications.

Dynamic and static quenching of 2-aminopurine fluorescence by the natural DNA nucleotides in solution

2-aminopurine (2AP) is a responsive fluorescent base analogue that is used widely as a probe of the local molecular environment in DNA. The ability of 2AP to report changes in local conformation and base-stacking interactions arises from the efficient quenching of its fluorescence by the natural DNA bases. However, the mechanism of this inter-base quenching remains imperfectly understood. Two previous studies of the collisional quenching of 2AP by the natural bases, in different buffer solutions, showed that dynamic quenching efficiency depends on the identity of the natural base, but disagreed on the relative quenching efficiencies of the bases. We report a comprehensive investigation of inter-base quenching of 2AP by the natural nucleoside monophosphates (NMPs), replicating the buffer conditions used in the previous studies. Using time-resolved fluorescence measurements to distinguish between dynamic and static quenching, we find that the dynamic quenching rate constants of the different bases show a consistent trend across both buffers, and this is in line with a charge-transfer mechanism. Time-resolved measurements also provide insight into static quenching, revealing formation of 2AP-NMP ground-state complexes in which 2AP displays a very short fluorescence lifetime, comparable to that seen in oligonucleotides. In these complexes, the dependence of the rate of quenching on the partner base also supports a charge-transfer mechanism.

Single-Stranded DNA Oligonucleotides Retain Rise Coordinates Characteristic of Double Helices

The structures of single-stranded DNA oligonucleotides from dimeric to hexameric sequences have been thoroughly investigated. Computations performed at the density functional level of theory including dispersion forces and solvation show that single-stranded helices adopt conformations very close to crystallographic B-DNA, with rise coordinates amounting up to 3.3 Å. Previous results, suggesting that single strands should be shorter than double helices, largely originated from the incompleteness of the adopted basis set. Although sensible deviations with respect to standard B-DNA are predicted, computations indicate that sequences rich in stacked adenines are the most ordered ones, favoring the B-DNA pattern and inducing regular arrangements also on flanking nucleobases. Several structural properties of double helices rich in adenine are indeed already reflected by the corresponding single strands.

Proton Transfer and Tautomerism in 2-Aminopurine-Thymine and Pyrrolocytosine-Guanine Base Pairs

Pyrrolocytosine (PC) and 2-aminopurine (2AP) are fluorescent nucleobase analogues of the DNA nucleobases cytosine and adenine, respectively, and form base pairs with guanine and thymine. Both fluorescent nucleobases are used extensively as probes for local structure in nucleic acids as the fluorescence properties of PC and 2AP are very sensitive to changes such as helix formation, although the reasons for this sensitivity are not clear. To address this question, ab initio calculations have been used to calculate energies, at the MP2 and CIS level, of three different tautomer pairings of PC-G, and two of 2AP-T, which can potentially be interconverted by double proton transfer between the bases. Potential energy curves linking the different tautomer pairs have been calculated. For both PC-G and 2AP-T, the most stable tautomer pair in the electronic ground state is that analogous to the natural C-G and A-T base pair. In the case of 2AP-T, an alternative, stable, tautomer base pair was located in the first electronically excited state; however, it lies higher in energy than the tautomer pair analogous to A-T, making conversion to the alternative form unlikely. In contrast, in the case of PC-G, an alternative tautomer base pair is found to be the most stable form in the first electronically excited state, and this form is accessible following initial excitation from the ground state tautomer pair, thus suggesting an alternative deactivation route via double proton transfer may be possible when PC is involved in hydrogen bonding, such as occurs in helical conformations.

Efficient intersystem crossing in 2-aminopurine riboside probed by femtosecond time-resolved transient vibrational absorption spectroscopy

The photophysical dynamics of 2-aminopurine, a fluorescent analogue of the canonical nucleobase adenine, has been studied by femtosecond transient vibrational absorption spectroscopy.

Comparative Assessment of Different RNA Tetranucleotides from the DFT-D3 and Force Field Perspective

Classical force field (FF) molecular dynamics (MD) simulations of RNA tetranucleotides have substantial problems in reproducing conformer populations indicated by NMR experiments. To provide more information about the possible sources of errors, we performed quantum mechanical (QM, TPSS-D3/def2-TZVP) and molecular mechanics (MM, AMBER parm99bsc0+χ_OL3) calculations of different r(CCCC), r(GACC), and r(UUUU) conformers obtained from explicit solvent MD simulations. Solvent effects in the static QM and MM calculations were mimicked using implicit solvent models (COSMO and Poisson-Boltzmann, respectively). The comparison of QM and MM geometries and energies revealed that the two methodologies provide qualitatively consistent results in most of the cases. Even though we found some differences, these were insufficient to indicate any systematic corrections of the RNA FF terms that could improve the performance of classical MD in simulating tetranucleotides. On the basis of these findings, we inferred that the overpopulation of intercalated conformers in the MD simulations of RNA tetramers, which were not observed experimentally, might be predominantly caused by imbalanced water-solvent and water-water interactions. Apart from the large-scale QM calculations performed to assess the performance of the AMBER FF, a representative spectrum of faster QM methods was tested.