Effect of backbone chemistry on hybridization thermodynamics of oligonucleic acids: a coarse-grained molecular dynamics simulation study

A unified view of polymer, dumbbell, and oligonucleotide nearest-neighbor (NN) thermodynamics is presented. DNA NN DeltaG degrees 37 parameters from seven laboratories are presented in the same format so that careful comparisons can be made. The seven studies used data from natural polymers, synthetic polymers, oligonucleotide dumbbells, and oligonucleotide duplexes to derive NN parameters; used different methods of data analysis; used different salt concentrations; and presented the NN thermodynamics in different formats. As a result of these differences, there has been much confusion regarding the NN thermodynamics of DNA polymers and oligomers. Herein I show that six of the studies are actually in remarkable agreement with one another and explanations are provided in cases where discrepancies remain. Further, a single set of parameters, derived from 108 oligonucleotide duplexes, adequately describes polymer and oligomer thermodynamics. Empirical salt dependencies are also derived for oligonucleotides and polymers.

DNA secondary structure plays an important role in biology, genotyping diagnostics, a variety of molecular biology techniques, in vitro-selected DNA catalysts, nanotechnology, and DNA-based computing. Accurate prediction of DNA secondary structure and hybridization using dynamic programming algorithms requires a database of thermodynamic parameters for several motifs including Watson-Crick base pairs, internal mismatches, terminal mismatches, terminal dangling ends, hairpins, bulges, internal loops, and multibranched loops. To make the database useful for predictions under a variety of salt conditions, empirical equations for monovalent and magnesium dependence of thermodynamics have been developed. Bimolecular hybridization is often inhibited by competing unimolecular folding of a target or probe DNA. Powerful numerical methods have been developed to solve multistate-coupled equilibria in bimolecular and higher-order complexes. This review presents the current parameter set available for making accurate DNA structure predictions and also points to future directions for improvement.

We developed OligoCalc as a web-accessible, client-based computational engine for reporting DNA and RNA single-stranded and double-stranded properties, including molecular weight, solution concentration, melting temperature, estimated absorbance coefficients, inter-molecular self-complementarity estimation and intra-molecular hairpin loop formation. OligoCalc has a familiar ‘calculator’ look and feel, making it readily understandable and usable. OligoCalc incorporates three common methods for calculating oligonucleotide-melting temperatures, including a nearest-neighbor thermodynamic model for melting temperature. Since it first came online in 1997, there have been more than 900 000 accesses of OligoCalc from nearly 200 000 distinct hosts, excluding search engines. OligoCalc is available at http://basic.northwestern.edu/biotools/OligoCalc.html, with links to the full source code, usage patterns and statistics at that link as well.