Thermodynamics of RNA duplexes with tandem mismatches containing a uracil-uracil pair flanked by C.G/G.C or G.C/A.U closing base pairs

The thermodynamics governing the denaturation of RNA duplexes containing 8 bp and a central tandem mismatch or 10 bp were evaluated using UV absorbance melting curves. Each of the eight tandem mismatches that were examined had one U-U pair adjacent to another noncanonical base pair. They were examined in two different RNA duplex environments, one with the tandem mismatch closed by G.C base pairs and the other with G.C and A.U closing base pairs. The free energy increments (Delta Gdegrees(loop)) of the 2 x 2 loops were positive, and showed relatively small differences between the two closing base pair environments. Assuming temperature-independent enthalpy changes for the transitions, (Delta Gdegrees(loop)) for the 2 x 2 loops varied from 0.9 to 1.9 kcal/mol in 1 M Na(+) at 37 degrees C. Most values were within 0.8 kcal/mol of previously estimated values; however, a few sequences differed by 1.2-2.0 kcal/mol. Single strands employed to form the RNA duplexes exhibited small noncooperative absorbance increases with temperature or transitions indicative of partial self-complementary duplexes. One strand formed a partial self-complementary duplex that was more stable than the tandem mismatch duplexes it formed. Transitions of the RNA duplexes were analyzed using equations that included the coupled equilibrium of self-complementary duplex and non-self-complementary duplex denaturation. The average heat capacity change (DeltaC(p)) associated with the transitions of two RNA duplexes was estimated by plotting DeltaH degrees and DeltaS degrees evaluated at different strand concentrations as a function of T(m) and ln T(m), respectively. The average DeltaC(p) was 70 +/- 5 cal K(-)(1) (mol of base pairs)(-)(1). Consideration of this heat capacity change reduced the free energy of formation at 37 degrees C of the 10 bp control RNA duplexes by 0.3-0.6 kcal/mol, which may increase Delta Gdegrees(loop) values by similar amounts.     

Thermodynamics of the melting of PNA(2)/DNA triple helices.

Equilibrium melting curves were obtained for triplexes, formed by single stranded DNA containing an A10 target with bis-PNA consisting of two T10 decamers. Thermodynamic parameters of melting were determined for Na(+) concentrations 50, 200 and 600mM by two methods. The melting enthalpy Delta H degrees was evaluated from the width of the differential melting curves and from the concentration dependence of the melting temperature. The latter method allowed also evaluating the melting entropy Delta S degrees. The following results were obtained: Delta H degrees = - 137 kcal/M, Delta S degrees = - 368 cal/M.K, Delta G degrees (298)= - 27 kcal/M. No dependence of Delta H degrees, Delta S degrees and Delta G degrees (298) was observed upon ionic strength within the accuracy of the experiment (+/- 10%). The absolute values of Delta H degrees, Delta S degrees and Delta G degrees(298) are 2 to 3 times higher than the published values of corresponding melting parameters for decameric PNA/DNA duplexes of various nucleic base sequences. The origin of the extremely high stability of the triplexes is discussed.     

High‐resolution calorimetric and optical melting profiles of DNA plasmids: Resolving contributions from intrinsic melting domains and specifically designed inserts

We demonstrate that differential scanning calorimetry (DSC) can be used to yield high‐resolution melting profiles for DNA plasmids that agree in all major features with the corresponding plasmid melting profiles derived using more traditional optical techniques. We further demonstrate that by combining information derived from both calorimetric and optical melting profiles one can glean insights that are unavailable from either melting curve alone. By using both optical and calorimetric observables, we show how one can resolve, identify, and measure the thermodynamic properties of particular sequences/domains of interest within a plasmid. We also show that complementary DSC and optical melting studies on plasmids with and without specifically designed inserts can provide fundamental advantages over the corresponding melting studies on other model system constructs for thermodynamically characterizing nucleic acid sequences/structures. © 1999 John Wiley & Sons, Inc. Biopoly 50: 303–318, 1999     

Solution conformation of an RNA hairpin loop.

The hairpin conformation adopted by the RNA sequence 5'GCGAUUUCUGACCGCC3' has been studied by one- and two-dimensional NMR spectroscopy. Exchangeable imino spectra in 60 mM Na+ indicate that the hairpin has a stem of six base pairs (indicated by boldface type) and a loop of three nucleotides. NOESY spectra of nonexchangeable protons confirm the formation of the stem region. The duplex has an A-conformation and contains an A.C apposition; a G.U base pair closes the loop region. The stem nucleotides have C3'-endo sugar conformations, as expected of an A-form duplex, whereas the three loop nucleotides adopt C2'-endo sugar puckers. Stacking within the loop, C8 upon the sugar of U7, stabilizes the structure. The pH dependence of both the exchangeable and nonexchangeable NMR spectra is consistent with the formation of an A+.C base pair, protonated at the N1 position of adenine. The stability of the hairpin was probed by using absorbance melting curves. The hairpin structure with the A+.C base pair is about +2 kcal/mol less stable in free energy at 37 degrees C than the hairpin formed with an A.U pair replacing the A+.C pair.     

Thermodynamic stability of DNA tandem mismatches

The thermodynamics of nine hairpin DNAs were evaluated using UV-monitored melting curves and differential scanning calorimetry (DSC). Each DNA has the same five-base loop and a stem with 8-10 base pairs. Five of the DNAs have a tandem mismatch in the stem, while four have all base pairs. The tandem mismatches examined (ga/ga, aa/gc, ca/gc, ta/ac, and tc/tc) spanned the range of stability observed for this motif in a previous study of 28 tandem mismatches. UV-monitored melting curves were obtained in 1.0 M Na(+), 0.1 M Na(+), and 0.1 M Na(+) with 5 mM Mg(2+). DSC studies were conducted in 0.1 M Na(+). Transition T(m) values were unchanged over a 50-fold range of strand concentration. Model-independent enthalpy changes (DeltaH degrees ) evaluated by DSC were in good agreement (+/-8%) with enthalpy values determined by van't Hoff analyses of the melting curves in 0.1 M Na(+). The average heat capacity change (DeltaC(p)) associated with the hairpin to single strands transitions was estimated from plots of DeltaH degrees and DeltaS degrees with T(m) and ln T(m), respectively, and from profiles of DSC curves. The average DeltaC(p) values (113 +/- 9 and 42 +/- 27 cal x K(-1) x mol(-1) of bp), were in the range of values reported in previous studies. Consideration of DeltaC(p) produced large changes in DeltaH degrees and DeltaS degrees extrapolated from the transition region to 37 degrees C and smaller but significant changes to free energies. The loop free energy of the five tandem mismatches at 37 degrees C varied over a range of approximately 4 kcal x mol(-1) for each solvent.     

Sequence and temperature effect on hydrogen bond disruption in DNA determined by a statistical analysis

We use the modified self-consistent phonon approximation theory to calculate temperature dependent interbase hydrogen bond disruption profiles for a number of six base pair repeating sequence infinite B-DNA polymers with various guanine-cytosine/adenine-thymine ratios. For comparison we also include results we have obtained in our earlier work on several B-DNA homopolymers, copolymers and a four-base-pair repeating sequence polymer. Our theory gives a statistical estimate of thermal fluctuational disruption probability of individual hydrogen bonds in individual base pairs in DNA as a function of temperature. The calculated probabilities show no sequence dependence at premelting temperatures, in agreement with proton exchange measurements. These probabilities however become very sensitive to base sequence at temperatures close to the observed melting temperatures. Multi-phasic critical transitions are found in which a portion of base pairs are disrupted at temperatures below the final disruption temperature. These transitions include localized as well as non-localized base pair opening. The localized transitions involve disruption of a few base-pairs at every other location without large scale base unstacking, and they may not appear in the observed UV curves with current resolution. On the other hand the overall disruption behavior is consistent with observations. The midpoint transition temperatures are close to the observed melting temperatures and these temperatures show the observed linear dependence on guanine-cytosine content. Our calculations indicate that our theory can be used effectively to calculate H-bond disruption behavior of different DNA sequences. Received: 20 February 1996 / Accepted: 2 May 1996     

Spectral analysis for base composition of DNA undergoing melting

A microcomputer-controlled spectrophotometer is described for obtaining the base composition of melting domains in DNA from derivative melting curves. Values have been determined for the differential molar extinction coefficients for the A-T and G-C base pair at the three wavelengths most useful for spectral analysis of base composition, 260, 270 and 282 nm. The average RMS error for these values was 29 l(mol X cm)-1 for the melting of 14 DNA specimens ranging in base composition from 0-0.72 F(G + C). A precision of approximately 1% in base composition of domains is possible. Such analysis is useful for confirming or establishing assignments of domains to particular subtransitional features in high resolution melting curves.     

Stacking energies in DNA

Variations in base mono- and dipoles result in variations in stacking energies for the 10 unique neighbor pairs in DNA. Stacking energies for pair M on N, expressed as TMN, were derived by matrix decomposition of a large set of linear algebraic expressions relating the measured Tm for subtransitions emanating from large polymeric DNAs, and the fractional neighbor frequencies, fMN, for the domains responsible for the transitions, Tm = sigma fMNTMN. Tm were determined for subtransitions that dissociate in approximately all-or-none fashion in high resolution melting profiles of partially deleted and recombinant forms of pBR322 DNA. Three different analytical maneuvers were undertaken to resolve subtransitions: site-specific cleavage of domains; deletion of domains; and addition of domains. Three dozen domains of widely divergent, quasi-random neighbor frequencies were identified and assigned, resulting in a unique set of values for TMN with standard deviation, sigma = +/- 0.23 degree C. The average difference between calculated and experimental Tm for domains is only +/- 0.17 degree C, indicating that the thermodynamic properties of these domains are not in any way unusual. Assuming delta S to be constant for all pairs, the corresponding delta HMN are found to have a precision of +/- 10 calories.mol-1 and an accuracy of +/- 606 calories.mol-1. TMN used to calculate melting curves by statistical mechanical analysis of sequences of the different plasmid specimens in this study were in quantitative agreement with observed curves for most sequences. These TMN differ significantly from those determined previously and also correlate poorly with values determined by quantum chemical analysis. Stabilities of neighbor pairs, expressed as the difference in free energy between that for a given pair (MN) and that for the average of like pairs (M, N), depend on the relationship of stacked purines and pyrimidines as follows. delta delta Gpu-py(-466 cal) greater than delta delta Gpu-pu(+52 cal) greater than delta delta Gpy-pu(+335 cal) Differences between experimental Tm and Tm calculated with TMN for the isolated neighbor pairs in the B-conformation are useful in the identification of altered structures and unusual modes of dissociation of helixes. A significantly higher Tm is observed for the highly biased repeated sequence synthetic helixes dA.dT, d(AGC).d(GCT), and d(GAT).d(ATC), reflecting auxiliary sources of stability such as bifurcated hydrogen bonds and/or altered structures for these helixes.     

Role of DNA regions flanking the tryptophan promoter of Escherichia coli. I. Insertion of synthetic oligonucleotides

To examine the effect of altering the nucleotide sequence near the promoter on its activity, pKO-1 vector derivatives have been constructed which allow insertion of DNA fragments at specified sites upstream or downstream from the trp promoter. Oligonucleotides that might be expected to alter the melting properties, or have a tendency to form a distinctive nonstandard structure were introduced. These oligonucleotides had the repeating dinucleotide sequences GC, AT or AG. Sequence analysis of the inserts and studies of the relative galactokinase expression from the altered plasmids indicated that changes upstream from the trp promoter at -39 or beyond had little effect on trp promoter activity, whereas changes at +2 or farther downstream produced up to two-fold increases in gene expression, as compared to the control plasmid.     

Melting profile and temperature dependent binding constant of an anticancer drug daunomycin-DNA complex

We calculate thermal fluctuational base pair opening probability and the drug binding constant of a daunomycin-bound Poly d(CGTA) · Poly d(TACG) at temperatures from room temperature to its melting temperature. For comparison we also carry out a calculation on a drug-free DNA with the same sequence. Our calculations are carried out by means of a statistical approach using microscopic structures and established force fields and with cooperative effects incorporated into the algorithm. Both hydrogen bond disruption probabilities and drug unstacking probability are determined self-consistently. These probabilities are then used to determine temperature dependent base pair opening probabilities and the drug binding constant. The calculated base pair opening probabilities and drug binding constant are found to be in fair agreement with experiments carried out at room temperature. Our calculation shows cooperative base pair disruption and drug dissociation at certain critical temperatures close to the observed melting temperatures for similar helices. We find that the temperature dependence of the drug binding constant fits well to the van't Hoff relation, in agreement with observations. Our calculation indicates the occurrence of a premelting transition in the drug-bound DNA helix. Some comments are made about this premelting transition.     

Loop energy in DNA

The loop function ?(N), representing the statistical weight of N complementary residues in a closed ring, has been determined by analysis of high-resolution melting curves of a series of recombinant homopoly(A · T)_N inserts in pBR322 DNA, where 150 > N > 50 base pairs (bp). Loops are found more stable and therefore presumably less elastic than expected for an ideal, freely jointed chain. A value of 97 ± 2 bp is obtained for the empirical orientation-stiffness parameter D in the expression for nonideal chains, ?(N) = (N + D)^?1.7. The 10% increase in apparent stiffness over that of an ideal chain closely coincides with the extent of residual stacking in the loop. It is thus concluded that the more favorable loop energy, such as expected of smaller loops, is due to the incipient helical orientation of some residues, predisposing the loop to reclosure. A quantitative loop function is essential for the prediction and assignment of domains in DNA.     

Intramolecular DNA triplexes in supercoiled plasmids. I. Effect of loop size on formation and stability

The capacity of four oligopurine.oligopyrimidine (pur.pyr) sequences with different lengths of interruptions in the center [GAA)4(N)n(GAA)4G) (n = 3, 5, 7, and 9) to adopt intramolecular DNA triplexes was evaluated in recombinant plasmids. The hyperreactive patterns of the pur.pyr inserts to specific chemical probes (OsO4, diethyl pyrocarbonate, and dimethyl sulfate) at the base pair level demonstrate that intramolecular triplexes with identical 12-base triads in the stem but with different loop sizes (4, 6, 8, and 10 bases) can form in supercoiled plasmids. Furthermore, the extent of OsO4 modification was measured as a function of temperature and of average negative supercoil density. In addition, the transition free energy of B-DNA to triplexes at pH 4.5 was determined by two-dimensional electrophoresis. These comparative studies show that longer loops require more supercoil energy for triplex formation and are less thermostable than triplexes with shorter loops. Also, it may be that not only the loop size but the base composition of the loop region affects the structural transition and triplex stability. Thus, these results significantly broaden the range of natural pur.pyr sequences that may adopt triplexes.     

Electrostatic forces at helix-coil boundaries in DNA

The Tm of internal loop-forming (dA.dT)N domains in pBR322 DNA has been measured over a tenfold range of [Na+]. The slopes SN = dTm/d log [Na+] are linear and decrease in magnitude with decreasing loop size N, signaling a reduction in Na+ released during the transition of these domains to the coil state. Values of SN decrease linearly with increasing N-1 in accordance with the expectation of a simple model for the occurrence of a gradient of long-range electrostatic forces at helix-coil boundaries, and extrapolate almost precisely to the value of S infinity observed for (dA.dT) infinity. These results indicate (1) less counterion is released per phosphate residue from the finite loop than from the infinite-sized loop, and (2) the difference in binding is constant for each boundary formed and independent of the size of the loop within the range examined: approximately 350 base pair (bp) greater than N greater than 71 bp. The slope of the dependence of SN on N-1 indicates the region of higher charge density at the boundary extends at least 18 A into the coil and probably 40-50 A before dropping to a value characteristic of the unperturbed coil. The free energy for excess counterion binding at boundaries can be expressed by -delta G/RT = 10.47 log[Na+] + 5.234 When the loop entropy function in a statistical mechanical algorithm for the dissociation of DNA is weighted by this quantity, calculated Tm are seen to vary by only +/- 0.09 degrees C from observed.     

DNA curvature influences the internal motions of supercoiled DNA.

We present evidence that short curved DNA segments can act as mediators for the ordering of large domains in superhelical DNA. Using a non-invasive solution method (dynamic light scattering), we investigated the effect of permanently curved inserts on the solution structure and on the internal motions of superhelical plasmid DNA. We find that the dynamics of superhelical DNA are strongly influenced by sequence- or protein-induced bending: in superhelical plasmids containing curved inserts the amplitude of the internal motion is lower than that of non-curved controls. Furthermore, the relative arrangement of curved sequences in the plasmids can influence the overall shape of the superhelical DNA. On linearized forms of the plasmids, these effects are not observed.     

New and versatile cloning vectors with kanamycin-resistance marker

Described here is a pair of small multi-copy kanamycin-resistance plasmids, containing the pUC lacZ alpha-complementation peptide and the pUC18 and pUC19 multiple cloning site. These plasmids and their derivatives allow simple and rapid transfer of inserts from one replicon to another without the necessity of purifying the insert from vector.     

Cloning in Escherichia coli K-12 of a Na+-dependent transport system from a marine bacterium.

The transport of D-alanine by Escherichia coli K-12 neither requires nor is stimulated by Na+. The transport of D-alanine by the marine bacterium Alteromonas haloplanktis 214 requires Na+ specifically. Mutants of E. coli which were unable to transport D-alanine were isolated by enrichment for D-cycloserine resistance. One of the mutants was transformed with a gene bank of A. haloplanktis chromosomal DNA. Two transformants, E. coli RM1(pPM1) and E. coli RM1(pPM2) were able to transport D-alanine by a Na+-dependent mechanism. Li+ and K+ were unable to replace Na+. Both transformants contained chimeric plasmids with inserts which hybridized with A. haloplanktis but not E. coli chromosomal DNA or each other. Despite the lack of homology between the inserts, Na+-dependent D-alanine transport in the two transformants could not be distinguished either by kinetic studies or by differences in the capacity of various amino acids to compete for D-alanine uptake.     

Betaine structure and the presence of hydroxyl groups alters the effects on DNA melting temperatures

Betaine lowers the melting temperature of deoxyribonucleic acid (DNA) and decreases its dependence on base composition. The effects of synthetic betaine analogs on the melting of DNA samples with different GC content were measured. Since many polyhydroxy compounds also lower DNA melting temperatures, hydroxyl-substituted betaine analogs were included. Some synthetic sulfonate analogs of betaine lowered the DNA melting temperatures by twice as much at the same molar concentration. They were up to twice as effective at decreasing the base pair dependence. Some carboxylate homologs of betaine, substituted with hydroxyl groups, increased the melting temperature. This effect was greater with low GC content DNA. Sulfonate analogs of betaine with hydroxyl groups usually destabilize the DNA, while their carboxylate analogs stabilize the DNA. Distances between the charges of these synthetic zwitterionic solutes influence the effect on DNA, with the optimum separation being two or three methylene groups. A betaine with two hydroxyl groups on one N-alkyl group had a greater effect than an isomer with two hydroxyl groups on separate N-alkyl substituents. We suggest that the effect of these solutes depends on structuring the hydration water of DNA, as well as interactions with the DNA structure itself.     

The nonequilibrium character of DNA melting: effects of the heating rate on the fine structure of melting curves.

To investigate the effects of heating rate on the DNA melting profile and to test the predictions of the theory of slow relaxation processes in DNA melting (1) concerning these effects, we obtained differential melting curves for the Bsp I C1 fragment of T7 DNA (1461 bp) and the Sma I-Eco RI fragment of Col E1 DNA (1291 bp) at heating rates of 0.05 and 0.5 deg/min. At low ionic strength (0.02 M Na+) the heating rate has been shown to affect the position of the third peak in melting curve for C1 fragment. According to the melting maps (2), this peak corresponds to the unwinding of the section between the end of the molecule and the region already melted. At high ionic strength (0.2 M Na+), when the melting of this DNA is reversible (3), the position of the peaks does not depend on the heating rate. In the case of the Col E1 DNA fragment the heating rate affects, as might be expected from the melting maps (4), only the last peak, as the melting of the last section is always nonequilibrium. The results of the study are in good qualitative agreement and in satisfactory quantitative agreement with the theoretical predictions (1).     

Studies of DNA dumbbells. III. Theoretical analysis of optical melting curves of dumbbells with a 16 base-pair duplex stem and Tn end loops (n = 2, 3, 4, 6, 8, 10, 14).

Optical melting curves of seven DNA dumbbells with the 16 base-pair duplex sequence 5'G-C-A-T-A-G-A-T-G-A-G-A-A-T-G-C3' linked on both ends by Tn (n = 2, 3, 4, 6, 8, 10, and 14) loops measured in 30, 70, and 120 mM Na+ are analyzed in terms of the numerically exact statistical thermodynamic model of DNA melting. The construction and characterization of these molecules were described in the previous paper (Amaratunga et al., 1992). As was recently reported for hairpins (T. M. Paner, M. Amaratunga, M. J. Doktycz, and A. S. Benight, 1990, Biopolymers, Vol. 29, pp. 1715-1734) theoretically calculated melting curves were fitted to experimental curves by simultaneously adjusting the parameters representing loop and circle formation to optimize the fits. The systematically determined empirical parameters provide evaluations of the free energies of hairpin loop formation delta Gloop (n) and single-strand circles delta Gcircle (N), as a function of end loop size, n = 2-14, and circle size, N = 32 + 2n. The dependence of these quantities on solvent ionic strength over the range from 30 to 120 mM Na+ was evaluated. An approximately analytical expression for the partition function Q(T) of the dumbbells was formulated that allowed a means for determining the transition enthalpy delta H degrees and entropy delta S degrees for every dumbbell, revealing the dependence of these quantities on loop size. In this multistate approach a manifold of partially melted intermediate microstates are considered and therefore no assumptions regarding the nature of the melting transitions (that they are two-state) are required. The transition thermodynamic parameters were also determined from a van't Hoff analysis of the melting curves. Comparisons between the results of the multistate analysis and the two-state van't Hoff analysis revealed significant differences for the dumbbells with larger end loops, indicating that the melting transitions of the larger looped dumbbells deviate considerably from two-state behavior. Results are then compared with published melting studies of much larger DNA dumbbells (D. B. Naritsin and Y. L. Lyubchenko, 1990, Journal of Biomolecular Structure and Dynamics, Vol. 8, pp. 1-13), of small hairpins (Paner et al., 1990; M. J. Doktycz, T. M. Paner, M. Amaratunga and A. S. Benight, 1990, Biopolymers, Vol. 30, pp. 829-845) and another dumbbell (A. S. Benight, J. M. Schurr, P. F. Flynn, B. R. Reid, and D. E. Wemmer, 1988) Journal of Molecular Biology, Vol. 200, pp. 377-399).(ABSTRACT TRUNCATED AT 400 WORDS)