PubMed HealthSearch

PubMed · 3459152

Predicting DNA duplex stability from the base sequence.

Abstract

We report the complete thermodynamic library of all 10 Watson-Crick DNA nearest-neighbor interactions. We obtained the relevant thermodynamic data from calorimetric studies on 19 DNA oligomers and 9 DNA polymers. We show how these thermodynamic data can be used to calculate the stability and predict the temperature-dependent behavior of any DNA duplex structure from knowledge of its base sequence. We illustrate our method of calculation by using the nearest-neighbor data to predict transition enthalpies and free energies for a series of DNA oligomers. These predicted values are in excellent agreement with the corresponding values determined experimentally. This agreement demonstrates that a DNA duplex structure thermodynamically can be considered to be the sum of its nearest-neighbor interactions. Armed with this knowledge and the nearest-neighbor thermodynamic data reported here, scientists now will be able to predict the stability (delta G degree) and the melting behavior (delta H degree) of any DNA duplex structure from inspection of its primary sequence. This capability should prove valuable in numerous applications, such as predicting the stability of a probe-gene complex; selecting optimal conditions for a hybridization experiment; deciding on the minimum length of a probe; predicting the influence of a specific transversion or transition on the stability of an affected DNA region; and predicting the relative stabilities of local domains within a DNA duplex.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

K J Breslauer, R Frank, H Blöcker, L A Marky. 1986. Predicting DNA duplex stability from the base sequence.. https://doi.org/10.1073/pnas.83.11.3746

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

Crystal and molecular structure of d(GTGCGCAC): investigation of the effects of base sequence on the conformation of octamer duplexes.

The structure of the self-complementary deoxyoctanucleotide d(GTGCGCAC), which crystallized as an A-type helix in the space group P4(3)2(1)2, with one strand in the crystallographic asymmetric unit has been determined and refined to a final R-value of 0.154 using 1.64-A diffraction data collected on an area detector. In contrast to the closely related sequence d(GTGTACAC)tet, there was no evidence for an ordered spermine molecule in the major groove of this octamer. Ordered water is found associated with almost all the exposed hydrogen bonding groups of the octamer. A pentagonal ring of water molecules is hydrogen bonded to O6 and N7 of G3 and the N4 and O6 of the C4.G13 base pair. A detailed comparison of the local helical parameters of d(GTGCGCAC) and d(GTGTACAC)tet is presented. The base sequence change at the center of the octamers affects several of the local helical parameters, via both intra- and interduplex interactions within the crystal.

Base Composition

Base pairing geometry in GA mismatches depends entirely on the neighboring sequence.

We have synthesized nine self-complementary DNA oligomers containing different flanking sequences adjacent to a pair of contiguous GA mismatches, and have used high resolution nuclear magnetic resonance (n.m.r.) to investigate the GpA phosphodiester backbone conformation and mismatch pairing schemes in these duplexes. We found dramatic effects of the flanking base pair on the hydrogen bonding and backbone conformation, which appear to be coupled. Thus the Ganti-Aanti base pairing scheme in a NAGATN sequence switches to a more stable sheared GA base pairing scheme in a NCGAGN or NTGAAN context, while no duplex is formed (or only GA bulges occur) when NAGATN is changed to NGGACN. Furthermore, the more stable sheared GA pairing in NPyGAPuN sequences is associated with a BII rather than BI backbone conformation for the phosphodiester between the adjacent mismatched GA pairs. The overall stability of these adjacent GA mismatches as measured by imino proton n.m.r. studies is Py-GA-Pu > A-GA-T > G-GA-C.

Base Composition

Structure-related properties of the mutagenic lesion 6-O-methylguanine in DNA.

Chemical probing of the structures of a few very similar 30 base-pair duplexes containing a 6-O-methylguanine (meG) residue at the 16th position reveals that the modified base simultaneously perturbs the helical structure in two ways; it preferentially unstacks the 3' neighbouring base residue (thymine in this study) on the same strand and it unstacks the pyrimidine to which it is base-paired. Depending on its neighbouring 5' base residue and the base-pairing pyrimidine, this perturbation can extend to a few base-pairs in both 3' and 5' directions from the abnormal base-pair. These perturbations can be detected by cleavage at the site for the restriction enzyme MaeII. The unstaking of the C in the meG.C and A.C base-pairs may explain the de novo methylation of these helices by the human DNA-(cytosine-5-)methyltransferase. Interestingly, the kinetics of repair of the 6-O-methylguanine-containing dinucleotides by the cloned human methylguanine methyltransferase appears to be largely determined by the strength of the stacking interaction between the 6-O-methylguanine and the 5' neighbouring base.

Base Composition