Comparison of thermodynamic stabilities between PNA (peptide nucleic acid)/DNA hybrid duplexes and DNA/DNA duplexes.

We have examined quantitatively stabilities of PNA/DNA hybrid duplexes with identical nearest-neighbor base pairs and compared stabilities between PNA/DNA and DNA/DNA. The average difference of stabilization energy of the short PNA/DNA was 0.9 kcal mol(-1), which suggests that the stability of the hybrids with identical nearest-neighbor base pairs can be predicted with the nearest-neighbor model as well as those of nucleic acid duplexes.     

Composites of peptide nucleic acids with titanium dioxide nanoparticles. III. Kinetics of PNA dissociation from nanocomposites containing DNA/PNA duplexes

When delivering peptide nucleic acids (PNA) into cells in the TiO₂ · PL · DNA/PNA nanocomposites consisting of titanium dioxide nanoparticles coated with polylysine (PL) and immobilized DNA/PNA duplexes, it is important to control the rate of the release of PNA from the carrier due to dissociation of the immobilized DNA/PNA duplex, followed by the desorption of PNA to solution while the DNA remains on the carrier. It was found that the rate constant of dissociation of the DNA/PNA duplex in the TiO₂ · PL · DNA/PNA nanocomposites depended on the number of complementary bases in the duplex. The half-retention time values for PNA in the studied nanocomposites containing the duplexes with 10, 12, 14, and 16 overlapping complementary base pairs were 10, 14, 22, and 70 min, respectively. Thus, it was shown that the rate of the release of PNA from the proposed nanocomposites can be controlled by varying the number of overlapping complementary base pairs in the immobilized DNA/PNA duplex. The method of the PNA immobilization may be used for designing nanocomposites having the optimum time value of the PNA release. The proposed TiO₂ · PL · DNA/PNA nanocomposites can be used to efficiently deliver therapeutically significant PNA drugs for their selective effect on pathogenic nucleic acids in cells.     

DNA sequence analysis by hybridization with oligonucleotide microchips: MALDI mass spectrometry identification of 5mers contiguously stacked to microchip oligonucleotides

Matrix-assisted laser desorption ionization mass spectrometry (MALDI MS) has been applied to increase the informational output from DNA sequence analysis. It has been used to analyze DNA by hybridization with microarrays of gel-immobilized oligonucleotides extended with stacked 5mers. In model experiments, a 28 nt long DNA fragment was hybridized with 10 immobilized, overlapping 8mers. Then, in a second round of hybridization DNA–8mer duplexes were hybridized with a mixture of 10 5mers. The stability of the 5mer complex with DNA was increased to raise the melting temperature of the duplex by 10–15°C as a result of stacking interaction with 8mers. Contiguous 13 bp duplexes containing an internal break were formed. MALDI MS identified one or, in some cases, two 5mers contiguously stacked to each DNA–8mer duplex formed on the microchip. Incorporating a mass label into 5mers optimized MALDI MS monitoring. This procedure enabled us to reconstitute the sequence of a model DNA fragment and identify polymorphic nucleotides. The application of MALDI MS identification of contiguously stacked 5mers to increase the length of DNA for sequence analysis is discussed.     

Demonstration of a word design strategy for DNA computing on surfaces.

A strategy for DNA computing on surfaces using linked sets of 'DNA words' that are short oligonucleotides (16mers) is proposed. The 16mer words have the format 5'-FFFFvvvvvvvvFFFF-3' in which 4-8 bits of data are stored in 8 variable ('v') base locations, and the remaining fixed ('F') base locations are used as a word label. Using a template and map strategy, a set of 108 8mers each of which possesses at least a 4 base mismatch with the complements to all the other members of the set (4bm complements) are identified for use as a variable base sequence set. In addition, sets of 4 and 12 word labels of the form ABCD....DCBA that are respectively 8bm and 6bm complements with each other are identified. The 16mers are chosen to have a G/C content of 50% in order to make the thermodynamic stability of the perfectly matched hybridized DNA duplexes similar; a simple pairwise additive method is used to estimate the perfect match and mismatch hybridization thermodynamics. A series of preliminary experiments are presented that use small arrays of 16mers attached to chemically modified gold surfaces and fluorescently labeled complements to study the hybridization adsorption and enzymatic manipulation of the oligonucleotides.     

Unique properties of purine/pyrimidine asymmetric PNA.DNA duplexes: differential stabilization of PNA.DNA duplexes by purines in the PNA strand

PNA.DNA duplexes are significantly stabilized by purine nucleobases in the PNA strand. To elucidate and understand the effect of switching the backbone in a nucleic acid duplex, we now report a thermodynamics study along with a solution conformations study of two purine/pyrimidine strand asymmetric duplexes and a strand symmetrical control by comparing the behavior of all four possible PNA/DNA combinations. In essence, we are comparing an identical basepair stack connected by either an aminoethyl glycine PNA or a deoxyribose DNA backbone. We show that the PNA.DNA duplexes containing purine-rich PNA strands are stabilized with regard to the thermal melting temperature and free energy as well as enthalpy (and concomitantly relatively less entropically disfavored). Based on our data, we find it unlikely that differences in counterion binding (identical ionic-strength dependence was observed), hydration (identical and insignificant water release was observed), or single-strand conformation can be responsible for the difference in duplex stability. The only consistent difference observed between the purine-rich PNA versus the pyrimidine-rich PNA in isosequential PNA.DNA duplexes is the significant increase in both binding enthalpy and entropy for the PNA.DNA duplexes containing pyrimidine-rich PNA in organic solvent, which would indicate that these duplexes are relatively enthalpically disfavored in water. Although our results so far do not allow us to identify the origin of the different stabilities of homopurine/homopyrimidine PNA.DNA duplexes, the evidence does point to a significant structural component, which involves enthalpic contributions both within the duplex structure and also from bound water molecules.     

Molecular dynamics structures of peptide nucleic acid x DNA hybrid in the wild-type and mutated alleles of Ki-ras proto-oncogene--stereochemical rationale for the low affinity of PNA in the presence of an AC mismatch

The low affinity of peptide nucleic acid (PNA) to hybridize with DNA in the presence of a mismatch endows PNA with a high degree of discriminatory capacity that has been exploited in therapeutics for the selective inhibition of the expression of point-mutated genes. To obtain a structural basis for this intriguing property, molecular dynamics simulations are carried out on PNA x DNA duplexes formed at the Ki-ras proto-oncogene, comprising the point-mutated (GAT), and the corresponding wild-type (GGT) codon 12. The designed PNA forms an A...C mismatch with the wild-type sequence and a perfect A...T pair with the point mutated sequence. Results show that large movements in the pyrimidine base of the A...C mismatch cause loss of stacking, especially with its penultimate base, concomitant with a variable mismatch hydrogen bond, including its occasional absence. These, in turn, bring about dynamic water interactions in the vicinity of the mismatch. Enthalpy loss and the disproportionate entropy gain associated with these are implicated as the factors contributing to the increase in free energy and diminished stability of PNA x DNA duplex with the A...C mismatch. Absence of these in the isosequential DNA duplex, notwithstanding the A...C mismatch, is attributed to the differences in topology of PNA x DNA vis-à-vis DNA duplexes. It is speculated that similar effects might be responsible for the reduced stability observed in PNA x DNA duplexes containing other base pair mismatches, and also in mismatch containing PNA x DNA duplexes.     

Non-nearest-neighbor dependence of the stability for RNA bulge loops based on the complete set of group I single-nucleotide bulge loops

Fifty-nine RNA duplexes containing single-nucleotide bulge loops were optically melted in 1 M NaCl, and the thermodynamic parameters DeltaH degrees, DeltaS degrees, DeltaG 37 degrees, and TM for each sequence were determined. Sequences from this study were combined with sequences from previous studies [Longfellow, C. E., et al. (1990) Biochemistry 29, 278-285; Znosko, B. M., et al. (2002) Biochemistry 41, 10406-10417], thus examining all possible group I single-nucleotide bulge loop and nearest-neighbor sequence combinations. The free energy increments at 37 degrees C for the introduction of a group I single-nucleotide bulge loop range between 1.3 and 5.2 kcal/mol. The combined data were used to develop a model for predicting the free energy of a RNA duplex containing a single-nucleotide bulge. For bulge loops with adjacent Watson-Crick base pairs, neither the identity of the bulge nor the nearest-neighbor base pairs had an effect on the influence of the bulge loop on duplex stability. The proposed model for prediction of the stability of a duplex containing a bulged nucleotide was primarily affected by non-nearest-neighbor interactions. The destabilization of the duplex by the bulge was related to the stability of the stems adjacent to the bulge. Specifically, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. The stability of a duplex containing a bulged nucleotide adjacent to a wobble base pair also was primarily affected by non-nearest-neighbor interactions. Again, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. However, when one or both of the bulge nearest neighbors was a wobble base pair, the free energy increment for insertion of a bulge loop is dependent upon the position and orientation of the wobble base pair relative the bulged nucleotide. Bulge sequences of the type ((5'UBX)(3'GY)), ((5'GBG)(3'UU)) and ((5'UBU)(3'GG)) are less destabilizing by 0.6 kcal/mol, and bulge sequences of the type ((5'GBX)(3'UY)) and ((5'XBU)(3'YG)) are more destabilizing by 0.4 kcal/mol than bulge loops adjacent to Watson-Crick base pairs.     

Site-specifically located 8-amino-2'-deoxyguanosine: thermodynamic stability and mutagenic properties in Escherichia coli

2-Nitropropane (2-NP), an important industrial solvent and a component of cigarette smoke, is mutagenic in bacteria and carcinogenic in rats. 8-Amino-2'-deoxyguanosine (8-amino-dG) is one of the types of DNA damage found in liver, the target organ in 2-NP-treated rats. To investigate the thermodynamic properties of 8-amino-dG opposite each of the four DNA bases, we have synthesized an 11mer, d(CCATCG*CTACC), in which G* represents the modified base. By annealing a complementary DNA strand to this modified 11mer, four sets of duplexes were generated each containing one of the four DNA bases opposite the lesion. Circular dichroism studies indicated that 8-amino-dG did not alter the global helical properties of natural right-handed B-DNA. The thermal stability of each duplex was examined by UV melting measurements and compared with its unmodified counterpart. For the unmodified 11mer, the relative stability of the complementary DNA bases opposite G was in the order C > T > G > A, as determined from their -DeltaG degrees values. The free energy change of each modified duplex was lower than its unmodified counterpart, except for the G*:G pair that exhibited a higher melting transition and a larger -DeltaG degrees than the G:G duplex. Nevertheless, the stability of the modified 11mer duplex also followed the order C > T > G > A when placed opposite 8-amino-dG. To explore if 8-amino-dG opposite another 8-amino-dG has any advantage in base pairing, a G*:G* duplex was evaluated, which showed that the stability of this duplex was similar to the G*:G duplex. Mutagenesis of 8-amino-dG in this sequence context was studied in Escherichia coli, which showed that the lesion is weakly mutagenic (mutation frequency approximately 10(-3)) but still can induce a variety of targeted and semi-targeted mutations.     

Comparative NMR study of A(n)-bulge loops in DNA duplexes: intrahelical stacking of A, A-A, and A-A-A bulge loops.

We have prepared a series of deoxyoligonucleotide duplexes of the sequence d(G-C-A-T-C-G-X-G-C-T-A-C-G).d(C-G-T-A-G-C-C-G-A-T-G-C), in which X represents either one (A), two (A-A), or three (A-A-A) unpaired adenine basis. Using two-dimensional proton and phosphorus NMR spectroscopy, we have characterized conformational features of these bulge-loop duplexes in solution. We find that Watson-Crick hydrogen bonding is intact for all 12 base pairs, including the GC bases that flank the bulge loop. Observation of NOE connectivities in both H2O and D2O allows us to unambiguously localize all of the bulged adenine residues to intrahelical positions within the duplex. This is in contrast to an earlier model for multiple-base bulge loops in DNA [Bhattacharyya, A., & Lilley, D. M. J. (1989) Nucleic Acids Res. 17, 6821-6840], in which all but the most 5' bulged base are looped out into solution. We find that insertion of two or three bases into the duplex results in the disruption of specific sequential NOEs for the base step across from the bulge loop site on the opposite strand. This disruption is characterized by a partial shearing apart of these bases, such that certain sequential NOEs for this base step are preserved. We observe a downfield-shifted phosphorus resonance, which we assign in the A-A-A bulge duplex to the 3' side of the last bulged adenine residue. Proton and phosphorus chemical shift trends within the An-bulge duplex series indicate that there is an additive effect on the structural perturbations caused by additional unpaired bases within the bulge loop. This finding parallels previous observations [Bhattacharyya, A., & Lilley, D. M. J. (1989) Nucleic Acids Res. 17, 6821-6840; Hsieh, C.-H., & Griffith, J. D. (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 4833-4837] on the magnitude of the induced bending of DNA duplexes by multiple-base bulge loops.     

Thermodynamic parameters for an expanded nearest-neighbor model for the formation of RNA duplexes with single nucleotide bulges

Thirty-four RNA duplexes containing single nucleotide bulges were optically melted, and the thermodynamic parameters deltaH degrees, deltaS degrees, deltaG degrees (37), and T(M) for each sequence were determined. Data from this study were combined with data from previous thermodynamic data [Longfellow, C. E., Kierzek, R., and Turner, D. H. (1990) Biochemistry 29, 278-85] to develop a model that will more accurately predict the free energy of an RNA duplex containing a single nucleotide bulge. Differences between purine and pyrimidine bulges as well as differences between Group I duplexes, those in which the bulge is not identical to either neighboring nucleotide, and Group II duplexes, those in which the bulge is identical to at least one neighboring nucleotide, were considered. The length of the duplex, non-nearest-neighbor effects, and bulge location were also examined. A model was developed which divides sequences into two groups: those with pyrimidine bulges and those with purine bulges. The proposed model for pyrimidine bulges predicts deltaG degrees (37,bulge) = 3.9 kcal/mol + 0.10deltaG degrees (37,nn) + beta, while the model for purine bulges predicts deltaG degrees (37,bulge) = 3.3 kcal/mol - 0.30deltaG degrees (37,nn) + beta, where beta has a value of 0.0 and -0.8 kcal/mol for Group I and Group II sequences, respectively, and deltaG degrees (37,nn) is the nearest-neighbor free energy of the base pairs surrounding the bulge. The conformation of bulge loops present in rRNA was examined. Three distinct families of structures were identified. The bulge loop was either extrahelical, intercalated, or in a "side-step" conformation.     

Stability of 3' double nucleotide overhangs that model the 3' ends of siRNA

Thermodynamic parameters are reported for duplex formation in 1 M NaCl for 16 RNA sequences, each containing a core tetramer duplex, GGCC, and a 3' overhang consisting of two bases. The results indicate additional double-helical stability is conferred by the double 3' terminal overhang relative to the single 3' terminal overhang. A nearest-neighbor analysis of the data indicates that the free energy contribution at 37 degrees C of the second base in the double 3' terminal overhang varies from 0 to 0.7 kcal/mol. The second base in the 3' double overhang can contribute nearly the same stability to a duplex as a base pair or a 3' dangling overhang. Stability contribution of a dangling base, two nucleotides removed from the 3' end of a duplex, is dependent upon both the identity of the base as well as that of the dangling base that it neighbors. A second dangling base only increases the stability of the duplex when it is neighboring a 3' purine dangling nucleotide. Furthermore, a second dangling pyrimidine provides a greater contribution to duplex stability than a purine. A nearest-neighbor model was developed to predict the influence of 3' double overhang on the stability of duplex formation. The model improves the prediction of free energy and melting temperature when tested against six sequences with different core duplexes.     

Interactions of oligonucleotide analogs containing methylphosphonate internucleotide linkages and 2'-O-methylribonucleosides.

The interactions of oligonucleotide analogs, 12-mers, which contain deoxyribo- or 2'-O-methylribose sugars and methylphosphonate internucleotide linkages with complementary 12-mer DNA and RNA targets and the effect of chirality of the methylphosphonate linkage on oligomer-target interactions was studied. Oligomers containing a single Rp or Sp methylphosphonate linkage (type 1) or oligomers containing a single phosphodiester linkage at the 5'-end followed by 10 contiguous methylphosphonate linkages of random chirality (type 2) were prepared. The deoxyribo- and 2'-O-methylribo- type 1 12-mers formed stable duplexes with both the RNA and DNA as determined by UV melting experiments. The melting temperatures, Tms, of the 2'-O-methylribo-12-mer/RNA duplexes (49-53 degrees C) were higher than those of the deoxyribo-12mer/RNA duplexes (31-36 degrees C). The Tms of the duplexes formed by the Rp isomers of these oligomers were approximately 3-5 degrees C higher than those formed by the corresponding Sp isomers. The deoxyribo type 2 12-mer formed a stable duplex, Tm 34 degrees C, with the DNA target and a much less stable duplex with the RNA target, Tm < 5 degrees C. In contrast, the 2'-O-methylribo type 2 12-mer formed a stable duplex with the RNA target, Tm 20 degrees C, and a duplex of lower stability with the DNA target, Tm < 5 degrees C. These results show that the previously observed greater stability of oligo-2'-O-methylribonucleotide/RNA duplexes versus oligodeoxyribonucleotide/RNA duplexes extends to oligomers containing methylphosphonate linkages and that the configuration of the methylphosphonate linkage strongly influences the stability of the duplexes.     

Thermodynamic and spectroscopic study of bulge loops in oligoribonucleotides

Thermodynamic parameters for bulge loops of one to three nucleotides in oligoribonucleotide duplexes have been measured by optical melting. The results indicate bulges Bn of An and Un have similar stabilities in the duplexes, GCGBmGCG + CGCCGC. The stability increment for a bulge depends on more than its adjacent base pairs. For example, the stability increment for a bulge is affected more than 1 kcal/mol by changing two nonadjacent base pairs or by adding terminal unpaired nucleotides (dangling ends) three base pairs away. Thus a nearest-neighbor approximation for helixes with bulges is oversimplified. Many of the non-self-complementary strands used in this study were observed to form homoduplexes. Such duplexes with GA mismatches were particularly stable.     

Peptide nucleic acid (PNA)/DNA hybrid duplexes: intercalation by an internally linked anthraquinone.

Peptide nucleic acids (PNA) mimic DNA and RNA by forming complementary duplex structures following Watson-Crick base pairing. A set of reporter compounds that bind to DNA by intercalation are known, but these compounds do not intercalate in PNA/DNA hybrid duplexes. Analysis of the hybrid PNA duplexes requires development of reporter compounds that probe their chemical and physical properties. We prepared a series of anthraquinone (AQ) derivatives that are linked to internal positions of a PNA oligomer. These are the first non-nucleobase functional groups that have been incorporated into a PNA. The resulting PNA(AQ) conjugates form stable hybrids with complementary DNA oligomers. We find that when the AQ groups are covalently bound to PNA that they stabilize the hybrid duplex and are, at least partially, intercalated.     

Oligodeoxynucleotides containing conformationally constrained abasic sites: a UV and fluorescence spectroscopic investigation on duplex stability and structure

The synthesis and incorporation into oligodeoxy­nucleotides of two novel, conformationally restricted abasic (AB) site analogs are described. The stability of oligonucleotide 18mer duplexes containing one such AB site opposite any of the four natural DNA bases was investigated by UV melting curve analysis and compared to that of duplexes containing a conformationally flexible propanediol unit 1 or a tetrahydrofuran unit 2 as an AB site analog. No major differences in the melting temperatures (ΔT_m 0–3°C) between the different abasic duplexes were observed. All AB duplexes were found to have T_ms that were lower by 9–15°C relative to a fully matched 18mer control duplex, and by 4–10°C relative to the corresponding 19mer duplexes in which the AB site is replaced by a mismatched nucleobase. Thus we conclude that the loss of stability of a duplex that is encountered by removal of a nucleobase from the stack cannot be compensated with conformational restriction of the AB site. From the van’t Hoff transition enthalpies obtained from the melting curves, it appears that melting cooperativity is higher for the duplexes containing the conformationally rigid AB sites. Fluorescence quenching experiments with duplexes containing the fluorescent base 2-amino­purine (2AP) opposite the AB sites showed a weak tendency towards more efficient stacking of this base in duplexes containing the conformationally constrained AB sites. Thus, such AB sites may structurally stabilize the cavity formed by the removal of a base. Potential applications emerging from the properties of such conformationally constrained AB sites in DNA diagnostics are discussed.     

Base pair opening kinetics study of the aegPNA:DNA hydrid duplex containing a site-specific GNA-like chiral PNA monomer

Peptide nucleic acids (PNA) are one of the most widely used synthetic DNA mimics where the four bases are attached to a N-(2-aminoethyl)glycine (aeg) backbone instead of the negative-charged phosphate backbone in DNA. We have developed a chimeric PNA (chiPNA), in which a chiral GNA-like γ(3)T monomer is incorporated into aegPNA backbone. The base pair opening kinetics of the aegPNA:DNA and chiPNA:DNA hybrid duplexes were studied by NMR hydrogen exchange experiments. This study revealed that the aegPNA:DNA hybrid is much more stable duplex and is less dynamic compared to DNA duplex, meaning that base pairs are opened and reclosed much more slowly. The site-specific incorporation of γ(3)T monomer in the aegPNA:DNA hybrid can destabilize a specific base pair and its neighbors, maintaining the thermal stabilities and dynamic properties of other base pairs. Our hydrogen exchange study firstly revealed the unique kinetic features of base pairs in the aegPNA:DNA and chiPNA:DNA hybrids, which will provide an insight into the development of methodology for specific DNA recognition using PNA fragments.     

Crystal structure of a 14 bp RNA duplex with non-symmetrical tandem GxU wobble base pairs

Adjacent GxU wobble base pairs are frequently found in rRNA. Atomic structures of small RNA motifs help to provide a better understanding of the effects of various tandem mismatches on duplex structure and stability, thereby providing better rules for RNA structure prediction and validation. The crystal structure of an RNA duplex containing the sequence r(GGUAUUGC-GGUACC)2 has been solved at 2.1 A resolution using experimental phases. Novel refinement strategies were needed for building the correct solvent model. At present, this is the only short RNA duplex structure containing 5'-U-U-3'/3'-G-G-5' non-symmetric tandem GxU wobble base pairs. In the 14mer duplex, the six central base pairs are all displaced away from the helix axis, yielding significant changes in local backbone conformation, helix parameters and charge distribution that may provide specific recognition sites for biologically relevant ligand binding. The greatest deviations from A-form helix occur where the guanine of a wobble base pair stacks over a purine from the opposite strand. In this vicinity, the intra-strand phosphate distances increase significantly, and the major groove width increases up to 3 A. Structural comparisons with other short duplexes containing symmetrical tandem GxU or GxT wobble base pairs show that nearest-neighbor sequence dependencies govern helical twist and the occurrence of cross-strand purine stacks.     

Improved thermodynamic parameters and helix initiation factor to predict stability of DNA duplexes.

To improve the previous DNA/DNA nearest-neighbor parameters, thermodynamic parameters (deltaH degrees, deltaS degrees and deltaG degrees) of 50 DNA/DNA duplexes were measured. Enthalpy change of a helix initiation factor is also considered though the parameters reported recently did not contain the factor. A helix initiation factor for DNA/DNA duplex determined here was the same as that of RNA/RNA duplex (deltaG degrees(37) = 3.4 kcal/mol). The improved nearest-neighbor parameters reproduced not only these 50 experimental values used here but also 15 other experimental values obtained in different studies. Comparing deltaG degrees(37) values of DNA/DNA nearest-neighbor parameters obtained here with those of RNA/RNA and RNA/DNA, RNA/RNA duplex was generally the most stable of the three kinds of duplexes with the same nearest-neighbor sequences. Which is more stable between DNA/DNA and RNA/DNA duplexes is sequence dependent.     

Massive parallel analysis of DNA-Hoechst 33258 binding specificity with a generic oligodeoxyribonucleotide microchip.

A generic oligodeoxyribonucleotide microchip was used to determine the sequence specificity of Hoechst 33258 binding to double-stranded DNA. The generic microchip contained 4096 oxctadeoxynucleo-tides in which all possible 4(6)= 4096 hexadeoxy-nucleotide sequences are flanked on both the 3'- and 5'-ends with equimolar mixtures of four bases. The microchip was manufactured by chemical immobilization of presynthesized 8mers within polyacrylamide gel pads. A selected set of immobilized 8mers was converted to double-stranded form by hybridization with a mixture of fluorescently labeled complementary 8mers. Massive parallel measurements of melting curves were carried out for the majority of 2080 6mer duplexes, in both the absence and presence of the Hoechst dye. The sequence-specific affinity for Hoechst 33258 was calculated as the increase in melting temperature caused by ligand binding. The dye exhibited specificity for A:T but not G:C base pairs. The affinity is low for two A:T base pairs, increases significantly for three, and reaches a plateau for four A:T base pairs. The relative ligand affinity for all trinucleotide and tetranucleotide sequences (A/T)(3)and (A/T)(4)was estimated. The free energy of dye binding to several duplexes was calculated from the equilibrium melting curves of the duplexes formed on the oligonucleotide microchips. This method can be used as a general approach for massive screening of the sequence specificity of DNA-binding compounds.