Imperfect DNA lesion repair in the semiconservative quasispecies model: Derivation of the Hamming class equations and solution of the single-fitness peak landscape

Abstract

This paper develops a Hamming class formalism for the semiconservative quasispecies equations with imperfect lesion repair, first presented and analytically solved in Y. Brumer and E.I. Shakhnovich (q-bio.GN/0403018, 2004). Starting from the quasispecies dynamics over the space of genomes, we derive an equivalent dynamics over the space of ordered sequence pairs. From this set of equations, we are able to derive the infinite sequence length form of the dynamics for a class of fitness landscapes defined by a master genome. We use these equations to solve for a generalized single-fitness-peak landscape, where the master genome can sustain a maximum number of lesions and remain viable. We determine the mean equilibrium fitness and error threshold for this class of landscapes, and show that when lesion repair is imperfect, semiconservative replication displays characteristics from both conservative replication and semiconservative replication with perfect lesion repair. The work presented here provides a formulation of the model which greatly facilitates the analysis of a relatively broad class of fitness landscapes, and thus serves as a convenient springboard into biological applications of imperfect lesion repair.