Growth states of catalytic reaction networks exhibiting energy metabolism

All cells derive nutrition by absorbing some chemical and energy resources from the environment; these resources are used by the cells to reproduce the chemicals within them, which in turn leads to an increase in their volume. In this study, we introduce a protocell model exhibiting catalytic reaction dynamics, energy metabolism, and cell growth. Results of extensive simulations of this model show the existence of four phases with regard to the rates of both the influx of resources and the cell growth. These phases include an active phase with high influx and high growth rates, an inefficient phase with high influx but low growth rates, a quasi-static phase with low influx and low growth rates, and a death phase with negative growth rate. A mean field model well explains the transition among these phases as bifurcations. The statistical distribution of the active phase is characterized by a power law and that of the inefficient phase is characterized by a nearly equilibrium distribution. We also discuss the relevance of the results of this study to distinct states in the existing cells.Comment: 21 pages, 5 figure

A stochastic model of catalytic reaction networks in protocells

Protocells are supposed to have played a key role in the self-organizing processes leading to the emergence of life. Existing models either (i) describe protocell architecture and dynamics, given the existence of sets of collectively self-replicating molecules for granted, or (ii) describe the emergence of the aforementioned sets from an ensemble of random molecules in a simple experimental setting (e.g. a closed system or a steady-state flow reactor) that does not properly describe a protocell. In this paper we present a model that goes beyond these limitations by describing the dynamics of sets of replicating molecules within a lipid vesicle. We adopt the simplest possible protocell architecture, by considering a semi-permeable membrane that selects the molecular types that are allowed to enter or exit the protocell and by assuming that the reactions take place in the aqueous phase in the internal compartment. As a first approximation, we ignore the protocell growth and division dynamics. The behavior of catalytic reaction networks is then simulated by means of a stochastic model that accounts for the creation and the extinction of species and reactions. While this is not yet an exhaustive protocell model, it already provides clues regarding some processes that are relevant for understanding the conditions that can enable a population of protocells to undergo evolution and selection.Comment: 20 pages, 5 figure

recent theoretical approaches to minimal artificial cells

Minimal artificial cells (MACs) are self-assembled chemical systems able to mimic the behavior of living cells at a minimal level, i.e. to exhibit self-maintenance, self-reproduction and the capability of evolution. The bottom-up approach to the construction of MACs is mainly based on the encapsulation of chemical reacting systems inside lipid vesicles, i.e. chemical systems enclosed (compartmentalized) by a double-layered lipid membrane. Several researchers are currently interested in synthesizing such simple cellular models for biotechnological purposes or for investigating origin of life scenarios. Within this context, the properties of lipid vesicles (e.g., their stability, permeability, growth dynamics, potential to host reactions or undergo division processes…) play a central role, in combination with the dynamics of the encapsulated chemical or biochemical networks. Thus, from a theoretical standpoint, it is very important to develop kinetic equations in order to explore first—and specify later—the conditions that allow the robust implementation of these complex chemically reacting systems, as well as their controlled reproduction. Due to being compartmentalized in small volumes, the population of reacting molecules can be very low in terms of the number of molecules and therefore their behavior becomes highly affected by stochastic effects both in the time course of reactions and in occupancy distribution among the vesicle population. In this short review we report our mathematical approaches to model artificial cell systems in this complex scenario by giving a summary of three recent simulations studies on the topic of primitive cell (protocell) systems

Necessary and sufficient conditions for protocell growth

International audienceWe consider a generic protocell model consisting of any conservative chemical reaction network embedded within a membrane. The membrane results from the self-assembly of a membrane precursor and is semi-permeable to some nutrients. Nutrients are metabolized into all other species including the membrane precursor, and the membrane grows in area and the protocell in volume. Faithful replication through cell growth and division requires a doubling of both cell volume and surface area every division time (thus leading to a periodic surface area-to-volume ratio) and also requires periodic concentrations of the cell constituents. Building upon these basic considerations, we prove necessary and sufficient conditions pertaining to the chemical reaction network for such a regime to be met. A simple necessary condition is that every moiety must be fed. A stronger necessary condition implies that every siphon must be either fed, or connected to species outside the siphon through a pass reaction capable of transferring net positive mass into the siphon. And in the case of nutrient uptake through passive diffusion and of constant surface area-to-volume ratio, a sufficient condition for the existence of a fixed point is that every siphon be fed. These necessary and sufficient conditions hold for any chemical reaction kinetics, membrane parameters or nutrient flux diffusion constants

Stochastic simulations of minimal cells: the Ribocell model

<p>Abstract</p> <p>Background</p> <p>Over the last two decades, lipid compartments (liposomes, lipid-coated droplets) have been extensively used as in vitro "minimal" cell models. In particular, simple and complex biomolecular reactions have been carried out inside these self-assembled micro- and nano-sized compartments, leading to the synthesis of RNA and functional proteins inside liposomes. Despite this experimental progress, a detailed physical understanding of the underlying dynamics is missing. In particular, the combination of solute compartmentalization, reactivity and stochastic effects has not yet been clarified. A combination of experimental and computational approaches can reveal interesting mechanisms governing the behavior of micro compartmentalized systems, in particular by highlighting the intrinsic stochastic diversity within a population of "synthetic cells".</p> <p>Methods</p> <p>In this context, we have developed a computational platform called ENVIRONMENT suitable for studying the stochastic time evolution of reacting lipid compartments. This software - which implements a Gillespie Algorithm - is an improvement over a previous program that simulated the stochastic time evolution of homogeneous, fixed-volume, chemically reacting systems, extending it to more general conditions in which a collection of similar such systems interact and change over the course of time. In particular, our approach is focused on elucidating the role of randomness in the time behavior of chemically reacting lipid compartments, such as micelles, vesicles or micro emulsions, in regimes where random fluctuations due to the stochastic nature of reacting events can lead an open system towards unexpected time evolutions.</p> <p>Results</p> <p>This paper analyses the so-called Ribocell (RNA-based cell) model. It consists in a hypothetical minimal cell based on a self-replicating minimum RNA genome coupled with a self-reproducing lipid vesicle compartment. This model assumes the existence of two ribozymes, one able to catalyze the conversion of molecular precursors into lipids and the second able to replicate RNA strands. The aim of this contribution is to explore the feasibility of this hypothetical minimal cell. By deterministic kinetic analysis, the best external conditions to observe synchronization between genome self-replication and vesicle membrane reproduction are determined, while its robustness to random fluctuations is investigated using stochastic simulations, and then discussed.</p

A spatial model of autocatalytic reactions

Biological cells with all of their surface structure and complex interior stripped away are essentially vesicles - membranes composed of lipid bilayers which form closed sacs. Vesicles are thought to be relevant as models of primitive protocells, and they could have provided the ideal environment for pre-biotic reactions to occur. In this paper, we investigate the stochastic dynamics of a set of autocatalytic reactions, within a spatially bounded domain, so as to mimic a primordial cell. The discreteness of the constituents of the autocatalytic reactions gives rise to large sustained oscillations, even when the number of constituents is quite large. These oscillations are spatio-temporal in nature, unlike those found in previous studies, which consisted only of temporal oscillations. We speculate that these oscillations may have a role in seeding membrane instabilities which lead to vesicle division. In this way synchronization could be achieved between protocell growth and the reproduction rate of the constituents (the protogenetic material) in simple protocells.Comment: Submitted to Phys. Rev.