Team:ESBS-Strasbourg/Modeling
From 2008.igem.org
Modeling
From the very beginning of the development of our construct, modeling played a critical role to help us identifying necessary characteristics for a working system. We used Simulink that is very intuitive to handle (thanks to the Uni Karlsruhe for the support). The model never was sophisticated enough to reflect all aspects of a living system, but nevertheless should be sufficient to help finding rational solutions for basic problems. In the following, there is a short presentation of problems we solved thanks to information received from simulation assays. The tasks are presented in a chronological way, despite the fact that it was impossible to just concentrate on one problem isolated from other factors. The real process therefore was iterative for each question to deliver the best possible result in the context of the overall system. There are two main problems that are tightly linked to one another, but that shall be presented in a separate manner for the sake of clarity:
Given the basic design with self-activator (sA), repressor (R) and cross-activator (cA), there are always two transcription factor species regulating one promoter. The ratio of activator to repressor shall be critical for either activation or repression of the transcriptional activity of the module. It would be ideal to receive a switch-like behaviour with the system being either silent or showing maximal expression with a steep transition between the two states. This can be achieved by taking advantage of cooperativity in eucaryotic expression systems in general with the additional benefits of several operator repetitions.
Two different approaches came up by either using just one DNA-binding domain for activator and repressor acting on the same operator or using two distinct DNA-binding domains acting on separate operators. In the following we refer to these two approaches as One- and Two-operator system respectively.
Characteristics with advantages and drawbacks:
One-operator system:
- activator and repressor bind to the same operator -> same DNA-binding domains; simplified promoter construction
- passive repression by competition for operator binding site
- active repression by repression domain
- problem of heterodimer formation as dimerisation domain is part of DNA-binding domain
- -> competitive inhibition equation for passive repression effect, active repression effect and heterodimer influence primarily neglected:
Two-operator system:
- Activator and repressor bind to distinct operators -> different DNA-binding domains
- active repression by repression domain, but no passive effect
- system repressed at high repressor concentration independently of activator/repressor ratio
- only what is activated can also be repressed
- -> activation equation as maximal expression rate within the repression equation:
The One-operator system is the simpler construction and also delivers better results with a broader range of TF concentrations for a successfully switching system. Nevertheless, the equation for the one-operator system needs further adaptations concerning the active repression effect and the heterodimer influence to describe the gene expression. It is unclear if the repressing effect of the repressor domain is stronger than the activating effect of the activator domain. There are few examples in literature where the two domains act in different concentrations on the same promoter. Furthermore, the impact of the heterodimers remains unclear as there might be preferences in homodimer formation. The position of the attached activator or repressor or homodimer in regard to the TATA-box might also play a critical role. Simplifying these problems, we considered an equation that is derived from the specific activation equation of allosteric enzymes as possible compromise:
Problem N° 2: required differences in activator/repressor impact
Expression of R and cA shall result in repression of the competing module, that means the repressor needs to be stronger than the cross-activator. This problem can be solved in several ways or through combinations of them:
- reducing the affinity of the DNA-binding domain of cA compared to R
- limiting the expression amount of cA compared to R
- taking advantage of the localisation effects with cell-cycle dependent NLS
Depending on the choice of expression equation, the effect of reducing the affinity of the DNA-binding domain of cA shows the same or similar result as a lower expression rate. The lower expression rate bears the advantage of a better standardisation between the two modules and also for further bits. Affinity mutations are different for each DNA-binding domain whereas differences in expression amount should always deliver the same concentration ratios. There are again two different ways to obtain different expression levels: one could either express the two genes from one single promoter with an IRES or from two different promoters comprising the same operator motif but different TATA box mutations. Taking advantage of cell-cycle dependent NLS would also be a very elegant methods, but is complicated to implement and also shows a narrow regulation window in several simulation assays.