Approximate analytic solution of a simple system of kinetic equations describing the enzymatic synthesis of nucleic acids.

Macroscopic kinetic models describing the process of polymerase chain reaction (PCR) are currently solved only by numerical methods, which hampers the development of effective software algorithms for processing the results of the reaction. This paper considers the application of the homotopy perturbation method for obtaining approximate analytical solution of the simplest system of enzymatic kinetic equations describing the synthesis of nucleic acid molecules during PCR. The resulting approximate analytic solution with high accuracy reproduces the results of a numerical solution of the system in a wide range of ratios of enzyme and substrate concentrations both for the case of a large excess of the substrate over the enzyme and vice versa.