Hydrogen bonding versus stacking stabilization by modified nucleobases incorporated in PNA.DNA duplexes

Mechanical characterization of base analogue modified nucleic acids by force spectroscopy

We use mechanical unfolding of single DNA hairpins with modified bases to accurately assess intra- and intermolecular forces in nucleic acids. As expected, the modification stabilizes the hybridized hairpin, but we also observe intriguing stacking interactions in the unfolded hairpin. Our study highlights the benefit of using base-modified nucleic acids in force-spectroscopy.

Nanomechanics of Diaminopurine-Substituted DNA

2,6-diaminopurine (DAP) is a nucleobase analog of adenine. When incorporated into double-stranded DNA (dsDNA), it forms three hydrogen bonds with thymine. Rare in nature, DAP substitution alters the physical characteristics of a DNA molecule without sacrificing sequence specificity. Here, we show that in addition to stabilizing double-strand hybridization, DAP substitution also changes the mechanical and conformational properties of dsDNA. Thermal melting experiments reveal that DAP substitution raises melting temperatures without diminishing sequence-dependent effects. Using a combination of atomic force microscopy (AFM), magnetic tweezer (MT) nanomechanical assays, and circular dichroism spectroscopy, we demonstrate that DAP substitution increases the flexural rigidity of dsDNA yet also facilitates conformational shifts, which manifest as changes in molecule length. DAP substitution increases both the static and dynamic persistence length of DNA (measured by AFM and MT, respectively). In the static case (AFM), in which tension is not applied to the molecule, the contour length of DAP-DNA appears shorter than wild-type (WT)-DNA; under tension (MT), they have similar dynamic contour lengths. At tensions above 60 pN, WT-DNA undergoes characteristic overstretching because of strand separation (tension-induced melting) and spontaneous adoption of a conformation termed S-DNA. Cyclic overstretching and relaxation of WT-DNA at near-zero loading rates typically yields hysteresis, indicative of tension-induced melting; conversely, cyclic stretching of DAP-DNA showed little or no hysteresis, consistent with the adoption of the S-form, similar to what has been reported for GC-rich sequences. However, DAP-DNA overstretching is distinct from GC-rich overstretching in that it happens at a significantly lower tension. In physiological salt conditions, evenly mixed AT/GC DNA typically overstretches around 60 pN. GC-rich sequences overstretch at similar if not slightly higher tensions. Here, we show that DAP-DNA overstretches at 52 pN. In summary, DAP substitution decreases the overall stability of the B-form double helix, biasing toward non-B-form DNA helix conformations at zero tension and facilitating the B-to-S transition at high tension.

Electronic properties of DNA: Description of weak interactions in TATA-box-like chains

DNA is one of the most important biomolecules since it contains all the genetic information about an organism. The tridimensional structure of DNA is a determinant factor that influences the physiological and biochemical mechanisms by which this molecule carries out its biological functions. It is believed that hydrogen bonds and π-π stacking are the most relevant non-covalent interactions regarding DNA stability. Due to its importance, several theoretical works have been made to describe these interactions, however, most of them often consider only the presence of two nitrogenous bases, having a limited overview of the participation of these in B-DNA stabilization. Furthermore, due to the complexity of the system, there are discrepancies between which involved interaction is more important in duplex stability. Therefore, in this project we describe these interactions considering the effect of chain length on the energy related to both hydrogen bonds and π-π stacking, using as model TATA-box-like chains with n base pairs (n=1 to 14) and taking into consideration two different models: ideal and optimized B-DNA. We have found that there is a cooperative effect on hydrogen bond and π-π stacking mean energies when the presence of other base pairs is considered. In addition, it was found that hydrogen bonds contribute more importantly than π-π stacking to B-DNA stability; nevertheless, the participation of π-π stacking is not negligible: when B-DNA looks for a conformation of lower energy, π-π stacking interaction are the first to be optimized. All work was realized under the framework of DFT using the DMol³ code (M06-L/DNP).

Focus on PNA Flexibility and RNA Binding using Molecular Dynamics and Metadynamics

Peptide Nucleic Acids (PNAs) can efficiently target DNA or RNA acting as chemical tools for gene regulation. Their backbone modification and functionalization is often used to increase the affinity for a particular sequence improving selectivity. The understanding of the trading forces that lead the single strand PNA to bind the DNA or RNA sequence is preparatory for any further rational design, but a clear and unique description of this process is still not complete. In this paper we report further insights into this subject, by a computational investigation aiming at the characterization of the conformations of a single strand PNA and how these can be correlated to its capability in binding DNA/RNA. Employing Metadynamics we were able to better define conformational pre-organizations of the single strand PNA and γ-modified PNA otherwise unrevealed through classical molecular dynamics. Our simulations driven on backbone modified PNAs lead to the conclusion that this γ-functionalization affects the single strand preorganization and targeting properties to the DNA/RNA, in agreement with circular dichroism (CD) spectra obtained for this class of compounds. MD simulations on PNA:RNA dissociation and association mechanisms allowed to reveal the critical role of central bases and preorganization in the binding process.

Influence of PNA containing 8-aza-7-deazaadenine on structure stability and binding affinity of PNA·DNA duplex: insights from thermodynamics, counter ion, hydration and molecular dynamics analysis

This paper describes the synthesis of a novel 8-aza-7-deazapurin-2,6-diamine (DPP)-containing peptide nucleic acid (PNA) monomer and Boc protecting group-based oligomerization of PNA, replacing adenine (A) with DPP monomers in the PNA strand. The PNA oligomers were synthesized against the biologically relevant SV40 promoter region (2494-AATTTTTTTTATTTA-2508) of pEGFP-N3 plasmid. The DPP-PNA·DNA duplex showed enhanced stability as compared to normal duplex (A-PNA·DNA). The electronic distribution of DPP monomer suggested that DPP had better electron donor properties over 2,6-diamino purine. UV melting and thermodynamic analysis revealed that the PNA oligomer containing a diaminopyrazolo(3,4-d)pyrimidine moiety (DPP) stabilized the PNA·DNA hybrids compared to A-PNA·DNA. DPP-PNA·DNA duplex showed higher water activity (Δnw = 38.5) in comparison to A-PNA·DNA duplex (Δnw = 14.5). The 50 ns molecular dynamics simulations of PNA·DNA duplex containing DPP or unmodified nucleobase-A showed average H-bond distances in the DPP-dT base pair of 2.90 Å (OH-N bond) and 2.91 Å (NH-N bond), which were comparably shorter than in the A-dT base pair, in which the average distances were 3.18 Å (OH-N bond) and 2.97 Å (NH-N bond), and there was one additional H-bond in the DPP-dT base pair of around 2.98 Å (O2H-N2 bond), supporting the higher stability of DPP-PNA·DNA. The analysis of molecular dynamics simulation data showed that the system binding free energy increased at a rate of approximately -4.5 kcal mol(-1) per DPP base of the PNA·DNA duplex. In summary, increased thermal stability, stronger hydrogen bonding and more stable conformation in the DPP-PNA·DNA duplex make it a better candidate as antisense/antigene therapeutic agents.

Structural perturbations induced by the alpha-anomer of the aflatoxin B(1) formamidopyrimidine adduct in duplex and single-strand DNA

The guanine N7 adduct of aflatoxin B₁exo-8,9-epoxide hydrolyzes to form the formamidopyrimidine (AFB-FAPY) adduct, which interconverts between α and β anomers. The β anomer is highly mutagenic in Escherichia coli, producing G → T transversions; it thermally stabilizes the DNA duplex. The AFB-α-FAPY adduct blocks replication; it destabilizes the DNA duplex. Herein, the structure of the AFB-α-FAPY adduct has been elucidated in 5′-d(C¹T²A³T⁴X⁵A⁶T⁷T⁸C⁹A¹⁰)-3′·5′-d(T¹¹G¹²A¹³A¹⁴T¹⁵C¹⁶A¹⁷T¹⁸A¹⁹G²⁰)-3′ (X = AFB-α-FAPY) using molecular dynamics calculations restrained by NMR-derived distances and torsion angles. The AFB moiety intercalates on the 5′ face of the pyrimidine moiety at the damaged nucleotide between base pairs T⁴·A¹⁷ and X⁵·C¹⁶, placing the FAPY C5−N⁵ bond in the R_a axial conformation. Large perturbations of the ε and ζ backbone torsion angles are observed, and the base stacking register of the duplex is perturbed. The deoxyribose orientation shifts to become parallel to the FAPY base and displaced toward the minor groove. Intrastrand stacking between the AFB moiety and the 5′ neighbor thymine remains, but strong interstrand stacking is not observed. A hydrogen bond between the formyl group and the exocyclic amine of the 3′-neighbor adenine stabilizes the E conformation of the formamide moiety. NMR studies reveal a similar 5′-intercalation of the AFB moiety for the AFB-α-FAPY adduct in the tetramer 5′-d(C¹T²X³A⁴)-3′, involving the R_a axial conformation of the FAPY C5−N⁵ bond and the E conformation of the formamide moiety. Since in duplex DNA the AFB moiety of the AFB-β-FAPY adduct also intercalates on the 5′ side of the pyrimidine moiety at the damaged nucleotide, we conclude that favorable 5′-stacking leads to the R_a conformational preference about the C5−N⁵ bond; the same conformational preference about this bond is also observed at the nucleoside and base levels. The structural distortions and the less favorable stacking interactions induced by the AFB-α-FAPY adduct explain its lower stability as compared to the AFB-β-FAPY adduct in duplex DNA. In this DNA sequence, hydrogen bonding between the formyl oxygen and the exocyclic amine of the 3′-neighboring adenine stabilizing the E configuration of the formamide moiety is also observed for the AFB-β-FAPY adduct, and suggests that the identity of the 3′-neighbor nucleotide modulates the stability and biological processing of AFB adducts.