Melting of DNA Oligomers: Dynamical Models and Comparison with Experimental Results

We compare experimental melting curves of short heterogeneous DNA oligomers with theoretical curves derived from statistical mechanics. Partition functions are computed with the one-dimensional Peyrard-Bishop (PB) Hamiltonian, already used in the study of the melting of long DNA chains. Working with short chains we take into account, in the computations, not only the breaking of the interstrand hydrogen bonds, but also the complete dissociation of the double helix into separate single strands. Since this dissociation equilibrium is of general relevance, independent of the particular microscopic model, we give some details of its treatment. We discuss how the non bonded three-dimensional interactions, not explicitly considered in the one-dimensional PB model, are taken into account through the treatment of the dissociation equilibrium. We also evaluate the relevance of the dissociation as a function of the chain length.     

A model for parallel triple helix formation by RecA: single-single association with a homologous duplex via the minor groove

The nucleoproteic filaments of RecA polymerized on single stranded DNA are able to integrate double stranded DNA in a coaxial arrangement (with DNA stretched by a factor 1.5), to recognize homologous sequences in the duplex and to perform strand exchange between the single stranded and double stranded molecules. While experimental results favor the hypothesis of an invasion of the minor groove of the duplex by the single strand, parallel minor groove triple helices have never been isolated or even modeled, the minor groove offering little space for a third strand to interact. Based on an internal coordinate modeling study, we show here that such a structure is perfectly conceivable when the two interacting oligomers are stretched by a factor 1.5, in order to open the minor groove of the duplex. The model helix presents characteristics that coincide with known experimental data on unwinding, base pair inclination and inter-proton distances. Moreover, we show that extension and unwinding stabilize the triple helix. New patterns of triplet interaction via the minor groove are presented.     

A Model for Parallel Triple Helix Formation by Ree A: Single-Strand Association with a Homologous Duple via the Minor Groove

Abstract

The nucleoproteic filaments of RecA polymerized on single stranded DNA are able to integrate double stranded DNA in a coaxial arrangement (with DNA stretched by a factor 1.5), to recognize homologous sequences in the duplex and to perform strand exchange between the single stranded and double stranded molecules. While experimental results favor the hypothesis of an invasion of the minor groove of the duplex by the single strand, parallel minor groove triple helices have never been isolated or even modeled, the minor groove offering little space for a third strand to interact. Based on an internal coordinate modeling study, we show here that such a structure is perfectly conceivable when the two interacting oligomers are stretched by a factor 1.5, in order to open the minor groove of the duplex. The model helix presents characteristics that coincide with known experimental data on unwinding, base pair inclination and inter-proton distances. Moreover, we show that extension and unwinding stabilize the triple helix. New patterns of triplet interaction via the minor groove are presented.     

Thermodynamics of alpha- and beta-structure formation in proteins

An atomic protein model with a minimalistic potential is developed and then tested on an alpha-helix and a beta-hairpin, using exactly the same parameters for both peptides. We find that melting curves for these sequences to a good approximation can be described by a simple two-state model, with parameters that are in reasonable quantitative agreement with experimental data. Despite the apparent two-state character of the melting curves, the energy distributions are found to lack a clear bimodal shape, which is discussed in some detail. We also perform a Monte Carlo-based kinetic study and find, in accord with experimental data, that the alpha-helix forms faster than the beta-hairpin.     

Mismatches and bubbles in DNA

Single mismatches in the DNA double helix form nucleation sites for bubbles. Although the overall melting temperature of the duplex is affected to different degrees depending on the probe length, the statistical weights of the bubble states around the defect are always strongly affected. Here we show experimentally that a single mismatch has indeed a dramatic effect on the distribution of intermediate (bubble) states in the melting transition of DNA oligomers. For probe lengths in the range 20-40 bases, the mismatch transforms a transition with many intermediates into a nearly two-state transition. One surprising consequence is the existence of a regime where the sensitivity of a mismatch detection assay based on monitoring intermediate states would increase with probe length. Our results provide experimental constraints on how mismatches should be implemented in models of DNA melting, such as the widely used thermodynamic nearest neighbor model, to which we compare our data.     

Effect of thiazole orange doubly labeled thymidine on DNA duplex formation

Nucleic acid oligonucleotides are widely used in hybridization experiments for specific detection of complementary nucleic acid sequences. For design and application of oligonucleotides, an understanding of their thermodynamic properties is essential. Recently, exciton-controlled hybridization-sensitive fluorescent oligonucleotides (ECHOs) were developed as uniquely labeled DNA oligomers containing commonly one thymidine having two covalently linked thiazole orange dye moieties. The fluorescent signal of an ECHO is strictly hybridization-controlled, where the dye moieties have to intercalate into double-stranded DNA for signal generation. Here we analyzed the hybridization thermodynamics of ECHO/DNA duplexes, and thermodynamic parameters were obtained from melting curves of 64 ECHO/DNA duplexes measured by ultraviolet absorbance and fluorescence. Both methods demonstrated a substantial increase in duplex stability (ΔΔG°(37) ～ -2.6 ± 0.7 kcal mol(-1)) compared to that of DNA/DNA duplexes of the same sequence. With the exception of T·G mismatches, this increased stability was mostly unaffected by other mismatches in the position opposite the labeled nucleotide. A nearest neighbor model was constructed for predicting thermodynamic parameters for duplex stability. Evaluation of the nearest neighbor parameters by cross validation tests showed higher predictive reliability for the fluorescence-based than the absorbance-based parameters. Using our experimental data, a tool for predicting the thermodynamics of formation of ECHO/DNA duplexes was developed that is freely available at http://genome.gsc.riken.jp/echo/thermodynamics/ . It provides reliable thermodynamic data for using the unique features of ECHOs in fluorescence-based experiments.     

The thermodynamic advantage of DNA oligonucleotide 'stacking hybridization' reactions: energetics of a DNA nick.

'Stacking hybridization reactions' wherein two or more short DNA oligomers hybridize in a contiguous tandem orientation onto a longer complementary DNA single strand have been employed to enhance a variety of analytical oligonucleotide hybridization schemes. If the short oligomers anneal in perfect head-to-tail register the resulting duplex contains a nick at every boundary between hybridized oligomers. Alternatively, if the short oligomers do not hybridize precisely in register, i.e. single strand regions on the longer strand are left unbound, gaps are formed between regions where short oligomers bind. The resulting gapped DNA duplexes are considerably less stable than their nicked duplex analogs. Formation of base pair stacking interactions between neighboring oligomers at the nicks that do not occur in gapped duplexes has been proposed as the source of the observed added stability. However, quantitative evidence supporting this hypothesis for DNA has not been reported. Until now, a direct comparison of the thermodynamics of DNA nicks versus DNA gaps has not been performed. In this communication we report such a comparison. Analysis of optical melting experiments in a well defined molecular context enabled quantitative evaluations of the relative thermodynamic difference between nicked and gapped DNA duplexes. Results of the analysis reveal that a nick may be energetically favored over a gap by at least 1.4 kcal/mol and perhaps as much as 2.4 kcal/mol. The presence of a 5'phosphate at a nick or gap fails to significantly affect their stabilities.     

A Three-State Model with Loop Entropy for the Overstretching Transition of DNA

We introduce a three-state model for a single DNA chain under tension that distinguishes among B-DNA, S-DNA, and M (molten or denatured) segments and at the same time correctly accounts for the entropy of molten loops, characterized by the exponent c in the asymptotic expression S ∼ -c ln n for the entropy of a loop of length n. Force extension curves are derived exactly by employing a generalized Poland-Scheraga approach and then compared to experimental data. Simultaneous fitting to force-extension data at room temperature and to the denaturation phase transition at zero force is possible and allows us to establish a global phase diagram in the force-temperature plane. Under a stretching force, the effects of the stacking energy (entering as a domain-wall energy between paired and unpaired bases) and the loop entropy are separated. Therefore, we can estimate the loop exponent c independently from the precise value of the stacking energy. The fitted value for c is small, suggesting that nicks dominate the experimental force extension traces of natural DNA.     

Thermodynamics of polymolecular duplexes between phosphate-methylated DNA and natural DNA

Phosphate-methylated (P.M.) DNA possesses a very high affinity for complementary natural DNA, as a result of the absence of interstrand electrostatic repulsions. In this study, a model system phosphate-methylated d[Cn] with natural d(Gk) (n less than k) is chosen for an investigation of the thermodynamic properties that determine duplex stability. The enthalpy change of a melting transition is shown to be considerably larger than is observed for corresponding natural DNA duplexes. It is found that delta Hn0 of GG/CC nearest neighbor pairwise interaction equals -15.6 kcal/mol, compared to -11.0 kcal/mol for the natural analog. The entropy change is strongly dependent on the length of the natural DNA strand and the number of phosphate-methylated DNA oligomers hybridized. The results are explained by means of a model in which a cooperative effect for subsequent hybridizations of phosphate-methylated DNA oligomers is assumed, thus giving additional stability.     

Formation of joint DNA molecules by two eukaryotic strand exchange proteins does not require melting of a DNA duplex

We have examined whether DNA strand exchange activities from nuclear extracts of HeLa cells or Drosophila melanogaster embryos have detectable helicase or melting activities. The partially purified recombinases have been shown to recognize homologous single strand and double strand DNA molecules and form joint molecules in a DNA strand exchange reaction. The joint molecule product consists of a linear duplex joined at one end by a region of DNA heteroduplex to a homologous single strand circular DNA. Using two different partially duplex helicase substrates, we are unable to detect any melting of duplex regions under conditions that promote joint molecule formation. One substrate consists of a 32P-labeled oligonucleotide 20 or 30 bases long annealed to M13mp18 circular single strand DNA. The second substrate consists of a linear single strand region flanked at each end by short duplex regions. We observe that even in the presence of excess recombinase protein or after prolonged incubation no helicase activity is apparent. Control experiments rule out the possibility that a helicase is masked by reannealing of displaced single strand fragments. Based on these findings and other data, we conclude that the human and D. melanogaster recombinases recognize and pair homologous sequences without significant melting of duplex DNA prior to strand exchange.     

Mesoscopic model parametrization of hydrogen bonds and stacking interactions of RNA from melting temperatures

Information about molecular interactions in DNA can be obtained from experimental melting temperature data by using mesoscopic statistical physics models. Here, we extend the technique to RNA and show that the new parameters correctly reproduce known properties such as the stronger hydrogen bonds of AU base pairs. We also were able to calculate a complete set of elastic constants for all 10 irreducible combinations of nearest neighbours (NNs). We believe that this is particularly useful as experimentally derived information about RNA elasticity is relatively scarce. The melting temperature prediction using the present model improves over those from traditional NN model, providing thus an alternative way to calculate these temperatures for RNA. Additionally, we calculated the site-dependent base pair oscillation to explain why RNA shows larger oscillation amplitudes despite having stronger AU hydrogen bonds.     

Melting of a self-complementary DNA minicircle. Comparison of optical melting theory with exchange broadening of the nuclear magnetic resonance spectrum

Melting curves are calculated for the 16-base-pair duplex DNA sequence 5' GTATCCGTACGGATAC 3' linked on the ends by TTTT single-strand loops. The equilibrium statistical thermodynamic theory of DNA melting is modified to include effects of end-loops on the melting transition. An excellent fit of the experimental melting curve in 0.2 M-NaCl is obtained using two adjustable parameters, one for end-loop formation and the other for formation of the complete 40-base single-strand loop. The best-fit calculated melting curve permits evaluation of these parameters. The free energy to close a TTTT end-loop is 2.12 kcal/mol (1 cal = 4.184 J). A TTTT end-loop or hairpin loop is significantly more stable than an internal loop of comparable size sandwiched between two helical regions, even after allowing for the different stacking contributions. Reasons for this increased stability are presented. The loop free energy of the 40-base single-strand open minicircle is evaluated to be +1.27 kcal/mol, thus favoring the melting of two end-loops into the large open minicircle. The present results are compared with those of others for d(T-A) oligomers. The sequence TTTT forms a more stable end-loop, or hairpin, than TATA by about 2.0 kcal/mol. Theoretical rate constants for the proton-transfer step in the standard hydrogen-exchange model are calculated by extending the theory of diffusion-controlled reactions to take account of the electrostatic potential of the DNA. The predicted ratios of rate constants for different pairs of catalysts exchanging an A.T proton agree satisfactorily with the available experimental data for a 14-base-pair linear duplex, which confirms the diffusion-control of the proton-transfer step. Data presented here for the 16 base-pair duplex of the minicircle are consistent with catalysis-limited exchange in which the proton-transfer step is likewise diffusion-controlled. Under catalysis-limited conditions, the imino proton exchange rates are predicted from the catalytic rate constants, prevailing buffer catalyst concentrations, and the equilibrium constants to form the unstacked open state of optical melting theory. The observed exchange rates of the A.T base-pairs show no sign of the strong predicted end-melting trend, and exceed the predicted values by factors of 10 to 400. Moreover, the succession of "melting" in the nuclear magnetic resonance line-broadening deviates from that predicted by optical melting theory.(ABSTRACT TRUNCATED AT 400 WORDS)     

Influence of base-pair changes and cooperativity parameters on the melting curves of short DNAs

Theoretical melting curves were calculated for four DNA restriction fragments, 157–257 base pairs (bp), and a series of hypothetical block DNAs with sequences d(C_2xA_xC_2x). d(C_2xT_xG_2x), 5 ? x ? 40. These DNAs provided a mixture of A·T/G·C sequence distributions with which to investigate the effects of parameters and base-pair changes on the melting of short DNAs. The sensitivity of DNA melting curves to changes in internal loop melting parameters σ and κ was examined. As Expected, theoretical melting curves of short DNAs with a quasirandom base-pair sequence vary little with changes in internal loop parameters. End melting dominates the transition behaviour of these moleucles. This was also observed for the block DNAs up to x = 22. Beyond this length, melting curves are highly sensitive to the internal loop parameters. Sensitivity is also predicted for a 157-bp fragment with a block distribution of A·T and G·C pairs. These results indicate that accurate evaluation of internal loop parameters is possible with short DNAs (100–200 bp) containing a G·C/A·T/G·C block distribution with at least 22 bp in each block. Duplex-to-single-strands dissociation parameters were reevaluated form experimental melting curve data of eight DNA fragments using a least squares fit approach. This analysis confirmed parameter values previously found with a simplified dissociation model. A Priori predictions are made on the effects of base-pair changes on the melting curves of three characterized DNA restriction fragments. Single base-pair changes are predicted to induce small but measurable changes in the melting curves. The characteristics of the altered melting curves depend on the location of the base-pair change.     

Force-induced melting of the DNA double helix. 2. Effect of solution conditions

In this paper, we consider the implications of the general theory developed in the accompanying paper, to interpret experiments on DNA overstretching that involve variables such as solution temperature, pH, and ionic strength. We find the DNA helix-coil phase boundary in the force-temperature space. At temperatures significantly below the regular (zero force) DNA melting temperature, the overstretching force, f(ov)(T), is predicted to decrease nearly linearly with temperature. We calculate the slope of this dependence as a function of entropy and heat-capacity changes upon DNA melting. Fitting of the experimental f(ov)(T) dependence allows determination of both of these quantities in very good agreement with their calorimetric values. At temperatures slightly above the regular DNA melting temperature, we predict stabilization of dsDNA by moderate forces, and destabilization by higher forces. Thus the DNA stretching curves, f(b), should exhibit two rather than one overstretching transitions: from single stranded (ss) to double stranded (ds) and then back at the higher force. We also predict that any change in DNA solution conditions that affects its melting temperature should have a similar effect on DNA overstretching force. This result is used to calculate the dependence of DNA overstretching force on solution pH, f(ov)(pH), from the known dependence of DNA melting temperature on pH. The calculated f(ov)(pH) is in excellent agreement with its experimental determination (M. C. Williams, J. R. Wenner, I. Rouzina, and V. A. Bloomfield, Biophys. J., accepted for publication). Finally, we quantitatively explain the measured dependence of DNA overstretching force on solution ionic strength for crosslinked and noncrosslinked DNA. The much stronger salt dependence of f(ov) in noncrosslinked DNA results from its lower linear charge density in the melted state, compared to crosslinked or double-stranded overstretched S-DNA.     

Rad54 oligomers translocate and cross-bridge double-stranded DNA to stimulate synapsis

Rad54 is a key component of the eukaryotic recombination machinery. Its presence in DNA strand-exchange reactions in vitro results in a significant stimulation of the overall reaction rate. Using untagged Rad54, we show that this stimulation can be attributed to enhancement of the formation of a key reaction intermediate known as DNA networks. Using a novel, single DNA molecule, dual-optical tweezers approach we show how Rad54 stimulates DNA network formation. We discovered that Rad54 oligomers possess a unique ability to cross-bridge or bind double-stranded DNA molecules positioned in close proximity. Further, Rad54 oligomers rapidly translocate double-stranded DNA while simultaneously inducing topological loops in the DNA at the locus of the oligomer. The combination of the cross-bridging and double-stranded DNA translocation activities of Rad54 stimulates the formation of DNA networks, leading to rapid and efficient DNA strand exchange by Rad51.     

Spectroscopic studies on the interaction of aristololactam-beta-D-glucoside with DNA and RNA double and triple helices: A comparative study.

The interaction of aristololactam-beta-D-glucoside (ADG), a DNA intercalating alkaloid, with the DNA triplexes, poly(dT). poly(dA)xpoly(dT) and poly(dC).poly(dG)xpoly(dC+), and the RNA triplex poly(rU).poly(rA)xpoly(rU) was investigated by circular dichroic, UV melting profile, spectrophotometric, and spectrofluorimetric techniques. Comparative interaction with the corresponding Watson-Crick duplexes has also been examined under identical experimental conditions. Triplex formation has been confirmed from biphasic thermal melting profiles and analysis of temperature-dependent circular dichroic measurements. The binding of ADG to triplexes and duplexes is characterized by the typical hypochromic and bathochromic effects in the absorption spectrum, quenching of steady-state fluorescence intensity, a decrease in fluorescence quantum yield, an increase or decrease of thermal melting temperatures, and perturbation in the circular dichroic spectrum. Scatchard analysis indicates that ADG binds both to the triplexes and the duplexes in a noncooperative manner. Binding parameters obtained from spectrophotometric measurements are best fit by the neighbor exclusion model. The binding affinity of ADG to the DNA triplexes is substantially stronger than to the RNA triplex. Thermal melting study further indicates that ADG stabilizes the Hoogsteen base-paired third strand of the DNA triplexes whereas it destabilizes the same strand of RNA triplex but stabilizes its Watson-Crick strands. Comparative data reveal that ADG exhibits a stronger binding to the triple helical structures than to the respective double helical structures.