Solution of the Crow-Kimura and Eigen models for alphabets of arbitrary size by Schwinger spin coherent states

Abstract

To represent the evolution of nucleic acid and protein sequence, we express the parallel and Eigen models for molecular evolution in terms of a functional integral representation with an $h$-letter alphabet, lifting the two-state, purine/pyrimidine assumption often made in quasi-species theory. For arbitrary $h$ and a general mutation scheme, we obtain the solution of this model in terms of a maximum principle. Euler's theorem for homogeneous functions is used to derive this `thermodynamic' formulation of evolution. The general result for the parallel model reduces to known results for the purine/pyrimidine $h=2$ alphabet and the nucleic acid $h=4$ alphabet for the Kimura 3 ST mutation scheme. Examples are presented for the $h=4$ and $h=20$ cases. We derive the maximum principle for the Eigen model for general $h$. The general result for the Eigen model reduces to a known result for $h=2$. Examples are presented for the nucleic acid $h=4$ and the amino acid $h=20$ alphabet. An error catastrophe phase transition occurs in these models, and the order of the phase transition changes from second to first order for smooth fitness functions when the alphabet size is increased beyond two letters to the generic case. As examples, we analyze the general analytic solution for sharp peak, linear, quadratic, and quartic fitness functions.

50 pages, 8 figures, to appear in J. Stat. Phys; some typos fixed

References in corpus (7)

• Quasispecies Theory for Multiple-Peak Fitness Landscapes
• Eigen model as a quantum spin chain: exact dynamics
• Quasispecies Theory for Horizontal Gene Transfer and Recombination
• Four-state quantum chain as a model of sequence evolution
• An asymptotic maximum principle for essentially linear evolution models
• A maximum principle for the mutation-selection equilibrium of nucleotide sequences
• Maximum principle and mutation thresholds for four-letter sequence evolution

Cited by in corpus (4)

• Quasispecies theory for finite populations
• Modularity Enhances the Rate of Evolution in a Rugged Fitness Landscape
• Semiconservative quasispecies equations for polysomic genomes: The general case
• Diagrammatic expressions for steady-state distribution and static responses in population dynamics