Modeling the polymerase chain reaction.
We introduce a mathematical model to treat the polymerase chain reaction (PCR), where we regard the accumulation of new molecules during a PCR cycle as a randomly bifurcating tree. This model enables us to compute an approximate formula for the distribution of the number of replications that have occurred between a pair of molecules, which depends on the efficiency lambda of the reaction, the number N0 of template molecules at the beginning of the PCR and the number c of PCR cycles. The reliability of the approximation is tested by computer simulations. Finally, to model the effect of the intrinsic error rate of the polymerase, we superimpose a substitution process on the tree. The resulting closed formula for the distribution of pairwise differences of sequences as a function of error rate mu and efficiency lambda can be used to estimate the error rate, if lambda is known.
D001483:Base Sequence; D004247:DNA / Q000737:chemistry; D003663:Decision Trees; D008962:Models, Theoretical; D016133:Polymerase Chain Reaction; D011336:Probability; D015203:Reproducibility of Results; D013698:Templates, Genetic
D004247:DNA