A diffusion-reaction model for DNA microarray assays

DNA microarrays are extensively used for the quantification of the degree of differential mRNA expression. The assay involves hybridization of mobile DNA strands with immobilized complementary DNA strands to form duplexes. The overall duplex formation rate depends on the rate of transport of strands in solution to the corresponding spot on the surface, and the rate of the hybridization reaction. We present a theoretical model that incorporates both kinetics of the reversible hybridization reaction and diffusional transport of the labeled strands, and analyze DNA microarray hybridization using this model. Simulations are carried out in a geometrically realistic domain for labeled DNA concentrations corresponding to rare and abundant transcripts for typical assay conditions. The rate of strand diffusion in solution is shown to strongly affect the overall hybridization rate. We compute the minimum inter-spot spacing for replicate spots to enhance sensitivity. We also determine the hybridization time for which reliable estimates of the relative mRNA abundance of two species can be obtained using total fluorescence intensities. An analytical solution for the concentration distribution of mobile strands at intermediate hybridization times provides a convenient tool to calculate the mobile strand concentration profiles. This model provides a framework for the process analysis of all microarray assays currently used for genomic transcriptional analysis.