Mathematical model of real-time PCR kinetics

Several real-time PCR (rtPCR) quantification techniques are currently used to determine the expression levels of individual genes from rtPCR data in the form of fluorescence intensities. In most of these quantification techniques, it is assumed that the efficiency of rtPCR is constant. Our analysis of rtPCR data shows, however, that even during the exponential phase of rtPCR, the efficiency of the reaction is not constant, but is instead a function of cycle number. In order to understand better the mechanisms belying this behavior, we have developed a mathematical model of the annealing and extension phases of the PCR process. Using the model, we can simulate the PCR process over a series of reaction cycles. The model thus allows us to predict the efficiency of rtPCR at any cycle number, given a set of initial conditions and parameter values, which can mostly be estimated from biophysical data. The model predicts a precipitous decrease in cycle efficiency when the product concentration reaches a sufficient level for template-template re-annealing to compete with primer-template annealing; this behavior is consistent with available experimental data. The quantitative understanding of rtPCR provided by this model can allow us to develop more accurate methods to quantify gene expression levels from rtPCR data.