Limits...
Fuzzy logic analysis of kinase pathway crosstalk in TNF/EGF/insulin-induced signaling.

Aldridge BB, Saez-Rodriguez J, Muhlich JL, Sorger PK, Lauffenburger DA - PLoS Comput. Biol. (2009)

Bottom Line: Simulations based on fuzzy logic recapitulate most features of the data and generate several predictions involving pathway crosstalk and regulation.We also find unexpected inhibition of IKK following EGF treatment, possibly due to down-regulation of autocrine signaling.More generally, fuzzy logic models are flexible, able to incorporate qualitative and noisy data, and powerful enough to produce quantitative predictions and new biological insights about the operation of signaling networks.

View Article: PubMed Central - PubMed

Affiliation: Center for Cell Decision Processes, Cambridge, Massachusetts, United States of America.

ABSTRACT
When modeling cell signaling networks, a balance must be struck between mechanistic detail and ease of interpretation. In this paper we apply a fuzzy logic framework to the analysis of a large, systematic dataset describing the dynamics of cell signaling downstream of TNF, EGF, and insulin receptors in human colon carcinoma cells. Simulations based on fuzzy logic recapitulate most features of the data and generate several predictions involving pathway crosstalk and regulation. We uncover a relationship between MK2 and ERK pathways that might account for the previously identified pro-survival influence of MK2. We also find unexpected inhibition of IKK following EGF treatment, possibly due to down-regulation of autocrine signaling. More generally, fuzzy logic models are flexible, able to incorporate qualitative and noisy data, and powerful enough to produce quantitative predictions and new biological insights about the operation of signaling networks.

Show MeSH

Related in: MedlinePlus

Fitting MK2 rule weights.(A) A heatmap depicts the data, untrained model (Figure 3B), and trained model timecourses for MK2. (B) The regressed rule weights are plotted for the 12candidate rules. The rules are indicated in tabular format; the firsttwo rows describe the state of the inputs, TNF and time, and the lastrow is the output MK2 state. L and H represent low and high states, andE is the state describing the early response lag. Symbols above the plotshow whether the rules were present (✓) or not applicable () in the untrained model.
© Copyright Policy
Related In: Results  -  Collection

License
getmorefigures.php?uid=PMC2663056&req=5

pcbi-1000340-g007: Fitting MK2 rule weights.(A) A heatmap depicts the data, untrained model (Figure 3B), and trained model timecourses for MK2. (B) The regressed rule weights are plotted for the 12candidate rules. The rules are indicated in tabular format; the firsttwo rows describe the state of the inputs, TNF and time, and the lastrow is the output MK2 state. L and H represent low and high states, andE is the state describing the early response lag. Symbols above the plotshow whether the rules were present (✓) or not applicable () in the untrained model.

Mentions: In the work described above, logic rules and membership functions for each gatewere established manually. A better approach is to use training to optimize theweights of all possible rules in a gate by minimizing the sum of the squareddifferences between the experimental data and local model output (see Methods). Following optimization, logic rulesthat are supported by the data should have weights near 1, whilepoorly-supported rules should have weights near 0. We tested the fittingalgorithm on the MK2 gate. For such a gate, which has two MFs each for the twoinputs (TNF and time) and the output (MK2 activity),23 = 8 explicit rules are possible.MK2 data from the 10 cytokine treatment conditions were used to optimize avector containing the 8 rule weights. Our initial optimization attempt failedbecause time-dependent MFs were not parameterized so as to capture rapidincreases in signals following cytokine treatment. We had implicitly ignoredthis discrepancy when fitting the model by hand. To improve the automatedfitting procedure, an additional MF for time was included to representimmediate-early responses, increasing the number of candidate rules to 12.Optimization yielded a gate with a good fit to data using only six rules withweights near one (Figure7A). These six rules were identical to those assembled manually with theexception of the new rule needed to represent immediate early signaling (Figure 7B). To test FL gateregression with more rules, we applied the algorithm to the same MK2 data usingone additional membership function (for medium activity levels) and compared itto an untrained model using the same MFs. The training process created severalrules that were nearly identical to those introduced manually as well as severalnew ones (FigureS1). The MK2 test case suggests that it is possible to optimize ruleweights as a means to fit logic rules without bias and is a first step towards amore rigorous approach to logic-based modeling.


Fuzzy logic analysis of kinase pathway crosstalk in TNF/EGF/insulin-induced signaling.

Aldridge BB, Saez-Rodriguez J, Muhlich JL, Sorger PK, Lauffenburger DA - PLoS Comput. Biol. (2009)

Fitting MK2 rule weights.(A) A heatmap depicts the data, untrained model (Figure 3B), and trained model timecourses for MK2. (B) The regressed rule weights are plotted for the 12candidate rules. The rules are indicated in tabular format; the firsttwo rows describe the state of the inputs, TNF and time, and the lastrow is the output MK2 state. L and H represent low and high states, andE is the state describing the early response lag. Symbols above the plotshow whether the rules were present (✓) or not applicable () in the untrained model.
© Copyright Policy
Related In: Results  -  Collection

License
Show All Figures
getmorefigures.php?uid=PMC2663056&req=5

pcbi-1000340-g007: Fitting MK2 rule weights.(A) A heatmap depicts the data, untrained model (Figure 3B), and trained model timecourses for MK2. (B) The regressed rule weights are plotted for the 12candidate rules. The rules are indicated in tabular format; the firsttwo rows describe the state of the inputs, TNF and time, and the lastrow is the output MK2 state. L and H represent low and high states, andE is the state describing the early response lag. Symbols above the plotshow whether the rules were present (✓) or not applicable () in the untrained model.
Mentions: In the work described above, logic rules and membership functions for each gatewere established manually. A better approach is to use training to optimize theweights of all possible rules in a gate by minimizing the sum of the squareddifferences between the experimental data and local model output (see Methods). Following optimization, logic rulesthat are supported by the data should have weights near 1, whilepoorly-supported rules should have weights near 0. We tested the fittingalgorithm on the MK2 gate. For such a gate, which has two MFs each for the twoinputs (TNF and time) and the output (MK2 activity),23 = 8 explicit rules are possible.MK2 data from the 10 cytokine treatment conditions were used to optimize avector containing the 8 rule weights. Our initial optimization attempt failedbecause time-dependent MFs were not parameterized so as to capture rapidincreases in signals following cytokine treatment. We had implicitly ignoredthis discrepancy when fitting the model by hand. To improve the automatedfitting procedure, an additional MF for time was included to representimmediate-early responses, increasing the number of candidate rules to 12.Optimization yielded a gate with a good fit to data using only six rules withweights near one (Figure7A). These six rules were identical to those assembled manually with theexception of the new rule needed to represent immediate early signaling (Figure 7B). To test FL gateregression with more rules, we applied the algorithm to the same MK2 data usingone additional membership function (for medium activity levels) and compared itto an untrained model using the same MFs. The training process created severalrules that were nearly identical to those introduced manually as well as severalnew ones (FigureS1). The MK2 test case suggests that it is possible to optimize ruleweights as a means to fit logic rules without bias and is a first step towards amore rigorous approach to logic-based modeling.

Bottom Line: Simulations based on fuzzy logic recapitulate most features of the data and generate several predictions involving pathway crosstalk and regulation.We also find unexpected inhibition of IKK following EGF treatment, possibly due to down-regulation of autocrine signaling.More generally, fuzzy logic models are flexible, able to incorporate qualitative and noisy data, and powerful enough to produce quantitative predictions and new biological insights about the operation of signaling networks.

View Article: PubMed Central - PubMed

Affiliation: Center for Cell Decision Processes, Cambridge, Massachusetts, United States of America.

ABSTRACT
When modeling cell signaling networks, a balance must be struck between mechanistic detail and ease of interpretation. In this paper we apply a fuzzy logic framework to the analysis of a large, systematic dataset describing the dynamics of cell signaling downstream of TNF, EGF, and insulin receptors in human colon carcinoma cells. Simulations based on fuzzy logic recapitulate most features of the data and generate several predictions involving pathway crosstalk and regulation. We uncover a relationship between MK2 and ERK pathways that might account for the previously identified pro-survival influence of MK2. We also find unexpected inhibition of IKK following EGF treatment, possibly due to down-regulation of autocrine signaling. More generally, fuzzy logic models are flexible, able to incorporate qualitative and noisy data, and powerful enough to produce quantitative predictions and new biological insights about the operation of signaling networks.

Show MeSH
Related in: MedlinePlus