Dataset Information

Simulations of stochastic biological phenomena.

ABSTRACT: This Teaching Resource provides lecture notes, slides, and a student assignment for a two-part lecture that introduces stochastic modeling of biological systems. The first lecture uses biological examples to present the concept of cell-to-cell variability and makes the connection between the variability of single-cell measurements and concepts from statistical mechanics and probability theory. This section makes the point that for low copy number of a species, the usual differential equation formalism is no longer applicable and needs to be replaced by a probabilistic approach based on the so-called Master Equation. As an example, a simple model of gene transcription is discussed in detail, the different contributions to the relevant Master Equation are highlighted, and the equation itself is derived. The second lecture describes how, for more complex and biologically interesting applications, direct solution of the Master Equation becomes difficult. Gillespie's algorithm, which is used in most cases of biological interest, is then introduced as a practical alternative. The lecture delves into the crux of Gillespie's algorithm, which entails the drawing of two random numbers at each time step. It establishes the corresponding formalism, details the connection between chemical rate constants and Gillespie rates, and culminates in a description and explanation of a core MATLAB program for the transcriptional model considered in the first lecture. This core program, written for a single cell, is expanded by the students in the homework assignment to consider both transcription and translation.

SUBMITTER: Hayot F

PROVIDER: S-EPMC3225182 | biostudies-literature | 2011 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Simulations of stochastic biological phenomena.

Hayot Fernand F

Science signaling 20110920 192

This Teaching Resource provides lecture notes, slides, and a student assignment for a two-part lecture that introduces stochastic modeling of biological systems. The first lecture uses biological examples to present the concept of cell-to-cell variability and makes the connection between the variability of single-cell measurements and concepts from statistical mechanics and probability theory. This section makes the point that for low copy number of a species, the usual differential equation for ...[more]

PMID: 21954294

Similar Datasets

Project description:BackgroundAnti-tumor therapies aim at reducing to zero the number of tumor cells in a host within their end or, at least, aim at leaving the patient with a sufficiently small number of tumor cells so that the residual tumor can be eradicated by the immune system. Besides severe side-effects, a key problem of such therapies is finding a suitable scheduling of their administration to the patients. In this paper we study the effect of varying therapy-related parameters on the final outcome of the interplay between a tumor and the immune system.ResultsThis work generalizes our previous study on hybrid models of such an interplay where interleukins are modeled as a continuous variable, and the tumor and the immune system as a discrete-state continuous-time stochastic process. The hybrid model we use is obtained by modifying the corresponding deterministic model, originally proposed by Kirschner and Panetta. We consider Adoptive Cellular Immunotherapies and Interleukin-based therapies, as well as their combination. By asymptotic and transitory analyses of the corresponding deterministic model we find conditions guaranteeing tumor eradication, and we tune the parameters of the hybrid model accordingly. We then perform stochastic simulations of the hybrid model under various therapeutic settings: constant, piece-wise constant or impulsive infusion and daily or weekly delivery schedules.ConclusionsResults suggest that, in some cases, the delivery schedule may deeply impact on the therapy-induced tumor eradication time. Indeed, our model suggests that Interleukin-based therapies may not be effective for every patient, and that the piece-wise constant is the most effective delivery to stimulate the immune-response. For Adoptive Cellular Immunotherapies a metronomic delivery seems more effective, as it happens for other anti-angiogenesis therapies and chemotherapies, and the impulsive delivery seems more effective than the piece-wise constant. The expected synergistic effects have been observed when the therapies are combined.

Project description:BackgroundOver 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".MethodsIn 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.ResultsThis 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.

Project description:BackgroundThe field of synthetic biology has greatly evolved and numerous functions can now be implemented by artificially engineered cells carrying the appropriate genetic information. However, in order for the cells to robustly perform complex or multiple tasks, co-operation between them may be necessary. Therefore, various synthetic biological systems whose functionality requires cell-cell communication are being designed. These systems, microbial consortia, are composed of engineered cells and exhibit a wide range of behaviors. These include yeast cells whose growth is dependent on one another, or bacteria that kill or rescue each other, synchronize, behave as predator-prey ecosystems or invade cancer cells.ResultsIn this paper, we study a synthetic ecosystem comprising of bacteria and yeast that communicate with and benefit from each other using small diffusible molecules. We explore the behavior of this heterogeneous microbial consortium, composed of Saccharomyces cerevisiae and Escherichia coli cells, using stochastic modeling. The stochastic model captures the relevant intra-cellular and inter-cellular interactions taking place in and between the eukaryotic and prokaryotic cells. Integration of well-characterized molecular regulatory elements into these two microbes allows for communication through quorum sensing. A gene controlling growth in yeast is induced by bacteria via chemical signals and vice versa. Interesting dynamics that are common in natural ecosystems, such as obligatory and facultative mutualism, extinction, commensalism and predator-prey like dynamics are observed. We investigate and report on the conditions under which the two species can successfully communicate and rescue each other.ConclusionsThis study explores the various behaviors exhibited by the cohabitation of engineered yeast and bacterial cells. The way that the model is built allows for studying the dynamics of any system consisting of two species communicating with one another via chemical signals. Therefore, key information acquired by our model may potentially drive the experimental design of various synthetic heterogeneous ecosystems.

Dataset Information

Simulations of stochastic biological phenomena.

Publications

Simulations of stochastic biological phenomena.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets