Integrative Multiscale Model

Molecular model

Motivation

The growth of a virus in its host cell is a complex and highly orchestrated multimolecular process; this complex molecular phenomenon gains importance when you want to study the behavior of designed systems that impact intracellular scales at first stages.

Based on previous works from Drew Endy et al. 1996 [reference number] we model chemical-kinetic systems of the life cycle of phage T7 in E. coli using a stochastic framework, one of the systems corresponds to the Wild-Type Model (WTM) of the life cycle of phage T7, in order to simulate the performance of our designed kamikaze system on the life cycle of phage T7 we made a second model or the KamikaZe Model (KZM), KZM integrates contribution of total number of ribosomes to overall translation rates and attack of toxins over bacterial translation machinery.

Data published over the last 40 years on various fields about phage T7 is integrated here [reference number maybe]. These stochastic models incorporate entry of T7 genome into the host, transcription and expression of T7 genes, translation from phage polycistronic mRNAs to proteins, replication of T7 DNA, assembly of T7 procapsids, packaging of T7 DNA, degradation of mRNAs, T7 final particle assembly and in the case of the KZM, contribution of ribosomes to global translation rates, inactivation of ribosomes and decay of the global translation rates.

With a fairly big number of runs of both models we can assemble probability distributions for each molecular species as a function of time. For this [in this intracellular model] multi-scale model we are particularly interested in the burst size Distribution (BSD); burst size is the number of phages an infected cell produce.

Once we have the BSD we are ready for the Spatial Population Model and finally for the whole integrative multi-scale model of our project.

Results and achievements

We decide to model the performance of the kamikaze system using a [http://en.wikipedia.org/wiki/Stochastic_simulation stochastic simulation algorithm]. Dan Gillespie [reference number] first implemented stochastic algorithms and methods to simulate the time evolution for chemically reacting systems. One of the most important advantages of the stochastic framework over the deterministic methods is that stochastic fluctuations of the reactions can play an important role in single cell/system behavior, which is hidden when you use deterministic procedures. We assembled WTM and KZM based on existing deterministic chemical-kinetic models [referencia a Endy y a Lingchong You] and resolved them using stochastic solvers from Simbiolgy. Both models integrate various events and chemical flows that determine the development of life cycle of phage T7 without and with kamikaze system.

In WTM we simulate different phenomena of the development of phage, from entry of phage’s DNA to final particle assembly and lysis, the importance of this model is that we reproduce the wild type life cycle of phage T7 and through an ensemble of runs we build a wild type BSD.

On the other hand, KZM tested the performance of our kamikaze. Apart from all wild-type processes, KZM incorporates the contribution of the ribosomes to the translation rates; this has been done in order to simulate the ribosome decay by the ribosome-inactivating toxin from the kamikaze construction. New BSD will result from the ensemble of runs of KZM, which will be used to study the behavior of our construction at the population scale.

Population model