In this thesis, I propose a general stochastic simulation algorithm for particle based lattice model using the concepts of Gillespie's stochastic simulation algorithm, which was originally designed for well-stirred systems. I describe the details about this method and analyze its complexity compared with the StochSim algorithm, another simulation algorithm originally proposed to simulate stochastic lattice model. I compare the performance of both algorithms with application to two different examples: the May-Leonard model and Ziff-Gulari-Barshad model. Comparison between the simulation results from both algorithms has validate our claim that our new proposed algorithm is comparable to the StochSim in simulation accuracy. I also compare the efficiency of both algorithms using the CPU cost of each code and conclude that the new algorithm is as efficient as the StochSim in most test cases, while performing even better for certain specific cases. / Computer simulation has been developed for almost one century. Stochastic lattice model, which follows the physics concept of lattice, is defined as a kind of system in which individual entities live on grids and demonstrate certain random behaviors according to certain specific rules. It is mainly studied using computer simulations. The most widely used simulation method to for stochastic lattice systems is the StochSim algorithm, which just randomly pick an entity and then determine its behavior based on a set of specific random rules. Our goal is to develop new simulation methods so that it is more convenient to simulate and analyze stochastic lattice system. In this thesis I propose another type of simulation methods for the stochastic lattice model using totally different concepts and procedures. I developed a simulation package and applied it to two different examples using both methods, and then conducted a series of numerical experiment to compare their performance. I conclude that they are roughly equivalent and our new method performs better than the old one in certain special cases.

The reaction-diffusion master equation (RDME) is a stochastic model for spatially heterogeneous chemical systems. Stochastic models have proved to be useful for problems from molecular biology since copy numbers of participating chemical species often are small, which gives a stochastic behaviour. The RDME is a discrete space model, in contrast to spatially continuous models based on Brownian motion. In this thesis two accuracy issues of the RDME on unstructured meshes are studied. The first concerns the rates of diffusion events. Errors due to previously used rates are evaluated, and a second order accurate finite volume method, not previously used in this context, is implemented. The new discretisation improves the accuracy considerably, but unfortunately it puts constraints on the mesh, limiting its current usability. The second issue concerns the rates of bimolecular reactions. Using the macroscopic reaction coefficients these rates become too low when the spatial resolution is high. Recently, two methods to overcome this problem by calculating mesoscopic reaction rates for Cartesian meshes have been proposed. The methods are compared and evaluated, and are found to work remarkably well. Their possible extension to unstructured meshes is discussed.

Gillespie's Stochastic Simulation Algorithm (SSA) provides an elegant simulation approach for simulating models composed of coupled chemical reactions. Although this approach can be used to describe a wide variety biological, chemical, and ecological systems, often systems have external behaviors that are difficult or impossible to characterize using chemical reactions alone. This work extends the applicability of the SSA by adding mechanisms for the inclusion of external events and external triggers. We define events as changes that occur in the system at a specified time while triggers are defined as changes that occur to the system when a particular condition is fulfilled. We further extend the SSA with the efficient implementation of these model parameters. This work allows numerous systems that would have previously been impossible or impractical to model using the SSA to take advantage of this powerful simulation technique.

O comportamento de sistemas epidêmicos é frequentemente descrito de maneira determinística, através do emprego de equações diferenciais ordinárias. Este trabalho visa fornecer uma visão estocástica do problema, traçando um paralelo entre o encontro de indivíduos em uma população e o choque entre partículas de uma reação química. Através dessa abordagem é apresentado o algoritmo de Gillespie, que fornece uma forma simples de simular a evolução de um sistema epidêmico. Fundamentos de processos estocásticos são apresentados para fundamentar uma técnica para a estimação de parâmetros através de dados reais. Apresentamos ainda o modelo de Tau-leaping e o modelo difusivo elaborados através de equações diferenciais estocásticas que são aproximações do modelo proposto por Gillespie. A aplicação dos modelos apresentados é exemplificada através do estudo de dados reais da epidemia de dengue ocorrida no estado do Rio de Janeiro entre os anos de 2012 e 2013. / The behavior of epidemic systems is often described in a deterministic way, through the use of ordinary differential equations. This paper aims to provide a stochastic view of the problem, drawing a parallel between the encounter between individuals in a population and the clash between particles of a chemical reaction. Through this approach is presented the Gillespie algorithm, which provides a simple way to simulate the evolution of an epidemic system. Fundamentals of stochastic process theory are presented to support a technique for estimating parameters through real data. We present the model of Tau-leaping and the diffusive model elaborated by stochastic differential equations that are approximations of the model proposed by Gillespie. The application of the presented models is exemplified through the study of real data of the dengue epidemic occurred in the state of Rio de Janeiro between the years of 2012 and 2013.

Cette thèse a pour sujet la biologie cellulaire, et plus particulièrement l'organisation interne des cellules eucaryotes. Bien que les différents acteurs régissant cette organisation aient été en grande partie identifiées, on ignore encore comment une architecture si complexe et dynamique peut émerger de simples interactions entres molécules. Un des objectifs des différentes études présentées dans cette thèse est de construire un cadre théorique permettant d'appréhender cette auto-organisation. Pour cela, nous étudions des problèmes spécifiques à différentes échelles allant du nanomètre (dynamique des hétérogénéités dans les membranes biologiques) au micromètre (organisation des organelles cellulaires), en utilisant des simulations numériques stochastiques et des méthodes analytiques. Le texte est organisé pour présenter les résultats des plus petites au plus grandes échelles. Dans le premier chapitre, nous étudions l'organisation de la membrane d'un seul compartiment en modélisant la dynamique d'hétérogénéités membranaires. Dans le second chapitre, nous étudions la dynamique d'un compartiment unique échangeant des vésicules avec le milieu extérieur. Nous étudions également comment deux compartiments différents peuvent être générés par les mêmes mécanismes d'échanges de vésicules. Enfin, dans le troisième chapitre, nous développons un modèle global de la dynamique des organelles cellulaires, dans le contexte particulier de la biogenèse de l'appareil de Golgi. / This thesis deals with cell biology, and particularly with the internal organization of eukaryotic cells. Although many of the molecular players contributing to the intra-cellular organization have been identified, we are still far from understanding how the complex and dynamical intra-cellular architecture emerges from the self-organization of individual molecules. One of the goals of the different studies presented in this thesis is to provide a theoretical framework to understand such self-organization. We cover specific problems at different scales, ranging from membrane organization at the nanometer scale to whole organelle structure at the micron scale, using analytical work and stochastic simulation algorithms. The text is organized to present the results from the smallest to the largest scales. In the first chapter, we study the membrane organization of a single compartment by modeling the dynamics of membrane heterogeneities. In the second chapter we study the dynamics of one membrane-bound compartment exchanging vesicles with the external medium. Still in the same chapter, we investigate the mechanisms by which two different compartments can be generated by vesicular sorting. Finally in the third chapter, we develop a global model of organelle biogenesis and dynamics in the specific context of the Golgi apparatus

Toxoplasma gondii (T. gondii) is an intracellular protozoan parasite. The parasite can infect all warm-blooded vertebrates. Up to 30% of the world's human population carry a Toxoplasma infection. However, the transmission dynamics of T. gondii has not been well understood, although a lot of mathematical models have been built. In this thesis, we adopt a complex life cycle model developed by Turner et al. and extend their work to include diffusion of hosts. Most of researches focus on the deterministic models. However, some scientists have reported that deterministic models sometimes are inaccurate or even inapplicable to describe reaction-diffusion systems, such as gene expression. In this case stochastic models might have qualitatively different properties than its deterministic limit. Consequently, the transmission pathways of T. gondii and potential control mechanisms are investigated by both deterministic and stochastic model by us. A stochastic algorithm due to Gillespie, based on the chemical master equation, is introduced. A compartment-based model and a Smoluchowski equation model are described to simulate the diffusion of hosts. The parameter analyses are conducted based on the reproduction number. The analyses based on the deterministic model are verified by stochastic simulation near the thresholds of the parameters. / Master of Science

No description available.

Biochemical reactions make up most of the activity in a cell. There is inherent stochasticity in the kinetic behavior of biochemical reactions which in turn governs the fate of various cellular processes. In this work, the precision of a method for dimensionality reduction for stochastic modeling of biochemical reactions is evaluated. Further, a method of stochastic simulation of reaction kinetics is implemented in case of a specific biochemical network involved in maintenance of long-term potentiation (LTP), the basic substrate for learning and memory formation. The dimensionality reduction method diverges significantly from a full stochastic model in prediction the variance of the fluctuations. The application of the stochastic simulation method to LTP modeling was used to find qualitative dependence of stochastic fluctuations on reaction volume and model parameters. / text

The finite state projection algorithm provides modelers a new way of directly solving the chemical master equation. The algorithm utilizes the matrix exponential function, and so the algorithm’s performance suffers when it is applied to large problems. Other work has been done to reduce the size of the exponentiation through mathematical simplifications, but efficiently exponentiating a large matrix has not been explored. This work explores implementing the finite state projection algorithm on several different high-performance computing platforms as a means of efficiently calculating the matrix exponential function for large systems. This work finds that general purpose graphics processing can accelerate the finite state projection algorithm by several orders of magnitude. Specific biological models and modeling techniques are discussed as a demonstration of the algorithm implemented on a general purpose graphics processor. The results of this work show that general purpose graphics processing will be a key factor in modeling more complex biological systems.

Nous nous proposons d'étudier la gestion du bruit stochastique d'expression génique. On s'intéresse plus particulièrement à la dynamique du bruit lors de la réponse cellulaire. Comment évolue le bruit? Quels sont les mécanismes en jeux? Quelle est l'importance du bruit dans le fonctionnement cellulaire? Pour répondre à ces questions, nous nous appuyons sur le réseau de régulation génétique qui gère la réponse au stress nutritionnel chez E. Coli. L'étude du comportement dynamique de ce réseau, au niveau d'une population de bactéries, a été initiée et est portée par la forte collaboration de deux équipes de la région : une de bio-informaticiens (l'équipe de Hidde de Jong de l'INRIA Rhône-Alpes) et la deuxième de biologistes (l'équipe de Hans Geiselmann, Laboratoire d'Adaptation et Pathogénie des Micro-organismes). En profitant donc de l'expérience et de la compréhension acquise par ces équipes, nous étudions les réponses individuelles de chaque bactérie lors de la transition entre état de stress nutritionnel, et état de croissance exponentielle. Le bruit d'expression génique est quantifié dans des nœuds clés du réseau de régulation. Pour ce faire, les bactéries sont suivies individuellement par microscopie de fluorescence sur plusieurs générations. Les données de fluorescence collectées sur cellules uniques permettent d'étudier la variabilité inter-cellulaire. Cette variabilité est quantifiée tout le long de la réponse: à chaque instant, on connaît la distribution des densités de fluorescence cellulaire dans la population de cellules. Et le suivi des lignées individuelles permet de travailler sur une population de cellules saines: les individus malades ou morts qui ne se divisent pas, sont écartés. En réduisant ainsi les phénomènes cellulaires en jeux, on réduit le nombre de paramètres. Les sources de bruit sont moins nombreuses, et il est plus facile de comprendre les mécanismes en jeux. Les informations de lignage cellulaire permettent aussi d'étudier la variabilité introduite par la phase du cycle cellulaire: les événements de division cellulaire peut être artificiellement synchronisés. Le bruit est alors étudié sur une population en phase lors de la division. Cette étude montre que le bruit sondé n'est pas dominé par les différences dans la phase du cycle cellulaire. On peut donc modéliser nos cellules sans tenir compte des différences introduites par le cycle cellulaire. Le modèle testé est simplifié aux étapes de transcription-traduction-maturation. Les paramètres du modèle sont inférés de nos données expérimentales, et le modèle est testé à travers des simulations. / We aim to investigate the management of the stochastic noise in gene expression and more precisely the study of noise in dynamical cellular responses. How the noise varies following a perturbation? What mechanisms are at play? How important is noise in the cellular function? To answer these questions, we are interested in the genetic regulatory network that handles the nutritional stress response in E. Coli. The noise of gene expression is quantified in a key node of the network control. For that bacteria are followed individually by fluorescence and phase contrast microscopy over several generations. This variability between cells is quantified throughout the response to the nutritional perturbation: at every moment, we know the density distribution of cellular fluorescence in the cell population. And monitoring of individual lines allows us to take into account only the population of healthy cells: individuals that do not divide neither grow, are discarded. Thereby reducing other sources of variability (e.g. cellular phenomena) we reduce the number of parameters. Noise sources are less numerous, and it is easier to understand the mechanisms at play. Also the information on cell lineage allow to study the variability introduced by the phase of the cell cycle: the events of cell division can be artificially synchronized. This study shows that the noise sounded is not dominated by differences in the phase of the cell cycle. We can therefore model our cells regardless of the differences introduced by the cell cycle. The tested model is simplified to the steps of transcription-translation-maturation. The model parameters are inferred from our experimental data and the model is tested through simulations.