Modeling of DNA hybridization kinetics for spatially resolved biochips.

The marriage of microfluidics with detection technologies that rely on highly selective nucleic acid hybridization will provide improvements in bioanalytical methods for purposes such as detection of pathogens or mutations and drug screening. The capability to deliver samples in a controlled manner across a two-dimensional hybridization detection platform represents a substantial technical challenge in the development of quantitative and reusable biochips. General theoretical and numerical models of heterogeneous hybridization kinetics are required in order to design and optimize such biochips and to develop a quantitative method for online interpretation of experimental results. In this work we propose a general kinetic model of heterogeneous hybridization and develop a technique for estimating the kinetic coefficients for the case of well-spaced, noninteracting surface-bound probes. The experimentally verified model is then incorporated into the BLOCS (biolab-on-a-chip simulation) 3D microfluidics finite element code and used to model the dynamic hybridization on a biochip surface in the presence of a temperature gradient. These simulations demonstrate how such a device can be used to discriminate between fully complementary and single-base-pair mismatched hybridization using fluorescence detection by interpretation of the unique spatially resolved intensity pattern. It is also shown how the dynamic transport of the targets is likely to affect the rate and location of hybridization as well as that, although nonspecific hybridization is present, the change in the concentration of hybridized targets over the sensor platform is sufficiently high to determine if a fully complementary match is present. Practical design information such as the optimum transport speed, target concentration, and channel height is presented. The results presented here will aid in the interpretation of results obtained with such a temperature-gradient biochip.