Thermodynamic stability of the 5' dangling-ended DNA hairpins formed from sequences 5'-(XY)2GGATAC(T)4GTATCC-3', where X, Y = A, T, G, C

Expressions for the partition function Q (T) of DNA hairpins are presented. Calculations of Q (T), in conjunction with our previously reported numerically exact algorithm [T. M. Paner, M. Amaratunga, M. J. Doktycz, and A. S. Benight (1990) Biopolymers, 29, 1715-1734], yield a numerical method to evaluate the temperature dependence of the transition enthalpy, entropy, and free energy of a DNA hairpin directly from its optical melting curve. No prior assumptions that the short hairpins melt in a two-state manner are required. This method is then applied in a systematic manner to investigate the stability of the six basepair duplex stem 5'-GGATAC-3' having four-base dangling single-strand ends with the sequences (XY)2, where X, Y = A, T, G, C, on the 5' end and a T4 loop on the 3' end. Results show that all dangling ends of the sample set stabilize the hairpin against melting. Increases in transition temperatures as great as 4.0 degrees C above the blunt-ended control hairpin were observed. The hierarchy of the hairpin transition temperatures is dictated by the identity of the first base of the dangling end adjoining the duplex in the order: purine greater than T greater than C. Calculated melting curves of every hairpin were fit to experimental curves by adjustment of a single parameter in the numerically exact theoretical algorithm. Exact fits were obtained in all cases. Experimental melting curves were also calculated assuming a two-state melting process. Equally accurate fits of all dangling-ended hairpin melting curves were obtained with the two-state model calculation. This was not the case for the melting curve of the blunt-ended hairpin, indicating the presence of a four-base dangling-end drives hairpin melting to a two-state process. Q (T) was calculated as a function of temperature for each hairpin using the theoretical parameters that provided calculated curves in exact agreement with the experimentally obtained optical melting curves. From Q (T), the temperature dependence of the transition enthalpy delta H, entropy delta S, and free energy delta G were calculated for every hairpin providing a quantitative assessment of the effects of dangling ends on hairpin thermodynamics. Comparisons of our results are made with those of the Breslauer group [M. Senior, R. A. Jones, and K. J. Breslauer (1988) Biochemistry 27, 3879-3885] on the T2 5' dangling-ended d(GC)3 duplexes.(ABSTRACT TRUNCATED AT 400 WORDS)     

Solving traveling salesman problems with DNA molecules encoding numerical values

We introduce a DNA encoding method to represent numerical values and a biased molecular algorithm based on the thermodynamic properties of DNA. DNA strands are designed to encode real values by variation of their melting temperatures. The thermodynamic properties of DNA are used for effective local search of optimal solutions using biochemical techniques, such as denaturation temperature gradient polymerase chain reaction and temperature gradient gel electrophoresis. The proposed method was successfully applied to the traveling salesman problem, an instance of optimization problems on weighted graphs. This work extends the capability of DNA computing to solving numerical optimization problems, which is contrasted with other DNA computing methods focusing on logical problem solving.     

The energetic basis of the DNA double helix: a combined microcalorimetric approach

Microcalorimetric studies of DNA duplexes and their component single strands showed that association enthalpies of unfolded complementary strands into completely folded duplexes increase linearly with temperature and do not depend on salt concentration, i.e. duplex formation results in a constant heat capacity decrement, identical for CG and AT pairs. Although duplex thermostability increases with CG content, the enthalpic and entropic contributions of an AT pair to duplex formation exceed that of a CG pair when compared at the same temperature. The reduced contribution of AT pairs to duplex stabilization comes not from their lower enthalpy, as previously supposed, but from their larger entropy contribution. This larger enthalpy and particularly the greater entropy results from water fixed by the AT pair in the minor groove. As the increased entropy of an AT pair exceeds that of melting ice, the water molecule fixed by this pair must affect those of its neighbors. Water in the minor groove is, thus, orchestrated by the arrangement of AT groups, i.e. is context dependent. In contrast, water hydrating exposed nonpolar surfaces of bases is responsible for the heat capacity increment on dissociation and, therefore, for the temperature dependence of all thermodynamic characteristics of the double helix.     

Stereodefined phosphorothioate analogues of DNA: relative thermodynamic stability of the model PS-DNA/DNA and PS-DNA/RNA complexes

Thermodynamic data regarding the influence of P-chirality on stability of duplexes formed between phosphorothioate DNA oligonucleotides (of either stereo-defined all-R(P) or all-S(P) or random configuration at the P atoms) and complementary DNA or RNA strands are presented. Measured melting temperatures and calculated DeltaG(37)(o) values showed that duplexes formed by PS-DNA oligomers with DNA strands are less stable than their unmodified counterparts. However, relative stability of the duplexes ([all-R(P)]-PS-DNA/DNA vs [all-S(P)]-PS-DNA/DNA) depends on their sequential composition rather than on the absolute configuration of PS-oligos, contrary to the results of theoretical considerations and molecular modeling reported in the literature. On the other hand, for all six analyzed pairs of diastereomers, the [all-R(P)]-PS isomers form more stable duplexes with RNA templates, but the origin of stereodifferentiation depends on the sequence with more favorable entropy and enthalpy factors which correlated with dT-rich and dA/dG-rich PS-oligomers, respectively.     

On the application of polyelectrolyte “limiting laws” to the helix-coil transition of DNA. I. Excess univalent cations

A general theory of polyelectrolyte solutions is here used to calculate the differences in Gibbs free energy, enthalpy, and entropy between the coil and helix forms of DNA at any temperature and salt concentration. The salt has univalent cations and is assumed present in excess over the base concentration. The results are restricted to sufficiently dilute solutions. It is shown that the salt concentrations effect is entirely entropic in origin. When applied to the melting temperature, the calculations yield a relation between the enthalpy difference at the melting temperature and the slope of the plot of melting temperature vs. the logarithm of the salt concentration. In accord with observation, both the Gibbs free energy difference at any fixed temperature and the melting temperature are predicted to be linear functions of the log of the salt concentration. However, the theory is not in quantitative agreement with enthalpy data. Data on various colligative and transport properties of both helix and coil forms are reviewed in the text and in Appendix B, and good agreement is found with theory for both forms. No attempt is made to explain why the theory is quantitative for these properties but not for heat measurements. Finally, in Appendix A, an approximate calculation is made of the free energy contributions due to ionic effects not associated with the salt concentration.     

Equilibrium structure and folding of a helix-forming peptide: circular dichroism measurements and replica-exchange molecular dynamics simulations

We have performed experimental measurements and computer simulations of the equilibrium structure and folding of a 21-residue alpha-helical heteropeptide. Far ultraviolet circular dichroism spectroscopy is used to identify the presence of helical structure and to measure the thermal unfolding curve. The observed melting temperature is 296 K, with a folding enthalpy of -11.6 kcal/mol and entropy of -39.6 cal/(mol K). Our simulations involve 45 ns of replica-exchange molecular dynamics of the peptide, using eight replicas at temperatures between 280 and 450 K, and the program CHARMM with a continuum solvent model. In a 30-ns simulation started from a helical structure, conformational equilibrium at all temperatures was reached after 15 ns. This simulation was used to calculate the peptide melting curve, predicting a folding transition with a melting temperature in the 330-350 K range, enthalpy change of -10 kcal/mol, and entropy change of -30 cal/(mol K). The simulation results were also used to analyze the peptide structural fluctuations and the free-energy surface of helix unfolding. In a separate 15-ns replica-exchange molecular dynamics simulation started from the extended structure, the helical conformation was first attained after approximately 2.8 ns, and equilibrium was reached after 10 ns of simulation. These results showed a sequential folding process with a systematic increase in the number of hydrogen bonds until the helical state is reached, and confirmed that the alpha-helical state is the global free-energy minimum for the peptide at low temperatures.     

Spectroscopic and thermodynamic studies on the binding of sanguinarine and berberine to triple and double helical DNA and RNA structures

A comparative study on the interaction of sanguinarine and berberine with DNA and RNA triplexes and their parent duplexes was performed, by using a combination of spectrophotometric, UV thermal melting, circular dichroic and thermodynamic techniques. Formation of the DNA and RNA triplexes was confirmed from UV-melting and circular dichroic measurements. The interaction process was characterized by increase of thermal melting temperature, perturbation in circular dichroic spectrum and the typical hypochromic and bathochromic effects in the absorption spectrum. Scatchard analysis indicated that both the alkaloids bound to the triplex and duplex structures in a non-cooperative manner and the binding was stronger to triplexes than to parent duplexes. Thermal melting studies further indicated that sanguinarine stabilized the Hoogsteen base paired third strand of both DNA and RNA triplexes more tightly compared to their Watson-Crick strands, while berberine stabilized the third strand only without affecting the Watson-Crick strand. However, sanguinarine stabilized the parent duplexes while no stabilization was observed with berberine under identical conditions. Circular dichroic studies were also consistent with the observation that perturbations of DNA and RNA triplexes were more compared to their parent duplexes in presence of the alkaloids. Thermodynamic data revealed that binding of sanguinarine and berberine to triplexes (T.AxT and U.AxU) and duplexes (A.T and A.U) showed negative enthalpy changes and positive entropy changes but that of sanguinarine to C.GxC(+) triplex and G.C duplex exhibited negative enthalpy and negative entropy changes. Taken together, these results suggest that both sanguinarine and berberine can bind and stabilize the DNA and RNA triplexes more strongly than their respective parent duplexes.     

热变性前后水合醛缩酶中水状态的热力学研究

 用差示扫描量热计,测定了不同水合态(水含量h从0.08g/g至2.50g /g)冻结的水合醛缩酶样品中冰的熔化温度和熔化焓,并计算出熔化熵。实验结果表明,变性前,水合醛缩酶中的水转化成四种状状:(1)不冻结水,(2)熔化焓和熔化温度均低了普通水的可冻结水,(3)熔化焓与普通水相同而熔化温度低于普通水的可冻结水,(4)熔化焓和熔化温度均与普通水相同的容积水。变性后,水合醛缩酶中的水可处于五种状态,即不冻结水,熔化焓低于普通水而熔化温度与普通水相同的可冻结水,以及熔化温度均低于普通水且彼此不相同的三种状态的可冻结水。实验还观察到,水合醛缩酶在热变性前后,水在这些态中的量是不同的。     

Prediction of hybridization and melting for double-stranded nucleic acids

This article presents a general statistical mechanical approach to describe self-folding together with the hybridization between a pair of finite length DNA or RNA molecules. The model takes into account the entire ensemble of single- and double-stranded species in solution and their mole fractions at different temperatures. The folding and hybridization models deal with matched pairs, mismatches, symmetric and asymmetric interior loops, bulges, and single-base stacking that might exist at duplex ends or at the ends of helices. All possible conformations of the single- and double-stranded species are explored. Only intermolecular basepairs are considered in duplexes at this stage.In particular we focus on the role of stacking between neighboring nucleotide residues of single unfolded strands as an important source of enthalpy change on helix formation which has not been modeled computationally thus far. Changes in the states of the single strands with temperature are shown to lead to a larger heat effect at higher temperature. An important consequence of this is that predictions of enthalpies, which are based on databases of nearest-neighbor energy parameters determined for molecules or duplexes with lower melting temperatures compared with the melting temperatures of the oligos for which they are used as a predictive tool, will be underestimated.     

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.     

On combinatorial DNA word design.

We consider the problem of designing DNA codes, namely sets of equi-length words over the alphabet [A, C, G, T] that satisfy certain combinatorial constraints. This problem is motivated by the task of reliably storing and retrieving information in synthetic DNA strands for use in DNA computing or as molecular bar codes in chemical libraries. The primary constraints that we consider, defined with respect to a parameter d, are as follows: for every pair of words w, x in a code, there are at least d mismatches between w and x if w not equal x and also between the reverse of w and the Watson-Crick complement of x. Extending classical results from coding theory, we present several upper and lower bounds on the maximum size of such DNA codes and give methods for constructing such codes. An additional constraint that is relevant to the design of DNA codes is that the free energies and enthalpies of the code words, and thus the melting temperatures, be similar. We describe dynamic programming algorithms that can (a) calculate the total number of words of length n whose free energy value, as approximated by a formula of Breslauer et al. (1986) falls in a given range, and (b) output a random such word. These algorithms are intended for use in heuristic algorithms for constructing DNA codes.     

Populous: a tool for building OWL ontologies from templates

Background

Computing accurate nucleic acid melting temperatures has become a crucial step for the efficiency and the optimisation of numerous molecular biology techniques such as in situ hybridization, PCR, antigene targeting, and microarrays. MELTING is a free open source software which computes the enthalpy, entropy and melting temperature of nucleic acids. MELTING 4.2 was able to handle several types of hybridization such as DNA/DNA, RNA/RNA, DNA/RNA and provided corrections to melting temperatures due to the presence of sodium. The program can use either an approximative approach or a more accurate Nearest-Neighbor approach.

Results

Two new versions of the MELTING software have been released. MELTING 4.3 is a direct update of version 4.2, integrating newly available thermodynamic parameters for inosine, a modified adenine base with an universal base capacity, and incorporates a correction for magnesium. MELTING 5 is a complete reimplementation which allows much greater flexibility and extensibility. It incorporates all the thermodynamic parameters and corrections provided in MELTING 4.x and introduces a large set of thermodynamic formulae and parameters, to facilitate the calculation of melting temperatures for perfectly matching sequences, mismatches, bulge loops, CNG repeats, dangling ends, inosines, locked nucleic acids, 2-hydroxyadenines and azobenzenes. It also includes temperature corrections for monovalent ions (sodium, potassium, Tris), magnesium ions and commonly used denaturing agents such as formamide and DMSO.

Conclusions

MELTING is a useful and very flexible tool for predicting melting temperatures using approximative formulae or Nearest-Neighbor approaches, where one can select different sets of Nearest-Neighbor parameters, corrections and formulae. Both versions are freely available at http://sourceforge.net/projects/melting/and at http://www.ebi.ac.uk/compneur-srv/melting/under the terms of the GPL license.     

Base-stacking and base-pairing contributions into thermal stability of the DNA double helix

Two factors are mainly responsible for the stability of the DNA double helix: base pairing between complementary strands and stacking between adjacent bases. By studying DNA molecules with solitary nicks and gaps we measure temperature and salt dependence of the stacking free energy of the DNA double helix. For the first time, DNA stacking parameters are obtained directly (without extrapolation) for temperatures from below room temperature to close to melting temperature. We also obtain DNA stacking parameters for different salt concentrations ranging from 15 to 100 mM Na⁺. From stacking parameters of individual contacts, we calculate base-stacking contribution to the stability of A•T- and G•C-containing DNA polymers. We find that temperature and salt dependences of the stacking term fully determine the temperature and the salt dependence of DNA stability parameters. For all temperatures and salt concentrations employed in present study, base-stacking is the main stabilizing factor in the DNA double helix. A•T pairing is always destabilizing and G•C pairing contributes almost no stabilization. Base-stacking interaction dominates not only in the duplex overall stability but also significantly contributes into the dependence of the duplex stability on its sequence.     

Thermodynamics of DNA branching.

Branched DNA molecules arise transiently as intermediates in genetic recombination or on extrusion of cruciforms from covalent circular DNA duplexes that contain palindromic sequences. The free energy of these structures relative to normal DNA duplexes is of interest both physically and biologically. Oligonucleotide complexes that can form stable branched structures, DNA junctions, have made it possible to model normally unstable branched states of DNA such as Holliday recombinational intermediates. We present here an evaluation of the free energy of creating four-arm branch points in duplex DNA, using a system of two complementary junctions and four DNA duplexes formed from different combinations of the same set of eight 16-mer strands. The thermodynamics of formation of each branched structure from the matching pair of intact duplexes have been estimated in two experiments. In the first, labeled strands are allowed to partition between duplexes and junctions in a competition assay on polyacrylamide gels. In the second, the heats of forming branched or linear molecules from the component strands have been determined by titration microcalorimetry at several temperatures. Taken together these measurements allow us to determine the standard thermodynamic parameters for the process of creating a branch in an otherwise normal DNA duplex. The free energy for reacting two 16-mer duplexes to yield a four-arm junction in which the branch site is incapable of migrating is + 1.1 (+/- 0.4) kcal mol-1 (at 18 degrees C, 10 mM-Mg2+). Analysis of the distribution of duplex and tetramer products by electrophoresis confirms that the free energy difference between the four duplexes and two junctions is small at this temperature. The associated enthalpy change at 18 degrees C is +27.1 (+/- 1.3) kcal mol-1, while the entropy is +89 (+/- 30) cal K-1 mol-1. The free energy for branching is temperature dependent, with a large unfavorable enthalpy change compensated by a favorable entropy term. Since forming one four-stranded complex from two duplexes should be an entropically unfavorable process, branch formation is likely to be accompanied by significant changes in hydration and ion binding. A significant apparent delta Cp is also observed for the formation of one mole of junction, +0.97 (+/-0.05) kcal deg-1 mol-1.     

Melting of Cross-Linked DNA IV. Methods for Computer Modeling of Total Influence on DNA Melting of Monofunctional Adducts,Intrastrand and Interstrand Cross-Links Formed by Molecules of an Antitumor Drug

Abstract

A theoretical method is developed for calculation of melting curves of covalent complexes of DNA with antitumor drugs. The method takes into account all the types of chemical modifications of the double helix caused by platinum compounds and DNA alkylating agents: 1) monofunctional adducts bound to one nucleotide; 2) intrastrand cross-links which appear due to bidentate binding of a drug molecule to two nucleotides that are included into the same DNA strand; 3) interstrand cross-links caused by bidentate binding of a molecule to two nucleotides of different strands. The developed calculation method takes into account the following double helix alterations at sites of chemical modifications: 1) a change in stability of chemically modified base pairs and neighboring ones, that is caused by all the types of chemical modifications; 2) a change in the energy of boundaries between helical and melted regions at sites of chemical modification (local alteration of the factor of cooperativity of DNA melting), that is caused by all the types of chemical modifications, too; 3) a change in the loop entropy factor of melted regions that include interstrand cross-links; 4) the prohibition of divergence of DNA strands in completely melted DNA molecules, which is caused by interstrand cross-links only. General equations are derived, and three calculation methods are proposed to calculate DNA melting curves and the parameters that characterize the helix-coil transition.     

Melting of cross-linked DNA IV. Methods for computer modeling of total influence on DNA melting of monofunctional adducts, intrastrand and interstrand cross-links formed by molecules of an antitumor drug

A theoretical method is developed for calculation of melting curves of covalent complexes of DNA with antitumor drugs. The method takes into account all the types of chemical modifications of the double helix caused by platinum compounds and DNA alkylating agents: 1) monofunctional adducts bound to one nucleotide; 2) intrastrand cross-links which appear due to bidentate binding of a drug molecule to two nucleotides that are included into the same DNA strand; 3) interstrand cross-links caused by bidentate binding of a molecule to two nucleotides of different strands. The developed calculation method takes into account the following double helix alterations at sites of chemical modifications: 1) a change in stability of chemically modified base pairs and neighboring ones, that is caused by all the types of chemical modifications; 2) a change in the energy of boundaries between helical and melted regions at sites of chemical modification (local alteration of the factor of cooperativity of DNA melting), that is caused by all the types of chemical modifications, too; 3) a change in the loop entropy factor of melted regions that include interstrand cross-links; 4) the prohibition of divergence of DNA strands in completely melted DNA molecules, which is caused by interstrand cross-links only. General equations are derived, and three calculation methods are proposed to calculate DNA melting curves and the parameters that characterize the helix-coil transition.     

ssDNA Pairing Accuracy Increases When Abasic Sites Divide Nucleotides into Small Groups

Accurate sequence dependent pairing of single-stranded DNA (ssDNA) molecules plays an important role in gene chips, DNA origami, and polymerase chain reactions. In many assays accurate pairing depends on mismatched sequences melting at lower temperatures than matched sequences; however, for sequences longer than ~10 nucleotides, single mismatches and correct matches have melting temperature differences of less than 3°C. We demonstrate that appropriately grouping of 35 bases in ssDNA using abasic sites increases the difference between the melting temperature of correct bases and the melting temperature of mismatched base pairings. Importantly, in the presence of appropriately spaced abasic sites mismatches near one end of a long dsDNA destabilize the annealing at the other end much more effectively than in systems without the abasic sites, suggesting that the dsDNA melts more uniformly in the presence of appropriately spaced abasic sites. In sum, the presence of appropriately spaced abasic sites allows temperature to more accurately discriminate correct base pairings from incorrect ones.     

Temperature dependence of the gel electrophoretic mobility of superhelical DNA.

We have determined the gel electrophoretic behavior of closed circular plasmid pSM1 DNA (5420 bp) as a function of both temperature and of linking number (Lk). At temperatures below 37 degrees, the electrophoretic mobility first increases, then becomes constant as Lk is decreased below that of the relaxed closed DNA. As the temperature is increased above 37 degrees the electrophoretic mobility first increases as Lk decreases and then varies in a cyclic manner with further decreases in Lk. As the temperature is increased over the range 37 degrees - 65 degrees the cyclic behavior is manifested at progressively smaller decreases in Lk and the amplitude of the cycles increases. We interpret the results in terms of the early melting of superhelical DNA, in which the free energy associated with superhelix formation is progressively transferred to local denaturation. Using a two state approximation, we estimate the free energy change in the first cyclic transition to be 35 Kcal/mole DNA at 37 degrees and to decrease linearly with temperature. The free energy becomes equal to zero at a temperature of 71.6 degrees, which lies within 3 degrees of the melting temperature for the corresponding nicked circular DNA. From the slope of this relationship we estimate the apparent entropy and enthalpy of the first mobility transition to be 6.0 Kcal/mole base pair and 17.3 cal/mole base pair/degree, values consistent with duplex melting.     

DNA sequencing and helix–coil transition. II. Loop entropy and DNA melting

An explicit analytical theory of DNA melting is constructed. It accounts for the loop entropy and the elasticity of DNA strands. Explicit analytical formulas are presented for the melting curves of natural DNA and periodic polymers. The nature of the DNA helix–coil transition is investigated, and it is found to crucially depend on the nucleotide sequence.