Developing a computational model that accurately reproduces the structural features of a dinucleoside monophosphate unit within B-DNA

How proton transfer affects the helical parameters in DNA:DNA microhelices

Proton transfer reactions are a widespread phenomenon in many areas of the life sciences and it is one of the origins of the spontaneous point mutations during DNA replication. Because of its importance, many studies have been reported on these reactions. However, the present work is the first one focused on the structural geometrical changes by double proton transfer (DPT). Thus, different Watson-Crick (WC) pairs were optimized first in a simple model with one nucleoside base pair, and in a microhelix form with three nucleoside base pairs. The canonical and few tautomeric forms were considered in DNA:DNA microhelices with A-type and B-type helical forms. The stability of these structures and how the DPT process affects the main geometrical parameters was analyzed, in particular the deformation of the helical parameters. The M06-2X DFT method was used for this purpose. The purine/pyrimidine ring in the keto form appears easier to be deformed than when it is in the enol form. The weaker WC base pair formed with mixed microhelices than with nucleobases alone and the significant deformation of the helical and backbone parameters with the DPT appears to complicate this process in microhelices.Communicated by Ramaswamy H. Sarma.

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.

Base pairs with 4-amino-3-nitrobenzonitrile: comparison with the natural WC pairs. Dimer and tetramer forms, Infrared and Raman spectra, and several proposed antiviral modified nucleosides

Base pairs of 4-amino-3-nitrobenzonitrile (4A-3NBN) molecule with uracil, thymine and cytosine nucleobases were optimized and compared to natural Watson-Crick (WC) pairs. The slightly greater flexibility of the -NO₂ group of 4A-3NBN than the N3-H group of the natural nucleobases together with a noticeable higher dipole moment of its pairs can facilitate disruption of the DNA/RNA helix formation. Several new mutagenic modified nucleosides with 4A-3NBN and 3-amino-2-nitrobenzonitrile (3A-2NBN) were proposed as antiviral prodrugs and their base pairs optimized. The special characteristics of these prodrugs appear appropriated for their clinical use. The counterpoise (CP) corrected interaction energies of the base pairs were calculated and compared to the natural ones. The M06-2X DFT method was used for this purpose. The molecular structure of 4A-3NBN was analyzed in detail and the crystal unit cell was simulated by a tetramer form and eight dimer forms. The performance of the B3LYP, X3LYP and M06-2X methods was tested on the vibrational wavenumbers in the monomer, dimer and tetramer forms of 4A-3NBN. The observed IR and Raman bands were assigned according to the optimum dimer II form determined by B3LYP and by the tetramer form calculated by M06-2X, which is the expected unit cell that forms the crystal net. The two best scaling procedures were used.Communicated by Ramaswamy H. Sarma.

Improving structural stability and anticoagulant activity of a thrombin binding aptamer by aromatic modifications

The thrombin binding aptamer (TBA) is a 15-mer DNA oligonucleotide (5'-GGTTGGTGTGGTTGG-3'), that can form a stable intramolecular antiparallel chair-like G-quadruplex structure. This aptamer shows anticoagulant properties by interacting with one of the two anion binding sites of thrombin, namely the fibrinogen-recognition exosite. Here, we demonstrate that terminal modification of TBA with aromatic fragments such as coumarin, pyrene and perylene diimide (PDI), improves the G-quadruplex stability. The large aromatic surface of these dyes can π-π stack to the G-quadruplex or to each other, thereby stabilizing the aptamer. With respect to the original TBA, monoPDI-functionalized TBA exhibited the most remarkable improvement in melting temperature (ΔT m ≈ +18 °C) and displayed enhanced anticoagulant activity.

Theoretical Investigation of the 4,5-Dibromorodamine Methyl Ester (TH9402) Photosensitizer Used in Photodynamic Therapy: Photophysics, Reactions in the Excited State, and Interactions with DNA

Photosensitizer (PS) molecules play a critical role in photodynamic therapy of cancer and the understanding of the molecular mechanism involved in the photophysics of these compounds, and their reactions in the excited state are, therefore, of great interest for the development of this technique. In this article, the photophysics of the cationic PS 4,5-dibromorodamine methyl ester (TH9402), its electron- and energy-transfer reactions in the excited triplet state, with molecular oxygen, nitric oxide, guanosine-5'-monophosphate (GMP), and guanine, and the interaction with DNA were evaluated. Time-dependent density functional theory calculations at the TPSSh/Def2-TZVP//B3LYP/Def2-TZVP level of theory in water solution reveals that the PS has a bright S₁ state 2.33 eV above the ground state that produces a fluorescent rate constant of 5.40 × 10⁷ s^-1, calculated using Fermi's golden rule within a path integral formalism. Once excited to the bright state, the main intersystem crossing (ISC) channel involves the coupling with the T₂ state just below S₁ (S₁ → T₂ → T₁) with an overall ISC rate constant of 10.1 × 10⁷ s^-1, in good agreement with the experimental data. Excited-state reaction thermodynamics, computed at the M06-2X/Def2-TZVP//B3LYP/Def2-TZVP level of theory in water, showed that from all the excited-state electron-transfer reactions studied, only the transfer from GMP to the PS is thermodynamically favorable, independent of the protonation state of guanosine, which indicates a possible DNA photo-oxidation mechanism for the PS. Triplet-triplet energy-transfer reactions from TH9402 to molecular oxygen, producing reactive singlet oxygen, and to the deprotonated guanosine, producing ³GMP^2-, are also thermodynamically favorable, with ΔG = -2.0 and -24.0 kcal//mol, respectively. However, the energy transfer to the monoprotonated guanosine is not favorable, (ΔG = 36.1), suggesting that in the DNA double-strand environment, this energy-transfer process may not be observed. The results show that the PS can act through electron transfer and triplet-triplet energy-transfer reactions involved in mechanism types I and II in photodynamic therapy. Interactions of TH9402 with the d(AGACGTCT)₂ octanucleotide revealed that the PS can intercalate between the d(GpC)-d(CpG) base pairs in three different orientations and, upon intercalation, the π → π* transition of the PS shows a bathochromic shift up to 90 nm and up to 60% decrease in intensity. Interactions through groove binding showed a smaller bathochromic shift of 52.2 nm and a 56% decrease in intensity of the main transition band.

Theoretical studies on the purine radical induced purine-purine type intrastrand cross-links

At the density functional theory (DFT) level, addition reactions between the guanine-8-yl radical and its 3'/5' neighboring purine deoxynucleosides forming the purine-purine type intrastrand cross-links were studied. It is found that addition of the guanine-8-yl radical to the C8 site of its 5' neighboring deoxyguanosine or deoxyadenosine is a two-step reaction consisting of a structurally relatively unfavourable conformational transformation step, while the corresponding 3' C8 addition is straightforward and kinetically more efficient. The 3' C8 preference of the guanine-8-yl radical additions indicates the existence of an obvious sequence effect, which is completely opposite to that observed in the formation of pyrimidine radicals induced DNA intrastrand cross-links. The detrimental effects from steric hindrance and stabilizing weak interactions make these addition reactions markedly suppressed in double stranded DNA. This work broadens our knowledge about the possible types of DNA intrastrand cross-links.

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.

DNA intrastrand cross-links induced by the purine-type deoxyguanosine-8-yl radical: a DFT study

Currently, all known DNA intrastrand cross-links are found to be induced by pyrimidine-type radicals; however, whether or not purine-type radicals are able to cause DNA intrastrand cross-links remains unclear. In the present study, probable additions of the highly reactive deoxyguanosine-8-yl radical to its 3'/5' neighboring pyrimidine nucleotides in four model compounds, 5'-G˙T-3', 5'-TG˙-3', 5'-G˙C-3', and 5'-CG˙-3', were studied using density functional theory (DFT) methods. In single-stranded DNA, the deoxyguanosine-8-yl radical is preferred to efficiently attack the C5 site of its 3' neighboring deoxythymidine or deoxycytidine, forming the G[8-5]T or G[8-5]C intrastrand cross-link rather than the C6 site forming the G[8-6]T or G[8-6]C intrastrand cross-link. The four corresponding sequence isomers, namely T[5-8]G, T[6-8]G, C[5-8]G, and C[6-8]G, formed by additions of deoxyguanosine-8-yl radical to its 5' neighboring pyrimidine nucleotides are predicted to be formed inefficiently. In double-stranded DNA, considering the detrimental effects of stabilizing weak interactions on related structural adjustments required in each addition reaction path, relatively lower reaction yields are suggested for the G[8-5]T and G[8-5]C intrastrand cross-links, while the formation of the other six intrastrand cross-links becomes quite difficult. All calculations definitely demonstrate that, in addition to pyrimidine-type radicals, the purine-type deoxyguanosine-8-yl radical is able to attack its 3'/5' neighboring pyrimidine nucleotides forming several DNA intrastrand cross-links.

Formation of pyrimidine-pyrimidine type DNA intrastrand cross-links: a theoretical verification

Pyrimidine-type radicals have been demonstrated to be able to attack their 3' or 5' neighboring purine nucleotides forming diverse DNA intrastrand cross-links, but whether or not these radicals can attack their surrounding pyrimidine nucleotides forming pyrimidine-pyrimidine type DNA intrastrand cross-links remains unclear. To resolve this question, probable additions of the uracil-5-methyl (˙U_CH2) radical to the C₅[double bond, length as m-dash]C₆ double bond of its 3'/5' neighboring pyrimidine nucleotides in the four models, 5'-T(˙U_CH2)-3', 5'-C(˙U_CH2)-3', 5'-(˙U_CH2)T-3', and 5'-(˙U_CH2)C-3', are explored in the present work employing density functional theory (DFT) methods. The C₆ site of its 5' neighboring thymidine is the preferred target for ˙U_CH2 radical addition, while additions of the ˙U_CH2 radical to the C₆ and C₅ sites of its 5' neighboring deoxycytidine are found to be competitive reactions. The ˙U_CH2 radical can react with both the C₆ and C₅ sites of its 3' neighboring pyrimidine nucleotides, but the efficiencies of these reactions are predicted to be much lower than those of the corresponding addition reactions to its 5' neighboring pyrimidine nucleotides, indicating the existence of an obvious sequence effect. All the addition products could be finally transformed into closed-shell intrastrand cross-links, the molecular masses of which are found to be exactly the same as certain MS values determined in a recent study of an X-irradiated deoxygenated aqueous solution of calf thymus DNA. The present study thus not only definitely corroborates the fact that the reactive ˙U_CH2 radical can attack its 3'/5' neighboring pyrimidine nucleotides forming several pyrimidine-pyrimidine type DNA intrastrand cross-links, but also provides a plausible explanation for the identities of these structurally unknown intrastrand cross-links.

Interaction of osmium(ii) redox probes with DNA: insights from theory

In the course of developing ultrasensitive and quantitative electrochemical point-of-care analytical tools for genetic detection of infectious diseases, osmium(ii) metallointercalators were revealed to be suitable and efficient redox probes to monitor the in vitro DNA amplification [Defever etal, Anal. Chem., 2011, 83, 1815-1821]. In this work, we thus propose a complete computational protocol in order to evaluate the affinity between Os(ii) complexes with double-stranded DNA. This protocol is based on molecular dynamics, with the parametrization of the GAFF force field for the Os(ii) complexes presenting an octahedral environment with polypyridine ligands, and QM/QM' calculations to evaluate the binding energy. For three Os(ii) probes and different binding sites, molecular dynamics simulations and interaction energies calculated at the QM/QM' level are successively discussed and compared to experimental data in order to identify the most stable binding sites. The computational protocol we propose should then be used to design more efficient Os(ii) metallointercalators.

Mechanism studies of addition reactions between the pyrimidine type radicals and their 3′/5′ neighboring deoxyguanosines

To clarify the biologically significant sequence effect existing in the formation of the pyrimidine-type radicals induced DNA intrastrand cross-links, addition mechanisms between the uridine-5-methyl (˙U_CH2), 6-hydroxy-5,6-dihydrothymidine-5-yl (˙T_6OH), and 6-hydroxy-5,6-dihydrocytidine-5-yl (˙C_6OH) radicals and their 3′/5′ neighboring deoxyguanosines (dG) are explored in the present study employing the model 5′-G(˙U_CH2)-3′, 5′-(˙U_CH2)G-3′, 5′-G(˙T_6OH)-3′, 5′-(˙T_6OH)G-3′, 5′-G(˙C_6OH)-3′, and 5′-(˙C_6OH)G-3′ sequences. It is found that the 5′ G/C₈ additions of the three radicals are all simple direct one-step reactions inducing only relatively small structural changes, while a conformational adjustment involving orientation transitions of both nucleobase moieties and twisting of the DNA backbone is indispensable for each 3′ G/C₈ addition. Furthermore, markedly positive reaction free energy requirements are estimated for these conformational transformations making the 3′ G/C₈ additions of the three radicals thermodynamically much more unfavorable than the corresponding 5′ G/C₈ additions. Such essential conformational adjustments along the 3′ G/C₈ addition paths that structurally greatly influence the local DNA structures and thermodynamically substantially reduce the addition efficiencies may be the reasons responsible for the differences in the formation yields and biological consequences of the pyrimidine-type radicals induced DNA intrastrand cross-link lesions.

For each radical, the 5′ G/C₈ addition is a simple direct one-step reaction, while a structurally significant and thermodynamically markedly unfavorable conformational adjustment is indispensable for the 3′ G/C₈ addition.

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

The fluorescence properties of dinucleotides incorporating 2-aminopurine (2AP) suggest that the simplest oligonucleotides adopt conformations similar to those found in duplex DNA. However, there is a lack of structural data for these systems. We report a density functional theory (DFT) study of the structures of 2AP-containing dinucleotides (deoxydinucleoside monophosphates), including full geometry optimisation of the sugar-phosphate backbone. Our DFT calculations employ the M06-2X functional for reliable treatment of dispersion interactions and include implicit aqueous solvation. Dinucleotides with 2AP in the 5'-position and each of the natural bases in the 3'-position are examined, together with the analogous 5'-adenine-containing systems. Computed structures are compared in detail with typical B-DNA base-step parameters, backbone torsional angles and sugar pucker, derived from crystallographic data. We find that 2AP-containing dinucleotides adopt structures that closely conform to B-DNA in all characteristic parameters. The structures of 2AP-containing dinucleotides closely resemble those of their adenine-containing counterparts, demonstrating the fidelity of 2AP as a mimic of the natural base. As a first step towards exploring the conformational heterogeneity of dinucleotides, we also characterise an imperfectly stacked conformation and one in which the bases are completely unstacked.

Role of 6-Mercaptopurine in the potential therapeutic targets DNA base pairs and G-quadruplex DNA: insights from quantum chemical and molecular dynamics simulations

The theoretical studies on DNA with the anticancer drug 6-Mercaptopurine (6-MP) are investigated using theoretical methods to shed light on drug designing. Among the DNA base pairs considered, 6-MP is stacked with GC with the highest interaction energy of -46.19 kcal/mol. Structural parameters revealed that structure of the DNA base pairs is deviated from the planarity of the equilibrium position due to the formation of hydrogen bonds and stacking interactions with 6-MP. These deviations are verified through the systematic comparison between X-H bond contraction and elongation and the associated blue shift and red shift values by both NBO analysis and vibrational analysis. Bent's rule is verified for the C-H bond contraction in the 6-MP interacted base pairs. The AIM results disclose that the higher values of electron density (ρ) and Laplacian of electron density (∇²ρ) indicate the increased overlap between the orbitals that represent the strong interaction and positive values of the total electron density show the closed-shell interaction. The relative sensitivity of the chemical shift values for the DNA base pairs with 6-MP is investigated to confirm the hydrogen bond strength. Molecular dynamics simulation studies of G-quadruplex DNA d(TGGGGT)₄ with 6-MP revealed that the incorporation of 6-MP appears to cause local distortions and destabilize the G-quadruplex DNA.

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.

π-Cooperativity effect on the base stacking interactions in DNA: is there a novel stabilization factor coupled with base pairing H-bonds?

The results from absolutely localized molecular orbital (ALMO)-energy decomposition analysis (EDA) and ALMO-charge transfer analysis (CTA) at M06-2X/cc-pVTZ level reveal that double-proton transfer (DPT) reactions through base pairing H-bonds have nonignorable effects on the stacking energies of dinucleotide steps, which introduces us to a novel stabilization (or destabilization) factor in the DNA duplex. Thus, intra- and inter-strand base stacking interactions are coalesced with each other mediated by H-bridged quasirings between base pairs. Changes in stacking energies of dinucleotide steps depending on the positions of H atoms are due to variations in local aromaticities of individual nucleobases, manifesting π-cooperativity effects. CT analyses show that dispersion forces in dinucleotide steps can lead to radical changes in the redox properties of nucleobases, in particular those of adenine and guanine stacked dimers in a strand. Besides Watson-Crick rules, novel base pairing rules were propounded by considering CT results. According to these, additional base pairing through π-stacks of nucleobases in dinucleotide steps does not cause any intrinsic oxidative damage to the associated nucleobases throughout DPT.

Quantum chemical benchmark study on 46 RNA backbone families using a dinucleotide unit

We have created a benchmark set of quantum chemical structure-energy data denoted as UpU46, which consists of 46 uracil dinucleotides (UpU), representing all known 46 RNA backbone conformational families. Penalty-function-based restrained optimizations with COSMO TPSS-D3/def2-TZVP ensure a balance between keeping the target conformation and geometry relaxation. The backbone geometries are close to the clustering-means of their respective RNA bioinformatics family classification. High-level wave function methods (DLPNO-CCSD(T) as reference) and a wide-range of dispersion-corrected or inclusive DFT methods (DFT-D3, VV10, LC-BOP-LRD, M06-2X, M11, and more) are used to evaluate the conformational energies. The results are compared to the Amber RNA bsc0χOL3 force field. Most dispersion-corrected DFT methods surpass the Amber force field significantly in accuracy and yield mean absolute deviations (MADs) for relative conformational energies of ∼0.4-0.6 kcal/mol. Double-hybrid density functionals represent the most accurate class of density functionals. Low-cost quantum chemical methods such as PM6-D3H+, HF-3c, DFTB3-D3, as well as small basis set calculations corrected for basis set superposition errors (BSSEs) by the gCP procedure are also tested. Unfortunately, the presently available low-cost methods are struggling to describe the UpU conformational energies with satisfactory accuracy. The UpU46 benchmark is an ideal test for benchmarking and development of fast methods to describe nucleic acids, including force fields.

A thiazole coumarin (TC) turn-on fluorescence probe for AT-base pair detection and multipurpose applications in different biological systems

Sequence-specific recognition of DNA by small turn-on fluorescence probes is a promising tool for bioimaging, bioanalytical and biomedical applications. Here, the authors report a novel cell-permeable and red fluorescent hemicyanine-based thiazole coumarin (TC) probe for DNA recognition, nuclear staining and cell cycle analysis. TC exhibited strong fluorescence enhancement in the presence of DNA containing AT-base pairs, but did not fluoresce with GC sequences, single-stranded DNA, RNA and proteins. The fluorescence staining of HeLa S3 and HEK 293 cells by TC followed by DNase and RNase digestion studies depicted the selective staining of DNA in the nucleus over the cytoplasmic region. Fluorescence-activated cell sorting (FACS) analysis by flow cytometry demonstrated the potential application of TC in cell cycle analysis in HEK 293 cells. Metaphase chromosome and malaria parasite DNA imaging studies further confirmed the in vivo diagnostic and therapeutic applications of probe TC. Probe TC may find multiple applications in fluorescence spectroscopy, diagnostics, bioimaging and molecular and cell biology.

B-DNA characteristics are preserved in double stranded d(A)3·d(T)3 and d(G)3·d(C)3 mini-helixes: conclusions from DFT/M06-2X study

We report the results of the first comprehensive DFT study on the d(A)3·d(T)3 and d(G)3·d(C)3 nucleic acid duplexes. The ability of mini-helixes to preserve the conformation of B-DNA in the gas phase and under the influence of such factors as: solvent, uncompensated charge, and counter-ions was evaluated using M06-2X functional with 6-31G(d,p) basis set. The accuracy of the models was ascertained based on their ability to reproduce key structural features of natural B-DNA. Analysis of the helicity suggests that the helical conformations adopt geometrical parameters which are close to those of the B-DNA form. The torsion angles fall somewhere between the values observed for BI/BII conformational classes. The comparative analysis of parameters of isolated Watson-Crick base pairs versus B-DNA-like conformations indicates the same tendency of base-pair polarization and hydration. Specifically, effects of polarization of nucleobases in continuum type dielectric medium mimicking water are stronger than those caused by the presence of backbone. Polar environment as well as the presence of counterions stabilizes duplexes, facilitating helix formation. Substantial conformational changes of nucleotides upon duplex formation decrease the binding energy. In spite of structural and energetic changes, the placement of a mini-helix into the gas phase does not lead to significant disruption of the structure. On the contrary, the duplex preserves its helicity and the strands remain bound.

Structural and magnetic properties of a variety of transition metal incorporated DNA double helices

By using density functional theory calculations, the structural, energetic, magnetic, and optical properties for a variety of transition metal (M = Mn, Fe, Co, Ni and Cu) ions incorporated modified-DNA (M-DNA) double helices has been investigated. The DNA is modified with either hydroxypyridone (H) or bis(salicylaldehyde)ethylenediamine (S-en) metalated bases. We find the formation of extended M-O network leading to the ferromagnetic interactions for the case of H-DNA for all the metal ions. More ordered stacking arrangement was found for S-en-DNA. We calculate the exchange coupling constant (J) considering Heisenberg Hamiltonian for quantitative description of magnetic interactions. The ferromagnetic and antiferromagnetic interactions are obtained by varying different transition metal ions. The extent of the magnetic interaction depends on the number of transition metal ions. Optical profiles show peaks below 2 eV, a clear signature of spin-spin coupling.

Energies and 2'-Hydroxyl Group Orientations of RNA Backbone Conformations. Benchmark CCSD(T)/CBS Database, Electronic Analysis, and Assessment of DFT Methods and MD Simulations

Sugar-phosphate backbone is an electronically complex molecular segment imparting RNA molecules high flexibility and architectonic heterogeneity necessary for their biological functions. The structural variability of RNA molecules is amplified by the presence of the 2'-hydroxyl group, capable of forming multitude of intra- and intermolecular interactions. Bioinformatics studies based on X-ray structure database revealed that RNA backbone samples at least 46 substates known as rotameric families. The present study provides a comprehensive analysis of RNA backbone conformational preferences and 2'-hydroxyl group orientations. First, we create a benchmark database of estimated CCSD(T)/CBS relative energies of all rotameric families and test performance of dispersion-corrected DFT-D3 methods and molecular mechanics in vacuum and in continuum solvent. The performance of the DFT-D3 methods is in general quite satisfactory. The B-LYP-D3 method provides the best trade-off between accuracy and computational demands. B3-LYP-D3 slightly outperforms the new PW6B95-D3 and MPW1B95-D3 and is the second most accurate density functional of the study. The best agreement with CCSD(T)/CBS is provided by DSD-B-LYP-D3 double-hybrid functional, although its large-scale applications may be limited by high computational costs. Molecular mechanics does not reproduce the fine energy differences between the RNA backbone substates. We also demonstrate that the differences in the magnitude of the hyperconjugation effect do not correlate with the energy ranking of the backbone conformations. Further, we investigated the 2'-hydroxyl group orientation preferences. For all families, we conducted a QM and MM hydroxyl group rigid scan in gas phase and solvent. We then carried out set of explicit solvent MD simulations of folded RNAs and analyze 2'-hydroxyl group orientations of different backbone families in MD. The solvent energy profiles determined primarily by the sugar pucker match well with the distribution data derived from the simulations. The QM and MM energy profiles predict the same 2'-hydroxyl group orientation preferences. Finally, we demonstrate that the high energy of unfavorable and rarely sampled 2'-hydroxyl group orientations can be attributed to clashes between occupied orbitals.

B-DNA Structure and Stability as Function of Nucleic Acid Composition: Dispersion-Corrected DFT Study of Dinucleoside Monophosphate Single and Double Strands

We have computationally investigated the structure and stability of all 16 combinations of two out of the four natural DNA bases A, T, G and C in a di-2′-deoxyribonucleoside-monophosphate model DNA strand as well as in 10 double-strand model complexes thereof, using dispersion-corrected density functional theory (DFT-D). Optimized geometries with B-DNA conformation were obtained through the inclusion of implicit water solvent and, in the DNA models, of sodium counterions, to neutralize the negative charge of the phosphate groups. The results obtained allowed us to compare the relative stability of isomeric single and double strands. Moreover, the energy of the Watson–Crick pairing of complementary single strands to form double-helical structures was calculated. The latter furnished the following increasing stability trend of the double-helix formation energy: d(TpA)₂ <d(CpA)₂ <d(ApT)₂ <d(ApA)₂ <d(GpT)₂ <d(GpA)₂ <d(ApG)₂ <d(CpG)₂ <d(GpG)₂ <d(GpC)₂, where the energy differences between the last four dimers, d(ApG)₂, d(CpG)₂, d(GpG)₂ and d(GpC)₂, is within 4.0 kcal mol⁻¹, and the energy between the most and the least stable isomers is 13.4 kcal mol⁻¹. This trend shows that the formation energy essentially increases with the number of hydrogen bonds per base pair, that is two between A and T and three between G and C. Superimposed on this main trend are more subtle effects that depend on the order in which bases occur within a strand from the 5’- to the 3’-end.

Photophysics and photochemistry of thymine deoxy-dinucleotide in water: a PCM/TD-DFT quantum mechanical study

We here report a fully quantum mechanical study of the main photochemical and photophysical decay routes in aqueous solution of thymine deoxy-dinucleotide (TpT(-) and TpTNa) and of its analogue locked in C3-endo puckering, characterizing five different representative backbone conformers and discussing the chemical physical effects modulating the yield of the different photoproducts. Our approach is based on time-dependent DFT calculations, using the last generation M052X functional, whereas solvent effects are included by means of the polarizable continuum model. Especially when at least one of the sugars adopts C3-endo puckering, a barrierless path on the bright ππ* excitons leads to the S(1)/S(0) crossing region corresponding to the formation of cyclobutane pyrimidine dimer. Charge transfer excited states involving the transfer of an electron from the 5' Thy toward the 3' Thy are involved in the formation of the oxetane intermediate in the path leading to 6-4 pyrimidine pyrimidinone adducts. A non-negligible energy barrier is associated with this latter pathway, which is possible only when one of the two nucleotides adopts C2-endo puckering. Monomer-like decay pathways, involving ππ* or nπ* excited states localized on a single base, are shown to be operative also for loosely stacked bases.

Glycosidic bond cleavage in deoxynucleotides: effects of solvent and the DNA phosphate backbone in the computational model

Density functional theory (B3LYP) was employed to examine the hydrolysis of the canonical 2'-deoxynucleotides in varied environments (gas phase or water) using different computational models for the sugar residue (methyl or phosphate group at C5') and nucleophile (water activated through full or partial proton abstraction). Regardless of the degree of nucleophile activation, our results show that key geometrical parameters along the reaction pathway are notably altered upon direct inclusion of solvent effects in the optimization routine, which leads to significant changes in the reaction energetics and better agreement with experiment. Therefore, despite the wide use of gas-phase calculations in the literature, small model computational work, as well as large-scale enzyme models, that strive to understand nucleotide deglycosylation must adequately describe the environment. Alternatively, although inclusion of the phosphate group at C5' also affects the geometries of important stationary points, the effects cancel to yield unchanged deglycosylation barriers, and therefore smaller computational models can be used to estimate the energy associated with nucleotide deglycosylation, with the 5' phosphate group included if full (geometric) details of the reaction are desired. Hydrogen-bonding interactions with the nucleobase can significantly reduce the barrier to deglycosylation, which supports suggestions that discrete hydrogen-bonding interactions with active-site amino acid residues can play a significant role in enzyme-catalyzed nucleobase excision. Taken together with previous studies, the present work provides vital clues about the components that must be included in future studies of the deglycosylation of isolated noncanonical nucleotides, as well as the corresponding enzyme-catalyzed reactions.

The DNA and RNA sugar-phosphate backbone emerges as the key player. An overview of quantum-chemical, structural biology and simulation studies

Knowledge of geometrical and physico-chemical properties of the sugar-phosphate backbone substantially contributes to the comprehension of the structural dynamics, function and evolution of nucleic acids. We provide a side by side overview of structural biology/bioinformatics, quantum chemical and molecular mechanical/simulation studies of the nucleic acids backbone. We highlight main features, advantages and limitations of these techniques, with a special emphasis given to their synergy. The present status of the research is then illustrated by selected examples which include classification of DNA and RNA backbone families, benchmark structure-energy quantum chemical calculations, parameterization of the dihedral space of simulation force fields, incorporation of arsenate into DNA, sugar-phosphate backbone self-cleavage in small RNA enzymes, and intricate geometries of the backbone in recurrent RNA building blocks. Although not apparent from the current literature showing limited overlaps between the QM, simulation and bioinformatics studies of the nucleic acids backbone, there in fact should be a major cooperative interaction between these three approaches in studies of the sugar-phosphate backbone.

Mechanistic and conformational flexibility of the covalent linkage formed during β-lyase activity on an AP-site: application to hOgg1

The β/δ-lyase activity of bifunctional glycosylases on damaged nucleotides in DNA involves the formation of a covalent linkage between the protein (lysine or N-terminal proline) and DNA (C1' of the damaged nucleotide). In the present study, the conformational and mechanistic flexibility of the cross-link is examined. Repair of 8-oxoguanine damage by hOgg1 is considered as a representative system, and the glycosylase through β-lyase steps are investigated using density functional theory. (PCM/SMD)-M06-2X/6-311+G(2df,2p)//PCM-B3LYP/6-31G(d) energetics were determined for eight unique mechanisms differing in the conformation of the imine linkage (E/Z), the proton (pro-S/R) abstracted during elimination, and whether the ring-opening step is base catalyzed. This initial study used a model system limited to the damaged nucleoside 3'-monophosphate and a model nucleophile to investigate this series of complex reaction steps. The great flexibility exhibited by the linkage and clustered β-elimination energetics indicate sterics will play a large role in predicting the preferred lyase mechanism for a given enzyme. The stationary points identified herein can be overlaid into a protein structure to assist in generating initial guesses for large model systems. By comparing the characterized geometries and enzyme active sites, methods for catalysis of the various chemical steps can be identified, and these possibilities are discussed in detail for hOgg1. Interestingly, the most stable structure on the potential energy surface occurs before elimination of the 3'-phosphate. Hydrolysis of the protein-DNA cross-link at this point would yield an AP-site, which provides support for the recently observed monofunctional activity of hOgg1.

Understanding the Sequence Preference of Recurrent RNA Building Blocks using Quantum Chemistry: The Intrastrand RNA Dinucleotide Platform

Folded RNA molecules are shaped by an astonishing variety of highly conserved noncanonical molecular interactions and backbone topologies. The dinucleotide platform is a widespread recurrent RNA modular building submotif formed by the side-by-side pairing of bases from two consecutive nucleotides within a single strand, with highly specific sequence preferences. This unique arrangement of bases is cemented by an intricate network of noncanonical hydrogen bonds and facilitated by a distinctive backbone topology. The present study investigates the gas-phase intrinsic stabilities of the three most common RNA dinucleotide platforms - 5'-GpU-3', ApA, and UpC - via state-of-the-art quantum-chemical (QM) techniques. The mean stability of base-base interactions decreases with sequence in the order GpU > ApA > UpC. Bader's atoms-in-molecules analysis reveals that the N2(G)…O4(U) hydrogen bond of the GpU platform is stronger than the corresponding hydrogen bonds in the other two platforms. The mixed-pucker sugar-phosphate backbone conformation found in most GpU platforms, in which the 5'-ribose sugar (G) is in the C2'-endo form and the 3'-sugar (U) in the C3'-endo form, is intrinsically more stable than the standard A-RNA backbone arrangement, partially as a result of a favorable O2'…O2P intra-platform interaction. Our results thus validate the hypothesis of Lu et al. (Lu Xiang-Jun, et al. Nucleic Acids Res. 2010, 38, 4868-4876), that the superior stability of GpU platforms is partially mediated by the strong O2'…O2P hydrogen bond. In contrast, ApA and especially UpC platform-compatible backbone conformations are rather diverse and do not display any characteristic structural features. The average stabilities of ApA and UpC derived backbone conformers are also lower than those of GpU platforms. Thus, the observed structural and evolutionary patterns of the dinucleotide platforms can be accounted for, to a large extent, by their intrinsic properties as described by modern QM calculations. In contrast, we show that the dinucleotide platform is not properly described in the course of atomistic explicit-solvent simulations. Our work also gives methodological insights into QM calculations of experimental RNA backbone geometries. Such calculations are inherently complicated by rather large data and refinement uncertainties in the available RNA experimental structures, which often preclude reliable energy computations.

Formation mechanism and structure of a guanine-uracil DNA intrastrand cross-link

The formation and structure of the 5'-G[8-5]U-3' intrastrand cross-link are studied using density functional theory and molecular dynamics simulations due to the potential role of this lesion in the activity of 5-halouracils in antitumor therapies. Upon UV irradiation of 5-halouracil-containing DNA, a guanine radical cation reacts with the uracil radical to form the cross-link, which involves phosphorescence or an intersystem crossing and a rate-determining step of bond formation. Following ionizing radiation, guanine and the uracil radical react, with a rate-limiting step involving hydrogen atom removal. Although cross-link formation from UV radiation is favored, comparison of calculated reaction thermokinetics with that for related experimentally observed purine-pyrimidine cross-links suggests this lesion is also likely to form from ionizing radiation. For the first time, the structure of 5'-G[8-5]U-3' within DNA is identified by molecular dynamics simulations. Furthermore, three conformations of cross-linked DNA are revealed, which differ in the configuration of the complementary bases. Distortions, such as unwinding, are localized to the cross-linked dinucleotide and complementary nucleotides, with minimal changes to the flanking bases. Global changes to the helix, such as bending and groove alterations, parallel cisplatin-induced distortions, which indicate 5'-G[8-5]U-3', may contribute to the cytotoxicity of halouracils in tumor cell DNA using similar mechanisms.