Solution equilibria of cytosine- and guanine-rich sequences near the promoter region of the n-myc gene that contain stable hairpins within lateral loops

Background

Cytosine- and guanine-rich regions of DNA are capable of forming complex structures named i-motifs and G-quadruplexes, respectively. In the present study the solution equilibria at nearly physiological conditions of a 34-base long cytosine-rich sequence and its complementary guanine-rich strand corresponding to the first intron of the n-myc gene were studied. Both sequences, not yet studied, contain a 12-base tract capable of forming stable hairpins inside the i-motif and G-quadruplex structures, respectively.

Methods

Spectroscopic, mass spectrometry and separation techniques, as well as multivariate data analysis methods, were used to unravel the species and conformations present.

Results

The cytosine-rich sequence forms two i-motifs that differ in the protonation of bases located in the loops. A stable Watson–Crick hairpin is formed by the bases in the first loop, stabilizing the i-motif structure. The guanine-rich sequence adopts a parallel G-quadruplex structure that is stable throughout the pH range 3–7, despite the protonation of cytosine and adenine bases at lower pH values. The presence of G-quadruplex aggregates was confirmed using separation techniques. When mixed, G-quadruplex and i-motif coexist with the Watson–Crick duplex across a pH range from approximately 3.0 to 6.5.

Conclusions

Two cytosine- and guanine-rich sequences in n-myc gene may form stable i-motif and G-quadruplex structures even in the presence of long loops. pH modulates the equilibria involving the intramolecular structures and the intermolecular Watson–Crick duplex.

General significance

Watson–Crick hairpins located in the intramolecular G-quadruplexes and i-motifs in the promoter regions of oncogenes could play a role in stabilizing these structures.     

A remarkably stable kissing-loop interaction defines substrate recognition by the Neurospora Varkud Satellite ribozyme

Kissing loops are tertiary structure elements that often play key roles in functional RNAs. In the Neurospora VS ribozyme, a kissing-loop interaction between the stem–loop I (SLI) substrate and stem–loop V (SLV) of the catalytic domain is known to play an important role in substrate recognition. In addition, this I/V kissing-loop interaction is associated with a helix shift in SLI that activates the substrate for catalysis. To better understand the role of this kissing-loop interaction in substrate recognition and activation by the VS ribozyme, we performed a thermodynamic characterization by isothermal titration calorimetry using isolated SLI and SLV stem–loops. We demonstrate that preshifted SLI variants have higher affinity for SLV than shiftable SLI variants, with an energetic cost of 1.8–3 kcal/mol for the helix shift in SLI. The affinity of the preshifted SLI for SLV is remarkably high, the interaction being more stable by 7–8 kcal/mol than predicted for a comparable duplex containing three Watson–Crick base pairs. The structural basis of this remarkable stability is discussed in light of previous NMR studies. Comparative thermodynamic studies reveal that kissing-loop complexes containing 6–7 Watson–Crick base pairs are as stable as predicted from comparable RNA duplexes; however, those with 2–3 Watson–Crick base pairs are more stable than predicted. Interestingly, the stability of SLI/ribozyme complexes is similar to that of SLI/SLV complexes. Thus, the I/V kissing loop interaction represents the predominant energetic contribution to substrate recognition by the trans-cleaving VS ribozyme.     

Non-nearest-neighbor dependence of stability for group III RNA single nucleotide bulge loops

Thirty-five RNA duplexes containing single nucleotide bulge loops were optically melted and the thermodynamic parameters for each duplex determined. The bulge loops were of the group III variety, where the bulged nucleotide is either a AG/U or CU/G, leading to ambiguity to the exact position and identity of the bulge. All possible group III bulge loops with Watson–Crick nearest-neighbors were examined. The data were used to develop a model to predict the free energy of an RNA duplex containing a group III single nucleotide bulge loop. The destabilization of the duplex by the group III bulge could be modeled so that the bulge nucleotide leads to the formation of the Watson–Crick base pair rather than the wobble base pair. The destabilization of an RNA duplex caused by the insertion of a group III bulge is primarily dependent upon non-nearest-neighbor interactions and was shown to be dependent upon the stability of second least stable stem of the duplex. In-line structure probing of group III bulge loops embedded in a hairpin indicated that the bulged nucleotide is the one positioned further from the hairpin loop irrespective of whether the resulting stem formed a Watson–Crick or wobble base pair. Fourteen RNA hairpins containing group III bulge loops, either 3′ or 5′ of the hairpin loop, were optically melted and the thermodynamic parameters determined. The model developed to predict the influence of group III bulge loops on the stability of duplex formation was extended to predict the influence of bulge loops on hairpin stability.     

A model capturing novel strand symmetries in bacterial DNA

Chargaff’s second parity rule for short oligonucleotides states that the frequency of any short nucleotide sequence on a strand is approximately equal to the frequency of its reverse complement on the same strand. Recent studies have shown that, with the exception of organellar DNA, this parity rule generally holds for double-stranded DNA genomes and fails to hold for single-stranded genomes. While Chargaff’s first parity rule is fully explained by the Watson–Crick pairing in the DNA double helix, a definitive explanation for the second parity rule has not yet been determined. In this work, we propose a model based on a hidden Markov process for approximating the distributional structure of primitive DNA sequences. Then, we use the model to provide another possible theoretical explanation for Chargaff’s second parity rule, and to predict novel distributional aspects of bacterial DNA sequences.     

A thermodynamic overview of naturally occurring intramolecular DNA quadruplexes

Loop length and its composition are important for the structural and functional versatility of quadruplexes. To date studies on the loops have mainly concerned model sequences compared with naturally occurring quadruplex sequences which have diverse loop lengths and compositions. Herein, we have characterized 36 quadruplex-forming sequences from the promoter regions of various proto-oncogenes using CD, UV and native gel electrophoresis. We examined folding topologies and determined the thermodynamic profile for quadruplexes varying in total loop length (5–18 bases) and composition. We found that naturally occurring quadruplexes have variable thermodynamic stabilities (ΔG₃₇) ranging from −1.7 to −15.6 kcal/mol. Overall, our results suggest that both loop length and its composition affect quadruplex structure and thermodynamics, thus making it difficult to draw generalized correlations between loop length and thermodynamic stability. Additionally, we compared the thermodynamic stability of quadruplexes and their respective duplexes to understand quadruplex–duplex competition. Our findings invoke a discussion on whether biological function is associated with quadruplexes with lower thermodynamic stability which undergo facile formation and disruption, or by quadruplexes with high thermodynamic stability.     

Isostericity and tautomerism of base pairs in nucleic acids

The natural bases of nucleic acids form a great variety of base pairs with at least two hydrogen bonds between them. They are classified in twelve main families, with the Watson–Crick family being one of them. In a given family, some of the base pairs are isosteric between them, meaning that the positions and the distances between the C1′ carbon atoms are very similar. The isostericity of Watson–Crick pairs between the complementary bases forms the basis of RNA helices and of the resulting RNA secondary structure. Several defined suites of non-Watson–Crick base pairs assemble into RNA modules that form recurrent, rather regular, building blocks of the tertiary architecture of folded RNAs. RNA modules are intrinsic to RNA architecture are therefore disconnected from a biological function specifically attached to a RNA sequence. RNA modules occur in all kingdoms of life and in structured RNAs with diverse functions. Because of chemical and geometrical constraints, isostericity between non-Watson–Crick pairs is restricted and this leads to higher sequence conservation in RNA modules with, consequently, greater difficulties in extracting 3D information from sequence analysis. Nucleic acid helices have to be recognised in several biological processes like replication or translational decoding. In polymerases and the ribosomal decoding site, the recognition occurs on the minor groove sides of the helical fragments. With the use of alternative conformations, protonated or tautomeric forms of the bases, some base pairs with Watson–Crick-like geometries can form and be stabilized. Several of these pairs with Watson–Crick-like geometries extend the concept of isostericity beyond the number of isosteric pairs formed between complementary bases. These observations set therefore limits and constraints to geometric selection in molecular recognition of complementary Watson–Crick pairs for fidelity in replication and translation processes.     

Molecular dynamics simulations and coupled nucleotide substitution experiments indicate the nature of A·A base pairing and a putative structure of the coralyne-induced homo-adenine duplex

Coralyne is an alkaloid drug that binds homo-adenine DNA (and RNA) oligonucleotides more tightly than it does Watson–Crick DNA. Hud’s laboratory has shown that poly(dA) in the presence of coralyne forms an anti-parallel duplex, however attempts to determine the structure by NMR spectroscopy and X-ray crystallography have been unsuccessful. Assuming adenine–adenine hydrogen bonding between the two poly(dA) strands, we constructed 40 hypothetical homo-(dA) anti-parallel duplexes and docked coralyne into the six most favorable duplex structures. The two most stable structures had trans glycosidic bonds, but distinct pairing geometries, i.e. either Watson–Crick Hoogsteen (transWH) or Watson–Crick Watson–Crick (transWW) with stability of transWH > transWW. To narrow down the possibilities, 7-deaza adenine base substitutions (dA→7) were engineered into homo-(dA) sequences. These substitutions significantly reduced the thermal stability of the coralyne-induced homo-(dA) structure. These experiments strongly suggest the involvement of N7 in the coralyne-induced A·A base pairs. Moreover, due to the differential effect on melting as a function of the location of the dA→7 mutations, these results are consistent with the N1–N7 base pairing of the transWH pairs. Together, the simulation and base substitution experiments predict that the coralyne-induced homo-(dA) duplex structure adopts the transWH geometry.     

Nearest neighbor parameters for Watson-Crick complementary heteroduplexes formed between 2'-O-methyl RNA and RNA oligonucleotides

Results from optical melting studies of Watson–Crick complementary heteroduplexes formed between 2′-O-methyl RNA and RNA oligonucleotides are used to determine nearest neighbor thermodynamic parameters for predicting the stabilities of such duplexes. The results are consistent with the physical model assumed by the individual nearest neighbor-hydrogen bonding model, which contains terms for helix initiation, base pair stacking and base pair composition. The sequence dependence is similar to that for Watson–Crick complementary RNA/RNA duplexes, which suggests that the sequence dependence may also be similar to that for other backbones that favor A-form RNA conformations.     

DNA multiplex hybridization on microarrays and thermodynamic stability in solution: a direct comparison

Hybridization intensities of 30 distinct short duplex DNAs measured on spotted microarrays, were directly compared with thermodynamic stabilities measured in solution. DNA sequences were designed to promote formation of perfect match, or hybrid duplexes containing tandem mismatches. Thermodynamic parameters ΔH°, ΔS° and ΔG° of melting transitions in solution were evaluated directly using differential scanning calorimetry. Quantitative comparison with results from 63 multiplex microarray hybridization experiments provided a linear relationship for perfect match and most mismatch duplexes. Examination of outliers suggests that both duplex length and relative position of tandem mismatches could be important factors contributing to observed deviations from linearity. A detailed comparison of measured thermodynamic parameters with those calculated using the nearest-neighbor model was performed. Analysis revealed the nearest-neighbor model generally predicts mismatch duplexes to be less stable than experimentally observed. Results also show the relative stability of a tandem mismatch is highly dependent on the identity of the flanking Watson–Crick (w/c) base pairs. Thus, specifying the stability contribution of a tandem mismatch requires consideration of the sequence identity of at least four base pair units (tandem mismatch and flanking w/c base pairs). These observations underscore the need for rigorous evaluation of thermodynamic parameters describing tandem mismatch stability.     

RNA structure and dynamics: A base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson–Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson–Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.     

Targeting the G-quadruplex-forming region near the P1 promoter in the human BCL-2 gene with the cationic porphyrin TMPyP4 and with the complementary C-rich strand

The B-cell lymphoma-2 (bcl-2) gene contains a region that has been implicated in the regulation of bcl-2 gene expression. This region can form G-quadruplex structures in solution [J.X. Dai, T.S. Dexheimer, D. Chen, M. Carver, A. Ambrus, R.A. Jones, D.Z. Yang, An intramolecular G-quadruplex structure with mixed parallel/antiparallel G-strands formed in the human BCL-2 promoter region in solution, J. Am. Chem. Soc. 128 (2006) 1096–1098.]. Here, we examined the acid–base and conformational equilibria of this G-quadruplex-forming region (BCL2G), as well as its interaction with both the porphyrin TMPyP4 and with the complementary C-rich strand. We used molecular absorption and circular dichroism techniques, in tandem with multivariate analysis tools. The results revealed the formation of an interaction complex BCL2G:TMPyP4 with a stoichiometry of 1:2 and an equilibrium constant equal to 5.0 (±2.3) × 10¹³ M⁻². Addition of the complementary C-rich strand to BCL2G induces the predominant formation of the Watson–Crick double-helix with an equilibrium constant equal to 10^7.7 M⁻¹ (at pH 7.1). Finally, the pH-induced formation of quadruplex structures from the Watson–Crick double-helix is characterized.     

NMR structures of double loops of an RNA aptamer against mammalian initiation factor 4A

A high affinity RNA aptamer (APT58, 58 nt long) against mammalian initiation factor 4A (eIF4A) requires nearly its entire nucleotide sequence for efficient binding. Since splitting either APT58 or eIF4A into two domains diminishes the affinity for each other, it is suggested that multiple interactions or a global interaction between the two molecules accounts for the high affinity. To understand the structural basis of APT58's global recognition of eIF4A, we determined the solution structure of two essential nucleotide loops (AUCGCA and ACAUAGA) within the aptamer using NMR spectroscopy. The AUCGCA loop is stabilized by a U-turn motif and contains a non-canonical A:A base pair (the single hydrogen bond mismatch: Hoogsteen/Sugar-edge). On the other hand, the ACAUAGA loop is stabilized by an AUA tri-nucleotide loop motif and contains the other type of A:A base pair (single hydrogen bond mismatch: Watson–Crick/Watson–Crick). Considering the known structural and functional properties of APT58, we propose that the AUCGCA loop is directly involved in the interaction with eIF4A, while the flexibility of the ACAUAGA loop is important to support this interaction. The Watson–Crick edges of C7 and C9 in the AUCGCA loop may directly interact with eIF4A.     

Insight into the structural organization of the omega leader of TMV RNA: the role of various regions of the sequence in the formation of a compact structure of the omega RNA

The 5′-untranslated sequence of tobacco mosaic virus RNA – the so called omega leader – is a well-known translational enhancer. The structure of the omega RNA has unusual features. Despite the absence of extensive secondary structure of the Watson–Crick type, the omega RNA possesses a stable compact conformation. The central part of the omega sequence contains many CAA repeats and is flanked by U-rich regions. In this work we synthesized the polyribonucleotides containing modified omega sequences, and studied them using analytical ultracentrifugation and thermal melting techniques. It was demonstrated that changes made in both the central and the 3′-proximal part of the sequence led to a strong destabilization of the omega RNA structure. We conclude that the regular (CAA)_n core region and the 3′-proximal AU-rich region of the omega RNA interact with each other and contribute together to the formation of a stable tertiary structure.     

Structural characterization of a new subclass of panicum mosaic virus-like 3′ cap-independent translation enhancer

Canonical eukaryotic mRNA translation requires 5′cap recognition by initiation factor 4E (eIF4E). In contrast, many positive-strand RNA virus genomes lack a 5′cap and promote translation by non-canonical mechanisms. Among plant viruses, PTEs are a major class of cap-independent translation enhancers located in/near the 3′UTR that recruit eIF4E to greatly enhance viral translation. Previous work proposed a single form of PTE characterized by a Y-shaped secondary structure with two terminal stem-loops (SL1 and SL2) atop a supporting stem containing a large, G-rich asymmetric loop that forms an essential pseudoknot (PK) involving C/U residues located between SL1 and SL2. We found that PTEs with less than three consecutive cytidylates available for PK formation have an upstream stem-loop that forms a kissing loop interaction with the apical loop of SL2, important for formation/stabilization of PK. PKs found in both subclasses of PTE assume a specific conformation with a hyperreactive guanylate (G*) in SHAPE structure probing, previously found critical for binding eIF4E. While PTE PKs were proposed to be formed by Watson–Crick base-pairing, alternative chemical probing and 3D modeling indicate that the Watson–Crick faces of G* and an adjacent guanylate have high solvent accessibilities. Thus, PTE PKs are likely composed primarily of non-canonical interactions.     

Integrity of duplex structures without hydrogen bonding: DNA with pyrene paired at abasic sites

DNA polymerases specifically insert the hydrophobic pyrene deoxynucleotide (P) opposite tetrahydrofuran (F), an stable abasic site analog, and DNA duplexes containing this non-hydrogen-bonded pair possess a high degree of thermodynamic stability. These observations support the hypothesis that steric complementarity and stacking interactions may be sufficient for maintaining stability of DNA structure and specificity of DNA replication, even in the absence of hydrogen bonds across the base pair. Here we report the NMR characterization and structure determination of two DNA molecules containing pyrene residues. The first is a 13mer duplex with a pyrene·tetrahydrofuran pair (P·F pair) at the ninth position and the second mimics a replication intermediate right after incorporation of a pyrene nucleoside opposite an abasic site. Our data indicate that both molecules adopt right-handed helical conformations with Watson– Crick alignments for all canonical base pairs. The pyrene ring stays inside the helix close to its baseless partner in both molecules. The single-stranded region of the replication intermediate folds back over the opposing strand, sheltering the hydrophobic pyrene moiety from water exposure. The results support the idea that the stability and replication of a P·F pair is due to its ability to mimic Watson–Crick structure.     

Recognition of T*G mismatched base pairs in DNA by stacked imidazole-containing polyamides: surface plasmon resonance and circular dichroism studies

An imidazole-containing polyamide trimer, f-ImImIm, where f is a formamido group, was recently found using NMR methods to recognize T·G mismatched base pairs. In order to characterize in detail the T·G recognition affinity and specificity of imidazole-containing polyamides, f-ImIm, f-ImImIm and f-PyImIm were synthesized. The kinetics and thermodynamics for the polyamides binding to Watson–Crick and mismatched (containing one or two T·G, A·G or G·G mismatched base pairs) hairpin oligonucleotides were determined by surface plasmon resonance and circular dichroism (CD) methods. f-ImImIm binds significantly more strongly to the T·G mismatch-containing oligonucleotides than to the sequences with other mismatched or with Watson–Crick base pairs. Compared with the Watson–Crick CCGG sequence, f-ImImIm associates more slowly with DNAs containing T·G mismatches in place of one or two C·G base pairs and, more importantly, the dissociation rate from the T·G oligonucleotides is very slow (small k_d). These results clearly demonstrate the binding selectivity and enhanced affinity of side-by-side imidazole/imidazole pairings for T·G mismatches and show that the affinity and specificity increase arise from much lower k_d values with the T·G mismatched duplexes. CD titration studies of f-ImImIm complexes with T·G mismatched sequences produce strong induced bands at ~330 nm with clear isodichroic points, in support of a single minor groove complex. CD DNA bands suggest that the complexes remain in the B conformation.     

Determination of thermodynamic parameters for HIV DIS type loop-loop kissing complexes

The HIV-1 type dimerization initiation signal (DIS) loop was used as a starting point for the analysis of the stability of Watson–Crick (WC) base pairs in a tertiary structure context. We used ultraviolet melting to determine thermodynamic parameters for loop–loop tertiary interactions and compared them with regular secondary structure RNA helices of the same sequences. In 1 M Na⁺ the loop–loop interaction of a HIV-1 DIS type pairing is 4 kcal/mol more stable than its sequence in an equivalent regular and isolated RNA helix. This difference is constant and sequence independent, suggesting that the rules governing the stability of WC base pairs in the secondary structure context are also valid for WC base pairs in the tertiary structure context. Moreover, the effect of ion concentration on the stability of loop–loop tertiary interactions differs considerably from that of regular RNA helices. The stabilization by Na⁺ and Mg²⁺ is significantly greater if the base pairing occurs within the context of a loop–loop interaction. The dependence of the structural stability on salt concentration was defined via the slope of a T_m/log [ion] plot. The short base-paired helices are stabilized by 8°C/log [Mg²⁺] or 11°C/log [Na⁺], whereas base-paired helices forming tertiary loop–loop interactions are stabilized by 16°C/log [Mg²⁺] and 26°C/log [Na⁺]. The different dependence on ionic strength that is observed might reflect the contribution of specific divalent ion binding to the preformation of the hairpin loops poised for the tertiary kissing loop–loop contacts.     

Triple helical DNA in a duplex context and base pair opening

It is fundamental to explore in atomic detail the behavior of DNA triple helices as a means to understand the role they might play in vivo and to better engineer their use in genetic technologies, such as antigene therapy. To this aim we have performed atomistic simulations of a purine-rich antiparallel triple helix stretch of 10 base triplets flanked by canonical Watson–Crick double helices. At the same time we have explored the thermodynamic behavior of a flipping Watson–Crick base pair in the context of the triple and double helix. The third strand can be accommodated in a B-like duplex conformation. Upon binding, the double helix changes shape, and becomes more rigid. The triple-helical region increases its major groove width mainly by oversliding in the negative direction. The resulting conformations are somewhere between the A and B conformations with base pairs remaining almost perpendicular to the helical axis. The neighboring duplex regions maintain a B DNA conformation. Base pair opening in the duplex regions is more probable than in the triplex and binding of the Hoogsteen strand does not influence base pair breathing in the neighboring duplex region.     

RNA aptamers that specifically bind to a 16S ribosomal RNA decoding region construct

RNA–RNA recognition is a critical process in controlling many key biological events, such as translation and ribozyme functions. The recognition process governing RNA–RNA interactions can involve complementary Watson–Crick (WC) base pair binding, or can involve binding through tertiary structural interaction. Hence, it is of interest to determine which of the RNA–RNA binding events might emerge through an in vitro selection process. The A-site of the 16S rRNA decoding region was chosen as the target, both because it possesses several different RNA structural motifs, and because it is the rRNA site where codon/anticodon recognition occurs requiring recognition of both mRNA and tRNA. It is shown here that a single family of RNA molecules can be readily selected from two different sizes of RNA library. The tightest binding aptamer to the A-site 16S rRNA construct, 109.2-3, has its consensus sequences confined to a stem–loop region, which contains three nucleotides complementary to three of the four nucleotides in the stem–loop region of the A-site 16S rRNA. Point mutations on each of the three nucleotides on the stem–loop of the aptamer abolish its binding capacity. These studies suggest that the RNA aptamer 109.2-3 interacts with the simple 27 nt A-site decoding region of 16S rRNA through their respective stem–loops. The most probable mode of interaction is through complementary WC base pairing, commonly referred to as a loop–loop ‘kissing’ motif. High affinity binding to the other structural motifs in the decoding region were not observed.