Mechanistic studies on the base-catalyzed transformation of neocarzinostatin chromophore: roles of bulged DNA

Nucleic acid bulges have been implicated in a number of biological processes and are specific cleavage targets for the enediyne antitumor antibiotic neocarzinostatin chromophore (NCS-chrom) in a base-catalyzed, radical-mediated reaction. Studies designed to elucidate the detailed mechanism of the base-catalyzed activation of NCS-chrom and to evaluate the roles of bulged DNA in its activation are described. They show that nucleobases in the DNA bulge are not required to form an effective bulge pocket but enhance the binding of the wedge-shaped activated drug molecule. Analysis of solvent deuterium isotope effects on NCS-chrom degradation and DNA cleavage efficiency experiments suggests that the spirolactone biradical 6 is a relatively stable species and that intramolecular quenching of the C2 radical of 6 to form the biologically active cyclospirolactone radical 7a occurs first (pathway a in Scheme 2), leaving the C6 radical to abstract the hydrogen atom from the DNA deoxyribose and to form the cyclospirolactone 8. Binding of the activated drug at the bulge site is required, but not sufficient, for efficient 8 formation, whereas cleavage of bulged DNA is not essential. Efficient generation of 8, but inefficient DNA damage generation, comes mainly from the likely high off-rate of 7a binding. The finding that thymidine 5'-carboxylic acid-ended oligonucleotide fragment can be formed in the reaction suggests that the process of DNA cleavage is rather slow and that sequential oxidations of the target 5'-carbon are possible. Study of the effect of solvent (methanol) concentration on NCS-chrom degradation indicates that bulged DNA acts to assist the intramolecular quenching of the radical at C2 by C8' ' of the naphthoate moiety by excluding solvent from the binding pocket, thus preventing the formation of spirolactones 9, and by blocking radical polymerization. Because in the absence or near absence of solvent methanol 8 formation does not reach even 10% that formed in the presence of bulged DNA, it is possible that the DNA bulge also induces a conformational change in the drug to promote the intramolecular reaction.     

New complex of post-activated neocarzinostatin chromophore with DNA: bulge DNA binding from the minor groove

Neocarzinostatin (NCS-chrom), a natural enediyne antitumor antibiotic, undergoes either thiol-dependent or thiol-independent activation, resulting in distinctly different DNA cleavage patterns. Structures of two different post-activated NCS-chrom complexes with DNA have been reported, revealing strikingly different binding modes that can be directly related to the specificity of DNA chain cleavage caused by NCS-chrom. The third structure described herein is based on recent studies demonstrating that glutathione (GSH) activated NCS-chrom efficiently cleaves DNA at specific single-base sites in sequences containing a putative single-base bulge. In this structure, the GSH post-activated NCS-chrom (NCSi-glu) binds to a decamer DNA, d(GCCAGAGAGC), from the minor groove. This binding triggers a conformational switch in DNA from a loose duplex in the free form to a single-strand, tightly folded hairpin containing a bulge adenosine embedded between a three base pair stem. The naphthoate aromatic moiety of NCSi-glu intercalates into a GG step flanked by the bulge site, and its substituent groups, the 2-N-methylfucosamine carbohydrate ring and the tetrahydroindacene, form a complementary minor groove binding surface, mostly interacting with the GCC strand in the duplex stem of DNA. The bulge site is stabilized by the interactions involving NCSi-glu naphthoate and GSH tripeptide. The positioning of NCSi-glu is such that only single-chain cleavage via hydrogen abstraction at the 5'-position of the third base C (which is opposite to the putative bulge base) in GCC is possible, explaining the observed single-base cleavage specificity. The reported structure of the NCSi-glu-bulge DNA complex reveals a third binding mode of the antibiotic and represents a new family of minor groove bulge DNA recognition structures. We predict analogue structures of NCSi-R (R = glu or other substituent groups) may be versatile probes for detecting the existence of various structures of nucleic acids. The NMR structure of this complex, in combination with the previously reported NCSi-gb-bulge DNA complex, offers models for specific recognition of DNA bulges of various sizes through binding to either the minor or the major groove and for single-chain cleavage of bulge DNA sequences.     

Molecular probes of DNA bulges: functional assay and spectroscopic analysis

Bulged structures in DNA and RNA have been linked to biomolecular processes involved in numerous diseases, thus probes with affinity for these nucleic acid targets would be of considerable utility to chemical biologists. Herein, we report guided chemical synthesis of small molecules capable of binding to DNA bulges by virtue of their unique (spirocyclic) geometry. The agents, modeled on a natural product congener, show pronounced selectivity for specific bulged motifs and are able to enhance slipped DNA synthesis, a hallmark functional assay of bulge binding. Significantly, bulge-agent complexes demonstrate characteristic fluorescent signatures depending on bulge and flanking sequence in the oligo. It is anticipated that these signature patterns can be harnessed as molecular probes of bulged hotspots in DNA and RNA.     

Designed DNA probes from the neocarzinostatin family: Impact of glycosyl linkage stereochemistry on bulge base binding

Bulged sites in DNA and RNA have become targets for rational drug design due to their suspected involvement in a number of key biomolecular processes. A lead compound, derived from the enediyne natural product NCS-chrom has been used to inform chemical synthesis of a family of designed probes of DNA bulges, one of which shows 80 nM affinity for a two base bulged target. Key contributors to binding of these spirocyclic compounds have been studied in order to correlate affinity and specificity with structural features. Herein, we demonstrate that the glycosyl linkage stereochemistry of the pendant aminofucosyl group plays a pivotal role in binding, and coupled with insight obtained with various bulged targets, will allow rational design of second generation ligands.     

Induced formation of a DNA bulge structure by a molecular wedge ligand-postactivated neocarzinostatin chromophore

Our previous structure elucidation of the complexes of DNA and postactivated neocarzinostatin chromophore (NCS-chrom) compounds revealed two distinctly different binding modes of this antitumor molecule. A thorough understanding of these results will provide the molecular basis for the binding and DNA chain cleavage properties of NCS-chrom. NCSi-gb is one of the postactivated mimics of NCS-chrom which is formed under thiol-free conditions and is able to bind to DNA. This report describes the structure refinement of the NCSi-gb-bulge-DNA complex [Stassinopoulos, A., Jie, J., Gao, X., and Goldberg, I. H. (1996) Science 272, 1943-1946] and the NMR characterization of the free bulge-DNA and free NCSi-gb. These results reveal that the formation of the complex involves conformational changes in both the DNA and the ligand molecule. Of mechanistic importance for the NCS-chrom-DNA interaction, the two ring systems of the drug are brought closer to each other in the complex. This conformation correlates well with the previously observed marked enhancement of the formation of a DNA bulge cleaving species in the presence of bulge-DNA sequences, due to the promotion of the intramolecular radical quenching of the activated NCS-chrom. Interestingly, the binding of NCSi-gb promotes the formation of a bulge binding pocket; this was not found in the unbound DNA. NCS-chrom is unique among the enediyne antibiotics in its ability to undergo two different mechanisms of activation to form two different DNA binding and cleaving species. The two corresponding DNA complexes are compared. One, the bulge-DNA binder NCSi-gb, involves the major groove, and the second, the duplex binder NCSi-glu which is generated by glutathione-induced activation, involves the minor groove. Since the two NCS-chrom-related ligand molecules contain some common chemical structural elements, such as the carbohydrate ring, the striking differences in their DNA recognition and chain cleavage specificity provide insights into the fundamental principles of DNA recognition and ligand design.     

Solution structure of a wedge-shaped synthetic molecule at a two-base bulge site in DNA

The solution structure of the complex formed between an oligonucleotide containing a two-base bulge (5'-CACGCAGTTCGGAC.5'-GTCCGATGCGTG) and DDI, a designed synthetic agent, has been elucidated using high-resolution NMR spectroscopy and restrained molecular dynamic simulation. DDI, which has been found to modulate DNA strand slippage synthesis by DNA polymerase I [Kappen, L. S., Xi, Z., Jones, G. B., and Goldberg, I. H. (2003) Biochemistry 42, 2166-2173], is a wedge-shaped spirocyclic molecule whose aglycone structure closely resembles that of the natural product, NCSi-gb, which strongly binds to an oligonucleotide containing a two-base bulge. Changes in chemical shifts of the DNA upon complex formation and intermolecular NOEs between DDI and the bulged DNA duplex indicate that agent specifically binds to the bulge site of DNA. The benzindanone moiety of DDI intercalates via the minor groove into the G7-T8-T9.A20 pocket, which consists of a helical base pair and two unpaired bulge bases, stacking with the G7 and A20 bases. On the other hand, the dihydronaphthalenone and aminoglycoside moieties are positioned in the minor groove. The aminoglycoside, which is attached to spirocyclic ring, aligns along the A20T21G22 sequence of the nonbulged strand, while the dihydronaphthalenone, which is restrained by the spirocyclic structure, is positioned near the G7-T8-T9 bulge site. The aminoglycoside is closely aligned with the dihydronaphthalenone, preventing its intercalation into the bulge site. In the complex, the unpaired purine (G7) is intrahelical and stacks with the intercalating moiety of DDI, whereas the unpaired pyrimidine (T8) is extrahelical. The structure of the complex formed by binding of the synthetic agent to the two-base bulged DNA reveals a binding mode that differs in important details from that of the natural product, explaining the different binding specificity for the bulge sites of DNA. The structure of the DDI-bulged DNA complex provides insight into the structure-binding affinity relationship, providing a rational basis for the design of specific, high-affinity probes of the role of bulged nucleic acid structures in various biological processes.     

Development of new simple molecular probes of DNA bulged structures

NCSi-gb is a neocarzinostatin chromophore (NCS-chrom) metabolite which binds strongly to certain two-base DNA bulges. Compared with previously reported NCSi-gb analogues, a new analogue with a different aminoglycoside position was synthesized, and it showed strong fluorescence and improved binding and sequence selectivity to DNA bulges. The N-dimethylated form of this analogue had a similar binding pattern, and it competitively inhibited bulge-specific cleavage by NCS-chrom.     

Stereochemical control of small molecule binding to bulged DNA: comparison of structures of spirocyclic enantiomer-bulged DNA complexes

The solution structure of the complex formed between an oligonucleotide containing a two-base bulge (5'-CACGCAGTTCGGAC.5'-GTCCGATGCGTG) and ent-DDI, a designed synthetic agent, has been elucidated using high-resolution NMR spectroscopy and restrained molecular dynamic simulation. Ent-DDI is a left-handed wedge-shaped spirocyclic molecule whose aglycone portion is an enantiomer of DDI, which mimics the spirocyclic geometry of the natural product, NCSi-gb, formed by base-catalyzed activation of the enediyne antibiotic neocarzinostatin. The benzindanone moiety of ent-DDI intercalates between the A6.T21 and the T9.A20 base pairs, overlapping with portions of the purine bases; the dihydronaphthalenone moiety is positioned in the minor groove along the G7-T8-T9 bulge sequence; and the aminoglycoside is in the middle of the minor groove, approaching A20 of the nonbulged strand. This alignment of ent-DDI along the DNA helical duplex is in the reverse direction to that of DDI. The aminoglycoside moiety of ent-DDI is positioned in the 3' direction from the bulge region, whereas that of the DDI is positioned in the 5' direction from the same site. This reverse binding orientation within the bulge site is the natural consequence of the opposite handedness imposed by the spirocyclic ring junction and permits the aromatic ring systems of the two spirocyclic enantiomers access to the bulge region. NMR and CD data indicate that the DNA in the DDI-bulged DNA complex undergoes a larger conformational change upon complex formation in comparison to the ent-DDI-bulged DNA, explaining the different binding affinities of the two drugs to the bulged DNA. In addition, there are different placements of the bulge bases in the helical duplex in the two complexes. One bulge base (G7) stacks inside the helix, and the other one (T8) is extrahelical in the DDI-bulged DNA complex, whereas both bulge bases in the ent-DDI-bulged DNA complex prefer extrahelical positions for drug binding. Elucidation of the detailed binding characteristics of the synthetic spirocyclic enantiomers provides a rational basis for the design of stereochemically controlled drugs for bulge binding sites.     

Specific binding of 2-amino-1,8-naphthyridine into a single guanine bulge as evidenced by photooxidation of GG doublet

Photoirradiation of 2-amino-1,8-naphthyridines in the presence of duplex DNA containing the GG doublet opposite a single bulge was examined. After hot piperidine treatment, DNA cleavage was observed preferentially at the GG opposite a single bulge. The cleavage efficiency was highly dependent on the nature of bulged base. The G cleavage at the GG opposite a single G bulge was exceptionally weak, suggesting an intercalative binding of 2-amino-1,8-naphthyridine chromophore into the GG step.     

Interaction of bulged DNA with leucine-containing mimics of NCS-chrom

Synthesis of chiral spirocyclic helical compounds containing leucine that mimic the molecular architecture of the potent DNA bulge binder obtained from the natural product metabolite NCSi-gb has been accomplished. The interaction between the compounds and DNA was studied by circular dichroism (CD) method. The results suggested that the two synthetic diastereoisomers specifically targeted the bulge site of DNA and induced conformational change of bulged DNA greatly.     

Solution structure of a designed spirocyclic helical ligand binding at a two-base bulge site in DNA

The solution structure of the complex formed between an oligodeoxynucleotide containing a two-base bulge (5'-CCATCGTCTACCTTTGGTAGGATGG) and SCA-alpha2, a designed spirocyclic helical molecule, has been elucidated. SCA-alpha2, a close mimic of the metabolite, NCSi-gb, of the DNA bulge-specific enediyne antibiotic neocarzinostatin, differs in possessing a more stable spirocyclic ring system and in lacking certain bulky groupings that compromise bulged DNA binding. This study provides a detailed comparison of the binding modes of the two complexes and provides new insights into the importance of shape and space, as opposed to simple nucleotide sequence, in complex formation at the bulge site. The two rigidly held aromatic rings of SCA-alpha2 form a right-handed helical molecular wedge that specifically penetrates the bulge-binding pocket and immobilizes the two bulge residues (GT), which point toward the minor groove, rather than the major groove as in the NCSi-gb.bulged DNA complex. The ligand aromatic ring systems stack on the DNA bulge-flanking base pairs that define the long sides of the triangular prism binding pocket. Like NCSi-gb, SCA-alpha2 possesses the natural N-methylfuranose moiety, alpha-linked to the benzindanol (BI) moiety. The amino sugar anchors in the major groove of the DNA and points toward the 3'-bulge-flanking base pair. Lacking the bulky cyclocarbonate of NCSi-gb, the SCA-alpha2.bulged DNA complex has a much less twisted and buckled 3'-bulge-flanking base pair (dG20.dC8), and the G20 residue stacks directly above the BI ring platform. Also, the absence of the methyl group and the free rotation of the methoxy group on the dihydronaphthanone (NA) moiety of SCA-alpha2 allow better stacking geometry of the NA ring above the 5'-bulge-flanking dG21.dC5 base pair. These and other considerations help to explain why NCSi-gb binds very poorly to bulged RNA and are consistent with the recent observation of good binding with SCA-alpha2. Thus, although the two complexes resemble each other closely, they differ in important local environmental details. SCA-alpha2 has a better hand-in-glove fit at the bulge site, making it an ideal platform for the placement of moieties that can react covalently with the DNA and for generating congeners specific for bulges in RNA.     

Recognition of bulged DNA by a neocarzinostatin product via an induced fit mechanism

The binding of the wedge-shaped isostructural analogue of the biradical species of the chromphore of antitumor antibiotic neocarzinostatin to sequence-specific bulged DNAs results in alterations in ellipticity of the DNAs. Circular dichroism (CD) spectroscopic results suggest that the drug specifically recognizes bulges of DNA via a combination of conformational selection and induced fit, not by binding to a preorganized site. Analysis of circular dichroism spectra indicates that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the bulge cavity. The effective recognition of the bulge site on duplex DNA appears to depend to a significant extent on the bent groove space being flexible enough to be able to adopt the geometrically optimal conformation compatible with the wedge-shaped drug molecule, rather than involving 'lock and key' recognition. The spectroscopic results indicate a change of DNA conformation, consistent with an allosteric binding model. Spectroscopic studies with various bulged DNAs also reveal that the binding strength directly correlates with the stability of the bulge structures.     

X-ray crystallographic observation of "in-line" and "adjacent" conformations in a bulged self-cleaving RNA/DNA hybrid

The RNA strand in an RNA/DNA duplex with unpaired ribonucleotides can undergo self-cleavage at bulge sites in the presence of a variety of divalent metal ions (Hüsken et al., Biochemistry, 1996, 35:16591-16600). Transesterification proceeds via an in-line mechanism, with the 2'-OH of the bulged nucleotide attacking the 3'-adjacent phosphate group. The site-specificity of the reaction is most likely a consequence of the greater local conformational freedom of the RNA backbone in the bulge region. A standard A-form backbone geometry prohibits formation of an in-line arrangement between 2'-oxygen and phosphate. However, the backbone in the region of an unpaired nucleotide appears to be conducive to an in-line approach. Therefore, the bulge-mediated phosphoryl transfer reaction represents one of the simplest RNA self-cleavage systems. Here we focus on the conformational features of the RNA that underlie site-specific cleavage. The structures of an RNA/DNA duplex with single ribo-adenosyl bulges were analyzed in two crystal forms, permitting observation of 10 individual conformations of the RNA bulge moiety. The bulge geometries cover a range of relative arrangements between the 2'-oxygen of the bulged nucleotide and the P-O5' bond (including adjacent and near in-line) and give a detailed picture of the conformational changes necessary to line up the 2'-OH nucleophile and scissile bond. Although metal ions are of crucial importance in the catalysis of analogous cleavage reactions by ribozymes, it is clear that local strain or conformational flexibility in the RNA also affect cleavage selectivity and rate (Soukup & Breaker, RNA, 1999, 5:1308-1325). The geometries of the RNA bulges frozen out in the crystals provide snapshots along the reaction pathway prior to the transition state of the phosphoryl transfer reaction.     

In vitro expansion of DNA triplet repeats with bulge binders and different DNA polymerases

The expansion of DNA repeat sequences is associated with many genetic diseases in humans. Simple bulge DNA structures have been implicated as intermediates in DNA slippage within the DNA repeat regions. To probe the possible role of bulged structures in DNA slippage, we designed and synthesized a pair of simple chiral spirocyclic compounds [Xi Z, Ouyang D & Mu HT (2006) Bioorg Med Chem Lett 16, 1180-1184], DDI-1A and DDI-1B, which mimic the molecular architecture of the enediyne antitumor antibiotic neocarzinostatin chromophore. Both compounds strongly stimulated slippage in various DNA repeats in vitro. Enhanced slippage synthesis was found to be synchronous for primer and template. CD spectra and UV thermal stability studies supported the idea that DDI-1A and DDI-1B exhibited selective binding to the DNA bulge and induced a significant conformational change in bulge DNA. The proposed mechanism for the observed in vitro expansion of long DNA is discussed.     

DNA bending by the bulge defect

Comparative gel electrophoresis measurements were used to characterize DNA bending in molecules containing an extra adenosine on one strand, the so-called bulge defect. We used oligomers containing A6 tracts separated from the bulged base by varying numbers of nucleotides to determine the direction and magnitude of the bulge bend. Helix unwinding by the bulge was determined from the electrophoretic anomaly as a function of the size of the repeated monomers. We conclude that the bulge bend is 21 degrees +/- 3 degrees, primarily in the direction of tilt away from the bulged base. The total helical advance of the DNA at the bulge site is smaller than would be the case if the complementary T were present, corresponding to an unwinding by 25 degrees +/- 6 degrees. These values are in good agreement with the results of NMR and energy minimization studies of the bulged base in double-helical deoxyoligonucleotides [Woodson, S. A., & Crothers, D.M. (1988) Biochemistry 27, 3130-3141]     

Probing DNA single strands for single-base bulges with neocarzinostatin chromophore.

Neocarzinostatin chromophore (NCS-Chrom) induces strong cleavage at a single site (C3) in the single-stranded and 5' (32)P-end-labeled 13-mer GCCAGATTTGAGC in a reaction dependent on a thiol. By contrast, in the duplex form of the same 13-mer, strand cleavage occurs only at the T and A residues, and C3 is not cleaved. To determine the minimal structural requirement(s) for C3 cleavage in the single-stranded oligomer, several deletions and mutations were made in the 13-mer. A 10-mer (GCCAGAGAGC) derived from the 13-mer by deletion of the three T residues was also cleaved exclusively at C3 by NCS-Chrom, generating fragments having 5' phosphate ends. That the cleavage at C3 is initiated by abstraction of its 5' hydrogen is confirmed in experiments using 3' (32)P-end-labeled 10-mer. The competent 13-mer and 10-mer were assigned hairpin structures with a stem loop and a single bulged out A base, placing C3 across from and 3' to the bulge. Removal of the bulged A base from the 13-mer and the 10-mer resulted in complete loss of cutting activity, proving that it is the essential determinant in competent substrates. Studies of thiol post-activated NCS-Chrom binding to the DNA oligomers show that the drug binds to the bulge-containing 13-mer (K(d) = 0.78 microM) and the 10-mer (K(d) = 1.11 microM), much more strongly than to the 12-mer (K(d) = 20 microM) and the 9-mer (K(d) = 41 microM), lacking the single-base bulge. A mutually induced-fit between NCS-Chrom and the oligomer resulting in optimal stabilization of the drug-DNA complex is proposed to account for the site-specific cleavage at C3. These studies establish the usefulness of NCS-Chrom as a probe for single-base bulges in DNA.     

DNA damage by thiol-activated neocarzinostatin chromophore at bulged sites

Bulge structures in nucleic acids are of general biological significance and are potential targets for therapeutic drugs. It has been shown in a previous study that thiol-activated neocarzinostatin chromophore is able to cleave duplex DNA selectively at a position opposite a single unpaired cytosine or thymine base on the 3' side. In this work, we studied in greater detail the nature of this type of cleavage and the basis for the selectivity of the bulge site cleavage over the usual strand cleavage at a T site in the duplex region by using duplexes containing an internal control and a bulge, which is composed of different types and number of bases. Experimental results indicated that the bulge-induced cleavage is initiated by 5' hydrogen abstraction and is greatly affected by the base composition of the bulge. A single-base bulge, especially when containing a purine, yields higher efficiency and greater selectivity for the bulge-induced cleavage. In particular, a single adenine base gives rise to the highest cleavage yield and provides over 20 times greater selectivity for cleavage at the bulge site compared with the internal control site in duplexes. The binding dissociation constants of postactivated drug for a stem-loop structure containing a one- or two-base bulge in the stem, measured by fluorescence quenching, show that the binding is about 3-4 times stronger for bulge-containing duplexes than for perfect hairpin duplexes. For RNA.DNA hybrid duplexes, where the DNA is the target strand and the RNA is the bulge-containing strand, bulge-induced cleavage was observed, although at low yield. On the other hand, when RNA is the nonbulge strand, no bulge-induced cleavage was found. When the reaction is performed in the absence of oxygen, the major product is a covalent adduct, and it is at the same location as the cleavage site under aerobic conditions.     

Selective interaction between tylophorine B and bulged DNA

Tylophorine B exhibits pronounced cytotoxicity and antitumor activity. In order to survey the structure selectivity to DNA afforded by tylophorine B, we have synthesized a variety of duplex, bulge- and hairpin-containing oligodeoxyribonucleotides. Their binding to tylophorine B has been assayed by fluorescence spectroscopy and thermal melting experiments. The results indicate that oligonucleotides interact with tylophorine B at submicromolar concentration, and the affinity for DNA bulge is optimal (with Kd of 0.018 microM). In addition, the bulged hairpin oligonucleotides are stabilized by binding to tylophorine B. These findings may shed some light on tylophorine B's mode of action in biological systems and result in the rational design of sequence-specific DNA binding molecules.     

Thermodynamic examination of trinucleotide bulged RNA in the context of HIV-1 TAR RNA

RNA structures contain many bulges and loops that are expected to be sites for inter- and intra-molecular interactions. Nucleotides in the bulge are expected to influence the structure and recognition of RNA. The same stability is assigned to all trinucleotide bulged RNA in the current secondary structure prediction models. In this study thermal denaturation experiments were performed on four trinucleotide bulged RNA, in the context of HIV-1 TAR RNA, to determine whether the bulge sequence affects RNA stability and its divalent ion interactions. Cytosine-rich bulged RNA were more stable than uracil-rich bulged RNA in 1 M KCl. Interactions of divalent ions were more favorable with uracil-rich bulged RNA by ~2 kcal/mol over cytosine-rich bulged RNA. The UCU-TAR RNA (wild type) is stabilized by 1.7 kcal/mol in 9.5 mM Ca²⁺ as compared with 1 M KCl, whereas no additional gain in stability is measured for CCC-TAR RNA. These results have implications for base substitution experiments traditionally employed to identify metal ion binding sites. To our knowledge, this is the first systematic study to quantify the effect of small sequence changes on RNA stability upon interactions with divalent ions.