Team:Edinburgh/Modelling/Genomic

From 2010.igem.org

Revision as of 08:58, 18 August 2010 by JRWK (Talk | contribs)







Overview: The Kappa modelling language


In order to model the core repressilator system and its attached signal transduction pathways, we have made use of a stochastic, agent- and rule-based language called Kappa.





In Kappa, biological entities such as proteins, DNA, and RNA are represented as agents, which are essentially named sets of sites that can be used to hold state or bind and interact with other agents. The example shown shows how a promoter BioBrick can be represented within Kappa – as a three-agent-long piece of DNA, connected via upstream and downstream sites, with binding sites for transcription factors and RNA polymerase, and a type site to keep track of its Registry code.





Interactions are represented by rules in the form of precondition and effect, with an associated rate of reaction that governs how frequently the interaction occurs. The example shown describes a repressor binding to an open binding site upon a promoter; note the preconditions that both promoter and repressor binding sites must be empty beforehand, and the effect that they are now bound together. In this case, the reaction is reversible; that is, there are both forward and backward reaction rates associated with binding and dissociation of the repressor upon the promoter. By combining agents with an appropriate set of rules and rates, a Kappa model can be used to simulate systems of varying complexity, from a simple MAPK cascade to the oscillating rhythm of a circadian clock.







The results of the simulation shown track the discrete counts of lambda-cI (red), TetR (blue), and LacI (yellow) in a slightly-modified version of the Elowitz repressilator, for one particular stochastic trajectory (obviously, different runs of the simulation will generate different results, some more variable than others). The modification made was the addition of a red luciferase BioBrick linked to a lacI promoter; high amounts of lacI repress the production of the red luciferase, but as soon as the concentration of lacI falls, the amount of red luciferase in the system rises, as expected.

Agent- and rule-based modelling allows for the circumvention of one of the primary problems with biological modelling: the fact that molecular entities existing under different conditions (phosphorylation, states of activation, etc.) can result in quantitative combinatorial explosion that greatly complicates traditional modelling methods such as differential equations. Kappa models can be unambiguously reduced to their ODE counterparts, but to do so would likely result in a system too large and too complicated to understand or create from scratch, as each rule in Kappa would equate to a potentially massive number of reactions in an ODE model.


For more information regarding Kappa and the basics of the language, please see the following resources: