Computer simulation study of probe-target hybridization in model DNA microarrays: effect of probe surface density and target concentration.

We use lattice Monte Carlo simulations to study the thermodynamics of hybridization of single-stranded "target" genes in solution with complementary "probe" DNA molecules immobilized on a microarray surface. The target molecules in our system contain 48 segments and the probes tethered on a hard surface contain 8-24 segments. The segments on the probe and target are distinct, with each segment representing a sequence of nucleotides that interacts exclusively with its unique complementary target segment with a single hybridization energy; all other interactions are zero. We examine how surface density (number of probes per unit surface area) and concentration of target molecules affect the extent of hybridization. For short probe lengths, as the surface density increases, the probability of binding long stretches of target segments increases at low surface density, reaches a maximum at an intermediate surface density, and then decreases at high surface density. Furthermore, as the surface density increases, the target is less likely to bind completely to one probe; instead, it binds simultaneously to multiple probes. At short probe lengths, as the target concentration increases, the fraction of targets binding completely to the probes (specificity) decreases. At long probe lengths, varying the target concentration does not affect the specificity. At all target concentrations as the probe length increases, the fraction of target molecules bound to the probes by at least one segment (sensitivity) increases while the fraction of target molecules completely bound to the probes (specificity) decreases. This work provides general guidelines to maximizing microarray sensitivity and specificity. Our results suggest that the sensitivity and specificity can be maximized by using probes 130-180 nucleotides long at a surface density in the range of 7 x 10(-5)- 3 x 10(-4) probe molecules per nm(2).