Team:Heidelberg/HEARTBEAT network/Results

{|
 * -valign="top" border="0" style="margin-left: 2px;"
 * width=650px" style="padding: 0 15px 15px 20px; background-color:#ede8e2"|

= HEARTBEAT Fuzzy Modeling - RESULTS =

Contents

 * Modeling inducible promoters
 * Modeling constitutively active promoters
 * Modeling p53 - a high sensitive promoter
 * Fuzzy Integrator: Final Model for Combined Activity


 * go back to HEARTBEAT fuzzy network OVERVIEW

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).

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.

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.

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.

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).

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.

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&kappa;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).

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&kappa;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&kappa;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&kappa;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&kappa;B action is required to obtain the highest promoter 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&kappa;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&kappa;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.

[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.

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&kappa;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&alpha; was held at a constant level. All promoters, including the constitutive ones, were stimulated. NF&kappa;B 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.

Stimulation with TNF&alpha; 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&alpha; as variable instead of CPT. A basal activity was exerted by the constitutive promoters as observed before. NF&kappa;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&alpha;, 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&kappa;B, and then the constitutive underwent an sharp increase. Interestingly, the slope of p53 sharply increased with a very little delay compared to NF&kappa;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.

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]