Position dependent mismatch discrimination on DNA microarrays – experiments and model
© Naiser et al; licensee BioMed Central Ltd. 2008
Received: 12 August 2008
Accepted: 01 December 2008
Published: 01 December 2008
The propensity of oligonucleotide strands to form stable duplexes with complementary sequences is fundamental to a variety of biological and biotechnological processes as various as microRNA signalling, microarray hybridization and PCR. Yet our understanding of oligonucleotide hybridization, in particular in presence of surfaces, is rather limited. Here we use oligonucleotide microarrays made in-house by optically controlled DNA synthesis to produce probe sets comprising all possible single base mismatches and base bulges for each of 20 sequence motifs under study.
We observe that mismatch discrimination is mostly determined by the defect position (relative to the duplex ends) as well as by the sequence context. We investigate the thermodynamics of the oligonucleotide duplexes on the basis of double-ended molecular zipper. Theoretical predictions of defect positional influence as well as long range sequence influence agree well with the experimental results.
Molecular zipping at thermodynamic equilibrium explains the binding affinity of mismatched DNA duplexes on microarrays well. The position dependent nearest neighbor model (PDNN) can be inferred from it. Quantitative understanding of microarray experiments from first principles is in reach.
The well-known double-helix structure of nucleic acids results from sequence-specific binding between complementary single strands. Sequential base pairing between A·T and C·G base pairs along the two complementary strands results in the formation of stable duplexes. This so called hybridization process is fundamental to many biological processes and biotechnologies. Microarrays consist of surface-tethered probe sequences, which act as specific scavengers for their respective complementary target sequence. The molecular recognition enables a highly parallel detection of nucleic acid sequences in complex target mixtures. Hybridization also occurs with single mismatched (MM) base pairs, however, these duplexes are significantly less stable than the corresponding perfect match (PM) [1, 2]. The single base pair mismatch-discrimination capability of short (~20 nt) oligonucleotide probes provides an important diagnostic tool for the detection of point-mutations and single nucleotide polymorphisms (SNPs) . DNA duplex stability arises from hydrogen bonding and base stacking interactions (the latter comprise van der Waals interactions, electrostatic and hydrophobic interactions between adjacent base pairs). According to the well-established nearest-neighbor model, thermodynamically a nucleic acid duplex can be considered the sum of these nearest-neighbor (NN) interactions [4–6]. The binding free energy of an oligonucleotide duplex can be predicted from the nearest-neighbor free energy parameters: The helix propagation parameters (one for each of the 10 possible base-pair doublets in case of a DNA/DNA duplex) account for the duplex sequence. Further parameters provide corrections for duplex initiation, A·T terminal pairs or a symmetry penalty in case of self-complementary sequences. The NN model adequately predicts oligonucleotide duplex melting temperatures T M in bulk solution . Datasets of Watson-Crick NN parameters  provide the basis for nucleic acid structure and melting temperature prediction software like the DINAMelt web server  (UNAFold), the HYTHER server and others. The NN model can be extended beyond the Watson-Crick pairs to include single base MM defects [7, 10].
In spite of good knowledge about nucleic acid hybridization in solution, the prediction of binding affinities on DNA microarrays remains empirical. Recent microarray studies [11–15] report, that the influence of even a point defect on hybridization signal intensity cannot be predicted easily. In particular the influence of defect position on the hybridization signal is stronger than the influence of MM-type [12, 14, 16].
Experiments show that the two-state nearest-neighbor (TSNN) approach , which has been very successful in predicting duplex stability in solution, does not appropriately describe MM binding affinities on DNA microarrays. The NN model does not account for the position of the individual NN pairs , except for the outermost ones. Based on microarray data, Zhang et al.  proposed a position dependent nearest-neighbor (PDNN) model. The model assumes that the duplex binding free energy can be expressed as a weighted sum of stacking energies with empirically derived positional weight parameters [17–21]. The purpose of this study is to investigate the influence of point defects on (surface bound) hybridization experimentally and theoretically. Previous studies investigate mismatch discrimination with samples of very different sequence motifs [11, 12]. However, other effects such as secondary structure formation or competitive binding may reduce the visibility of the impact of the MM-defect on the binding affinity. To avoid such complications we performed experiments with fixed sequence motifs: We focus on small variations of the probe sequences. We perform hybridization studies with home-made microarrays comprising sets of very similar probe sequences. We use a single target sequence in each hybridization assay in order to avoid inter-target binding as well as target competition of different sequences for one and the same probe sequence. In order to avoid excluded volume interactions or secondary structure we limit the length of the target sequence to be of the order of the probes. These simplifications (described in detail in ) enable a detailed investigation of the influences of defect type, defect position, flanking base pairs and the sequence motif on the binding affinity. The extensive set of hybridization affinities obtained from our experiments enables us to perform a very complete analysis. We compare the experimental data to theoretical modeling based on a double-ended molecular zipper approach (the double-ended nucleic acid zipper has been previously described by [22–26]). We find that in order to reproduce the microarray hybridization signal in our model, the heterogeneity of binding affinities – mostly owing to in situ synthesis-related probe defects (e.g. probe polydispersity) – needs to be taken into account. More than that, synthesis defects arise as useful for parallel detection of many different sequences.
DNA Microarray Hybridization Experiments
Protocols for the preparation of dendrimer-functionalized microarray substrates (adapted from ) and for the light-directed synthesis (based on NPPOC-phosphoramidites ), as well as details on the hybridization assay and on fluorescence microscopy based microarray analysis (Fig. 1) are provided in Naiser et al. [14, 15].
In each microarray hybridization assay a probe set of cognate probes with purposefully introduced point mutations – derived from a common probe sequence motif – is hybridized against a single target sequence, which perfectly matches the probe sequence motif. We systematically vary defect type and defect position to provide the complete "defect profile" of hybridization affinities with probe sets. We include not only all single base mismatches (MMs), but also, in order to investigate mismatch discrimination in a broader context of other sequence defects, we consider single base bulges (originating from insertions and deletions) as well as probes with multiple defects. Since the ca. 130 probes within each probe set differ only by single bases we are able to distinguish between defect-positional and sequence influence. In our experimental conditions hybridization equilibrium is reached after a few tens of minutes. Further details can be found in .
Position Dependent Influence of Single Base MMs and Bulges on Probe-Target Binding Affinity
As can directly be inferred from Fig. 2, defects in the middle of the probes are most destabilizing. In the center of a 16 mer duplex a single nucleotide MM typically reduces the hybridization signal to 0–40% of the corresponding PM duplex hybridization signal. Defect type and nearest-neighbor effects have less influence on the hybridization signal than defect position. Our experiments show a mostly monotonous decrease of hybridization signals over a range of typically 5–8 defect positions (for 16 mer probes and up to 14 positions for some 25 mer sequence motifs) from the duplex ends towards the center of the duplex. This is consistent with previous work [11, 12].
We also perform hybridization experiments on oligonucleotide duplexes with two single base deletion defects at varying positions x and y. The results show that the binding affinity depends also on the relative position of the defects (for details see Additional file 1, Fig. S5 and ). The hybridization signal is largest if each defect is located close to an end. Lowest binding affinities are observed for defect configurations which divide the sequence into three roughly equally long subsequences. Closely spaced defects (with a distance of less than four nucleotides) systematically increase their impact with distance.
In thermodynamic equilibrium duplex nucleation (determined by the slow nucleation rate k nuc ) is balanced by duplex dissociation with the dissociation rate k diss . The widely used two-state nearest-neighbor model (including mismatched NN-dimers as described by ) cannot provide an explanation for this positional influence, it does not account for the position of the individual nearest-neighbor dimers. We assume that the nucleation rates k nuc of very similar duplexes (differing by a single base pair, e.g. a PM duplex and a corresponding mismatched duplex) are virtually identical. Thus, the positional dependence observed experimentally can be expected to result from differences in k diss . In agreement with  we show that the positional influence originates from end-domain unzipping. Our experimental findings suggest a common mechanism for DPI, that is independent of the defect type. Further, the relatively long range of the DPI (Fig. 3A and 3B) suggests that molecular dynamics may well be a good candidate for an explanation. The symmetry of DPI (with respect to the duplex ends) and sequence-specific deviations from the symmetry indicate a zipping related mechanism. Thus, in order to account for partial denatured duplex states, we use a double-ended zipper model of the oligonucleotide duplex to determine mismatched oligonucleotide duplex stabilities as a function of defect position. We consider a situation in thermodynamic equilibrium.
Double-ended Zipper Model
k+ and k- are the fast zipping and unzipping rates determined by the nearest-neighbor propagation parameters of the individual base-pair doublets. The time-evolution of the oligonucleotide zipper can be considered a biased random walk with a finite probability for complete dissociation (described by the duplex dissociation rate k diss ). Since we consider thermodynamic equilibrium, we can use a partition function for fast numerics.
Partition Function Approach (PFA) to Investigate Oligonucleotide Duplex Thermodynamics
We use a partition function approach [22–25] and investigate if the double-ended zipper model can reproduce our experimental results. On the basis of unified NN-Parameters  we calculate statistical weights of partially denatured duplex states. The effect of partial binding with respect to microarray data was discussed earlier in [24–26].
Consideration of Point Defects
While defects near the duplex ends result in low mismatch discrimination only (i.e. small reduction of K with respect to the PM binding affinity) defects in the center result in higher MM discrimination as K then approaches the value of the two-state equilibrium constant. NN-pair free energy increments for single base MMs are in the range of 1 to 3 kcal/mol per NN-pair (derived from NN parameters [8, 10]). Employing these values in Eq. 9 for Δg° = -1.4 kcal/mol (Fig. 5A), DPI propagation is restricted to 3 or 6 NN-pairs, respectively. However, in subsequences with weakly bound NN-pairs (as demonstrated in Fig. 5B) the DPI can propagate further towards the middle of the duplex.
Relation Between the Hybridization Signal and the Binding Free Energy ΔG D
In order to compare our numerical analysis to the experimentally observed hybridization signals we need to understand how the hybridization signal (fluorescence intensity from hybridized targets) is linked to duplex stability. As detailed below the assumption of a single (homogeneous) binding affinity within a microarray feature of the Langmuir adsorption model does not describe the experimentally observed hybridization signal intensities well. In this section we account for the heterogeneity that is introduced by in situ synthesis related random mutations of the microarray probe sequences.
Assuming a stepwise error rate of 10%, more than 90% of the 25 mer duplexes contain at least one synthesis error . Since the number of synthesis errors per probe follows a binomial distribution, the majority of the strands contains between one and three single base defects.
We calculate binding constants K i of the individual, randomly "mutated" probe sequences on the basis of the zipper model. Using the approach of Forman et al.  we obtain the total hybridization signal by summing up over the distribution of probes, where the contribution of each individual mutated probe θ i is described by a Langmuir equation (Eq. 10) with the binding constant K i . Probe polydispersity (in length as well as in sequence) reproduces a "stretched isotherm"  (similar to a Sips isotherm), with a significantly broadened transition region. This explains our experimental results in Figs. 6 and 2 well. A simulation of the transfer function θ() for various error rates and a comparison between the experimental data in Fig. 6 and the corresponding simulation results are provided in Additional file 1, Fig. S6.
Numerical Analysis of Mismatched Duplex Stability-Comparison with Experimental Results
To model experimental results with the partition function approach we choose the NN free energy of the mismatched base pair as a free parameter. = 1 kcal/mol (at T = 325 K) describes our experimental observations (in particular the dominating positional influence with respect to defect type-related influences) best – see Fig. 3. This value is also in good agreement with bulk solution parameters . Results of the numerical simulation (in Fig. 3A) demonstrate that the shallower slope of the hybridization signal at the right duplex-end corresponds to a series of weak NN pairs (as anticipated by Eq. 9). The partition function Z(x) largely determines the positional influence. Additionally, as shown in Fig. 3B, defect-type related influence (the difference between MM and PM free energies δ Δg° affects the statistical weight of the completely dissociated duplex w D ) is reflected in the hybridization affinity K(x) and in the hybridization signal θ(x). In addition to single base pair defects our binding affinity model reproduces well our experimental results on the binding affinities of oligonucleotide duplexes with two single base deletion defects (for details see Additional file 1, Fig. S5).
Approximation of the Zipper Model with a Position Dependent Nearest-Neighbor (PDNN) Model
In order to investigate the generality of our finding, we investigate if PDNN models, which fit experimental data well, can be inferred from our model framework. We note that zippering has been previously proposed as the rationale behind the PDNN model in .
In the following we investigate the contribution of each base pair to duplex stability and ask if there is a position-dependent contribution of Watson-Crick NN pairs in the same way as for defects.
In this paper we studied, experimentally and theoretically, the stability of short (l < 26 bp) linear surface-bound oligonucleotide duplexes with single base defects. We demonstrated that the rationale behind positional dependent models of oligonucleotide duplex stability is the partial denaturation of the duplexes. We have shown, that the strong influence of the defect position on mismatch discrimination [11–14, 16, 49] and the influence of the sequence context – beyond nearest neighbors [14, 34] can be quantitatively inferred from a molecular zipper model. Partial (end-domain-)denaturation of the duplex as proposed by us in  as well as in [16, 24, 25] results in a positional influence that is entropic in nature. The zipping process is modulated by the sequential arrangement of the base pairs. The model confirms the observed influence of the sequence context beyond the nearest-neighbors. Further the zipper model provides a theoretical foundation to the positional dependent nearest-neighbor model of Zhang et al. .
In the commonly employed two-state nearest-neighbor model, nucleic acid duplex hybridzation/denaturation is considered to be an all-or-none process. According to literature indeed end-fraying effects are expected to be small beyond three bases , however, in our studied case, we conclude that end-fraying plays a non-negligible role. This is surprising since the dissociation probability of individual base pairs decreases towards the center of the duplex in an exponential fashion (see Additional file 1, Fig. S4) and remains very low for most NN-pairs.
We propose that the effect of the defect position on probe-target binding affinities becomes apparent in the hybridization signal intensities due to the unavoidable probe polydispersity of optical synthesis. It indeed appears that the positional dependence of single base MM discrimination is more commonly observed on photolithographically produced DNA oligonucleotide arrays [11–14] rather than (in large scale studies) on spotted microarrays [40, 50, 51] or in solution-phase experiments. We notice, however, that in small studies (investigating few sequences) a positional influence in solution  and on spotted microarrays  has been reported. The probe polydispersity in our experiments smoothes out the steep sigmoid relation between the hybridization intensity and binding free energy ΔG D that is expected for defect free probes, and explains why (within a relatively broad range of ≈ 20 kcal/mol) variations of the binding free energies – like for example the influence of the defect position – are reflected (by means of an approximately linear relation) in the hybridization signal intensities.
The authors thank Dr. Pramod Pullarkat for many helpful discussions and suggestions on this work. Our research was supported by the University of Bayreuth.
- Nelson JW, Martin FH, Tinoco I: DNA and RNA Oligomer Thermodynamics – The Effect of Mismatched Bases on Double-Helix Stability. Biopolymers 1981, 20(12):2509–2531. 10.1002/bip.1981.360201204View ArticlePubMed
- Patel DJ, Kozlowski SA, Marky LA, Rice JA, Broka C, Dallas J, Itakura K, Breslauer KJ: Structure, dynamics, and energetics of deoxyguanosine-thymidine wobble base pair formation in the self-complementary d(cgtgaattcgcg) duplex in solution. Biochemistry 1982, 21(3):437–444. 10.1021/bi00532a003View ArticlePubMed
- Conner BJ, Reyes AA, Morin C, Itakura K, Teplitz RL, Wallace RB: Detection of sickle-cell beta-s-globin allele by hybridization with synthetic oligonucleotides. Proceedings of the National Academy of Sciences of the United States of America 1983, 80: 278–282. 10.1073/pnas.80.1.278PubMed CentralView ArticlePubMed
- Devoe H, Tinoco I: Stability of helical polynucleotides – base contributions. Journal of Molecular Biology 1962, 4(5):500-&.View ArticlePubMed
- Borer PN, Dengler B, Tinoco I, Uhlenbeck OC: Stability of Ribonucleic-Acid Double-Stranded Helices. Journal of Molecular Biology 1974, 86(4):843–853. 10.1016/0022-2836(74)90357-XView ArticlePubMed
- Breslauer KJ, Frank R, Blocker H, Marky LA: Predicting DNA Duplex Stability from the Base Sequence. Proceedings of the National Academy of Sciences of the United States of America 1986, 83(11):3746–3750. 10.1073/pnas.83.11.3746PubMed CentralView ArticlePubMed
- SantaLucia J, Hicks D: The thermodynamics of DNA structural motifs. Annual Review of Biophysics and Biomolecular Structure 2004, 33: 415–440. 10.1146/annurev.biophys.32.110601.141800View ArticlePubMed
- SantaLucia J: A unified view of polymer, dumbbell, and oligonucleotide DNA nearest-neighbor thermodynamics. Proceedings of the National Academy of Sciences of the United States of America 1998, 95(4):1460–1465. 10.1073/pnas.95.4.1460PubMed CentralView ArticlePubMed
- Markham NR, Zuker M: DINAMelt web server for nucleic acid melting prediction. Nucleic Acids Research 2005, 33: W577-W581. 10.1093/nar/gki591PubMed CentralView ArticlePubMed
- Allawi HT, SantaLucia J: Thermodynamics and NMR of internal GT mismatches in DNA. Biochemistry 1997, 36(34):10581–10594. 10.1021/bi962590cView ArticlePubMed
- Wick LM, Rouillard JM, Whittam TS, Gulari E, Tiedje JM, Hashsham SA: On-chip non-equilibrium dissociation curves and dissociation rate constants as methods to assess specificity of oligonucleotide probes. Nucleic Acids Research 2006, 34(3):e26. 10.1093/nar/gnj024PubMed CentralView ArticlePubMed
- Pozhitkov A, Noble PA, Domazet-Loso T, Nolte AW, Sonnenberg R, Staehler P, Beier M, Tautz D: Tests of rRNA hybridization to microarrays suggest that hybridization characteristics of oligonucleotide probes for species discrimination cannot be predicted. Nucleic Acids Research 2006, 34(9):e66. 10.1093/nar/gkl133PubMed CentralView ArticlePubMed
- Naiser T, Mai T, Michel W, Ott A: Versatile maskless microscope projection photolithography system and its application in light-directed fabrication of DNA microarrays. Review of Scientific Instruments 2006, 77(6):063711. 10.1063/1.2213152View Article
- Naiser T, Kayser J, Mai T, Michel W, Ott A: Impact of point-mutations on the hybridization affinity of surface-bound DNA/DNA and RNA/DNA oligonucleotide duplexes: comparison of single base mismatches and base bulges. BMC Biotechnology 2008, 8: 48. 10.1186/1472-6750-8-48PubMed CentralView ArticlePubMed
- Naiser T: Characterization of Oligonucleotide Microarray Hybridization: Microarray Fabrication by Light-Directed in situ Synthesis – Development of an Automated DNA Microarray Synthesizer, Characterization of Single Base Mismatch Discrimination and the Position-Dependent Influence of Point Defects on Oligonucleotide Duplex Binding Affinities. PhD thesis, University of Bayreuth 2008. [http://opus.ub.uni-bayreuth.de/volltexte/2008/461/]
- Rennie C, Noyes HA, Kemp SJ, Hulme H, Brass A, Hoyle DC: Strong position-dependent effects of sequence mismatches on signal ratios measured using long oligonucleotide microarrays. BMC Genomics 2008, 9: 317. 10.1186/1471-2164-9-317PubMed CentralView ArticlePubMed
- Zhang L, Miles MF, Aldape KD: A model of molecular interactions on short oligonucleotide microarrays. Nature Biotechnology 2003, 21(7):818–821. 10.1038/nbt836View ArticlePubMed
- Naef F, Magnasco MO: Solving the riddle of the bright mismatches: Labeling and effective binding in oligonucleotide arrays. Physical Review E 2003, 68: 011906. 10.1103/PhysRevE.68.011906View Article
- Mei R, Hubbell E, Bekiranov S, Mittmann M, Christians FC, Shen MM, Lu G, Fang J, Liu WM, Ryder T, Kaplan P, Kulp D, Webster TA: Probe selection for high-density oligonucleotide arrays. Proceedings of the National Academy of Sciences of the United States of America 2003, 100(20):11237–11242. 10.1073/pnas.1534744100PubMed CentralView ArticlePubMed
- Binder H, Preibisch S, Kirsten T: Base pair interactions and hybridization isotherms of matched and mismatched oligonucleotide probes on microarrays. Langmuir 2005, 21(20):9287–9302. 10.1021/la051231sView ArticlePubMed
- Held GA, Grinstein G, Tu Y: Relationship between gene expression and observed intensities in DNA microarraysa modeling study. Nucleic Acids Research 2006, 34(9):e70. 10.1093/nar/gkl122PubMed CentralView ArticlePubMed
- Gibbs JH, Dimarzio EA: Statistical mechanics of helix-coil transitions in biological macromolecules. Journal of Chemical Physics 1959, 30: 271–282. 10.1063/1.1729886View Article
- Kittel C: Phase transition of a molecular zipper. American Journal of Physics 1969, 37(9):917-&. 10.1119/1.1975930View Article
- Deutsch JM, Liang S, Narayan O: Modelling of microarray data with zippering. Preprint q-bio.BM/0406039 v1 2004.
- Binder H: Thermodynamics of competitive surface adsorption on DNA microarrays. Journal of Physics-Condensed Matter 2006, 18(18):S491-S523. 10.1088/0953-8984/18/18/S02View Article
- Naiser T, Ehler O, Mai T, Michel W, Ott A: Hybridization to surface-bound oligonucleotide probes: Influence of point defects. 2006.
- Fodor SPA, Read JL, Pirrung MC, Stryer AT, Land Lu, Solas D: Light-directed, spatially addressable parallel chemical synthesis. Science 1991, 251(4995):767–773. 10.1126/science.1990438View ArticlePubMed
- McGall GH, Barone AD, Diggelmann M, Fodor SA, Gentalen E, Ngo N: The efficiency of light-directed synthesis of DNA arrays on glass substrates. Journal of the American Chemical Society 1997, 119(22):5081–5090. 10.1021/ja964427aView Article
- Singh-Gasson S, Green RD, Yue YJ, Nelson C, Blattner F, Sussman MR, Cerrina F: Maskless fabrication of light-directed oligonucleotide microarrays using a digital micromirror array. Nature Biotechnology 1999, 17(10):974–978. 10.1038/13664View ArticlePubMed
- Luebke KJ, Balog RP, Mittelman D, Garner HR: Digital optical chemistry: A novel system for the rapid fabrication of custom oligonucleotide arrays. In Microfabricated Sensors, Application of Optical Technology for DNA Analysis. Volume 815. Edited by: Richard Kordal, Usmani, Wai Tak Law. American Chemical Society Publications; 2002:87–106.View Article
- Nuwaysir EF, Huang W, Albert TJ, Singh J, Nuwaysir K, Pitas A, Richmond T, Gorski T, Berg JP, Ballin J, McCormick M, Norton J, Pollock T, Sumwalt T, Butcher L, Porter D, Molla M, Hall C, Blattner F, Sussman MR, Wallace RL, Cerrina F, Green RD: Gene expression analysis using oligonucleotide arrays produced by maskless photolithography. Genome Research 2002, 12(11):1749–1755. 10.1101/gr.362402PubMed CentralView ArticlePubMed
- Le Berre V, Trevisiol E, Dagkessamanskaia A, Sokol S, Caminade AM, Majoral JP, Meunier B, Francois J: Dendrimeric coating of glass slides for sensitive DNA microarrays analysis. Nucleic Acids Research 2003, 31(16):e88. 10.1093/nar/gng088PubMed CentralView ArticlePubMed
- Hasan A, Stengele KP, Giegrich H, Cornwell P, Isham KR, Sachleben RA, Pfleiderer W, Foote RS: Photolabile protecting groups for nucleosides: Synthesis and photodeprotection rates. Tetrahedron 1997, 53(12):4247–4264. 10.1016/S0040-4020(97)00154-3View Article
- Hall TS, Pancoska P, Riccelli PV, Mandell K, Benight AS: Sequence context and thermodynamic stability of a single base pair mismatch in short deoxyoligonucleotide duplexes. Journal of the American Chemical Society 2001, 123(47):11811–11812. 10.1021/ja016360jView ArticlePubMed
- Craig ME, Crothers DM, Doty P: Relaxation kinetics of dimer formation by self complementary oligonucleotides. Journal of Molecular Biology 1971, 62(2):383-&. 10.1016/0022-2836(71)90434-7View ArticlePubMed
- Halperin A, Buhot A, Zhulina EB: Sensitivity, specificity, and the hybridization isotherms of DNA chips. Biophysical Journal 2004, 86(2):718–730. 10.1016/S0006-3495(04)74150-8PubMed CentralView ArticlePubMed
- Burden CJ, Pittelkow YE, Wilson SR: Statistical Analysis of Adsorption Models for Oligonucleotide Microarrays. Statistical Applications in Genetics and Molecular Biology 2004, 3: 35. 10.2202/1544-6115.1095View Article
- Held GA, Grinstein G, Tu Y: Modeling of DNA microarray data by using physical properties of hybridization. Proceedings of the National Academy of Sciences of the United States of America 2003, 100(13):7575–7580. 10.1073/pnas.0832500100PubMed CentralView ArticlePubMed
- Binder H, Kirsten T, Loeffler M, Stadller PF: Sensitivity of microarray oligonucleotide probes: Variability and effect of base composition. Journal of Physical Chemistry B 2004, 108(46):18003–18014. 10.1021/jp049593gView Article
- Fish DJ, Horne MT, Brewood GP, Goodarzi JP, Alemayehu S, Bhandiwad A, Searles RP, Benight AS: DNA multiplex hybridization on microarrays and thermodynamic stability in solution: a direct comparison. Nucleic Acids Research 2007, 35(21):7197–7208. 10.1093/nar/gkm865PubMed CentralView ArticlePubMed
- Carlon E, Heim T: Thermodynamics of RNA/DNA hybridization in high-density oligonucleotide microarrays. Physica A-Statistical Mechanics and its Applications 2006, 362(2):433–449. 10.1016/j.physa.2005.09.067View Article
- Peterson AW, Wolf LK, Georgiadis RM: Hybridization of mismatched or partially matched DNA at surfaces. Journal of the American Chemical Society 2002, 124(49):14601–14607. 10.1021/ja0279996View ArticlePubMed
- Glazer M, Fidanza JA, McGall GH, Trulson MO, Forman JE, Suseno A, Frank CW: Kinetics of oligonucleotide hybridization to photolithographically patterned DNA arrays. Analytical Biochemistry 2006, 358(2):225–238. 10.1016/j.ab.2006.07.042View ArticlePubMed
- Vainrub A, Pettitt BM: Coulomb blockage of hybridization in two-dimensional DNA arrays. Physical Review E 2002, 66(4):041905. 10.1103/PhysRevE.66.041905View Article
- Halperin A, Buhot A, Zhulina EB: Brush effects on DNA chips: thermodynamics, kinetics, and design guidelines. Biophysical Journal 2005, 89(2):796–811. 10.1529/biophysj.105.063479PubMed CentralView ArticlePubMed
- Jobs M, Fredriksson S, Brookes AJ, Landegren U: Effect of oligonucleotide truncation on sing le-nucleotide distinction by solid-phase hybridization. Analytical Chemistry 2002, 74: 199–202. 10.1021/ac010555sView ArticlePubMed
- Halperin A, Buhot A, Zhulina EB: On the hybridization isotherms of DNA microarrays: the Langmuir model and its extensions. Journal of Physics-Condensed Matter 2006, 18(18):S463-S490. 10.1088/0953-8984/18/18/S01View Article
- Forman JE, Walton ID, Stern D, Rava RP, Trulson MO: Thermodynamics of duplex formation and mismatch discrimination on photolithographically synthesized oligonucleotide arrays. Molecular Modeling of Nucleic Acids 1998, 682: 206–228.View Article
- Dorris DR, Nguyen A, Gieser L, Lockner R, Lublinsky A, Patterson M, Touma E, Sendera TJ, Elghanian R, Mazumder A: Oligodeoxyribonucleotide probe accessibility on a three-dimensional DNA microarray surface and the effect of hybridization time on the accuracy of expression ratios. BMC Biotechnology 2003, 3: 6. 10.1186/1472-6750-3-6PubMed CentralView ArticlePubMed
- Urakawa H, Noble PA, El Fantroussi S, Kelly JJ, Stahl DA: Single-base-pair discrimination of terminal mismatches by using oligonucleotide microarrays and neural network analyses. Applied and Environmental Microbiology 2002, 68: 235–244. 10.1128/AEM.68.1.235-244.2002PubMed CentralView ArticlePubMed
- Weckx S, Carlon E, De Vuyst L, Van Hummelen P: Thermodynamic behavior of short oligonucleotides in microarray hybridizations can be described using Gibbs free energy in a nearest-neighbor model. Journal of Physical Chemistry B 2007, 111(48):13583–13590. 10.1021/jp075197xView Article
- Kierzek R, Burkard ME, Turner DH: Thermodynamics of single mismatches in RNA duplexes. Biochemistry 1999, 38(43):14214–14223. 10.1021/bi991186lView ArticlePubMed
- Zhu J, Wartell RM: The effect of base sequence on the stability of RNA and DNA single base bulges. Biochemistry 1999, 38(48):15986–15993. 10.1021/bi9916372View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.