Team:Heidelberg/HEARTBEAT network/Results

From 2009.igem.org

Revision as of 03:52, 22 October 2009 by Naoiwamoto (Talk | contribs)

HEARTBEAT Fuzzy Modeling - RESULTS

Contents


Modeling inducible promoters

Capturing the basic behaviour of an inducible promoter

Based on prior data analysis of data derived from flow cytometry experiments using inducible synthetic NFkB responsive promoters we developed a fuzzy controller which reflects the behavior of inducible promoters.

By using our model we were able to capture the basic behaviour of an inducible promoter (Fig.1).

Figure 1: Activity of an inducible promoter. Promoter of interest is only active if it is inducible and a stimulation is present.


Capturing the impact of one parameter

Once this was accomplished, we aimed for capturing the dependency of the activity of the promoter on more heterogeneous parameters, such as affinity score, position score and presence stimulation (Figs. 2a and 2b). As expected, promoter activity was observed only in the presence of a defined stimulation. Obviously, the affinity score produces a more dynamic, steep gradient in a certain range whereas the position score changes the promoter activity only slightly.

Figure 2a: Activity of an inducible promoter depending on the position score (PosScore)
Figure 2b: Activity of an inducible promoter depending on the affinity score (AfScore)


The combinatorial dependency of promoter activity upon stimulation

Finally, we tested the combinatorial dependency of promoter activity on both position and affinity for an inducible promoter upon stimulation (Fig. 3). According to our result, an optimal promoter activity is obtained only in case the promoter sequence harbours a high position as well as a high affinity score. As obtained above, a high affinity score raises the promoter activity to a certain level (and therefore can be considered as a more crucial parameter). A shift towards a high position score, however, is still needed to reach an maximal promoter activation.

Figure 3: Affinity Score is responsible for the switch to a high promoter activity level. A high position score is required though to reach the maximal activity.

Our final model of inducible promoters shows the relationship between stimuli, type of promoter, i.e. inducible, affinity score and position score. The model predicted a high dependency of the promoter activity on affinity score revealing its crucial role in TF binding and transcriptional activation. In addition, a certain level of the position score is required for the promoter sequence to augment its promoter activity and reach the maximum, thereby becoming an optimal inducible promoter

Importantly, our model captures the behaviour of an inducible promoter. We next targeted another type of promoters. We wondered whether we would be able to establish the same functional model for an non-inducible promoter (constitutively active, short: constitutive).

[TOP]

Modeling constitutively active promoters

In parallel to the rationally designed synthetic promoters our team created a set of randomly assembled synthetic promoters which are constitutively active by using RA-PCR. Based on data analysis of flow cytometry experiments we developed another fuzzy controller which reflects the behaviour of constitutively active promoters.

Capturing the basic behaviour of an constitutive promoter

By using our model we were able to capture the basic behaviour of an inducible promoter (Fig.4). According to our model, constitutive promoter is only active if it's type is selected to be "not inducible". We showed that the constitutive promoter acted independently from the presence of an stimulus.

Figure 4: Activity of an constitutive promoter


Capturing the impact of one parameter

By accomplishing this, we again aimed for capturing the effect of single parameters on the "constitutive system"(Figs. 5a and 5b). A gradient of transcriptional activity was observed when varying the position score over the whole range (Fig. 5a). No difference was shown between a invoked system (stimulus = "yes") and a non-invoked system (stimulus = "no"). The affinity score showed the same behaviour (data not shown). We also showed the effect of the sequence affinity on the outcome of the model (Fig. 5b). Interestingly, the system showed a three-step behaviour decribing a low (basal activity), an as well as a high activity state. Hereby we also took the promoter type into account. As expected, the promoter was only active if the promoter type was set to "constitutive". The same behaviour was documented for the position score (data not shown).

Figure 5a: Activity of an constitutively active promoter depending on the position score (PosScore)
Figure 5b: Activity of an constitutively active promoter depending on the affinity score (AfScore)


Effect of TFBS position and TF affinity on promoter activity

Fig. 3 shows the combinatorial dependency of constitutive promoter activity on both position and affinity. In the range of a high affinity score, the system behaves similar to our system describing inducible promoters. In general we can state that these results support our hypothesis that the affinity score has to be considered as a crucial parameter besides the position of a TFBS on a certain promoter sequence. Still, both scores hast to be in a high range in order to reach maximal promoter activation.

Figure 6: Both high affinity score as well as position score are required to reach maximal transcriptional activity.

Analogously to our previous model of inducible promoters, in here we were also able to capture the relationship between stimuli, affinity score and position score for an constitutively active promoter. It was able to reproduce the behaviour of the sequences derived by RA-PCR. It showed a three-step-dependency of the transcriptional activity on the affinity score which was in good accordance to prior observations.

Having accomplished the development of two primary fuzzy models, our next task was to integrate a highly senstive, inducible promoter (p53) which comprises not only the basic behaviour but also reflects more complex behaviour. Different stimulant concentration levels and a further pathway will be considered as auxiliary inputs.


[TOP]

Modeling p53 - a high sensitive promoter

It has be previously described that the crosstalk between the tumour supressor p53 and NFκB pathways regulate apoptosis [1]. As a further challange, we aimed at integrating this complex behaviour into our HEARTBEAT fuzzy network model. Herefore we used HEARTBEAT to evaluate eight p53 responsive sequences which were generated using our newly postulated RA-PCR. Subsequently calculating their respective respective positions and affinity scores, we would like to measure the activity of these sequence upon induction with CPT (camptothetin), as CPT has been established to promote apoptosis [2]. Thus, our model was designed to reflect the relationship between NFkB activity, increasing concentration of CPT as well as affinity and position score of the eight p53 responsive sequences to the respective inducible promoter activity.

Basal performance of the p53 promoter

Our model which has p53-responsive promters as a integrative system was able to capture the basic behaviour of an inducible promoter similar to Figs. 1-3. It was only active if it was invoked with a certain stimulation and each single parameter hat an enhancing impact on the transcriptional activity (data not shown).

Advanced performance of the p53 promoter

Differential stimulation resulted in several activity levels

We first characterized the impact of newly introduced parameter, the stimulation by camptothecin (CPT). Figs. 8a and 8b show the outcome of our model for variable CPT concentration depending on the position score as well as on the affinity score, respectively. As shown in Fig. 8, only a high CPT stimulation led to an enhanced promoter activity. Thereby the position score had an decisive role in jumping to a highly active state. A similar result was obtained when level of CPT stimulation was compared to the affinitiy score (data not shown).

The combinatorial effect of CPT and NFκB

Following the parameter analysis of variable stimulation levels, we characterized the dependency of the promoter activity on the second newly introduced parameter. The activity of NFκB-pathway was shown to activate the p53 [2] . A synergistic effect of these two gene regulatory pathways was reported by our model which succeeded in capturing of that behaviour as well. Increasing NFκB activity corresponds to a higher transcriptional activity when CPT stimulation is high. Both high stimulation by CPT as well as a high mode of NFκB action is required to obtain the highest promoter activity.

Figure 8: Differential CPT stimulation leads to a differential promoter activity. The system is risen up to a high level when a high position score is invoked by high CPT concentration.
Figure 9: A combined high induction is by CPT and NFκB is lifting the transcriptional activity..

An maximal set of parameter inputs is essential for optimal promoter activation

So far we our model characterized the behaviour of an inducible promoter depending on a single or two parameter/s. The overall combinatorial dependency of the p53 promoter system on all five inputs was determined (Fig. 10). In this surface diagram showing the activity among differing position and affinity scores all other model inputs were set to their maximum. Regarding this specific conditions, our model predicts that the system is able to reach its maximum only if all five parameters are reaching their maximum. This supports our hypothesis that the combinatorial inducement of the system by CPT and NFκB reveals a synergistic effect on the system. All in all, our p53 model was able to reproduce the performance of an inducible promoter (see inducible NFκB model) and furthermore demonstrated the importance of a differential excitation of a signal integrity circuit.


We successfully developed three individual fuzzy inference systems which captured the behaviour of different promoter types. In the next section, all three systems are combined to one large fuzzy network system.

Figure 10: All parameters need to be at their maximum to reach optimal transcriptional activity. For this surface plot on a p53 responsive, inducible promoter, NFκB and CPT were set to 4 and 1, respectively (maximal input).

[TOP]

Fuzzy Integrator: Final Model for Combined Activity

Can all these three single systems for evaluating promoter activity be combined and fuzzyfied together? For answering this important question we developed a final model which was capable of characterizing several possible input combinations. This global structure was then used for understanding the performance of combinatorial promoter activities.

Figure 11: Snapshot of our model integrating all three fuzzy inference systems. All possible inputs are integrated, evaluated and passed to our final integrator which then calculates the combinatorial activity of the model as a whole.

A model predicting the behaviour of a system with inducible promoters

We first focused our interest of investigation on a system with inducible promoters. All parameters except for the affinity score was set to the maximum so that we could study the whole activity of the model. Fig. 12 shows the behaviour of such a system. Both NFκB and p53 are stimulated and added up to the combinatorial activity. As expected, the constitutive promoters show only basal activity, not being affected by the changing affinity score at all. Also, the two inducible promoters behaved according to our expectation. The combined activity reached its maximum at around 6.5 Promoter Activity Units and reached its saturation.

Additional invoking of constitutive promoters show increased joined activity

The next step was to simulate the mode of action of a system combining inducible and constitutive promoters with a varying CPT input while TNFα was held at a constant level. All promoters, including the constitutive ones, were stimulated. NF&kappaB; does not show any changes in its activity over the whole range since it does not respond to capthotecin at all (Fig. 13). By simulating these conditions we accurately captured that the maximal activity is reached upon increasing induction by CPT, and that the dynamics of the combinatorial promoter activity is determined only by the dynamics of the p53 activity.

Figure 12: A model predicting the behaviour of a system with inducible promoters.
Figure 13: Additional invoking of constitutive promoters show increased joined activity.

Stimulation with TNFα unravels positive feedback on combined activity

In addition to our prior analysis we examined the same situation as in the model above except for we chose TNFα as variable instead of CPT. A basal activity was exerted by the constitutive promoters as observed before. NFκB is induced as expected and switches to a high activity level when the concentration of the stimulant reaches its high state (Fig. 14). Of particular interest is the dynamics of the p53 responsive promoters. As though not directly influenced by TNFα, p53 doubles its activity upon higher stimulation. This additional switch-like behaviour contributes to a dynamic jump of the integrative module by exhibiting a two-fold increase of the combined activity.

A final combined model uncovers global maximal behaviour

Another case in which all sequence types are considered to be evoked was studied at last. Fig. 15 shows the combined promoter activity depending on varying affinity scores. All three modules excessed only a low action when the activity is low. With increasing activity, first p53, then NFκB, and then the constitutive underwent an sharp increase. Interestingly, the slope of p53 sharply increased with a very little delay compared to NFκB which is another observation supporting the positive feedback between these two pathways. As the inducible systems reaches their saturation, the constitutive increases. This corresponds to the global peak at Affinity Score = 1 which is highly determined by the peak of the constitutive promoters. The overall decrease beyond reaching the global maximum is in good accordance with the decrease of the constitutive modules.

Figure 14: Stimulation with TNFα unravels positive feedback on combined activity.
Figure 15: A final combined model uncovers global maximal behaviour. The promoter activity is plotted against the affinity score.


All in all, we were able to create a dynamic model by joining our three single fuzzy logic submodules together to one whole global system. We showed several interesting characteristics of the system including a pathway crosstalk and a positive feedback. We thus conclude that our model is capable of reflecting not only basic gene regulatory modules but also complex behaviours which integrate a variety of gene regulatory elements.


[TOP]

References

[1] Fujioka, S. et al. Stabilization of p53 is a novel mechanism for proapoptotic function of NF-κB. The Journal of Biological Chemistry, 279, 27549-59 (2004).

[2] Lee, S.-J. et al. Combined application of camptohecin and the guanylate cyclase activator YC-1: Impact on cell death and apoptosis-related proteins in ovarian carcinoma cell lines. Chem.-Biol. Interactions, 181, 185-192 (2009).


[TOP]