"
Team:Groningen/modeling simulation engine
From 2011.igem.org
The Cumulus simulation engine
Cumulus uses an internal single cell simulation engine to evaluate parameter settings. The simulation engine is a simple temporal integration over a collection of differential equations. As an input it takes a model of the cell to simulate, some starting conditions, and the input data relating to controlled environmental factors. The timescale we use is seconds since this seems a good middle ground that allows us to characterize the processes we are interested in typically the effects of reporters such as GFP. All faster processes such as diffusion of compound through the cell or the binding of induction factors to operating regions are assumed to occur instantaneously between time steps. Slower processes such as transcription and translation delays are explicitly modeled. The reason we do not allow the use of variable timescales is that these can yield numerically different solutions for the same data. This would break the functionality of Cumulus to compare between models. Hence we decided to maintain consistency within the context of our modeling tool.
The model
For our model we use a deterministic set of differential equations (ODE's and PDE's). All these act on the current model state to compute the next state of the model. This state represents all the important characteristics in a cell which we refer to as factors. However, if need be one could access previous states to facilitate complex temporal progressions such as delays emerging from transcriptions and translations. Also we can save this entire history for any factor so examining in retrospect is possible.
We distinguish four types of factors:
- Promoter activation: Measured in polymerase bindings per second, for instance the activation of a LacI promoter
- Transcript concentrations: Measured in nanomole per liter, for instance the concentration of an RNA molecule coding for cI.
- Compound concentrations: Measured in nanomole per liter, for instance the concentration of an GFP molecule.
- Environmental factors: Measured in different units, for instance temperature or light exposure.
Using these four types of factors we can generate new models states. By saving every model state we can plot the behavior of our system in an easy and readable way.
Modularety
The Cumulus simulation engine is a system is of modular object oriened design. Because of this adding new components can be done in a matter of minutes. We only implemented the components we needed for our own system and and some test construcs, but the potential for extendabilety is very large. At the moment the simulators biggest contraines are a lack of noise modeling and multy cellular support.
- ModelComponent: The abstract base class of all model components. All ModelComponents have a name.
- Compound: The abstract base class of all the stuff floating around inside the cell. Compounts all have a concentration which is registered as nano moles for every timestep.
- Protein: The class we used for everything that is a compound but not a transcript. Proteins have a concentration decay rate over time, which covers degredation as well as the dilution because of cell growth. We used this class not only for proteins but also for smaller molecules that we assumed shared this behavoir (such as arabinose or IPTG).
- Transcript: The abstract base class for all transcripts. We treat everything that is a direct result of a transcription process. Every Transcript has a concentration decay rate. It has several subfunctionsclasses that express the different types of transcript we modeled.
- Regular Transcript: This is a transcript that moves on to translation. It has a Ribosome binding site strength and a halflife.
- taRNA Transcript: A small transcript that binds to a crRNA transcript top release its ribosome bindingsite. It has halflife but no ribosome binding afinity.
- crRNA Transcript: The crRNA is made by wrapping an other contruct in this modeling component and adding the refrence to a taRNA Transcript. The concetration of the wrapped contruct is increased by the crRNA binding to the taRNA transcript. The binding reducing both its own concentration and that of the taRNA in equal amounts and increasing the concentration of the internal transcript by that same amount. For this reaction we asume a mass action law.
- Gene: Genes are simple components that code for a single protein. They Exist as a
- Environmental: This abstract base class is the superclass of all modeling componentes that can not be modified by the cells. They typicaly contain input and output components.
- Luminescence: An environmental for converting fluorescent molecule concetrations into measurement data. It is initialised with a background an a multiplier but these are not parameters since they are not part of the biological functioning.
- Temperature: An environmental for making temprature input data, it contains function for making coldshock time series.
- Extracellular Concentration: An environmental for extracellular concentration such as arbinose or IPTG. We assume the slowly diffuse into
- Promoter: The abstract base class of all promoters it has an activation but nothing else.
- Additive Hill Promoter: The model of a promoter that can be induced via hill kinetics. It has four
- Jammed Promoter: Our way of modeling the jammers is creating a meta promoter that takes two other promoters as its input. The activety of the meta promoter is calculated by substracting the activety of the reverse promoter multiplied by some factor from the activety of the forward promoter. The multiplication factor involved is a paramete of the Jammed Promoter
