A particular case of the famous quasispecies model - the Crow-Kimura model with a permutation invariant fitness landscape - is investigated. Using the fact that the mutation matrix in the case of a permutation invariant fitness landscape has a special tridiagonal form, a change of the basis is suggested such that in the new coordinates a number of analytical results can be obtained. In particular, using the eigenvectors of the mutation matrix as the new basis, we show that the quasispecies distribution approaches a binomial one and give simple estimates for the speed of convergence. Another consequence of the suggested approach is a parametric solution to the system of equations determining the quasispecies. Using this parametric solution we show that our approach leads to exact asymptotic results in some cases, which are not covered by the existing methods. In particular, we are able to present not only the limit behavior of the leading eigenvalue (mean population fitness), but also the exact formulas for the limit quasispecies eigenvector for special cases. For instance, this eigenvector has a geometric distribution in the case of the classical single peaked fitness landscape. On the biological side, we propose a mathematical definition, based on the closeness of the quasispecies to the binomial distribution, which can be used as an operational definition of the notorious error threshold. Using this definition, we suggest two approximate formulas to estimate the critical mutation rate after which the quasispecies delocalization occurs.
Analytical analysis of spatially extended autocatalytic and hypercyclic systems is presented. It is shown that spatially explicit systems in the form of reaction-diffusion equations with global regulation possess the same major qualitative features as the corresponding local models. In particular, using the introduced notion of the stability in the mean integral sense we prove the competitive exclusion principle for the autocatalytic system and the permanence for the hypercycle system. Existence and stability of stationary solutions are studied. For some parameter values it is proved that stable spatially non-uniform solutions appear.
We study general properties of the leading eigenvalue w(q) of Eigen's evolutionary matrices depending on the probability q of faithful reproduction. This is a linear algebra problem that has various applications in theoretical biology, including such diverse fields as the origin of life, evolution of cancer progression, and virus evolution. We present the exact expressions for w(q), w ′ (q), w ′′ (q) for q = 0, 0.5, 1 and prove that the absolute minimum of w(q), which always exists, belongs to the interval [0, 0.5]. For the specific case of a single peaked landscape we also find lower and upper bounds on w(q), which are used to estimate the critical mutation rate, after which the distribution of the types of individuals in the population becomes almost uniform. This estimate is used as a starting point to conjecture another estimate, valid for any fitness landscape, and which is checked by numerical calculations. The last estimate stresses the fact that the inverse dependence of the critical mutation rate on the sequence length is not a generally valid fact. Therefore, the discussions of the error threshold applied to biological systems must take this fact into account.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.