The currently used microarray technology, where signal intensity after stringent washing is a measure for the number of hybridized targets is already in a mature state. There is still an ongoing discussion on quality improvements of microarray experiments. In general the measured intensity is
I = I_{Sp}(c) + I_{Bg}
where I_{Sp}(c) is the specific signal due to hybridization of the complementary probetarget duplex and I_{Bg} is the background due to nonspecific effects. I_{Sp} depends on the concentration c of the complementary strand in solution (target). For hybridization to surface immobilized nucleic acids there are several factors influencing the unspecific term I_{Bg} e.g. cross hybridization, optical effects, microarray surface composition and intra and interarray process variations. The discrimination of specific and nonspecific hybridization of similar targets is still a challenge in microarray probe design and data analysis. Further the development of advanced model based algorithms for background determination is still a matter of research.
Microarray based realtime data acquisition during hybridization and melting tries to overcome these problems by recording the change of specific intensity signal I_{Sp} during the hybridization and denaturation relative to the unspecific signal I_{Bg}. Thus it gives the possibility to monitor melting of nucleic acid duplexes within the environment of the application. This way the fact that hybridization kinetics can provide information for the discrimination between specific and unspecific hybridization and for quantification of targets can be exploited. The multiple measurements of intensity changes opens the door for advanced signal processing and noise reduction techniques, modelbased fitting can solve the inherent weakness of single sampling measurement.
Here we report on the development of a physical model based nonlinear regression algorithm for realtime microarray data analysis and the implementation within a Matlab programmed GUI. Based on the experimental data of 40 experiments we used equations of physicochemical systems with similar behavior as fit function for the nonlinear regression putting the focus on highdensity microarray applications.
Kinetic
Currently it is not possible to give a numerical description of the complex hydrodynamic transport processes. Thus nonlinear regression to physicochemical model systems was used.
The intensity of kinetic measurements in microarray experiments is not necessarily monotonically increasing. In case of competitive hybridization a target with a high reaction rate but lower binding energy is replaced during hybridization process with targets of a higher binding energy
The kinetic parameters are being calculated as follow
I(X) = 
a * k_{h} 
e^{ kh (X  x0)}  e^{ kc (X  x0)} 
k_{c}  k_{h} 
For competitive hybridization the two constants k_{h} and k_{c} give the reaction rates of the two consecutive steps of hybridization and competitive displacement.
In case of non competitive hybridization k_{c} = 0 the formula is still valid.
Melting
The fitting of the melting curves was based on a convolution of linear functions for the description of the temperature behavior of the dye and an exponential distribution function.
I(T) = 
k_{a }(k_{Cy } x) 
 (k_{HT} * T) + k_{const} 
e ^{kb (x  ktm)} +1 
with the intensity I(T) and the fitting parameters k_{tm}, k_{a}, k_{b}, k_{Cy}, k_{HT} and kconst. k_{a} is related to the amount of hybridized target. k_{b} is related to the slope of the melting curve at melting temperature (Tm). The larger kb the steeper is the melting transition. In general, k_{b} is larger for longer oligos. k_{Cy} considers the linear temperature behavior of the dye. The term k_{HT} * T + k_{const} is a first order correction considering the remaining fluorescent intensity at temperatures above the melting of the nucleic acids due to unspecific binding of fluorophores to the probe spots. The whole nonlinear behavior of the system is described by the exponential function with the two parameters k_{tm} and k_{b}. The parameter k_{tm} gives the melting temperature.
