Sequence-dependent oligonucleotide-target duplex stabilities: rules from empirical studies with a set of twenty-mers

We were interested in developing a better method to predict the thermal stability of specific oligonucleotide-target duplexes. Recognizing that the base sequence can have important effects, we investigated the use of a simple parameter based on nearest-neighbor stacking interactions, the mean stacking temperature. We took values for doublet stabilities from the literature and used a computer program to calculate mean stacking temperatures for all oligonucleotides of specified length and G + C content in the M13 phage genome. As expected, the program predicted a fairly broad range of stabilities for different sequences of equal G + C content. We selected 20-mer sequences representing the highest and lowest mean stacking temperatures at 25%, 50% and 75% G + C and synthesized them for use as probes against M13 DNA immobilized on filters. By hybridizing and washing at different temperatures, we demonstrated that mean stacking temperatures correlate well with observed stabilities. Relative stabilities of the six oligos were predicted correctly in every case. We used conditions appropriate to oligonucleotide probing and polymerase chain reaction and we were able to derive simple linear equations relating the empirical data and mean stacking temperature for both. These observations should be useful in planning experiments with oligonucleotides.