Tetramolecular G-quadruplex formation pathways studied by electrospray mass spectrometry

Ionic strength of electrospray droplets affects charging of DNA oligonucleotides

The fundamental aspects of charging in electrospray ionization (ESI) are hotly debated. In the present study, ESI charging of DNA oligonucleotides was explored in both positive (ESI+) and negative (ESI-) polarity using mass spectrometry detection. Single-stranded 12-mer CCCCAATTCCCC in buffer solution (aqueous NH4Ac, 100 mM) produced similar charge state distribution (CSD) in either ESI+ or ESI-. Similarity of CSD in ESI+ and ESI- was also observed for the double-stranded 12-mer CGCGAATTCGCG. By adding typical low-vapor reagents (e.g. m-nitro benzyl alcohol, m-NBA; sulfolane) into the same buffer solution (<0.5% w/v), both CCCCAATTCCCC and CGCGAATTCGCG revealed strong supercharging (SC) effect in ESI-, while very little or no SC effect was observed in ESI+. With either sulfolane or m-NBA, the CGCGAATTCGCG duplex dissociated into single strands in ESI-. No SC was observed in both ESI+ and ESI- for thermally denatured CGCGAATTCGCG duplex in NH4 Ac buffer without the reagents. These findings are difficult to reconcile with the earlier model, which attributes SC in aqueous buffer solution to the conformational changes of analytes. Our observations suggest that the ionic strength of ESI droplets strongly affects the CSD of biopolymers such as DNA oligonucleotides and that SC effect is related to the depletion of ionic strength during the ESI process.

Heparin decamer bridges a growth factor and an oligolysine by different charge-driven interactions

Full-length heparin is widely used in tissue engineering applications due its multiple protein-binding sites that allow it to retain growth factor affinity while associating with oligopeptide components of the tissue scaffold. However, the extent to which oligopeptide coupling interferes with cognate protein binding is difficult to predict. To investigate such simultaneous interactions, we examined a well-defined ternary system comprised of acidic fibroblast growth factor (FGF), tetralysine (K4), with a heparin decamer (dp10) acting as a noncovalent coupler. Electrospray ionization mass spectrometry was used to assess binding affinities and complex stoichiometries as a function of ionic strength for dp10·K4 and FGF·dp10. The ionic strength dependence of K4·dp10 formation is qualitatively consistent with binding driven by the release of condensed counterions previously suggested for native heparin with divalent oligopeptides (Mascotti, D. P.; Lohman, T. M. Biochemistry 1995, 34, 2908-2915). On the other hand, FGF binding displays more complex ionic strength dependence, with higher salt resistance. Remarkably, dp10 that can bind two FGF molecules can only bind one tetralysine. The limited binding of K4 to dp10 suggests that the tetralysine might not block growth factor binding, and the 1:1:1 ternary complex is indeed observed. The analysis of mass distribution of the bound dp10 chains in FGF·dp10, FGF2·dp10, and FGF·dp10·K4 complexes indicated that higher degrees of dp10 sulfation promote the formation of FGF2·dp10 and FGF·dp10·K4. Thus, the selectivity of appropriately chosen short heparin chains could be used to modulate growth factor sequestration and release in a way not feasible with heterogeneous native heparin. In support of this, human hepatocellular carcinoma cells (HEP3Bs) treated with FGF·dp10·K4 were found to exhibit biological activity similar to cells treated with FGF.

Sequence and solvent effects on telomeric DNA bimolecular G-quadruplex folding kinetics

Telomeric DNA sequences are particularly polymorphic: the adopted structure is exquisitely sensitive to the sequence and to the chemical environment, for example, solvation. Dehydrating conditions are known to stabilize G-quadruplex structures, but information on how solvation influences the individual rates of folding and unfolding of G-quadruplexes remains scarce. Here, we used electrospray mass spectrometry for the first time to monitor bimolecular G-quadruplex formation from 12-mer telomeric strands, in the presence of common organic cosolvents (methanol, ethanol, isopropanol, and acetonitrile). Based on the ammonium ion distribution, the total dimer signal was decomposed into contributions from the parallel and antiparallel structures to obtain individual reaction rates, and the antiparallel G-quadruplex structure was found to form faster than the parallel one. A dimeric reaction intermediate, in rapid equilibrium with the single strands, was also identified. Organic cosolvents increase the stability of the final structures mainly by increasing the folding rates. Our quantitative analysis of reaction rate dependence on cosolvent percentage shows that organic cosolvent molecules can be captured or released upon G-quadruplex formation, highlighting that they are not inert with DNA. In contrast to the folding rates, the G-quadruplex unfolding rates are almost insensitive to solvation effects, but are instead governed by the sequence and by the final structure: parallel dimers dissociate slower than antiparallel dimers only when thymine bases are present at the 5'-end. These results contribute unraveling the folding pathways of telomeric G-quadruplexes. The solvent effects revealed here enlighten that G-quadruplex structure in dehydrated, and molecularly crowded environments are modulated by the nature of cosolvent (e.g., methanol favors antiparallel structures) due to direct interactions, and by the time scale of the reaction, with >200-fold acceleration of bimolecular G-quadruplex formation in the presence of 60% cosolvent.

Structural dynamics of possible late-stage intermediates in folding of quadruplex DNA studied by molecular simulations

Explicit solvent molecular dynamics simulations have been used to complement preceding experimental and computational studies of folding of guanine quadruplexes (G-DNA). We initiate early stages of unfolding of several G-DNAs by simulating them under no-salt conditions and then try to fold them back using standard excess salt simulations. There is a significant difference between G-DNAs with all-anti parallel stranded stems and those with stems containing mixtures of syn and anti guanosines. The most natural rearrangement for all-anti stems is a vertical mutual slippage of the strands. This leads to stems with reduced numbers of tetrads during unfolding and a reduction of strand slippage during refolding. The presence of syn nucleotides prevents mutual strand slippage; therefore, the antiparallel and hybrid quadruplexes initiate unfolding via separation of the individual strands. The simulations confirm the capability of G-DNA molecules to adopt numerous stable locally and globally misfolded structures. The key point for a proper individual folding attempt appears to be correct prior distribution of syn and anti nucleotides in all four G-strands. The results suggest that at the level of individual molecules, G-DNA folding is an extremely multi-pathway process that is slowed by numerous misfolding arrangements stabilized on highly variable timescales.

Self-assembled G-quadruplex nanostructures: AFM and voltammetric characterization

G-rich oligodeoxynucleotides (ODNs) have great medical and nanotechnological potential, because they can self-assemble into G-quadruplexes and higher-order nanostructures. The folding properties of d(G)10, d(TG9) and d(TG8T) ODNs were studied using atomic force microscopy (AFM) and voltammetry at carbon electrodes. Single-stranded ODNs, in Na(+) containing solutions and for short incubation times, were detected using AFM as network films and polymeric structures and using voltammetry by the occurrence of only the guanine oxidation peak. G-quadruplexes, in Na(+) containing solutions and long incubation times, or in K(+) containing solutions, were detected using AFM as spherical aggregates and using voltammetry by the decrease of the guanine oxidation peak and the occurrence of the G-quartet oxidation peak. Concerning the self-assembling into higher-order nanostructures, d(G)10 was the only sequence forming G-nanowires observed using AFM, d(TG9) formed short G-based super-structures that adsorbed as rod-like shape aggregates, and d(TG8T) formed no nanostructures, due to the presence of thymine residues at both 5' and 3' ends.

Tri-G-quadruplex: controlled assembly of a G-quadruplex structure from three G-rich strands

In my (DNA) dreams: A tri-G-quadruplex was constructed from three strands (T1-T3) of DNA using duplex formation to guide the G-rich tracts into close proximity with the addition of Li(+) ions (see scheme). The defined G-quadruplex structure was formed upon addition of Na(+) ions and characterized by gel electrophoresis and spectroscopy.

Investigating the interactions between cations, peroxidation substrates and G-quadruplex topology in DNAzyme peroxidation reactions using statistical testing

Cations affect the topology and enzymatic proficiency of most macromolecular catalysts but the role of cations in DNAzyme peroxidation reactions remains unresolved. Herein, we use statistical methods (ANOVA, t-test and Wilcoxon Mann-Whitney non-parametric test) to demonstrate that there are strong associations between cations, DNAzyme topology, peroxidation substrate and peroxidation rates of G-quadruplex peroxidises. Ammonium cation was found to be superior to all tested cations, including potassium. A t-test indicated that NH(4)(+) was better than K(+) with a p-value=0.05. Interestingly, the nature of the peroxidation substrate employed affected the dependence of peroxidation rate on the cation present and of the three substrates tested, 2,2'-azino-bis (3-ethylbenzothiazoline-6-sulfonic acid (ABTS), tyramine and 3,3',5,5'-tetramethylbenzidine (TMB), ABTS was the most sensitive to the nature of cation present.

Molecular dynamics simulations of G-DNA and perspectives on the simulation of nucleic acid structures

The article reviews the application of biomolecular simulation methods to understand the structure, dynamics and interactions of nucleic acids with a focus on explicit solvent molecular dynamics simulations of guanine quadruplex (G-DNA and G-RNA) molecules. While primarily dealing with these exciting and highly relevant four-stranded systems, where recent and past simulations have provided several interesting results and novel insight into G-DNA structure, the review provides some general perspectives on the applicability of the simulation techniques to nucleic acids.

Isomer-selected photoelectron spectroscopy of isolated DNA oligonucleotides: phosphate and nucleobase deprotonation at high negative charge states

Fractionation according to ion mobility and mass-to-charge ratio has been used to select individual isomers of deprotonated DNA oligonucleotide multianions for subsequent isomer-resolved photoelectron spectroscopy (PES) in the gas phase. Isomer-resolved PE spectra have been recorded for tetranucleotides, pentanucleotides, and hexanucleotides. These were studied primarily in their highest accessible negative charge states (3-, 4-, and 5-, respectively), as provided by electrospraying from room temperature solutions. In particular, the PE spectra obtained for pentanucleotide tetraanions show evidence for two coexisting classes of gas-phase isomeric structures. We suggest that these two classes comprise: (i) species with excess electrons localized exclusively at deprotonated phosphate backbone sites and (ii) species with at least one deprotonated base (in addition to several deprotonated phosphates). By permuting the sequence of bases in various [A(5-x)T(x)](4-) and [GT(4)](4-) pentanucleotides, we have established that the second type of isomer is most likely to occur if the deprotonated base is located at the first or last position in the sequence. We have used a combination of molecular mechanics and semiempirical calculations together with a simple electrostatic model to explore the photodetachment mechanism underlying our photoelectron spectra. Comparison of predicted to measured photoelectron spectra suggests that a significant fraction of the detected electrons originates from the DNA bases (both deprotonated and neutral).

Interaction of pyrrolobenzodiazepine (PBD) ligands with parallel intermolecular G-quadruplex complex using spectroscopy and ESI-MS

Studies on ligand interaction with quadruplex DNA, and their role in stabilizing the complex at concentration prevailing under physiological condition, has attained high interest. Electrospray ionization mass spectrometry (ESI-MS) and spectroscopic studies in solution were used to evaluate the interaction of PBD and TMPyP4 ligands, stoichiometry and selectivity to G-quadruplex DNA. Two synthetic ligands from PBD family, namely pyrene-linked pyrrolo[2,1-c][1,4]benzodiazepine hybrid (PBD1), mixed imine-amide pyrrolobenzodiazepine dimer (PBD2) and 5,10,15,20-tetrakis(N-methyl-4-pyridyl)porphyrin (TMPyP4) were studied. G-rich single-stranded oligonucleotide d(5'GGGGTTGGGG3') designated as d(T(2)G(8)), from the telomeric region of Tetrahymena Glaucoma, was considered for the interaction with ligands. ESI-MS and spectroscopic methods viz., circular dichroism (CD), UV-Visible, and fluorescence were employed to investigate the G-quadruplex structures formed by d(T(2)G(8)) sequence and its interaction with PBD and TMPyP4 ligands. From ESI-MS spectra, it is evident that the majority of quadruplexes exist as d(T(2)G(8))(2) and d(T(2)G(8))(4) forms possessing two to ten cations in the centre, thereby stabilizing the complex. CD band of PBD1 and PBD2 showed hypo and hyperchromicity, on interaction with quadruplex DNA, indicating unfolding and stabilization of quadruplex DNA complex, respectively. UV-Visible and fluorescence experiments suggest that PBD1 bind externally where as PBD2 intercalate moderately and bind externally to G-quadruplex DNA. Further, melting experiments using SYBR Green indicate that PBD1 unfolds and PBD2 stabilizes the G-quadruplex complex. ITC experiments using d(T(2)G(8)) quadruplex with PBD ligands reveal that PBD1 and PBD2 prefer external/loop binding and external/intercalative binding to quadruplex DNA, respectively. From experimental results it is clear that the interaction of PBD2 and TMPyP4 impart higher stability to the quadruplex complex.

Mass spectrometry and ion mobility spectrometry of G-quadruplexes. A study of solvent effects on dimer formation and structural transitions in the telomeric DNA sequence d(TAGGGTTAGGGT)

We survey here state of the art mass spectrometry methodologies for investigating G-quadruplexes, and will illustrate them with a new study on a simple model system: the dimeric G-quadruplex of the 12-mer telomeric DNA sequence d(TAGGGTTAGGGT), which can adopt either a parallel or an antiparallel structure. We will discuss the solution conditions compatible with electrospray ionisation, the quantification of complexes using ESI-MS, the interpretation of ammonium ion preservation in the complexes in the gas phase, and the use of ion mobility spectrometry to resolve ambiguities regarding the strand stoichiometry, or separate and characterise different structural isomers. We also describe that adding electrospray-compatible organic co-solvents (methanol, ethanol, isopropanol or acetonitrile) to aqueous ammonium acetate increases the stability and rate of formation of dimeric G-quadruplexes, and causes structural transitions to parallel structures. Structural changes were probed by circular dichroism and ion mobility spectrometry, and the excellent correlation between the two techniques validates the use of ion mobility to investigate G-quadruplex folding. We also demonstrate that parallel G-quadruplex structures are easier to preserve in the gas phase than antiparallel structures.

Charge and substituent effects on the stability of porphyrin/G-quadruplex adducts

The adduct ions of two tetramolecular G-quadruplexes formed from the d(TGGGGT) and d(TTGGGGGT) single strands with a group of cationic porphyrins, with different charges and substituents, and one neutral porphyrin, were investigated by ESI-MS and ESI-MS/MS in the negative ion mode. Formation of [Q + nNH(4)(+)+P(p+)-(z + n + p)H(+)](z-) adduct ions (where Q = quadruplex, n = number of quartets minus 1, P = porphyrin and p(+) = 0,1,2,3,4) indicates that the porphyrins are bound outside the quadruplexes providing an additional stabilization to those structures. The fragmentation pathways of the [Q + nNH(4)(+)+P(p+)-(z + n + p)H(+)](z-) adduct ions depend on the number of positive charges (p(+)) of the porphyrins and on the overall complex charge (z(-)), but do not show a significant dependence on the type of the substituent groups in the porphyrins. Formation of the 'unfilled' ions [Q + P(p+)-(z + p)H(+)](z-) predominates for porphyrins with a higher number of positive charges. Strand separation with the formation of [T + P(p+)-(z-2 + p)H(+)]((z-2)-) and (SS-2H(+))(2-) ions, where T = [d(TG(4)T)](3) and [d(T(2)G(5)T)](3) and SS = d(TG(4)T) and d(T(2)G(5)T) is only observed for the complexes with a higher overall negative charge. Porphyrin loss with the formation of [Q + nNH(4)(+)-(z + n)H(+)](z-) ions occurs predominantly for the neutral and monocharged porphyrins. The predominant formation of the 'unfilled' ions, [Q + P(p+)-(z + n)H(+)](z-), for porphyrins with a higher number of charges shows that these porphyrins can prevent strand separation and preserve, at least partially, the quadruplex structure.

Conformational changes of non-B DNA

In contrast to B-DNA that has a right-handed double helical structure with Watson-Crick base pairing under the ordinary physiological conditions, repetitive DNA sequences under certain conditions have the potential to fold into non-B DNA structures such as hairpin, triplex, cruciform, left-handed Z-form, tetraplex, A-motif, etc. Since the non-B DNA-forming sequences induce the genetic instability and consequently can cause human diseases, the molecular mechanism for their genetic instability has been extensively investigated. On the contrary, non-B DNA can be widely used for application in biotechnology because many DNA breakage hotspots are mapped in or near the sequences that have the potential to adopt non-B DNA structures. In addition, they are regarded as a fascinating material for the nanotechnology using non-B DNAs because they do not produce any toxic byproducts and are robust enough for the repetitive working cycle. This being the case, an understanding on the mechanism and dynamics of their structural changes is important. In this critical review, we describe the latest studies on the conformational dynamics of non-B DNAs, with a focus on G-quadruplex, i-motif, Z-DNA, A-motif, hairpin and triplex (189 references).

Insight into G-DNA structural polymorphism and folding from sequence and loop connectivity through free energy analysis

The lengths of G-tracts and their connecting loop sequences determine G-quadruplex folding and stability. Complete understanding of the sequence–structure relationships remains elusive. Here, single-loop G-quadruplexes were investigated using explicit solvent molecular dynamics (MD) simulations to characterize the effect of loop length, loop sequence, and G-tract length on the folding topologies and stability of G-quadruplexes. Eight loop types, including different variants of lateral, diagonal, and propeller loops, and six different loop sequences [d0 (i.e., no intervening residues in the loop), dT, dT₂, dT₃, dTTA, and dT₄] were considered through MD simulation and free energy analysis. In most cases the free energetic estimates agree well with the experimental observations. The work also provides new insight into G-quadruplex folding and stability. This includes reporting the observed instability of the left propeller loop, which extends the rules for G-quadruplex folding. We also suggest a plausible explanation why human telomere sequences predominantly form hybrid-I and hybrid-II type structures in K⁺ solution. Overall, our calculation results indicate that short loops generally are less stable than longer loops, and we hypothesize that the extreme stability of sequences with very short loops could possibly derive from the formation of parallel multimers. The results suggest that free energy differences, estimated from MD and free energy analysis with current force fields and simulation protocols, are able to complement experiment and to help dissect and explain loop sequence, loop length, and G-tract length and orientation influences on G-quadruplex structure.

The application of DNA and RNA G-quadruplexes to therapeutic medicines

The intriguing structural diversity in folded topologies available to guanine-rich nucleic acid repeat sequences have made four-stranded G-quadruplex structures the focus of both basic and applied research, from cancer biology and novel therapeutics through to nanoelectronics. Distributed widely in the human genome as targets for regulating gene expression and chromosomal maintenance, they offer unique avenues for future cancer drug development. In particular, the recent advances in chemical and structural biology have enabled the construction of bespoke selective DNA based aptamers to be used as novel therapeutic agents and access to detailed structural models for structure based drug discovery. In this critical review, we will explore the important underlying characteristics of G-quadruplexes that make them functional, stable, and predictable nanoscaffolds. We will review the current structural database of folding topologies, molecular interfaces and novel interaction surfaces, with a consideration to their future exploitation in drug discovery, molecular biology, supermolecular assembly and aptamer design. In recent years the number of potential applications for G-quadruplex motifs has rapidly grown, so in this review we aim to explore the many future challenges and highlight where possible successes may lie. We will highlight the similarities and differences between DNA and RNA folded G-quadruplexes in terms of stability, distribution, and exploitability as small molecule targets. Finally, we will provide a detailed review of basic G-quadruplex geometry, experimental tools used, and a critical evaluation of the application of high-resolution structural biology and its ability to provide meaningful and valid models for future applications (255 references).

d(CGGTGGT) forms an octameric parallel G-quadruplex via stacking of unusual G(:C):G(:C):G(:C):G(:C) octads

Among non-canonical DNA secondary structures, G-quadruplexes are currently widely studied because of their probable involvement in many pivotal biological roles, and for their potential use in nanotechnology. The overall quadruplex scaffold can exhibit several morphologies through intramolecular or intermolecular organization of G-rich oligodeoxyribonucleic acid strands. In particular, several G-rich strands can form higher order assemblies by multimerization between several G-quadruplex units. Here, we report on the identification of a novel dimerization pathway. Our Nuclear magnetic resonance, circular dichroism, UV, gel electrophoresis and mass spectrometry studies on the DNA sequence dCGGTGGT demonstrate that this sequence forms an octamer when annealed in presence of K⁺ or NH₄⁺ ions, through the 5′-5′ stacking of two tetramolecular G-quadruplex subunits via unusual G(:C):G(:C):G(:C):G(:C) octads.

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.

Conformational changes in quadruplex oligonucleotide structures probed by Raman spectroscopy

Quadruplex structures are higher order structures formed by guanine-rich oligonucleotides. In the present study, temperature-induced conformational changes in the quadruplex structures of aptamers and other guanine-rich oligonucleotides are probed by Raman spectroscopy. In particular, dramatic changes in the fingerprint region are observed in the spectra of thrombin binding aptamer at higher temperatures. These changes are accompanied by a decrease in the intensity of the 1480 cm(-1) peak (attributed to C8 = N7-H2), which is diagnostic of the quadruplex structure. We also show that these changes can be reversed (to a certain extent) by addition of K(+) ions.

Kinetically governed polymorphism of d(G₄T₄G₃) quadruplexes in K+ solutions

It has been generally recognized that understanding the molecular basis of some important cellular processes is hampered by the lack of knowledge of forces that drive spontaneous formation/disruption of G-quadruplex structures in guanine-rich DNA sequences. According to numerous biophysical and structural studies G-quadruplexes may occur in the presence of K(+) and Na(+) ions as polymorphic structures formed in kinetically governed processes. The reported kinetic models suggested to describe this polymorphism should be considered inappropriate since, as a rule, they include bimolecular single-step associations characterized by negative activation energies. In contrast, our approach in studying polymorphic behavior of G-quadruplexes is based on model mechanisms that involve only elementary folding/unfolding transitions and structural conversion steps that are characterized by positive activation energies. Here, we are investigating a complex polymorphism of d(G(4)T(4)G(3)) quadruplexes in K(+) solutions. On the basis of DSC, circular dichroism and UV spectroscopy and polyacrylamide gel electrophoresis experiments we propose a kinetic model that successfully describes the observed thermally induced conformational transitions of d(G(4)T(4)G(3)) quadruplexes in terms of single-step reactions that involve besides single strands also one tetramolecular and three bimolecular quadruplex structures.