Alessandro Filisetti

The Stochastic Evolution of a Protocell: The Gillespie Algorithm in a Dynamically Varying Volume

We propose an improvement of the Gillespie algorithm allowing us to study the time evolution of an ensemble of chemical reactions occurring in a varying volume, whose growth is directly related to the amount of some specific molecules, belonging to the reactions set. This allows us to study the stochastic evolution of a protocell, whose volume increases because of the production of container molecules. Several protocell models are considered and compared with the deterministic models.

Sufficient conditions for emergent synchronization in protocell models.

In this paper, we study general protocell models aiming to understand the synchronization phenomenon of genetic material and container productions, a necessary condition to ensure sustainable growth in protocells and eventually leading to Darwinian evolution when applied to a population of protocells. Synchronization has been proved to be an emergent property in many relevant protocell models in the class of the so-called surface reaction models, assuming both linear- and non-linear dynamics for the involved chemical reactions. We here extend this analysis by introducing and studying a new class of models where the relevant chemical reactions are assumed to occur inside the protocell, in contrast with the former model where the reaction site was the external surface. While in our previous studies, the replicators were assumed to compete for resources, without any direct interaction among them, we here improve both models by allowing linear interaction between replicators: catalysis and/or inhibition. Extending some techniques previously introduced, we are able to give a quite general analytical answer about the synchronization phenomenon in this more general context. We also report on results of numerical simulations to support the theory, where applicable, and allow the investigation of cases which are not amenable to analytical calculations.

A stochastic model of the emergence of autocatalytic cycles

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.

The detection of intermediate-level emergent structures and patterns.

Artificial life is largely concerned with systems that exhibit different emergent phenomena; yet, the identification of emergent structures is frequently a difficult challenge. In this paper we introduced a system to identify candidate emergent mesolevel dynamical structures in dynamical networks. This method is based on an extension of a measure introduced for detecting clusters in biological neural networks; its main novelty in comparison to previous application of similar measures is that we used it to consider truly dynamical networks, and not only fluctuations around stable asymptotic states. The identified structures are clusters of elements that behave in a coherent and coordinated way and that loosely interact with the remainder of the system. We have evidence that our approach is able to identify these “emerging things” in some artificial network models and in more complex data coming from catalytic reaction networks and biological gene regulatory systems (A.thaliana). We think that this system could suggest interesting new ways in dealing with artificial and biological systems.

The search for candidate relevant subsets of variables in complex systems

We describe a method to identify relevant subsets of variables, useful to understand the organization of a dynamical system. The variables belonging to a relevant subset should have a strong integration with the other variables of the same relevant subset, and a much weaker interaction with the other system variables. On this basis, extending previous work on neural networks, an information-theoretic measure, the dynamical cluster index, is introduced in order to identify good candidate relevant subsets. The method does not require any previous knowledge of the relationships among the system variables, but relies on observations of their values over time. We show its usefulness in several application domains, including: (i) random Boolean networks, where the whole network is made of different subnetworks with different topological relationships (independent or interacting subnetworks); (ii) leader-follower dynamics, subject to noise and fluctuations; (iii) catalytic reaction networks in a flow reactor; (iv) the MAPK signaling pathway in eukaryotes. The validity of the method has been tested in cases where the data are generated by a known dynamical model and the dynamical cluster index is applied in order to uncover significant aspects of its organization; however, it is important that it can also be applied to time series coming from field data without any reference to a model. Given that it is based on relative frequencies of sets of values, the method could be applied also to cases where the data are not ordered in time. Several indications to improve the scope and effectiveness of the dynamical cluster index to analyze the organization of complex systems are finally given.

Growth and Division in a Dynamic Protocell Model

In this paper a new model of growing and dividing protocells is described, whose main features are (i) a lipid container that grows according to the composition of the molecular milieu (ii) a set of “genetic memory molecules” (GMMs) that undergo catalytic reactions in the internal aqueous phase and (iii) a set of stochastic kinetic equations for the GMMs. The mass exchange between the external environment and the internal phase is described by simulating a semipermeable membrane and a flow driven by the differences in chemical potentials, thereby avoiding to resort to sometimes misleading simplifications, e.g., that of a flow reactor. Under simple assumptions, it is shown that synchronization takes place between the rate of replication of the GMMs and that of the container, provided that the set of reactions hosts a so-called RAF (Reflexive Autocatalytic, Food-generated) set whose influence on synchronization is hereafter discussed. It is also shown that a slight modification of the basic model that takes into account a rate-limiting term, makes possible the growth of novelties, allowing in such a way suitable evolution: so the model represents an effective basis for understanding the main abstract properties of populations of protocells.

Exploring the organisation of complex systems through the dynamical interactions among their relevant subsets

Complex systems often show forms of organisation where a clear-cut hierarchy of levels with a well-defined direction of information flow cannot be found. In this paper we propose an information-theoretic method aimed at identifying the dynamically relevant parts of a system along with their relationships, interpreting in such a way the system’s dynamical organisation. The analysis is quite general and can be applied to many dynamical systems. We show here its application to two relevant biological examples, the case of mammalian cell cycle network and of Mitogen Activated Protein Kinase (MAPK) cascade. The result of our analysis shows that the elements of the mammalian cell cycle network act as a single compact group, whereas the MAPK system can be decomposed into two dynamically distinct parts, with asymmetric information flows.

SYNCHRONIZATION PHENOMENA IN PROTOCELL MODELS

This work aims to consider simplified models of protocells in order to describe their general behaviors. The advantage of the modelling approach is that the early protocells of life-forms on Earth are not reproduced in the present time. However, the problem is considered as a right track to understand the origin of life as well as to work with more objective synthesis of new drugs.

Research article: Parameter sensitivity analysis of stochastic models: Application to catalytic reaction networks

A general numerical methodology for parametric sensitivity analysis is proposed, which allows to determine the parameters exerting the greatest influence on the output of a stochastic computational model, especially when the knowledge about the actual value of a parameter is insufficient. An application of the procedure is performed on a model of protocell, in order to detect the kinetic rates mainly affecting the capability of a catalytic reaction network enclosed in a semi-permeable membrane to retain material from its environment and to generate a variety of molecular species within its boundaries. It is shown that the former capability is scarcely sensitive to variations in the model parameters, whereas a kinetic rate responsible for profound modifications of the latter can be identified and it depends on the specific reaction network. A faster uptaking of limited resources from the environment may have represented a significant advantage from an evolutionary point of view and this result is a first indication in order to decipher which kind of structures are more suitable to achieve a viable evolution.

Dynamical Criticality: Overview and Open Questions

Systems that exhibit complex behaviours are often found in a particular dynamical condition, poised between order and disorder. This observation is at the core of the so-called criticality hypothesis, which states that systems in a dynamical regime between order and disorder attain the highest level of computational capabilities and achieve an optimal trade-off between robustness and flexibility. Recent results in cellular and evolutionary biology, neuroscience and computer science have revitalised the interest in the criticality hypothesis, emphasising its role as a viable candidate general law in adaptive complex systems. This paper provides an overview of the works on dynamical criticality that are — To the best of our knowledge — Particularly relevant for the criticality hypothesis. The authors review the main contributions concerning dynamics and information processing at the edge of chaos, and illustrate the main achievements in the study of critical dynamics in biological systems. Finally, the authors discuss open questions and propose an agenda for future work.