Dynamics of Evolutionary Optimization

General criteria of selection are derived from the kinetic equations of polynucleotide replication. As an illustrative example we discuss replication in the continuously stirred tank reactor (CSTR). The total rate of RNA synthesis is optimized during selection. The conjecture that the rate of approach towards the stable steady state is a maximum can be easily disproved. It is possible, nevertheless, to derive a potential function for polynucleotide replication in the CSTR. Following a method first introduced by Shahshahani we define a non Euclidean metric on the space of polynucleotide concentrations. In this space with a Riemannian metric the systems follows the corresponding generalized gradient during the process of selection and, therefore, the rate of ascent is now maximum. Potential functions can be derived also for some second order autocatalytic systems which are of interest in evolution, for a multidimensional Schloegl model in the CSTR and, as originally has been shown by Shahshahani, for the Fisher-Haldane-Wright equation of population genetics. In the general case, however, second order autocatalysis is not compatible with the existence of a potential. The elementary hypercycle is discussed as one simple example of a reaction network whose dynamics cannot be described by means of a generalized gradient system. Finite population size introduces a stochastic element into the selection process. Under certain conditions fluctuations in particle numbers become extremely important for the dynamics of selection. Two examples of this kind are: kinetic degeneracy of rate constants and low accuracy of replication.

[1]  J. S. McCaskill,et al.  A localization threshold for macromolecular quasispecies from continuously distributed replication rates , 1984 .

[2]  P. Schuster,et al.  Random selection — a simple model based on linear birth and death processes , 1984 .

[3]  W. Ebeling,et al.  Stochastic Theory of Molecular Replication Processes with Selection Character , 1977 .

[4]  S. Shahshahani,et al.  A New Mathematical Framework for the Study of Linkage and Selection , 1979 .

[5]  P. Schuster,et al.  Self-replication with errors. A model for polynucleotide replication. , 1982, Biophysical chemistry.

[6]  H. Leung,et al.  Stochastic analysis of a nonlinear model for selection of biological macromolecules , 1981 .

[7]  F. Schlögl Chemical reaction models for non-equilibrium phase transitions , 1972 .

[8]  Colin J. Thompson,et al.  On Eigen's theory of the self-organization of matter and the evolution of biological macromolecules , 1974 .

[9]  P. Schuster Polynucleotide Replication and Biological Evolution , 1984 .

[10]  Selection under random mutations in stochastic Eigen model , 1982 .

[11]  Reinhart Heinrich,et al.  Analysis of the selection equations for a multivariable population model. Deterministic and stochastic solutions and discussion of the approach for populations of self-reproducing biochemical networks , 1981 .

[12]  M. Kimura The Neutral Theory of Molecular Evolution: Introduction , 1983 .

[13]  S. Spiegelman,et al.  An approach to the experimental analysis of precellular evolution , 1971, Quarterly Reviews of Biophysics.

[14]  P. Schuster,et al.  The Dynamics of Catalytic Hypercycles — A Stochastic Simulation , 1984 .

[15]  P. Schuster,et al.  Random Selection and the Neutral Theory — Sources of Stochasticity in Replication , 1984 .

[16]  Some principles governing selection in self-reproducing macromolecular systems , 1978, Journal of mathematical biology.

[17]  Manfred Eigen,et al.  Macromolecular Evolution: Dynamical Ordering in Sequence Space , 1985 .

[18]  Josef Hofbauer,et al.  Dynamics of Linear and Nonlinear Autocatalysis and Competition , 1984 .

[19]  Manfred Eigen,et al.  Kinetics of ribonucleic acid replication. , 1983 .

[20]  K Sigmund,et al.  Dynamical systems under constant organization I. Topological analysis of a family of non-linear differential equations--a model for catalytic hypercycles. , 1978, Bulletin of mathematical biology.