A significant amount of experimental data on the reaction kinetics for the mitochondrial DNA polymerase gamma exist, but interpreting that data is difficult due to the complex nature of the function of the polymerase. In order to model how these measured kinetics values for polymerase gamma affect the final function of the polymerase, the replication of an entire strand of mitochondrial DNA, we implement a stochastic simulation of the series of reaction events that the polymerase carries out. These reactions include the correct and incorrect polymerization events, exonuclease events which may remove both incorrectly and correctly matched base pairs, and the disassociation of the polymerase from the mitochondrial DNA template. We also describe other reactions which may be included, such as the addition of nucleoside analog tri-phosphates as substrates. The simulation analysis of the kinetics data is implemented through a standard Gillespie algorithm. We describe the methods necessary to define, code and test this algorithm, as well as describing the hardware and software options that are available.
Copyright (c) 2010 Elsevier Inc. All rights reserved.