On the Growth Rate of Non-Enzymatic Molecular Replicators

It is well known that non-enzymatic template directed molecular replicators X + nO ---> 2X exhibit parabolic growth d[X]/dt = k [X]^{1/2}. Here, we analyze the dependence of the effective replication rate constant k on hybridization energies, temperature, strand length, and sequence composition. First we derive analytical criteria for the replication rate k based on simple thermodynamic arguments. Second we present a Brownian dynamics model for oligonucleotides that allows us to simulate their diffusion and hybridization behavior. The simulation is used to generate and analyze the effect of strand length, temperature, and to some extent sequence composition, on the hybridization rates and the resulting optimal overall rate constant k. Combining the two approaches allows us to semi-analytically depict a fitness landscape for template directed replicators. The results indicate a clear replication advantage for longer strands at low temperatures.