Analysis of GAA/TTC DNA triplexes using nuclear magnetic resonance and electrospray ionization mass spectrometry

Native Mass Spectrometric Insights into the Formation and Stability of DNA Triplexes

Deoxyribonucleic acid is a genetic biomacromolecule that contains the inherited information required to build and maintain a living organism. While the canonical duplex DNA structure is rigorously characterized, the structure and function of higher order DNA structures such as DNA triplexes are comparatively poorly understood. Previous literature has shown that these triplexes can form under physiological conditions, and native mass spectrometry offers a useful platform to study their formation and stability. However, experimental conditions including buffer salt concentration, pH, and instrumentation parameters such as ion mode can confound analysis by impacting the amount of triplex observed by mass spectrometry. Using model 30mer Y-type triplex sequences, we demonstrate the influence a range of experimental variables have on the detection of DNA triplex structures, informing suitable conditions for the study. When carefully considered conditions are used, mass spectrometry offers a powerful complementary tool for the analysis of higher order DNA assemblies.

Mass Spectrometry of Nucleic Acid Noncovalent Complexes

Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.

DNA triplex structure, thermodynamics, and destabilisation: insight from molecular simulations

Molecular dynamics simulations are used to elucidate the structure and thermodynamics of DNA triplexes associated with the neurodegenerative disease Friedreich's ataxia (FRDA), as well as complexes of these triplexes with the small molecule netropsin, which is known to destabilise triplexes. The ability of molecular simulations in explicit solvent to accurately capture triplex thermodynamics is verified for the first time, with the free energy to dissociate a 15-base antiparallel purine triplex-forming oligomer (TFO) from the duplex found to be slightly higher than reported experimentally. The presence of netropsin in the minor groove destabilises the triplex as expected, reducing the dissociation free energy by approximately 50%. Netropsin binding is associated with localised narrowing of the minor groove near netropsin, an effect that has previously been under contention. This leads to localised widening of the major groove, weakening hydrogen bonds between the TFO and duplex. Consequently, destabilisation is found to be highly localised, occurring only when netropsin is bound directly opposite the TFO. The simulations also suggest that near saturation of the minor groove with ligand is required for complete triplex dissociation. A structural analysis of the DNA triplexes that can form with the FRDA-related duplex sequence indicates that the triplex with a parallel homopyrimidine TFO is likely to be more stable than the antiparallel homopurine-TFO triplex, which may have implications for disease onset and treatment.

Translocation and deletion breakpoints in cancer genomes are associated with potential non-B DNA-forming sequences

Gross chromosomal rearrangements (including translocations, deletions, insertions and duplications) are a hallmark of cancer genomes and often create oncogenic fusion genes. An obligate step in the generation of such gross rearrangements is the formation of DNA double-strand breaks (DSBs). Since the genomic distribution of rearrangement breakpoints is non-random, intrinsic cellular factors may predispose certain genomic regions to breakage. Notably, certain DNA sequences with the potential to fold into secondary structures [potential non-B DNA structures (PONDS); e.g. triplexes, quadruplexes, hairpin/cruciforms, Z-DNA and single-stranded looped-out structures with implications in DNA replication and transcription] can stimulate the formation of DNA DSBs. Here, we tested the postulate that these DNA sequences might be found at, or in close proximity to, rearrangement breakpoints. By analyzing the distribution of PONDS-forming sequences within ±500 bases of 19 947 translocation and 46 365 sequence-characterized deletion breakpoints in cancer genomes, we find significant association between PONDS-forming repeats and cancer breakpoints. Specifically, (AT)_n, (GAA)_n and (GAAA)_n constitute the most frequent repeats at translocation breakpoints, whereas A-tracts occur preferentially at deletion breakpoints. Translocation breakpoints near PONDS-forming repeats also recur in different individuals and patient tumor samples. Hence, PONDS-forming sequences represent an intrinsic risk factor for genomic rearrangements in cancer genomes.

Developments in Electrospray Ionization Mass Spectrometry of Non-Covalent DNA–Ligand Complexes

Many anti-cancer drugs function by binding non-covalently to double-stranded (ds) DNA. Electrospray ionization mass spectrometry (ESI-MS) has emerged over the past decade as a sensitive technique for the determination of stoichiometries and relative binding affinities of DNA–ligand interactions. The chromosome contains nucleotide sequences, for example, guanosine-rich regions, that predispose them to the formation of higher order structures such as quadruplex DNA (qDNA). Sequences that form qDNA are found in the telomeres. The proposal that ligands that stabilize qDNA might interfere with the activity of telomerase in cancer cells has stimulated the search for ligands that are selective for qDNA over dsDNA. The insights gained from the development of ESI-MS methods for analysis of non-covalent dsDNA–ligand complexes are now being applied in the search for qDNA-selective ligands. ESI-MS is a useful first-pass screening technique for qDNA-binding ligands. This short review describes some experimental considerations for ESI-MS analysis of DNA–ligand complexes, briefly addresses the question of whether non-covalent DNA–ligand complexes are faithfully transferred from solution to the gas phase, discusses ion mobility mass spectrometry as a technique for probing this issue, and highlights some recent ESI-MS studies of qDNA-selective ligands.

Evaluation of effects of bivalent cations on the formation of purine-rich triple-helix DNA by ESI-FT-MS

The GGA triplet repeats are widely dispersed throughout eukaryotic genomes. (GGA)n or (GGT)n oligonucleotides can interact with double-stranded DNA containing (GGA:CCT)n to form triple-stranded DNA. The effects of 8 divalent metal ions (3 alkaline-earth metals and 5 transition metals) on formation of these purine-rich triple-helix DNA were investigated by electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FT-MS). In the absence of metal ions, no triplex but single-strand, duplex, and purine homodimer ions were observed in mass spectra. The triple-helix DNA complexes were observed only in the presence of certain divalent ions. The effects of different divalent cations on the formation of purine-rich triplexes were compared. Transition-metal ions, especially Co(2+) and Ni(2+), significantly boost the formation of triple-helix DNA, whereas alkaline-earth metal ions have no positive effects on triplex formation. In addition, Ba(2+) is notably beneficial to the formation of homodimer instead of triplex.

Interactions of mitoxantrone with duplex and triplex DNA studied by electrospray ionization mass spectrometry

We have examined interactions between mitoxantrone (MXT) and DNA duplexes or triplexes with different base compositions by using electrospray ionization mass spectrometry (ESI-MS), respectively. MXT interacts preferentially with DNA duplexes compared to the triplexes. In the mass spectrum of the duplex-MXT mixture, the complex peaks dominated in the ratios of duplex/MXT of 1:1, 1:2 and 1:3, and the 1:2 duplex/MXT peak was the most abundant. In contrast, only 1:1 triplex-MXT complexes were observed in the mass spectrum of the triplex-MXT mixture, and the most intensive peak was a free triplex ion without MXT. Moreover, no sequence selectivity of MXT to different DNA duplexes was found while MXT showed greater affinity to the triplexes that have adjacent TAT or C(+)GC sequences. In the course of sustained off-resonance irradiation collision-induced dissociation (SORI-CID), the MXT-duplex complexes generated two separated strands, and the MXT remained on the purine strand side. UV/Vis spectra showed that MXT interacted with DNA by intercalation. Compared with emodin (a duplex intercalator) and napthylquinoline (a triplex binder), we found that the side chain of MXT might play a role in the binding of MXT to the duplexes and the triplexes. ESI-MS shows an advantage in speed and straightforwardness for the study of drug interactions with nucleic acids.

Studies of the intermolecular DNA triplexes of C+.GC and T.AT triplets by electrospray ionization Fourier-transform ion cyclotron resonance mass spectrometry

Formation and stabilities of four 14-mer intermolecular DNA triplexes, consisting of third strands with repeating sequence CTCT, CCTT, CTT, or TTT, were studied by electrospray ionization Fourier-transform ion cyclotron resonance mass spectrometry (ESI-FTICR-MS) in the gas phase. The gas-phase stabilities of the triplexes were compared with their CD spectra and melting behaviors in solution, and parallel correlation between two phases were obtained. In the presence of 20 mM NH(4) (+) (pH 5.5), the formation of the TTT triplex was not detected in both solution and the gas phase. Other triplexes showed the same order, CTCT > CCTT > CTT, of ion abundances in mass spectra and T(m) values in solution. The more stable triplexes are those that contained higher percentage of C(+).GC triplets and an alternating CT sequence. However, the CCTT with the same C(+).GC triplets as the CTCT showed a higher stability than the latter during the gas-phase dissociation. Furthermore, a biphasic triplex-to-duplex-to-single transition was detected in the gas phase, while a monophasic triplex-to-single dissociation was observed in solution. The present results reveal that hydrogen bonds and electrostatic interactions dominate in the gas phase, while base stacking and hydrophobic interactions dominate in solution to stabilize the triplexes. Moreover, weak acidic conditions (pH 5-6) promote the formation of the parallel triplexes.

Gas-phase DNA oligonucleotide structures. A QM/MM and atoms in molecules study

QM/MM calculations have been employed to investigate the role of hydrogen bonding and pi-stacking in single- and double-stranded DNA oligonucleotides. DFT calculations and Atoms in Molecules analysis on QM/MM-optimized structures allow characterization and estimation of the energies of pi-stacking and hydrogen-bond interactions. This shows that pi-stacking interactions depend on the number and the nature of the DNA bases for single-stranded nucleotides; for instance, guanines are found to be involved in strong hydrogen bonds, whereas adenines interact mainly via stacking interactions. The role of interbase hydrogen bonding was explored: the -NH2 groups of guanine, adenine, and cytosine participate in N-H...O and N-H...N interactions. These are much stronger in single-strand oligonucleotides, where the -NH2 groups are highly nonplanar. In double-stranded DNA, the strong base-pairing hydrogen bonds of complementary bases lead to more planar -NH2 groups, which tend to be involved in pi-stacking interactions rather than H-bonds. The use of AIM also allows us to evaluate the interplay of pi-stacking and H-bonding, suggesting that cooperativity does occur, but is generally limited to about 1-2 kcal/mol.

Sequence-specific conjugates of oligo(2'-O-methylribonucleotides) and hairpin oligocarboxamide minor-groove binders: design, synthesis, and binding studies with double-stranded DNA

New conjugates of triplex-forming pyrimidine oligo(2'-O-methylribonucleotides) with one or two 'head-to-head' hairpin oligo(N-methylpyrrole carboxamide) minor-groove binders (MGBs) attached to the terminal phosphate of the oligonucleotides with a oligo(ethylene glycol) linker were synthesized. It was demonstrated that, under appropriate conditions, the conjugates form stable complexes with double-stranded DNA (dsDNA) similarly to triplex-forming oligo(deoxyribonucleotide) (TFO) conjugates containing 5-methylated cytosines. Kinetic and thermodynamic parameters of the complex formation were evaluated by gel-shift assay and thermal denaturation. Higher melting temperatures (Tm), faster complex formation, and lower dissociation constants (Kd) of the triple helices (6-7 nM) were observed for complexes of MGB-oligo(2'-O-methylribonucleotide) conjugates with the target dsDNA compared to the nonconjugated individual components. Interaction of MGB moieties with the HIV proviral DNA fragment was indicated by UV/VIS absorption changes at 320 nm in the melting curves. The introduction of thymidine via a 3',3'-type 'inverted' phosphodiester linkage at the 3'-end of oligo(2'-O-methylribonucleotide) conjugates (3'-protection) had no strong influence on triplex formation, but slightly affected complex stability. At pH 6.0, when one or two hairpin MGBs were attached to the oligonucleotide, both triplex formation and minor-groove binding played important roles in complex formation. When two 'head-to-head' oligo(N-methylpyrrole) ligands were attached to the same terminal phosphate of the oligonucleotide or the linker, binding was observed at pH >7.5 and at high temperatures (up to 74 degrees). However, under these conditions, binding was retained only by the MGB part of the conjugate.

Sticky DNA: in vivo formation in E. coli and in vitro association of long GAA*TTC tracts to generate two independent supercoiled domains

The expanded GAA*TTC repeat sequence associated with Friedreich's ataxia (FRDA) adopts non-B DNA structures, (triplexes and sticky DNA). Sticky DNA is formed in plasmids by the association of two long GAA*TTC tracts at lengths that are found in the sequence of the frataxin gene in patients. Most FRDA patients have expanded GAA*TTC repeats (up to 1700 triplets), which inhibit the transcription of the gene, thus diminishing the synthesis of frataxin, a mitochondrial protein involved in iron-sulfur cluster biogenesis. Negative supercoiling and MgCl(2) (or MnCl(2)) are required to stabilize sticky DNA (a dumbbell-shaped structure) in plasmids with a pair of repeat tracts where n> or =60 in the direct repeat orientation in vitro. Since the triplet repeat sequences (TRS) were symmetrically positioned in the plasmids and because a number of unique restriction sites were present in the vector, studies were conducted to evaluate the influence of selectively linearizing one or the other supercoiled domains created by the DNA*DNA associated region, i.e. the stable complex at the pair of TRS's. The two domains behave independently, thus confirming the association of the two tracts and the dumbbell-shaped plasmid in our model for sticky DNA. Linking number investigations were performed on a family of plasmids harboring different lengths (30, 60, or 176 repeats), orientations and number of tracts (one or two) of a GAA*TTC repeat in Escherichia coli to evaluate the in vivo role, if any, of sticky DNA. Unexpectedly, this non-B DNA conformation elicited the formation of a TRS-length dependent change in the global topology of the plasmids, indicative of an apparent compression of the primary helices. Thus, linking number determinations confirm that sticky DNA has an important consequence in vivo.