Duplex Stem Research Articles

Inhibition of Ricin A-Chain with Pyrrolidine Mimics of the Oxacarbenium Ion Transition State

Ricin A-chain (RTA) catalyzes the hydrolytic depurination of a specific adenosine at position 4324 of 28S rRNA. Kinetic isotope effects on the hydrolysis of a small 10mer stem-tetraloop oligonucleotide substrate established the mechanism of the reaction as D(N)*A(N), involving an oxacarbenium ion intermediate in a highly dissociative transition state. An inhibitor with a protonated 1,4-dideoxy-1,4-imino-D-ribitol moiety, a 4-azasugar mimic, at the depurination site in the tetraloop of a 14mer oligonucleotide with a 5 bp duplex stem structure had previously been shown to bind to RTA with a K(d) of 480 nM, which improved to 12 nM upon addition of adenine. Second-generation stem-tetraloop inhibitors have been synthesized that incorporate a methylene bridge between the nitrogen of a 1-azasugar mimic, namely, (3S,4R)-3-hydroxy-4-(hydroxymethyl)pyrrolidine, and substituents, including phenyl, 8-aza-9-deazaadenyl, and 9-deazaadenyl groups, that mimic the activated leaving group at the transition state. The values for the dissociation constants (K(i)) for these were 99 nM for the phenyl 10mer, 163 and 94 nM for the 8-aza-9-deazaadenyl 10- and 14mers, respectively, and 280 nM for the 9-deazaadenyl 14mer. All of these compounds are among the tightest binding molecules known for RTA. A related phenyl-substituted inhibitor with a deoxyguanosine on the 5'-side of the depurination site was also synthesized on the basis of stem-loop substrate specificity studies. This molecule binds with a K(i) of 26 nM and is the tightest binding "one-piece" inhibitor. 8-Aza-9-deaza- and 9-deazaadenyl substituents provide an increased pK(a) at N7, a protonation site en route to the transition state. The binding of these inhibitors is not improved relative to the binding of their phenyl counterpart, however, suggesting that RTA might also employ protonation at N1 and N3 of the adenine moiety to activate the substrate during catalysis. Studies with methylated adenines support this argument. That the various stem-loop inhibitors have similar potencies suggests that an optimal one-piece inhibitor remains to be identified. The second-generation inhibitors described here incorporate ribose mimics missing the 2-hydroxy group. On the basis of inhibition data and substrate specificity studies, the 2'-hydroxyl group at the depurination site seems to be critical for recruitment as well as catalysis by RTA.

GA and AG Sequences of DNA React with Cisplatin at Comparable Rates.

The sequence selectivity of the antitumor drug cisplatin (cis-[PtCl(2)(NH(3))(2)] (1)) between the 5'-AG-3' and 5'-GA-3' sites of DNA has been a matter of discussion for more than twenty years. In this work, we compared the reactivity of GA and AG sequences of DNA towards the aquated forms of cisplatin (cis-[PtCl(NH(3))(2)(H(2)O)](+) (2), cis-[Pt(NH(3))(2)(H(2)O)(2)](2+) (3), and cis-[Pt(OH)(NH(3))(2)(H(2)O)](+) (4)) using two sets of experiments. In the first, we investigated a DNA hairpin, whose duplex stem contained a TGAT sequence as the single reactive site, and determined the individual rate constants of platination with 2 and 3 for G and A in acidic solution. The rate constants at 20 degrees C in 0.1M NaClO(4) at pH 4.5+/-0.1 were 0.09(4) M(-1)s(-1) (G) and 0.11(3) M(-1)s(-1) (A) for 2, and 9.6(1) M(-1)s(-1) (G) and 1.7(1) M(-1)s(-1) (A) for 3. These values are similar to those obtained previously for an analogous hairpin that contained a TAGT sequence. The monoadducts formed with 2 by both GA purines are extremely long-lived, partly as a result of the slow hydrolysis of the chloro monoadduct at A, and partly because of the very low chelation rate (1.4 x 10(-5)s(-1) at 20 degrees C) of the aqua monoadduct on the guanine. In the second set of experiments, we incubated pure or enriched samples of 1, 2, 3, or 4 for 18-64 h at 25 degrees C with a 19 base pair (bp) DNA duplex, whose radiolabeled top strand contained one GA and one AG sequence as the only reactive sites. Quantification of the number of GA and AG cross-links afforded a ratio of about two in favor of AG, irrespective of the nature of the leaving ligands. These results disagree with a previous NMR spectroscopy study, and indicate that GA sequences of DNA are substantially more susceptible to attack by cisplatin than previously thought.

Genetic analysis of the vaccinia virus I6 telomere-binding protein uncovers a key role in genome encapsidation.

The linear, double-stranded DNA genome of vaccinia virus contains covalently closed hairpin termini. These hairpin termini comprise a terminal loop and an A+T-rich duplex stem that has 12 extrahelical bases. DeMasi et al. have shown previously that proteins present in infected cells and in virions form distinct complexes with the telomeric hairpins and that these interactions require the extrahelical bases. The vaccinia virus I6 protein was identified as the protein showing the greatest specificity and affinity for interaction with the viral hairpins (J. DeMasi, S. Du, D. Lennon, and P. Traktman, J. Virol. 75:10090-10105, 2001). To gain insight into the role of I6 in vivo, we generated eight recombinant viruses bearing altered alleles of I6 in which clusters of charged amino acids were changed to alanine residues. One allele (temperature-sensitive I6-12 [tsI6-12]) conferred a tight ts phenotype and was used to examine the stage(s) of the viral life cycle that was affected at the nonpermissive temperature. Gene expression, DNA replication, and genome resolution proceeded normally in this mutant. However, proteolytic processing of structural proteins, which accompanies virus maturation, was incomplete. Electron microscopic studies confirmed a severe block in morphogenesis in which immature, but no mature, virions were observed. Instead, aberrant spherical virions and large crystalloids were seen. When purified, these aberrant virions were found to have normal protein content but to be devoid of viral DNA. We propose that the binding of I6 to viral telomeres directs genome encapsidation into the virus particle.

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.

Sequence requirements for viral RNA replication and VPg uridylylation directed by the internal cis-acting replication element (cre) of human rhinovirus type 14.

Until recently, the cis-acting signals required for replication of picornaviral RNAs were believed to be restricted to the 5' and 3' noncoding regions of the genome. However, an RNA stem-loop in the VP1-coding sequence of human rhinovirus type 14 (HRV-14) is essential for viral minus-strand RNA synthesis (K. L. McKnight and S. M. Lemon, RNA 4:1569-1584, 1998). The nucleotide sequence of the apical loop of this internal cis-acting replication element (cre) was critical for RNA synthesis, while secondary RNA structure, but not primary sequence, was shown to be important within the duplex stem. Similar cres have since been identified in other picornaviral genomes. These RNA segments appear to serve as template for the uridylylation of the genome-linked protein, VPg, providing the VPg-pUpU primer required for viral RNA transcription (A. V. Paul et al., J. Virol. 74:10359-10370, 2000). Here, we show that the minimal functional HRV-14 cre resides within a 33-nucleotide (nt) RNA segment that is predicted to form a simple stem-loop with a 14-nt loop sequence. An extensive mutational analysis involving every possible base substitution at each position within the loop segment defined the sequence that is required within this loop for efficient replication of subgenomic HRV-14 replicon RNAs. These results indicate that three consecutive adenosine residues (nt 2367 to 2369) within the 5' half of this loop are critically important for cre function and suggest that a common RNNNAARNNNNNNR loop motif exists among the cre sequences of enteroviruses and rhinoviruses. We found a direct, positive correlation between the capacity of mutated cres to support RNA replication and their ability to function as template in an in vitro VPg uridylylation reaction, suggesting that these functions are intimately linked. These data thus define more precisely the sequence and structural requirements of the HRV-14 cre and provide additional support for a model in which the role of the cre in RNA replication is to act as template for VPg uridylylation.

Binding of oligonucleotides to a viral hairpin forming RNA triplexes with parallel G*G*C triplets

Infrared and UV spectroscopies have been used to study the assembly of a hairpin nucleotide sequence (nucleotides 3-30) of the 5' non-coding region of the hepatitis C virus RNA (5'-GGCGGGGAUUAUCCCCGCUGUGAGGCGG-3') with a RNA 20mer ligand (5'-CCGCCUCACAAAGGUGGGGU-3') in the presence of magnesium ion and spermidine. The resulting complex involves two helical structural domains: the first one is an intermolecular duplex stem at the bottom of the target hairpin and the second one is a parallel triplex generated by the intramolecular hairpin duplex and the ligand. Infrared spectroscopy shows that N-type sugars are exclusively present in the complex. This is the first case of formation of a RNA parallel triplex with purine motif and shows that this type of targeting RNA strands to viral RNA duplexes can be used as an alternative to antisense oligonucleotides or ribozymes.

Structure of the 1,N(2)-propanodeoxyguanosine adduct in a three-base DNA hairpin loop derived from a palindrome in the Salmonella typhimurium hisD3052 gene.

The solution structure of the 1,N(2)-propanodeoxyguanosine (PdG) adduct was determined in a 3-base hairpin loop formed by d(CGCGGTXTCCGCG) (X = PdG). This sequence is contained within the Salmonella typhimurium hisD3052 gene, a hotspot for frameshift mutagenesis. PdG provides a structural model for the primary adduct induced in DNA by malondialdehyde, the 3-(2'-deoxy-beta-D-erythro-pentofuranosyl)pyrimido[1,2-a]-purin-10(3H)-one (M(1)G) lesion. The solution structure of the PdG-containing hairpin was refined by molecular dynamics calculations restrained by a combination of NMR-derived distances and dihedral angles, using a simulated annealing protocol. The structure of the PdG-modified hairpin consisted of a five-base-pair stem and a three-base loop consisting of T(6), X(7), and T(8). T(6) projected into the minor groove of the stem adjacent to G(4). The modified base X(7) stacked on top of the duplex stem and wedged between bases T(8) and C(9). The PdG moiety was oriented such that the imidazole proton was facing the minor groove of the stem and the exocyclic protons projected into the major groove. The structure of the adducted hairpin was compared with the structure of the corresponding unmodified oligodeoxynucleotide, and was found to be similar. There was a minor difference in the backbone angles of the G and PdG Hairpins at the phosphate linkage between G(5) and T(6) involving the G(5) epsilon angle and T(6) alpha and beta angles. The PdG-modified hairpin exhibited an increase in T(m) of approximately 2 degrees C compared to the unmodified hairpin. The structural and thermodynamic similarities suggested that PdG does not stabilize this hairpin and thus may not promote its extrusion in duplex DNA. The structural results are correlated with the results of site-specific mutagenesis experiments in the same sequence, which do not show evidence of frameshift mutations associated with hairpin loop formation. The geometry of this three-base loop is similar to that of other DNA hairpins containing three-base loops, and suggests a common motif for the folding of these loops.

Remarkable stability of hairpins containing 2',5'-linked RNA loops.

We report here the results of a comparative study of hairpin loops that differ in the connectivity of phosphodiester linkages (3',5'- versus 2',5'-linkages). In addition, we have studied the effect of changing the stem composition on the thermodynamic stability of hairpin loops. Specifically, we constructed hairpins containing one of six stem duplex combinations, i.e., DNA:DNA ("DD"), RNA:RNA ("RR"), DNA:RNA ("DR"), 2',5'-RNA:RNA ("RR"), 2',5'-RNA:DNA ("RD"), and 2',5'-RNA:2',5'-RNA ("RR"), and one of three tetraloop compositions, i.e., 2',5'-RNA ("R"), RNA ("R"), and DNA ("D"). All hairpins contained the conserved and well-studied loop sequence 5'-...C(UUCG)G...-3' [Cheong et al. Nature 1990, 346, 680-682]. We show that the 2',5'-linked loop C(UUCG)G, i.e.,...C(3'p5')U(2'p5')U(2'p5')C(2'p5')G(2'p5')G(3'p5')..., like its "normal" RNA counterpart, forms an unusually stable tetraloop structure. We also show that the stability imparted by 2',5'-RNA loops is dependent on base sequence, a property that is shared with the regioisomeric 3',5'-RNA loops. Remarkably, we find that the stability of the UUCG tetraloop is virtually independent of the hairpin stem composition (DD, RR, RR, etc.), whereas the native RNA tetraloop exerts extra stability only when the stem is duplex RNA (R:R). As a result, the relative stabilities of hairpins with a 2',5'-linked tetraloop, e.g. ggac(UUCG)gtcc (T(m) = 61.4 degrees C), are often superior to those with RNA tetraloops, e.g. ggac(UUCG)gtcc (T(m) = 54.6 degrees C). In fact, it has been possible to observe the formation of a 2',5'-RNA:DNA hybrid duplex by linking the hybrid's strands to a (UUCG) loop. These duplexes (RD), which are not stable enough to form in an intermolecular complex [Wasner et al. Biochemistry 1998, 37, 7478-7486], were stable at room temperature (T(m) approximately 50 degrees C). Thus, 2',5'-loops have potentially important implications in the study of nucleic acid complexes where structural data are not yet available. Furthermore, they may be particularly useful as structural motifs for synthetic ribozymes and nucleic acid "aptamers".

NMR structure of a ribosomal RNA hairpin containing a conserved CUCAA pentaloop.

The structure of a 23 nt RNA sequence, rGGACCCGGGCUCAACCUGGGUCC, was elucidated using homonuclear NMR, distance geometry and restrained molecular dynamics. This RNA is analogous to residues 612-628 of the Escherichia coli 16S rRNA. The structure of the RNA reveals the presence of a pentaloop closed by a duplex stem in typical A-form conformation. The loop does not form a U-turn motif, as previously predicted. A non-planar A.C.A triple base interaction (hydrogen bonds A13 NH6-C10 O2 and C10 N3-A14 NH6) stabilizing the loop structure is inferred from structure calculations. The CUCAA loop structure is asymmetrical, characterized by a reversal of the phosphodiester backbone at the UC step (hydrogen bond C12 NH4-C10 O2') and 3'-stacking within the CAA segment. Loop base U11 is oriented towards the major groove and the consecutive adenosines on the 3'-end of the loop are well stacked, exposing their reactive functional groups in the minor groove defined by the duplex stem. The solution structure of the loop resembles that seen in the 3.3 A X-ray structure of the entire 30S subunit, where the analogous loop interacts with a ribosomal protein and a receptor RNA helix.

High-Purity Preparation of a Large DNA Dumbbell

We report on the efficient biochemical synthesis of a large DNA dumbbell starting from a pair of short DNA hairpins with long single-stranded tails of arbitrary sequence. The DNA dumbbell is obtained by enzymatic ligation yielding a 94-bp duplex stem closed at both termini by single-stranded loops of 5 nt. Following ligation, all unligated precursors and monoligated by-products were multiply biotinylated via nick-translation or primer-extension or both. Thus, they could readily be removed from the DNA dumbbell preparation by a mild biomagnetic separation procedure. The closed conformation of the purified DNA dumbbell was verified by its altered gel mobility as compared with unligated or monoligated samples and by an exonuclease assay. Considering the promising therapeutic potential of DNA dumbbells, the developed biosynthetic approach could be used for high-purity preparation of longer, covalently closed DNA decoys.

Rapid discrimination among individual DNA hairpin molecules at single-nucleotide resolution using an ion channel.

RNA and DNA strands produce ionic current signatures when driven through an alpha-hemolysin channel by an applied voltage. Here we combine this nanopore detector with a support vector machine (SVM) to analyze DNA hairpin molecules on the millisecond time scale. Measurable properties include duplex stem length, base pair mismatches, and loop length. This nanopore instrument can discriminate between individual DNA hairpins that differ by one base pair or by one nucleotide.

Hybridization of single-stranded DNA targets to immobilized complementary DNA probes: comparison of hairpin versus linear capture probes.

A microtiter-based assay system is described in which DNA hairpin probes with dangling ends and single-stranded, linear DNA probes were immobilized and compared based on their ability to capture single-strand target DNA. Hairpin probes consisted of a 16 bp duplex stem, linked by a T(2)-biotin.dT-T(2) loop. The third base was a biotinylated uracil (U(B)) necessary for coupling to avidin coated microtiter wells. The capture region of the hairpin was a 3' dangling end composed of either 16 or 32 bases. Fundamental parameters of the system, such as probe density and avidin adsorption capacity of the plates were characterized. The target DNA consisted of 65 bases whose 3' end was complementary to the dangling end of the hairpin or to the linear probe sequence. The assay system was employed to measure the time dependence and thermodynamic stability of target hybridization with hairpin and linear probes. Target molecules were labeled with either a 5'-FITC, or radiolabeled with [gamma-(33)P]ATP and captured by either linear or hairpin probes affixed to the solid support. Over the range of target concentrations from 10 to 640 pmol hybridization rates increased with increasing target concentration, but varied for the different probes examined. Hairpin probes displayed higher rates of hybridization and larger equilibrium amounts of captured targets than linear probes. At 25 and 45 degrees C, rates of hybridization were better than twice as great for the hairpin compared with the linear capture probes. Hairpin-target complexes were also more thermodynamically stable. Binding free energies were evaluated from the observed equilibrium constants for complex formation. Results showed the order of stability of the probes to be: hairpins with 32 base dangling ends > hairpin probes with l6 base dangling ends > 16 base linear probes > 32 base linear probes. The physical characteristics of hairpins could offer substantial advantages as nucleic acid capture moieties in solid support based hybridization systems.

Actinomycin D Binds Strongly to d(CGACGACG) and d(CGTCGTCG)

Earlier calorimetric studies had indicated that despite the absence of a GpC sequence, the self-complementary octamer d(CGTCGACG) binds strongly to actinomycin D (ACTD) with high cooperativity and a 2:1 drug/duplex ratio. A subsequent optical spectral study with related oligomers led us to suggest that ACTD may likely stack at the G · C basepairs of the duplex termini. New findings are reported herein to indicate that despite the lack of complete self-complementarity, oligomers of d(CGXCGXCG) [X = A or T] motif exhibit unusually strong ACTD affinities with binding constants of roughly 2 × 10 7 M −1 and binding densities of 1 drug molecule per strand. The ACTD binding affinity for the corresponding heteroduplex obtained by annealing these two oligomers is, however, considerably reduced. Although spectroscopic results with related oligomers obtained by removing, replacing, or appending bases at the termini appear to be consistent with the end-stacking model, capillary electrophoretic (CE) evidence provides additional insights into the binding mode. CE experiments with the self-complementary oligomers d(CGAGCTCG) and d(CGTCGACG) revealed contrasting migration patterns in the presence of ACTD, with mobility retardation and acceleration exhibited by the GpC- and non-GpC-containing octamers, respectively, whereas the X/X-mismatched d(CGXCGXCG) experienced retardation. These results, along with those of related oligomers, suggest that ACTD may in fact stack at the duplex stem end of a monomeric hairpin or at the 3′-end of dG as a single strand. The seemingly cooperative ACTD binding and the curved Scatchard plot for the self-complementary d(CGTCGACG) may thus be attributed to the drug-induced duplex denaturation resulting from strong binding to single strands of d(CGXCGYCG) motif. Detailed structural information on the ACTD-DNA complexes, however, must await further NMR investigations.

Melting studies of short DNA hairpins containing the universal base 5-nitroindole

Effects of the universal base 5-nitroindole on the thermodynamic stability of DNA hairpins having a 6 bp stem and four base loops were investigated by optical absorbance and differential scanning calorimetry techniques. Melting studies were conducted in buffer containing 115 mM Na(+). Five different modified versions of DNA hairpins containing a 5-nitroindole base or bases substituted at different positions in the stem and loop regions were examined. Thermo-dynamic parameters of the melting transitions estimated from a two-state analysis of optical melting curves and measured directly by calorimetry revealed that the presence of 5-nitroindole bases in the duplex stem or loop regions of short DNA hairpins significantly affects both their enthalpic and entropic melting components in a compensating manner, while the transition free energy varies linearly with the transition temperature. The calorimetrically determined enthalpy and entropy values of the modified hairpins were considerably smaller (43-53%) than the two-state optical parameters, suggesting that solvent effects may be significant in the melting processes of these hairpins. Results of circular dichroism measurements also revealed slight differences between the modified hairpins and the control in both the duplex and melted states, suggesting subtle structural differences between the control and DNA hairpins containing a 5-nitroindole base or bases.

Spectral and physical characterization of the inverted terminal repeat DNA structure from adenoassociated virus 2.

An oligodeoxynucleotide (ODN) that includes elements found in secondary structures at the 5'- and 3'- ends of adenoassociated virus 2 virion DNA was synthesized by ligation of three overlapping ODNs. The most stable secondary structure was calculated to be branched, with a 61 bp duplex stem, terminating in a three-way junction with 9 bp arms. The electrophoretic mobility of the ODN is slower than expected for normal duplex DNA of the same size, suggesting a bent or branched conformation. CD spectra indicate that the ITR structure is largely B form DNA, although there is a slight blue shift compared to the spectra of the isolated stem and loop elements. Thermal melting experiments indicate that the hairpin is significantly more stable than the isolated stem and loop elements. Singular value decomposition of UV spectra obtained as a function of temperature indicates that four species contribute to changes in the spectra upon denaturation, indicating that the melting is not a simple two-state process. Characterization of the branched ODN by differential scanning calorimetry permits estimation of the enthalpy of melting by a model-free analysis, yielding DeltaHcal= 614 kcal mol-1. This value agrees with the enthalpy computed for the most stable secondary structure.

Efficient synthesis of DNA dumbbells using template-induced chemical ligation in double-stranded polynucleotides closed by minihairpin fragments.

The chemical ligation of 17 50-54-membered nicked DNA dumbbells with different closing fragments, nick positions, and nucleotides facing the nick were investigated. T4, T5, GTA4C, GCGA2GC, and GCGA3GC sequences were chosen as the closing fragments. The nicks were placed in the center of the duplex stem or were adjacent to the closing fragments. N-(3-dimethylaminopropyl)-N'-ethylcarbodiimide and cyanogen bromide were used as the condensing agents. We showed that the ligation efficiency is 10%-90% depending on the sequence of the closing fragments, nick position, and nucleotides facing the nick. Coupling yields of 80%-90% were observed when the nick was situated in the middle of the molecule between two T residues or was adjacent to GCGA2GC or GCGA3GC minihairpins. In the last case, the reacting 3'-phosphate and 5'-hydroxy groups were brought close together by only two base pair minihairpins. The coupling yields did not depend on the nature of the condensing agent. On the basis of the results obtained, we believe a rational design of nicked DNA dumbbells has been developed for efficient chemical synthesis of closed dumbbells.

Studies of DNA dumbbells. VI. Analysis of optical melting curves of dumbbells with a sixteen-base pair duplex stem and end-loops of variable size and sequence.

Optical melting curves of 22 DNA dumbbells with the 16-base pair duplex sequence 5'-G-C-A-T-C-A-T-C-G-A-T-G-A-T-G-C-3' linked on both ends by single-strand loops of A, or C, sequences (iota = 2, 3, 4, 6, 8, 10, 14). T sequences (iota = 2, 3, 4, 6, 8, 10), and G iota sequences (iota = 2, 4) were measured in phosphate buffered solvents containing 30, 70, and 120 mM Na+. For dumbbells with loops comprised of at least three nucleotides, stability is inversely proportional to end-loop size. Dumbbells with loops comprised of only two nucleotide bases generally have lower stabilities than dumbbells with three base nucleotide loops. Experimental melting curves were analyzed in terms of the numerically exact (multistate) statistical thermodynamic model of DNA dumbbell melting previously described (T. M. paner, M. Amaratunga & A. S. Benight (1992), Biopolymers 32, 881). Theoretically calculated melting curves were fitted to experimental curves by simultaneously adjusting model parameters representing statistical weights of intramolecular hairpin loop and single-strand circle states. The systematically determined empirical parameters provided evaluations of the energetics of hairpin loop formation as a function of loop size, sequence, and salt environment. Values of the free energies of hairpin loop formation delta Gloop(n > iota) and single-strand circles, delta Gcir(N) as a function of end-loop size, tau = 2-14, circle size, N = 32 + 2 iota, and loop sequence were obtained. These quantities were found to depend on end-loop size but not loop sequence. Their empirically determined values also varied with solvent ionic strength. Analytical expressions for the partition function Q(T) of the dumbbells were evaluated using the empirically evaluated best-fit loop parameters. From Q(T), the melting transition enthalpy delta H, entropy delta S, and free energy delta G, were evaluated for the dumbbells as a function of end-loop size, sequence, and [Na+]. Since the multistate analysis is based on the numerically exact model, and considers a statistically significant number of theoretically possible partially melted states, it does not require prior assumptions regarding the nature of the melting transition, i.e., whether or not it occurs in a two-state manner. For comparison with the multistate analysis, thermodynamic transition parameters were also evaluated directly from experimental melting curves assuming a two-state transition and using the graphical van't Hoff analysis. Comparisons between results of the multistate and two-state analyses suggested dumbbells with loops comprised of six or fewer residues melted in a two-state manner, while the melting processes for dumbbells with larger end-loops deviate from two-state behavior. Dependence of thermodynamic transition parameters on [Na+] as a function of loop size suggests single-strand end-loops have different counterion binding properties than the melted circle. Results are compared with those obtained in an earlier study of dumbbells with the slightly different stem sequence 5'-G-C-A-T-A-G-A-T-G-A-G-A-A-T-G-C-3' linked on the ends by T iota loops (iota = 2, 3, 4, 6, 8, 10, 14).

Studies of DNA dumbbells. VI. Analysis of optical melting curves of dumbbells with a sixteen-base pair duplex stem and end-loops of variable size and sequence

Optical melting curves of 22 DNA dumbbells with the 16-base pair duplex sequence 5'-G-C-A-T-C-A-T-C-G-A-T-G-A-T-G-C-3' linked on both ends by single-strand loops of A, or C, sequences (iota = 2, 3, 4, 6, 8, 10, 14). T sequences (iota = 2, 3, 4, 6, 8, 10), and G iota sequences (iota = 2, 4) were measured in phosphate buffered solvents containing 30, 70, and 120 mM Na+. For dumbbells with loops comprised of at least three nucleotides, stability is inversely proportional to end-loop size. Dumbbells with loops comprised of only two nucleotide bases generally have lower stabilities than dumbbells with three base nucleotide loops. Experimental melting curves were analyzed in terms of the numerically exact (multistate) statistical thermodynamic model of DNA dumbbell melting previously described (T. M. paner, M. Amaratunga & A. S. Benight (1992), Biopolymers 32, 881). Theoretically calculated melting curves were fitted to experimental curves by simultaneously adjusting model parameters representing statistical weights of intramolecular hairpin loop and single-strand circle states. The systematically determined empirical parameters provided evaluations of the energetics of hairpin loop formation as a function of loop size, sequence, and salt environment. Values of the free energies of hairpin loop formation delta Gloop(n > iota) and single-strand circles, delta Gcir(N) as a function of end-loop size, tau = 2-14, circle size, N = 32 + 2 iota, and loop sequence were obtained. These quantities were found to depend on end-loop size but not loop sequence. Their empirically determined values also varied with solvent ionic strength. Analytical expressions for the partition function Q(T) of the dumbbells were evaluated using the empirically evaluated best-fit loop parameters. From Q(T), the melting transition enthalpy delta H, entropy delta S, and free energy delta G, were evaluated for the dumbbells as a function of end-loop size, sequence, and [Na+]. Since the multistate analysis is based on the numerically exact model, and considers a statistically significant number of theoretically possible partially melted states, it does not require prior assumptions regarding the nature of the melting transition, i.e., whether or not it occurs in a two-state manner. For comparison with the multistate analysis, thermodynamic transition parameters were also evaluated directly from experimental melting curves assuming a two-state transition and using the graphical van't Hoff analysis. Comparisons between results of the multistate and two-state analyses suggested dumbbells with loops comprised of six or fewer residues melted in a two-state manner, while the melting processes for dumbbells with larger end-loops deviate from two-state behavior. Dependence of thermodynamic transition parameters on [Na+] as a function of loop size suggests single-strand end-loops have different counterion binding properties than the melted circle. Results are compared with those obtained in an earlier study of dumbbells with the slightly different stem sequence 5'-G-C-A-T-A-G-A-T-G-A-G-A-A-T-G-C-3' linked on the ends by T iota loops (iota = 2, 3, 4, 6, 8, 10, 14).

Structural similarities and differences between H1- and H2-family DNA minihairpin loops: NMR studies of octameric minihairpins

TheDNA sequences 5′-d(CGC-AC-GCG)-3′ (HPAC), 5′-d(CGC-AA-GCG)-3′ (HPAA), 5′-d(CGC-TC-GCG)-3′ (HPTC), and 5′-d(CGC-CT-GCG)-3′ (HPCT), were studied by means of nmr spectroscopy. At low DNA concentration and no added salt all four molecules adopt a minihairpin structure, containing three Watson–Crick base pairs and a two-residue loop. The structure of the HPAC hairpin is based on quantitative distance restraints, derived by a full relaxation matrix approach (iterative relaxation matrix approach), together with torsion angles obtained from coupling constant analysis. The loop folding is of the H1-family type, characterized by continuous 3′-5′ stacking of the loop bases on the duplex stem. The structure of the HPAA hairpin is similar to that of HPAC, but is more flexible and has a lower thermodynamic stability (Tm 326 K vs 320 K). According to “weakly” distance-constrained simulations in water on the HPAC minihairpin, the typical H1-family loop folding remains intact during the simulation. However, residue-based R factors of simulated nuclear Overhauser effect spectroscopy spectra, free molecular dynamics simulations in vacuo, and unusual chemical shift profiles indicate partial destacking of the loop bases at temperatures below the overall melting midpoint. The dynamic nature of the loop bases gives insight into the geometrical tolerances of stacking between bases in H1-family minihairpin loops. The HPTC and HPCT minihairpins, both containing a pyrimidine base at the first position in the loop, adopt a H2-family type folding, in which the first loop base is loosely bound in the minor groove and the second loop base is stacked upon the helix stem. The thermal stability for these two hairpins corresponds to 327–329 K, but depends on local base sequence. Preference for the type of folding depends on a single substitution from a pyrimidine (H2 family) to a purine (H1 family) at the first position of the miniloop and is explained by differences in base stacking energies, steric size, and the number of possible candidates for hydrogen bonds in the minor groove. In view of newly collected data, previous models of the H1-family and H2-family hairpins had to be revised and are now compatible with the reported HPTC and HPAC structures. The structural difference between the refined structure of HPAC and HPTC show that a conversion between H1-family and H2-family hairpins is geometrically possible by a simple pivot point rotation of 270° along two torsion angles, thereby swiveling the first loop base from a stacked position in a H1-family folding toward a position in the minor groove in a H2-family folding. The second loop residue subsequently shifts to the position of the first base in a concerted fashion. © 1998 John Wiley & Sons, Inc. Biopoly 46: 375–393, 1998

Structural and dynamic helix geometry alterations induced by mismatch base pairs in double-helical RNA

A ribooligonucleotide microhelix derived from the acceptor stem of Escherichia coli tRNA Ala having a C3-A70 mismatch in place of the G3-U70 wobble pair in the wild-type tRNA Ala, and a sequence variant with a regular U3-A70 base pair have been investigated by NMR. In vivo, suppressor tRNA Ala variants with C3-A70 (as well as several other) mismatch pairs are substrates for alanyl-tRNA synthetase (ARS), supporting the hypothesis of an `indirect' recognition of the identity element 3-70 mismatch pair via structural modifications caused by the mispair in comparison to canonical A-RNA helices. It is demonstrated that the C-A mismatch likewise induces helix geometry alterations, in particular with respect to base stacking in the vicinity of the mismatch. However, with reference to the `wild-type' G3-U70 microhelix, destacking in the C3-A70 acceptor stem duplex occurs in the opposite direction from the mismatch pair. Therefore it is concluded that the locally enhanced conformational flexibility or dynamics associated with the structural changes induced by the mismatch pairs could be an essential prerequisite for optimal adaptation of the tRNA Ala acceptor stem to the contact region of the ARS.