Search

Autocatalytic cycles are rather widespread in nature and in several theoretical models of catalytic reaction networks their emergence is hypothesized to be inevitable when the network is or becomes sufficiently complex. Nevertheless, the emergence of autocatalytic cycles has been never observed in wet laboratory experiments. Here, we present a novel model of catalytic reaction networks with the explicit goal of filling the gap between theoretical predictions and experimental findings. The model is based on previous study of Kauffman, with new features in the introduction of a stochastic algorithm to describe the dynamics and in the possibility to increase the number of elements and reactions according to the dynamical evolution of the system. Furthermore, the introduction of a temporal threshold allows the detection of cycles even in our context of a stochastic model with asynchronous update. In this study, we describe the model and present results concerning the effect on the overall dynamics of varying (a) the average residence time of the elements in the reactor, (b) both the composition of the firing disk and the concentration of the molecules belonging to it, (c) the composition of the incoming flux. © Springer-Verlag 2011.

Autocatalytic cycles are rather common in biological systems and they might have played a major role in the transition from non-living to living systems. Several theoretical models have been proposed to address the experimentalists during the investigation of this issue and most of them describe a phase transition depending upon the level of heterogeneity of the chemical soup. Nevertheless, it is well known that reproducing the emergence of autocatalytic sets in wet laboratories is a hard task. Understanding the rationale at the basis of such a mismatch between theoretical predictions and experimental observations is therefore of fundamental importance. We here introduce a novel stochastic model of catalytic reaction networks, in order to investigate the emergence of autocatalytic cycles, sensibly considering the importance of noise, of small-number effects and the possible growth of the number of different elements in the system. Furthermore, the introduction of a temporal threshold that defines how long a specific reaction is kept in the reaction graph allows to univocally define cycles also within an asynchronous framework. The foremost analyses have been focused on the study of the variation of the composition of the incoming flux. It was possible to show that the activity of the system is enhanced, with particular regard to the emergence of autocatalytic sets, if a larger number of different elements is present in the incoming flux, while the specific length of the species seems to entail minor effects on the overall dynamics. © 2011 Filisetti et al; licensee Chemistry Central Ltd.

Recently some methods have been presented to extract ordinary differential equations (ODE) directly from an experimental time series. Here, we introduce a new method to find an ODE which models both the short time and the long time dynamics. The experimental data are represented in a state space and the corresponding flow vectors are approximated by polynomials of the state vector components. We apply these methods both to simulated data and experimental data from human limb movements, which like many other biological systems can exhibit limit cycle dynamics. In systems with only one oscillator there is excellent agreement between the limit cycling displayed by the experimental system and the reconstructed model, even if the data are very noisy. Furthermore we study systems of two coupled limit cycle oscillators. There, a reconstruction was only successful for data with a sufficiently long transient trajectory and relatively low noise level.