Front. Physiol. Frontiers in Physiology Front. Physiol. 1664-042X Frontiers Media S.A. 10.3389/fphys.2013.00021 Physiology Original Research Research Topic: From structural to molecular systems biology: experimental and computational approaches to unravel mechanisms of kinase activity regulation in cancer and neurodegeneration Modeling Wnt/β-Catenin Target Gene Expression in APC and Wnt Gradients Under Wild Type and Mutant Conditions Benary Uwe 1 † Kofahl Bente 1 † Hecht Andreas 2 Wolf Jana 1 * 1Mathematical Modelling of Cellular Processes, Max Delbrück Center for Molecular Medicine Berlin-Buch Berlin, Germany 2Institute of Molecular Medicine and Cell Research, BIOSS Centre for Biological Signalling Studies, Albert-Ludwigs-Universität Freiburg Freiburg, Germany

Edited by: Matteo Barberis, Humboldt University Berlin, Germany; Max Planck Institute for Molecular Genetics, Berlin, Germany

Reviewed by: Thomas Hofer, German Cancer Research Center, Germany; Edda Klipp, Humboldt-Universität zu Berlin, Germany

*Correspondence: Jana Wolf, Mathematical Modelling of Cellular Processes, Max Delbrück Center for Molecular Medicine Berlin-Buch, Robert-Rössle-Str. 10, 13351 Berlin, Germany. e-mail: jana.wolf@mdc-berlin.de

†Uwe Benary and Bente Kofahl have contributed equally to this work.

This article was submitted to Frontiers in Systems Biology, a specialty of Frontiers in Physiology.

 25 02 2013 2013 4 21 28 09 2012 28 01 2013 Copyright © 2013 Benary, Kofahl, Hecht and Wolf. 2013

This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.

The Wnt/β-catenin pathway is involved in the regulation of a multitude of physiological processes by controlling the differential expression of target genes. In certain tissues such as the adult liver, the Wnt/β-catenin pathway can attain different levels of activity due to gradients of Wnt ligands and/or intracellular pathway components like APC. How graded pathway activity is converted into regionally distinct patterns of Wnt/β-catenin target gene expression is largely unknown. Here, we apply a mathematical modeling approach to investigate the impact of different regulatory mechanisms on target gene expression within Wnt or APC concentration gradients. We develop a minimal model of Wnt/β-catenin signal transduction and combine it with various mechanisms of target gene regulation. In particular, the effects of activation, inhibition, and an incoherent feedforward loop (iFFL) are compared. To specify activation kinetics, we analyze experimental data that quantify the response of β-catenin/TCF reporter constructs to different Wnt concentrations, and demonstrate that the induction of these constructs occurs in a cooperative manner with Hill coefficients between 2 and 5. In summary, our study shows that the combination of specific gene regulatory mechanisms with a time-independent gradient of Wnt or APC is sufficient to generate distinct target gene expression patterns as have been experimentally observed in liver. We find that cooperative gene activation in combination with a TCF feedback can establish sharp borders of target gene expression in Wnt or APC gradients. In contrast, the iFFL renders gene expression independent of gradients of the upstream signaling components. Our subsequent analysis of carcinogenic pathway mutations reveals that their impact on gene expression is determined by the gene regulatory mechanism and the APC concentration of the cell in which the mutation occurs.

 mathematical modeling gradients gene expression Wnt/β-catenin signaling liver zonation mutation feedback Introduction

Wnt ligands are secreted signaling molecules that can generate concentration gradients across tissues (Zecca et al., 1996; Whangbo and Kenyon, 1999). In adult liver, time-independent gradients in Wnt signaling components, notably Wnt and APC (adenomatous polyposis coli), have been observed along the periportal-pericentral axis. In the periportal region of the mouse liver, high APC concentrations have been detected that decrease toward the pericentral region (Benhamouche et al., 2006; Torre et al., 2010). These changes in APC levels are paralleled by alterations in the expression profiles of Wnt/β-catenin target genes such as GS (glutamine synthetase), Lect2 (leukocyte cell-derived chemotaxin 2), or Gls2 (glutaminase 2) (Benhamouche et al., 2006; Braeuning et al., 2006). In addition to APC levels, it has been suggested that the available Wnt concentration is involved in the formation of the expression profiles (Benhamouche et al., 2006). The spatial differences in Wnt/β-catenin target gene expression have been suggested to contribute to the zonation of the liver, which may provide the basis for the liver’s ability to fulfill different metabolic functions (Jungermann and Kietzmann, 1996; Burke and Tosh, 2006; Colnot and Perret, 2011). Examples are urea formation, which is mostly observed in the periportal region, and glutamine synthesis, which occurs in the pericentral area (Jungermann and Kietzmann, 1996; Torre et al., 2010; Colnot and Perret, 2011).

The Wnt/β-catenin pathway exerts its control on gene expression mainly by regulating the concentration and the activity of β-catenin (Kimelman and Xu, 2006; Cadigan and Peifer, 2009; MacDonald et al., 2009). In resting cells, β-catenin is constantly produced and degraded via a destruction-complex-dependent mechanism. The destruction complex, which contains the scaffold proteins APC and Axin, mediates the phosphorylation and ubiquitination of β-catenin. A Wnt stimulus leads to the inhibition of the destruction complex. In consequence, less β-catenin is degraded and more protein is able to enter the nucleus. There it regulates the expression of specific target genes. To fulfill its function as a transcriptional regulator, β-catenin interacts with members of the TCF family of transcription factors which comprises T-cell factors (TCF7, TCF7L1, and TCF7L2) and lymphoid enhancer-binding factor (LEF1).

A number of Wnt/β-catenin target genes have been identified (a regularly updated list can be found at http://www.stanford.edu/group/nusselab/cgi-bin/wnt/). However, detailed regulatory mechanisms of their expression have been suggested for only a few. Experiments have shown that the expression of E-cadherin is directly inhibited by the β-catenin/TCF complex (Jamora et al., 2003). A subgroup of target genes, e.g. patterning genes, has been proposed to be regulated by an incoherent feedforward loop (iFFL) motif (Goentoro and Kirschner, 2009). Different members of the TCF family have been reported to be negatively or positively regulated by Wnt/β-catenin signaling thereby giving rise to different types of feedback loops (Roose et al., 1999; Hovanes et al., 2001; Saegusa et al., 2005; Vadlamudi et al., 2005; Driskell et al., 2007). In addition to these β-catenin/TCF-dependent regulatory mechanisms, the action of co-regulators as well as the occurrence of posttranslational modifications have been suggested as contributors to the regulation of Wnt/β-catenin target gene expression (Hecht and Kemler, 2000; Archbold et al., 2012). Some of these co-regulators may regulate the switch of TCF between its repressing and activating functions (Archbold et al., 2012).

An aberrant activity of the Wnt/β-catenin pathway is associated with various diseases, most notably with different types of cancer (Clevers, 2006; Klaus and Birchmeier, 2008). Many mutations of pathway components lead to a constitutive activation of the Wnt/β-catenin pathway by interfering with the destruction-complex-dependent degradation of β-catenin and thereby changing target gene expression. For instance, mutations in APC have been detected in 80% of colorectal cancer. Mutations in β-catenin occur with a frequency of 20–30% in hepatocellular carcinoma and up to 65% in hepatoblastomas. Hepatic cancer is also associated with mutations in Axin (Giles et al., 2003).

In recent years, several mathematical models of the Wnt/β-catenin pathway have been developed, reviewed in Kofahl and Wolf (2010). They range from detailed models which take into account the molecular interactions of the pathway’s components (Lee et al., 2003), to very condensed descriptions, which aim to reproduce the general dynamics of the signaling pathway (Mirams et al., 2010). These models have been used to investigate the interplay between the adhesive and transcriptional functions of β-catenin (van Leeuwen et al., 2007) as well as the effects of carcinogenic mutations of pathway components (Cho et al., 2006). A multitude of modeling approaches in different biological contexts has been focused on gradients in order to understand how gradients are created, maintained, and interpreted by cells (Gurdon and Bourillot, 2001; Jaeger et al., 2008; Meinhardt, 2009; Nahmad and Lander, 2011; Wolpert, 2011). However, the impact of gradients of Wnt growth factors or downstream pathway components on Wnt/β-catenin target gene expression has not yet been addressed by modeling approaches.

In our theoretical approach, we aim to understand the contribution of different gene regulatory mechanisms to mRNA profiles of Wnt/β-catenin target genes in the presence of time-independent Wnt or APC concentration gradients as observed in adult liver. To this end, we develop a minimal mathematical model for Wnt/β-catenin signal transduction and extend it by different gene regulatory mechanisms that have been proposed in the literature for Wnt/β-catenin target genes. In particular, (i) an activation mechanism, (ii) an inhibition mechanism, and (iii) an iFFL are investigated. Moreover, we consider the impact of an additional possible feedback via TCF. By simulating the effects of pathway mutations, we furthermore compare target gene mRNA concentrations under wild type and carcinogenic conditions. Our study demonstrates that the combination of specific gene regulatory mechanisms and a gradient of either Wnt or APC is sufficient to generate the distinct target gene expression patterns that have been observed in the liver.

 Results Model development Minimal model of Wnt/β-catenin signal transduction and target gene expression

Various aspects of the Wnt/β-catenin pathway have been investigated in theoretical studies (Kofahl and Wolf, 2010). However, the available models either do not focus on possible regulation of downstream target gene expression and/or miss pathway components which are of importance in our analysis. We therefore develop a minimal model for the Wnt/β-catenin signaling pathway that enables us to easily combine signal transduction with different mechanisms of downstream target gene regulation. We consider individual target genes because gene regulatory networks related to Wnt/β-catenin signaling have not yet been described in the liver. Since we focus on transcriptional rather than translational regulation, only the mRNA of the target genes is included in the minimal model. Two schemes of the minimal model with identical Wnt/β-catenin signal transduction modules, but different downstream gene regulatory motifs, are shown in Figures 1A,B.

Two reaction schemes of the minimal model illustrating different mechanisms of Wnt/β-catenin target gene regulation. The signal transduction module is shared by (A,B) and highlighted by a light gray background. It consists of the reactions 1–9 that are explained in the text. The number next to the arrow denotes the particular reaction. One-headed arrows denote reactions taking place in the indicated direction. Double-headed arrows illustrate binding reactions. Dashed arrows represent activation steps and inverted, dashed “T”s denote inhibition. The gray, dashed arrow connecting the mRNA and reaction 6 indicates the possible feedback considered in Section “Transcriptional feedback regulation via TCF family members”. The two schemes differ in their downstream target gene regulation motif. (A) Direct regulatory mechanisms: the dashed line ending with a dot between β-catenin/TCF complex and reaction 10 indicates that the β-catenin/TCF complex regulates the mRNA production either by a direct activating or repressing mechanism (see Sections “Activation of Wnt/β-catenin target genes” and “Repression of Wnt/β-catenin target genes”, respectively). (B) Incoherent feedforward loop: β-catenin/TCF complexes induce the expression of a repressor of the mRNA production, which counteracts the simultaneous direct activation of mRNA expression by β-catenin/TCF complexes in reaction 10, see Section “Wnt/β-catenin target gene regulation by an incoherent feedforward loop motif”.

The Wnt/β-catenin signal transduction module of our minimal model consists of β-catenin, TCF, APC, the β-catenin/TCF complex, the β-catenin/APC complex as well as active and inactive Dishevelled (Dsha and Dshi, respectively). β-Catenin is constantly produced (reaction 1) and degraded in an APC-independent as well as an APC-dependent manner (reactions 2 and 3, respectively). β-Catenin can reversibly form complexes with APC (reaction 4) and TCF (reaction 5). TCF is produced and degraded (reactions 6 and 7, respectively) while the total concentration of APC is assumed to remain constant within the considered time period. The APC-dependent degradation of β-catenin is inhibited by activated Dishevelled. Dishevelled is activated by a Wnt stimulus (reaction 8) and can become deactivated again (reaction 9). Like APC, Dsh obeys a conservation relation.

The gene regulatory mechanisms investigated here can be grouped into two different representations illustrated in the two model schemes (Figures 1A,B). In Figure 1A, the β-catenin/TCF complex directly regulates the transcription of target gene mRNA (reaction 10). Both direct positive and direct negative regulatory mechanisms are considered in this study. In Figure 1B, the motif of the iFFL is coupled to the signal transduction module. In comparison to the scheme in Figure 1A, an additional repressor is involved in this mechanism. Its synthesis is induced by the β-catenin/TCF complex (reaction 12) and it is independently degraded via reaction 13. The repressor inhibits the production of target gene mRNA, thereby counteracting direct mRNA activation by the β-catenin/TCF complex (reaction 10). The mRNA is degraded in all cases via reaction 11 (Figures 1A,B). In Section “Transcriptional feedback regulation via TCF family members”, we investigate the implications of an additional feedback acting in combination with the regulatory mechanisms. This feedback is indicated by a gray dashed arrow from the mRNA to the TCF production (reaction 6) in Figures 1A,B.

Overall, the signal transduction module consists of five ordinary differential equations (ODEs) and two conservation relations (see Section “Model equations and parameters” in the Appendix, Eqs A1–A16). In the model presented in Figure 1A, one additional ODE is included to account for the mRNA dynamics. In the model given in Figure 1B, two additional ODEs, one accounting for the repressor dynamics and one describing the mRNA dynamics, are coupled to the signal transduction module. The kinetic parameters of the signal transduction module are set in such a way that a strong accordance between the minimal and the quantitative model by Lee et al. (2003) (hereafter referred to as the detailed model) regarding the β-catenin steady state and response to a constant Wnt stimulus is achieved, see Section “Comparison of the minimal and the detailed model”. Most of the kinetic parameters as well as the total concentrations of APC and Dsh (APCtot and Dshtot, respectively) are adopted from the detailed model. They are given in Tables A1 and A2 of Section “Model equations and parameters” in the Appendix and are used in all simulations. The expression of target gene mRNA downstream of the signaling module is described in a qualitative way since quantitative data are lacking so far. The steady state concentration of the mRNA is therefore set to an arbitrary value of about 0.68 nM under the condition of a total APC concentration of 100 nM and the absence of Wnt. This absolute mRNA concentration can easily be rescaled to more appropriate values once specific experimental data become available.

The analysis focuses on the consequences of gene regulatory mechanisms within a Wnt or APC gradient. To cover the conditions suggested to occur in the liver, the theoretical analyses are performed for a time-independent gradient of the Wnt and the total APC concentration. Quantitative data on the range and the spatial distribution of the concentration in these gradients are missing. To include all possible conditions, we vary the Wnt and APC concentrations over a wide range.

 Comparison of the minimal and the detailed model

First, we compare the dynamics of β-catenin between the minimal model (black, solid line) and the detailed model (gray, dashed line; Figure 2). The unstimulated and the Wnt-stimulated steady state levels of β-catenin are almost identical in the minimal and the detailed model; the deviation is less than 1.5% in both cases (Figure 2A). The minimal model, moreover, reproduces the dependence of the β-catenin steady state concentration on the Wnt or total APC concentration (Figures 2B,C). In general, it holds that the β-catenin steady state concentration is higher the higher the Wnt concentration or the lower the APC concentration. Note that the β-catenin concentration range is considerably broader in the APC gradient than in the Wnt gradient (around 4–1646 nM compared to 25–166 nM β-catenin, respectively), as both components act differently on the degradation of β-catenin (reaction 3). We also observe a strong agreement between the two models in their β-catenin time courses upon pathway activation by a constant Wnt stimulus. In both models, the constant Wnt stimulus is realized by setting Wnt from 0 to 1 nM at time point t = 0 min. The β-catenin dynamics in response to a transient Wnt stimulus cannot completely be captured (not shown). Taken together, the minimal model very well reproduces the dynamical properties of β-catenin that are relevant for our analysis.

Comparison of the dynamical properties of β-catenin in the minimal and the detailed model. (A) Time courses of β-catenin in response to a constant Wnt stimulation. (B,C) Steady states of β-catenin for different levels of Wnt or APC, respectively, within corresponding concentration gradients. The simulations of the minimal (black, solid line) and the detailed model (gray, dashed line) are compared.

In the detailed model, the TCF concentration obeys a conservation relation. Since here we also address a possible feedback mechanism via the transcriptional regulation of TCF, a turnover of TCF is included in our signaling module. This modification has implications for the steady states and dynamics of the β-catenin/TCF complex. While the unstimulated steady state of the β-catenin/TCF complex differs less than 1% between the minimal and the detailed model, the complex concentration increases more strongly upon Wnt stimulation in the minimal model, resulting in a higher stimulated steady state (Figure 3A). In both models the steady states of the β-catenin/TCF complex increase with an increasing Wnt or a decreasing APC concentration (Figures 3B,C). Overall, due to the restriction of the total TCF concentration in the detailed model, the concentration range of the β-catenin/TCF complex is much smaller compared to the minimal model.

Comparison of the dynamical properties of the β-catenin/TCF complex in the minimal and the detailed model. (A) Time courses of the β-catenin/TCF complex in response to constant Wnt stimulation, and (B,C) steady states of the β-catenin/TCF complex in gradients of Wnt and APC concentration, respectively, are compared between the minimal (black, solid line) and the detailed model (gray, dashed line).

Wnt and APC are differentially integrated by the signaling module. This has consequences for the signaling readout. First, the β-catenin as well as the β-catenin/TCF complex steady state concentration changes in opposite directions along the Wnt and the APC gradient. Secondly, the signal transduction module restricts the effective concentration ranges of the gradients differently, outside of which a further increase or decrease of the concentration of Wnt or APC does not lead to a significant change of the β-catenin concentration (see Figures 2B,C). A change in APC concentration can directly modulate β-catenin degradation (reaction 3; Eq. A10). Therefore, the β-catenin steady state is only limited by the ratio of its production to its alternative degradation rate constants (Eqs A8 and A9, respectively). In contrast, a change in the Wnt concentration is not linearly passed to the β-catenin degradation (reaction 3), but has an indirect effect via the activation of Dsh. Dsh obeys a conservation relation (Eq. A6) which leads to a limitation of the Wnt concentration range that influences the β-catenin degradation.

 Analysis of target gene expression Wnt/β-catenin target gene regulation by an incoherent feedforward loop motif

It has been proposed that a subgroup of Wnt target genes, in particular patterning genes, is regulated by an iFFL motif (Goentoro and Kirschner, 2009; Goentoro et al., 2009). We investigate the impact of this motif using our model (Figure 1B). We chose the parameters of the iFFL motif in a way that the motif reproduces all the properties described in Goentoro et al. (2009); parameters are given in Table A3 of Section “Model equations and parameters” in the Appendix. These are the ability to generate a perfect adaptation of the mRNA concentration and the independence of the mRNA response to the initial absolute concentration of β-catenin/TCF (Figure A1 in the Appendix). It has been shown that the iFFL motif meets these two requirements only for parameters chosen from a certain parameter range (Goentoro et al., 2009).

Having verified that the iFFL in the minimal model reproduces these characteristics, we next analyzed the consequences of this motif on the mRNA steady state within a Wnt or an APC gradient. Our simulations show that an increase in the Wnt concentration along a gradient results in an increase of the β-catenin/TCF level (Figure 4A). In contrast to the response of β-catenin/TCF, the mRNA steady state concentration remains almost constant along the Wnt gradient (Figure 4C). Similar to our observations with the Wnt gradient, levels of β-catenin/TCF change along the APC concentration gradient (Figure 4B). Again, the iFFL motif renders the mRNA steady state concentration robust over a wide APC concentration range (Figure 4D). Overall, the iFFL motif allows for a constant target gene expression if the parameters are chosen in the way described above.

Expression levels of differentially regulated target genes in a Wnt or an APC concentration gradient. (A,B) β-Catenin/TCF steady states within the Wnt or the APC concentration gradient. (C–H) Steady state concentrations of target genes induced by an iFFL motif (C,D), induced by cooperative activation with the Hill coefficient n = 3 (E,F), and upon inhibition (G,H) are shown within (C,E,G) a Wnt gradient or (D,F,H) an APC gradient.

 Activation of Wnt/β-catenin target genes

For many Wnt target genes, an induction of their expression upon activation of the Wnt/β-catenin pathway has been reported. Examples are Axin-2 as well as the liver-specific Wnt target genes OAT (ornithine aminotransferase) and Glt1 (glutamate transporter 1) (Cadoret et al., 2002; Jho et al., 2002). In the promoters of many Wnt target genes, multiple TCF/LEF binding elements (TBEs) have been detected. For instance, eight TBEs have been identified in the promoter and first intron of Axin-2 (Jho et al., 2002), and three TBEs have been reported for the liver-specific Wnt target gene Lect2 (Ovejero et al., 2004). Multiple TBEs open the possibility of cooperative target gene activation. However, detailed information on the activation of target gene expression is limited. To gain insights into possible activation mechanisms, we analyzed published experimental data of β-catenin reporter constructs (Biechele and Moon, 2008). In these experiments, the responses to different Wnt concentrations of reporter constructs with 3, 8, and 12 TBEs have been compared. In our approach, we considered three possible kinetics for the transcription rate: (i) linear mass action kinetics, (ii) Michaelis–Menten kinetics, and (iii) Hill-type activation. We fitted the kinetics (i), (ii), and (iii) to the data and compared the calculated values of the corrected Akaike information criterion (AICc; see Section “Fitting of reporter data” in the Appendix) (Hurvich and Tsai, 1995; Lukacs et al., 2010). The analysis revealed that the data for all four constructs are best explained by Hill-type functions (Figure 5A; Figure A2 in the Appendix), with Hill coefficients n higher than 1 (Figure 5B) suggesting a cooperative regulation of these reporters. We furthermore observed that the medians of the Hill coefficients for the construct with 8 TBEs and those with 12 do not vary significantly (3.9, 4.0, and 3.9, respectively) although the number of TBEs increases by 50%. The median of the construct with three TBEs is smaller than the other three values (2.6). This analysis suggests that the expression of the β-catenin reporter constructs occurs in a cooperative manner.

Fitting of different kinetic expressions to reporter construct data. (A) The readout intensity data of the reporter construct with three TBEs is plotted with standard deviations; data provided by Biechele and Moon (2008). The respective linear activation function (light gray), Michaelis–Menten activation function (dark gray), and Hill-type activation function (black) realizing the best fit (i.e., yielding the minimal log-likelihood function value) are shown. The fits of the linear and the Michaelis–Menten activation lie on top of each other. (B) The distributions of the fitted Hill coefficients are shown for the four reporter constructs. Only those Hill coefficients of fits are taken into account, whose log-likelihood function value deviates at most 20% from the log-likelihood function value of the best fit. The constructs with 3, 8, transiently transfected 12 TBEs, and stable integrated 12 TBEs are color-coded in black, dark gray, gray, and light gray bars, respectively.

We next investigated the consequences of cooperative activation on target gene mRNA expression within a Wnt or an APC concentration gradient. To this end, we used the overall minimal model including signal transduction and gene expression, with the mRNA production rate (v10) given by the Eq. A25 (Section “Model equations and parameters” in the Appendix). A Hill coefficient of n = 3 as derived from the analysis of the reporter data is used (see Figure 5B and Section “Fitting of reporter data” in the Appendix). All kinetic parameters are given in Table A4 of Section “Model equations and parameters” in the Appendix. Our simulations show that the target gene mRNA concentration increases with increasing Wnt or decreasing APC concentrations (Figures 4E,F). There exist distinct ranges of Wnt and APC concentrations in which the mRNA levels strongly change. Outside these ranges, changes in Wnt or APC levels do not significantly affect mRNA levels. Furthermore, we observe that the APC gradient creates a slightly larger mRNA concentration range than the Wnt gradient.

Next we compared the effects of the cooperative target gene activation with Hill coefficients n = 2, 3, or 5 to that of a linear transcription rate in the APC gradient. The linear mRNA production rate (v10) is described by Eq. A24 (Section “Model equations and parameters” in the Appendix). We observed that the possible range of the relative mRNA concentration within the APC gradient is narrower for the cooperative activation than for a linear activation (Figure 6A) due to the saturation of the Hill-type activation. Figure 6C shows that the APC concentration range in which the mRNA levels change is broad in the case of the mass action kinetics (Figure 6C, first line) whereas the range is narrow in the case of cooperative activation (Figure 6C, lines 2–4). The larger the Hill coefficient is, the smaller the APC concentration range in which the mRNA levels change. The absolute APC concentration for which the mRNA levels change decreases the larger the KM value is (Eq. A25), as shown in Figures 6B,D. If the KM value becomes larger than the maximal possible β-catenin/TCF complex concentration in our model (KM > 450 nM), KM additionally influences the width of the APC concentration range in which the mRNA levels change. Under these conditions, KM also reduces the possible relative mRNA concentration range.

Impact of the mRNA expression kinetics and their parameters on the mRNA concentration profile within an APC concentration gradient. (A) The steady state concentration of the target gene mRNA is compared assuming a mass action activation (dashed black line) and cooperative activation with Hill coefficients n = 2, 3, and 5 (black, dark gray, and light gray, respectively) and KM = 23 nM in all three cases. For each curve the steady state is normalized to the respective maximal mRNA concentration. (B) The steady state concentration of the target gene mRNA is compared assuming cooperative activation with KM values of 10, 23, 100, and 500 nM (black, dark gray, gray, and light gray, respectively). The Hill coefficient is set to n = 3 in all cases. For each curve the steady state is normalized to the respective maximal mRNA concentration. (C,D) The absolute value of the concentration changes along the curves of (A,B), respectively, are compared. The maximal calculated change in each panel is color-coded in black.

 Repression of Wnt/β-catenin target genes

E-cadherin is a well-known example of negatively regulated Wnt target genes (Huber et al., 1996; Jamora et al., 2003). Also several liver-specific Wnt target genes such as Gls2, Arg1 (arginase 1), and Cps1 (carbamoyl-phosphatase 1) have been shown to be repressed by Wnt/β-catenin signaling (Benhamouche et al., 2006). In most cases the mechanisms mediating the repression have not been characterized in detail. In the case of E-cadherin, a direct repression has been demonstrated (Huber et al., 1996; Jamora et al., 2003). To investigate a direct negative regulation mechanism, we integrated an inhibitory function in the rate of transcription (Eq. A28; the kinetic parameters are given in Table A5 of Section “Model equations and parameters” in the Appendix).

First, the consequences of a negative regulation of a target gene within the minimal model were analyzed within a Wnt gradient (Figure 4G). We observed a decreasing steady state concentration of the mRNA along the increasing Wnt gradient. However, the mRNA level is kept in a very restricted concentration range. The underlying reason is the limitation of the inhibition by the available β-catenin/TCF concentration in the gradient (Figure 4A).

In the APC gradient a critical concentration threshold of APC must be reached to induce a target gene response (Figure 4H). At APC concentrations above this threshold, the mRNA levels increase with the APC concentration until saturation is reached at very high APC levels. The range of the mRNA concentration covered by the APC gradient is larger than that by the Wnt gradient (Figures 4H,G). The inhibition constant Ki (Eq. A28) strongly determines the critical APC concentration above which the mRNA level react to changes in the APC concentration (Figure 7). The larger Ki is, the smaller the critical concentration threshold. Furthermore, Ki influences the possible relative mRNA concentration range; larger Ki values reduce this range.

Dependence of the critical APC concentration for target gene induction on the inhibition constant. (A) The steady state concentration of target gene mRNA is compared assuming different values of the inhibition constant Ki: 23, 100, and 500 nM (black, dark gray, and light gray, respectively). For each curve, the steady state is normalized to the respective maximal mRNA concentration. (B) The absolute values of the concentration changes along the curves of (A) are plotted. The maximal calculated change is color-coded in black.

 Transcriptional feedback regulation via TCF family members

The Wnt/β-catenin signaling pathway is regulated by several types of transcriptional feedback (Logan and Nusse, 2004). Cases of negative feedback via Axin-2 or Dickkopf (DKK) are prominent examples that have been investigated both experimentally (Leung et al., 2002; Lustig et al., 2002; Niida et al., 2004) and theoretically (Wawra et al., 2007; Goldbeter and Pourquie, 2008; Jensen et al., 2010). Here, we concentrate on possible transcriptional feedbacks via TCF family members (Roose et al., 1999; Hovanes et al., 2001; Saegusa et al., 2005; Vadlamudi et al., 2005; Driskell et al., 2007). To analyze the impact of this transcriptional feedback on Wnt/β-catenin target gene expression, we replaced the constant TCF production rate (Eq. A13) with an mRNA-dependent translation rate (Eq. A30; parameter given in Table A6 of Section “Model equations and parameters” in the Appendix). In the following we study the effect of the different target gene expression mechanisms in the presence of a TCF feedback, beginning with the iFFL motif.

As described in Section “Wnt/β-catenin target gene regulation by an incoherent feedforward loop motif” the iFFL motif renders the target gene expression level independent of the Wnt concentration present. The addition of a feedback (i.e., replacing Eq. A13 by Eq. A30) does not change this observation (Figure 8A). The comparison of Figures 4C and 8A shows that the target gene steady state levels are not changed by the addition of the feedback. These observations also hold for the mRNA level within the APC gradient (Figure 8B). The feedback leads to the occurrence of a second steady state solution at very low mRNA concentrations in both gradients. However, this steady state is instable and can therefore not be detected in experiments.

Steady state responses of differently regulated target genes in the presence of a TCF feedback. Steady state responses of target genes regulated by an iFFL motif (A,B), in a cooperative manner with the Hill coefficient n = 3 (C,D), and by inhibition (E,F) are shown within (A,C,E) a Wnt gradient or (B,D,F) an APC gradient. The color-code denotes stable steady states as black solid lines and instable steady states as gray dashed lines.

By combining the cooperative activation of mRNA production with the feedback mechanism, a discontinuous dependence of mRNA expression level on the Wnt or APC concentration can arise (Figures 8C,D). At low Wnt concentrations, mRNA levels remain low until a critical Wnt concentration is reached (Figure 8C). If the Wnt concentration exceeds this critical value, the mRNA concentration rises instantaneously to a high level. Beyond the critical Wnt concentration, a further increase of the Wnt concentration only leads to a minor increase of the mRNA expression level. A decrease of the Wnt concentration results in the drop of the mRNA concentration to the initial level. This drop also occurs instantaneously, but at a lower Wnt concentration than the critical concentration mentioned above. In the range between the two critical Wnt concentrations, cells of low and high mRNA expression levels may coexist. Such a bistability, that is, the coexistence of low and high mRNA-expressing cells, is also possible in an APC gradient (Figure 8D).

The comparison of the absolute mRNA concentrations in the models with and without the feedback (Figures 8C and 4E, respectively) reveals that within the Wnt gradient, the mRNA can reach higher maximal levels in the presence of the feedback compared to its absence. In the APC gradient, the mRNA concentration range does not significantly change due to the feedback (Figures 4F and 8D). The reason is that for low APC levels, the β-catenin/TCF concentration is already much larger than the KM value in the transcription rate (Eq. A25). A further increase of the β-catenin/TCF concentration due to the feedback does therefore hardly affect the transcription rate and consequently the mRNA concentration.

The combination of transcriptional target gene inhibition with a feedback mechanism creates a negative feedback loop. Under certain conditions, negative feedback loops are able to create limit cycle oscillations around an instable steady state. However, we detected only stable steady state solutions in our minimal model for the given kinetic parameters (Figures 8E,F). By comparing the model with and without the feedback we observe that the feedback reduces the concentration range of the target gene mRNA. In particular, at high Wnt concentrations, the mRNA concentration is higher in the presence of the feedback than in its absence (Figures 8E and 4G, respectively). This implies that the feedback mechanism overall weakens the inhibitory effect of β-catenin/TCF on the mRNA expression. A similar effect can be observed in the APC gradient. At low APC concentrations the presence of the feedback weakens the repressive effect of the β-catenin/TCF complex on the mRNA expression. This results in a higher mRNA concentration in the presence compared to the absence of the feedback mechanism (Figures 8F and 4H). Despite its impact on the mRNA concentration, the feedback mechanism hardly alters the concentration range along the APC gradient in which the mRNA concentration is sensitive to APC concentration changes (Figure 9A, lines 4 and 5).

Simulation of the distinct target gene expression patterns within an APC gradient. (A) The relative expression profiles along an APC gradient are simulated for different target genes each subject to a different regulatory mechanism as indicated. (B) Relative mRNA expression at two positions in the APC gradient (APC = 10 nM, left, and APC = 100 nM, right) for the six indicated gene regulatory mechanisms. The relative expression under wild type and mutated conditions is shown as well as the absolute difference (Δ) as a measure of the impact strength of the mutation (color-coded in blue and orange, respectively). The mutation exclusively affects the APC-dependent β-catenin degradation (reaction 3); the parameter k3 in Eq. A10 is set to 1% of its original value in the mutant case.

 The effects of carcinogenic mutations on target gene expression

Mutations in Axin or β-catenin frequently lead to a constitutive activation of Wnt/β-catenin signaling in hepatic cancer cells (Giles et al., 2003; Takigawa and Brown, 2008). These mutations often interfere with the APC-dependent degradation of β-catenin. In our model, they all affect reaction 3 and can be simulated by changing the rate constant of that reaction. To realize one representative mutation, we set the rate constant to 1% of its original value. In Figure 9B, the mRNA expression level of the mutant situation is simulated at two positions in the APC gradient and compared to the corresponding wild type situation.

For an APC concentration of 100 nM, the mRNA levels of all but the iFFL-regulated genes are strongly changed by the mutation. In contrast, with an APC concentration of 10 nM, the same mutation strongly affects only mass-action-regulated target genes. Expression levels of negatively regulated target genes are weakly changed by the mutation. The mutation does not significantly augment the mRNA levels of cooperatively activated target genes, since the critical APC concentration for target gene activation has been passed, resulting in already high mRNA expression. We observe that at 100 nM as well as 10 nM APC, iFFL-regulated gene expression is not influenced by the mutation. The results show that the impact of a mutation on the mRNA concentration depends on both the APC content of the cell and the specific regulatory mechanism of the target gene.

 Discussion

In this study, we qualitatively analyzed the implications of different gene regulatory mechanisms in time-independent gradients of concentrations of either Wnt or APC. To this end, we developed a minimal model of Wnt/β-catenin signal transduction and extended it to different regulatory mechanisms of target gene expression. Our analysis demonstrated that combining these mechanisms with a gradient in the concentration of either Wnt or APC is sufficient to generate spatially distinct target gene expression patterns as have been experimentally observed in liver (Benhamouche et al., 2006; Braeuning et al., 2006). The gene regulatory mechanism and the available concentration of APC also determine the impact of a mutation on the target gene expression.

The development of an integrative model including signal transduction and different gene regulatory mechanisms has allowed us to investigate how the concentration of Wnt or APC is translated via the signaling processes and transcriptional regulation into target gene expression. So far, most studies have focused either on the impact of a gradient in the context of the signal transduction (Melen et al., 2005; Murray et al., 2010) or in the context of gene regulatory networks (Balaskas et al., 2012) but did not interlace both aspects. Our analysis showed that the differential integration of the Wnt and the APC concentration into the signaling processes has two consequences. First, the signaling pathway transmits the Wnt and APC gradient into opposing β-catenin/TCF complex concentration profiles (Figures 3B,C). Second, the signaling processes limit the possible concentration range of the transcriptional active β-catenin/TCF complex along the gradients. By this means, the possible Wnt or APC concentration range is confined to a narrower effective range that influences target gene mRNA levels.

Applying our model, we investigated how the different gene regulatory mechanisms influence target gene expression within Wnt or APC concentration gradients. Our simulations showed that direct activation and inhibition mechanisms for target gene mRNA can generate opposite expression patterns (Figure 9A, lines 1 and 4, respectively, for an APC gradient). Such opposing expression profiles have also been detected by in situ hybridization experiments of liver sections. In regions with high expression of positively regulated target genes such as Lect2 and Axin-2, the expression of negatively regulated target genes such as Gls2 and Arg1 was low and vice versa (Colnot and Perret, 2011).

Until now, details on the molecular regulation of target gene expression are not available in the literature. For the characterization of the activation of gene expression, we investigated different possible activation kinetics using reporter construct data (Biechele and Moon, 2008). Our analysis indicated that the expression of these constructs occurs in a cooperative manner; the data were better approximated by Hill-type kinetics than by linear or Michaelis–Menten kinetics (Figure 5A; Figure A2 and Table A7 in Appendix). The Hill coefficients, which are a measure of cooperativity, were quantified to be in the range of 2–5. Interestingly, the Hill coefficients do not depend linearly on the number of TBEs in the constructs. While an increase in the number of TBEs from 3 to 8 is accompanied by a considerable increase in the Hill coefficient, an increase from 8 to 12 TBEs hardly changes the Hill coefficient. One might speculate that a critical number of TBEs is necessary to realize maximal cooperativity; any additional TBE may affect the readout intensity but not cooperativity. Whether these results also hold for natural Wnt/β-catenin target genes has to be verified experimentally. An indication for the cooperative activation of natural Wnt/β-catenin target genes is provided by site-directed mutation studies in the siamois promoter that harbors three TBEs. The elimination of one TBE weakly reduced the response of target gene expression to stimulation while the elimination of any pair of sites reduced most of the response (Brannon et al., 1997).

Taking cooperative activation of target gene expression into account, the simulations revealed that the Hill coefficient correlates with the steepness of the change in the mRNA level (Figures 6A,C). Furthermore, we demonstrated that a combination of a cooperative activation mechanism with n = 3 and a feedback via TCF is sufficient to create an all-or-none response for target gene expression (Figure 9A, line 3). Thus, cooperative activation with large Hill coefficients, or in combination with feedback, could create sharp borders in the zonation of target gene expression. Positive feedbacks and cooperativity have already been described to establish switch-like responses in other systems (Tyson et al., 2003; Melen et al., 2005). Whether the sharp expression profile of liver-specific Wnt target genes such as GS (Gebhardt et al., 2007; Colnot and Perret, 2011) is established by a cooperative regulation and to what extent additional mechanisms contribute remains to be investigated. Generally, several additional mechanisms are discussed in the literature as contributors to the zonation of target gene expression in the liver, e.g., (i) the integration of different signaling pathways (Hailfinger et al., 2006; Gebhardt and Hovhannisyan, 2010; Archbold et al., 2012), (ii) varying microenvironments (Jungermann and Kietzmann, 1996; Torre et al., 2010; Colnot and Perret, 2011), and (iii) differences in the activation of the Wnt/β-catenin pathway (Torre et al., 2010; Colnot and Perret, 2011). In the case of GS, a contribution of other signaling pathways induced by oxidative stress or growth hormones has been suggested (Gebhardt and Hovhannisyan, 2010).

Different target genes regulated by the same mechanism may still be expressed in different regions along the periportal-pericentral axis in the liver. Examples are the two negatively regulated target genes Gls2 and Arg1. While Gls2 is only expressed in the proximal periportal region, Arg1 is expressed in the proximal and distal periportal area (Benhamouche et al., 2006). Our theoretical analysis showed that the critical Wnt or APC concentration, at which changes in mRNA levels occur, strongly depends on the parameters Ki in target gene repression or KM in target gene activation. Different parameter values shift the APC concentration range in which changes of gene expression occurs (Figures 6B,D as well as Figure 7). Thus, differences in kinetic parameters might be an explanation for the different expression areas of Gls2 and Arg1.

While the activation and inhibition mechanisms discussed above cause changes in the mRNA level along a Wnt or an APC gradient, the iFFL motif renders the steady state of the target gene independent of the gradient (Figure 9A, line 6). The steady state concentration of the target gene throughout the gradient is also insensitive to the addition of the feedback via TCF (Figures 8A,B). Provided that the iFFL motif exists and operates in hepatic gene regulation, it might be a mechanism to neutralize the impact of the gradients in order to establish constant expression levels of genes involved in zonation-independent processes. In that respect, it is interesting to note that not all metabolic processes are locally confined in the liver; e.g. the synthesis of serum proteins (Colnot and Perret, 2011).

Signal transduction in the Wnt/β-catenin pathway is perturbed in several types of cancer (Giles et al., 2003). We analyzed mutations in Axin and β-catenin that are frequently detected in hepatic cancer. Our study demonstrated that the impact of the mutation on the mRNA concentration depends on the specific gene regulatory mechanism as well as the APC concentration in the mutant cell. This implies that the same mutation may result in different phenotypes (i.e., gene expression profiles) depending on the APC content of the cell in which the mutation occurs. A prominent example is the expression profile of genes regulated by cooperative activation. While the mutation has a strong impact on the mRNA steady state concentration at 100 nM APC, there is almost no effect if this mutation occurs at 10 nM APC (Figure 9B). The analysis shows that at both APC concentrations considered here, the expression of genes regulated by a mass action mechanism is strongly susceptible to mutations. In contrast, the steady state expression of iFFL-regulated genes is not changed by the mutation at either APC concentrations.

In our study we simplified the mutation mechanisms and only considered their potential to modify the APC-dependent degradation of β-catenin. It might be speculated that a mutation in β-catenin not only affects its stability but also its transcriptional activity. In an experimental study, the expression of target genes in hepatocellular carcinoma harboring a mutation in either Axin or β-catenin has been measured (Zucman-Rossi et al., 2007). This study has revealed that gene expression might be different depending on whether Axin or β-catenin is mutated. These differences cannot be reproduced by our minimal model due to the strong simplification made for the molecular mechanisms involved in APC-dependent β-catenin degradation and the simplified implementation of the mutation.

In our modeling approach, we consider individual cells independently of their neighbors. However, some Wnt/β-catenin target genes encode for proteins which affect intercellular signaling, the gradient and/or its effectiveness. For instance, DKK1 is a Wnt/β-catenin target gene and its product influences Wnt/β-catenin signaling of the surrounding cells by competing with Wnt ligands for receptor binding (Mao et al., 2001; Chamorro et al., 2005). The expression and cell surface representation of Wnt receptors LRP and Frizzled is also regulated by Wnt/β-catenin signaling (Hao et al., 2012). Target gene products that modulate the effectiveness of the gradient are discussed to increase the stability against fluctuations in the gradient (Jaeger et al., 2008). Other theoretical investigations have shown that the exchange of molecules between cells in a gradient may produce even more complex spatio-temporal behavior (Schütze and Wolf, 2010).

Our minimal signaling model reproduces the β-catenin steady states and dynamics upon constant Wnt stimulation of the detailed model and in fact extends the detailed model with respect to the representation of β-catenin/TCF complex concentration ranges and the inclusion of gene regulatory mechanisms. The parameters and total concentrations of the minimal model are predominantly derived from the detailed model (Lee et al., 2003). There are no comprehensive experimental time series data sets available to estimate the kinetic parameters of mammalian cell types. Recent experiments on the quantification of pathway components in resting mammalian cells show that key components such as Axin, APC, and β-catenin may be expressed in a highly cell-type specific manner (Chen et al., 2010; Schwanhausser et al., 2011; Tan et al., 2012). To what extent variations in the absolute concentrations have consequences for the signaling dynamics and downstream gene expression remains to be investigated. Our minimal model with its extension toward target gene regulation provides a useful and versatile basis for future modeling approaches of the Wnt/β-catenin pathway.

 Materials and Methods Numerical simulations

The changes of the concentrations of the system’s components are described by a set of ODEs and algebraic conservation equations. The particular equations and parameters are listed in Section “Model equations and parameters” in the Appendix. All calculations are done with Mathematica 8.0 (Wolfram Research, Inc.). Steady state solutions are numerically obtained.

 Fitting of reporter data

Published experimental data of β-catenin reporter constructs (Biechele and Moon, 2008) are analyzed to gain insights into possible mechanisms of target gene activation. Three possible kinetics are considered: (i) linear mass action kinetics, (ii) Michaelis–Menten kinetics, and (iii) Hill-type activation. In 500 independent fitting runs each, the kinetics (i), (ii), and (iii) are fitted to the data using a log-likelihood approach. The AICc values are calculated to identify the kinetics which explains the data best. Details on the fitting procedure are provided in Section “Fitting of reporter data” in the Appendix.

 Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Appendix Model equations and parameters The signal transduction module of the minimal model

Figure 1 shows the schematic representation of the minimal model illustrating two different mechanisms of Wnt/β-catenin target gene regulation. In both schemes, the signal transduction module (highlighted with a gray background) is identical. It consists of reactions 1–9 and is given by the Eqs A1–A16. In the schemes, the number next to an arrow denotes the number of the particular reaction. Components in a complex are separated by a slash.

Differential equations d(β_catenin)dt=v1-v2-v3-v4-v5 d(APC)dt=-v4 d(TCF)dt=-v5+v6-v7 d(β_catenin∕TCF)dt=v5 d(Dsha)dt=v8-v9

Conservation relations Dshtot=Dshi+Dsha APCtot=APC+(β_catenin∕APC)

Rate equations v1=const. v2=k2⋅β_catenin v3=k3⋅APC⋅β_cateninK+Dsha v4=k4⋅APC⋅β_catenin-k-4⋅(β_catenin∕APC) v5=k5⋅TCF⋅β_catenin-k-5⋅(β_catenin∕TCF) v6=const. v7=k7⋅TCF v8=k8⋅Dshi⋅Wnt v9=k9⋅Dsha

 Regulation of the target gene expression by the iFFL motif

The effects of the iFFL motif (Figure 1B) are investigated in Section “Wnt/β-catenin target gene regulation by an incoherent feedforward loop motif”. For the analyses, the signal transduction module (Eqs A1–A16) is extended by the Eqs A17–A22.

Additional differential equations of the iFFL motif d(repressor)dt=v12-v13 d(mRNA)dt=v10-v11

Additional rate equations of the iFFL motif v10=v10max⋅(β−catenin/TCF)K11+(β−catenin/TCF)K1+repressorK2+(β−catenin/TCF)⋅repressorK3 v11=k11⋅mRNA v12=k12⋅(β_catenin∕TCF) v13=k13⋅repressor v10max is the maximal transcription rate, K1, K2, and K3 are effective binding constants.

 Activation of the target gene expression

In Section “Activation of Wnt/β-catenin target genes”, the differences between linear and cooperative activation of the mRNA production by β-catenin/TCF complexes are analyzed. To this end, the signal transduction module is extended by either a linear or a Hill kinetics mRNA production term. All together, the model consists of Eqs A1–A16 (signal transduction module), Eq. A23 (ODE of the temporal change of the mRNA), Eq. A26 (degradation of the mRNA), and either Eq. A24 or Eq. A25 in the case of linear or Hill-type activation kinetics of the mRNA production, respectively.

Additional differential equations of the activation of the gene expression d(mRNA)dt=v10-v11

Additional rate equations of the activation of the gene expression v10=k10⋅(β_catenin∕TCF) v10=v10max⋅(β_catenin∕TCF)nKMn+(β_catenin/TCF)n+v10basal v11=k11⋅mRNA v10max is the maximal inducible transcription rate, v10basal is an additional basal transcription rate, KM describes the β-catenin/TCF complex concentration at half maximal inducible transcription rate, and n is the Hill coefficient, which represents a measure of the cooperativity between the individual β-catenin/TCF complexes.We use a Hill coefficient of three throughout the manuscript unless stated otherwise. Note that for n = 1, Eq. A25 describes a Michaelis–Menten kinetics.

 Inhibition of target gene expression by β-catenin/TCF complexes

The effects of repression of mRNA production by the β-catenin/TCF complex are investigated in Section “Repression of Wnt/β-catenin target genes”. For these analyses, the model consists of the Eqs A1–A16 (signal transduction module), Eq. A27 (ODE of the temporal change of the mRNA), Eq. A29 (mRNA degradation), and Eq. A28 which describes the inhibition of the mRNA production by the β-catenin/TCF complex.

Additional differential equations of the inhibition of the mRNA expression d(mRNA)dt=v10-v11

Additional rate equations of the inhibition of the mRNA expression v10=v10max⋅11+(β−catenin/TCF)Ki v11=k11⋅mRNA v10max is the maximal transcription rate; Ki is the inhibition constant.

Kinetic parameters of the signal transduction module.

Parameter Value
v1 0.423 nM min−1
k2 2.57 · 10−4 min−1
k3 3.08 · 10−3 min−1
K 18 nM
k4 105 nM−1 min−1
k−4 1.2 · 108 min−1
k5 0.0333 nM−1 min−1
k−5 1 min−1
v6 0.686 nM min−1
k7 0.084 min−1
k8 8 · 10−3 nM−1 min−1
k9 6.7 · 10−4 min−1

Total protein concentrations of proteins obeying a conservation relation.

Component Value
APCtot 100 nM
Dshtot 100 nM
Transcriptional feedback via TCF

In the Sections “Wnt/β-catenin target gene regulation by an incoherent feedforward loop motif” to “Repression of Wnt/β-catenin target genes”, the consequences of the three different regulatory mechanisms are analyzed under the assumption of a constant production of TCF (Eq. A13). In Section “Transcriptional feedback regulation via TCF family members”, the effect of the different regulatory motifs is investigated in the presence of a transcriptional feedback via TCF. In this case, the TCF synthesis is not constant but depends on the mRNA level. For the analyses of transcriptional feedback effects, we thus replace the Eq. A13 by Eq. A30 in all cases.

v6=k6⋅mRNA

Kinetic parameters of the iFFL motif.

Parameter Value
v10max 1368.18 nM min–1
k11 10–2 min–1
k12 4.44 · 10–3 min–1
k13 10–3 min–1
K1 1300 nM
K2 0.03 nM
K3 4 · 104 nM

Kinetic parameters of the linear and Hill-type activation of the gene expression.

Parameter Value
v10max 0.03 nM min–1
v10basal 6.1 · 10–3 nM min–1
n 3
k10 0.03 min–1
k11 10–2 min–1
KM 23 nM

Additional kinetic parameters of the repression motif.

Parameter Value
v10max 8.84 · 10–3 nM min–1
k11 10–2 min–1
Ki 23 nM

Kinetic parameter of the TCF feedback.

Parameter Value
k6 1 min–1
Fitting of reporter data

β-Catenin/TCF reporter experiments have been published in which the response of reporter constructs with 3, 8, and 12 TBEs to different Wnt concentrations has been compared (Biechele and Moon, 2008). The data that has been the basis of Figure 8.1.A in Biechele and Moon (2008) has been kindly provided. It tables the Wnt concentration as “percentage of Wnt conditioned medium” with the corresponding mean and standard deviation of three normalized reporter readout replicates.

To test for possible cooperativity between β-catenin/TCF complexes in the regulation of reporter expression, we first calculated the β-catenin/TCF complex concentrations for the given “percentage of Wnt conditioned medium” using our signal transduction module (Eqs A1–A16). For the calculations, we assumed that the “percentage of Wnt conditioned medium” equals the respective percentage of the simulated Wnt concentration of 1 nM. Moreover, following the basic assumption in reporter experiments, we consider the fold change of the reporter readout, i.e., Wnt-stimulated to control measurement, to be linearly correlated to the fold change of the reporter mRNA. The reporter mRNA concentration change over time can be described by the following ODE: d(mRNA)dt=vinduced+vbasal-kdegradation⋅mRNA in which vinduced denotes the β-catenin/TCF-complex-dependent transcriptional activation rate, vbasal is a basal transcription rate, and kdegradation denotes the degradation rate constant in the linear reporter mRNA degradation rate. As the reporter readout is measured at steady state condition, we derive the mRNA steady state concentration (mRNAstst) from Eq. A31: mRNAstst=vinduced+vbasalkdegradation

Corrected Akaike information criterion of the best fits.

Number of TBEs Linear kinetics Michaelis–Menten kinetics Hill-type kinetics
3 159.5 165.1 29.0
8 558.7 564.3 51.9
12a 7021.6 7027.2 66.1
12b 16728.3 16733.9 328.0

For all three kinetics, (i) linear mass action kinetics, (ii) Michaelis–Menten kinetics, and (iii) Hill-type activation, the calculated values of the corrected Akaike information criterion are given for the four data sets.

aTransiently transfected reporter construct with 12 TBEs; bstable integrated reporter construct with 12 TBEs.

In the absence of a Wnt stimulus, i.e., vinduced = 0, this term simplifies to the basal reporter mRNA steady state (mRNAbasal). Since the experimental data shows a standard deviation in the control data point, we take a small error ε in the basal reporter mRNA steady state into account.

mRNAbasal=vbasal±εkdegradation

The experiments report the fold change of the reporter mRNA, which gives the ratio of the mRNA steady state to the basal mRNA steady state. According to Eqs A32 and A33, the fold change is given by: foldchange=a⋅vinduced+c where a=1vbasal±ε is a scaling factor and c=vbasalvbasal±ε is the offset constant along the ordinate.

Next, we considered three possible transcriptional activation kinetics in vinduced in Eq. A34: (i) linear mass action kinetics, (ii) Michaelis–Menten kinetics, and (iii) Hill-type activation: (i)foldchange=a⋅(β_catenin/TCF-x0)+c (ii) fold change=a⋅(β−catenin/TCF−x0)KM+(β−catenin/TCF−x0)+c (iii) fold change=a⋅(β−catenin/TCF−x0)nKMn+(β−catenin/TCF−x0)n+c where a is the scaling constant, x0 and c are constants accounting for an offset along the abscissa and ordinate, respectively. KM describes the β-catenin/TCF complex concentration at half maximal transcriptional activation rate and the Hill coefficient n represents a measure of the cooperativity between the individual β-catenin/TCF complexes.

We fitted the functions (i), (ii), and (iii) to the data set of each reporter construct. In the fitting procedure the log-likelihood function, that is the weighted sum of squares of the residuals between the data and the simulations of the activation kinetics, was minimized. Each residual was weighted with respect to the corresponding standard deviation of the data point. We performed 500 independent fitting runs choosing different starting parameter values for each construct to obtain the global minimum of the log-likelihood function. Considering the best fit for each function we compared the calculated values of the AICc to identify the kinetics which explains the data best (Hurvich and Tsai, 1995; Lukacs et al., 2010). The AICc was calculated by: AICc=2⋅p-2⋅lnL+2⋅p⋅(p+1)d-p-1 where p is the number of parameters of the tested kinetics, d denotes the sample size (i.e., d = 8 data points), and ln[L] is the log-likelihood function. Our calculations show that the Hill-type activation kinetics (iii) yields the smallest AICc value and thus approximates the data best (Figure 5A; Figure A2) suggesting a cooperative regulation of the reporter expression.

Input (β-catenin/TCF) and output (mRNA) time courses of the iFFL motif. (A) A stepwise increase of the β-catenin/TCF concentration leads to (B) a perfect adaptation of the target gene mRNA. The dynamics of the mRNA is independent of the initial absolute concentration of β-catenin/TCF if the same fold change in the β-catenin/TCF level is induced.

Comparison of the best fit functions. The readout intensity of the reporter constructs with 8, 12 (transfected), and 12 (integrated) TBEs (A–C) are plotted with standard deviations. The linear activation function (light gray), the Michaelis–Menten activation function (dark gray), and the Hill-type activation function (black) realizing the best fit (i.e., yielding the minimal log-likelihood function values) are shown for each reporter construct. Note, that the fits of the linear and the Michaelis–Menten activation lie on top of each other. The experimental data was kindly provided by Biechele and Moon (2008).

Next, the group of best fits, i.e., those whose log-likelihood function value deviated at most 20% from the minimal value of the log-likelihood function, is used to gain information on the variability of the Hill coefficients. Their analysis revealed that the Hill coefficients n are all higher than 1 (Figure 5B). The values of the Hill coefficients span from 2.0 to 4.2 for the 3 TBEs construct. For the 8 TBEs and the two different 12 TBEs constructs, the Hill coefficients take values in a narrower range between 3.9 and 4.2 with medians of 3.9, 4.0, and 3.9, respectively. In contrast, the median of the construct with the 3 TBEs is lower (2.6) than that of the other three constructs. This may point to a limitation of the number of β-catenin/TCF complexes that can act cooperatively in the regulation of the reporter constructs.

We thank Randall T. Moon and Travis L. Biechele for kindly providing the experimental data on reporter constructs, Rolf Gebhardt for stimulating discussions, Katharina Baum and Manuela Benary for help in model fitting, and Russ Hodge for reading the manuscript. Studies in this project were funded by the ForSys-programme of the German Ministry of Education and Research (grant number 0315289 to Jana Wolf), the MSBN project within the Helmholtz Alliance on Systems Biology funded by the Initiative and Networking Fund of the Helmholtz Association (to Uwe Benary, Bente Kofahl, and Jana Wolf), the Deutsche Forschungsgemeinschaft (CRC-850 subproject B5 to Andreas Hecht), and the project Virtual Liver of the German Ministry of Education and Research (grant number 0315766 TP A2.5/B2.2 to Andreas Hecht). Bente Kofahl was partly supported by the project CancerSys of the European Union (grant number EU-FP7 HEALTH-F4-2008-223188).

 Abbreviations

AICc, corrected Akaike information criterion; APC, adenomatous polyposis coli; Arg1, arginase 1; Cps1, carbamoyl-phosphatase 1; DKK, Dickkopf; Dsh, Dishevelled; Gls2, glutaminase 2; Glt1, glutamate transporter 1; GS, glutamine synthetase; iFFL, incoherent feedforward loop; Lect2, leukocyte cell-derived chemotaxin 2; LEF, lymphoid enhancer-binding factor; OAT, ornithine aminotransferase; ODE, ordinary differential equation; TBE, TCF/LEF binding elements; TCF, T-cell factor.

 References Archbold H. C. Yang Y. X. Chen L. Cadigan K. M. (2012). How do they do Wnt they do?: regulation of transcription by the Wnt/beta-catenin pathway. Acta Physiol. (Oxf.) 204, 74–109.10.1111/j.1748-1716.2011.02293.x21624092 Balaskas N. Ribeiro A. Panovska J. Dessaud E. Sasai N. Page K. M. (2012). Gene regulatory logic for reading the Sonic Hedgehog signaling gradient in the vertebrate neural tube. Cell 148, 273–284.10.1016/j.cell.2011.10.04722265416 Benhamouche S. Decaens T. Godard C. Chambrey R. Rickman D. S. Moinard C. (2006). Apc tumor suppressor gene is the “zonation-keeper” of mouse liver. Dev. Cell 10, 759–770.10.1016/j.devcel.2006.03.01516740478 Biechele T. L. Moon R. T. (2008). Assaying beta-catenin/TCF transcription with beta-catenin/TCF transcription-based reporter constructs. Methods Mol. Biol. 468, 99–110.10.1007/978-1-59745-249-6_819099249 Braeuning A. Ittrich C. Kohle C. Hailfinger S. Bonin M. Buchmann A. (2006). Differential gene expression in periportal and perivenous mouse hepatocytes. FEBS J. 273, 5051–5061.10.1111/j.1742-4658.2006.05503.x17054714 Brannon M. Gomperts M. Sumoy L. Moon R. T. Kimelman D. (1997). A beta-catenin/XTcf-3 complex binds to the siamois promoter to regulate dorsal axis specification in Xenopus. Genes Dev. 11, 2359–2370.10.1101/gad.11.18.23599308964 Burke Z. D. Tosh D. (2006). The Wnt/beta-catenin pathway: master regulator of liver zonation? Bioessays 28, 1072–1077.10.1002/bies.2048517041892 Cadigan K. M. Peifer M. (2009). Wnt signaling from development to disease: insights from model systems. Cold Spring Harb. Perspect. Biol. 1, a002881.10.1101/cshperspect.a00288120066091 Cadoret A. Ovejero C. Terris B. Souil E. Levy L. Lamers W. H. (2002). New targets of beta-catenin signaling in the liver are involved in the glutamine metabolism. Oncogene 21, 8293–8301.10.1038/sj.onc.120611812447692 Chamorro M. N. Schwartz D. R. Vonica A. Brivanlou A. H. Cho K. R. Varmus H. E. (2005). FGF-20 and DKK1 are transcriptional targets of beta-catenin and FGF-20 is implicated in cancer and development. EMBO J. 24, 73–84.10.1038/sj.emboj.760046015592430 Chen Y. Gruidl M. Remily-Wood E. Liu R. Z. Eschrich S. Lloyd M. (2010). Quantification of beta-catenin signaling components in colon cancer cell lines, tissue sections, and microdissected tumor cells using reaction monitoring mass spectrometry. J. Proteome Res. 9, 4215–4227.10.1021/pr100205k20590165 Cho K. H. Baek S. Sung M. H. (2006). Wnt pathway mutations selected by optimal beta-catenin signaling for tumorigenesis. FEBS Lett. 580, 3665–3670.10.1016/j.febslet.2006.04.06216764864 Clevers H. (2006). Wnt/beta-catenin signaling in development and disease. Cell 127, 469–480.10.1016/j.cell.2006.10.01817081971 Colnot S. Perret C. (2011). “Liver zonation,” in Molecular Pathology of Liver Diseases, ed. Monga S. P. S. (Springer: Molecular Pathology Library), 7–16. Driskell R. R. Goodheart M. Neff T. Liu X. Luo M. Moothart C. (2007). Wnt3a regulates Lef-1 expression during airway submucosal gland morphogenesis. Dev. Biol. 305, 90–102.10.1016/j.ydbio.2007.01.03817335794 Gebhardt R. Baldysiak-Figiel A. Krugel V. Ueberham E. Gaunitz F. (2007). Hepatocellular expression of glutamine synthetase: an indicator of morphogen actions as master regulators of zonation in adult liver. Prog. Histochem. Cytochem. 41, 201–266.10.1016/j.proghi.2006.12.00117368308 Gebhardt R. Hovhannisyan A. (2010). Organ patterning in the adult stage: the role of Wnt/beta-catenin signaling in liver zonation and beyond. Dev. Dyn. 239, 45–55.19705440 Giles R. H. Van Es J. H. Clevers H. (2003). Caught up in a Wnt storm: Wnt signaling in cancer. Biochim. Biophys. Acta 1653, 1–24.12781368 Goentoro L. Kirschner M. W. (2009). Evidence that fold-change, and not absolute level, of beta-catenin dictates Wnt signaling. Mol. Cell 36, 872–884.10.1016/j.molcel.2009.11.01720005849 Goentoro L. Shoval O. Kirschner M. W. Alon U. (2009). The incoherent feedforward loop can provide fold-change detection in gene regulation. Mol. Cell 36, 894–899.10.1016/j.molcel.2009.11.01720005851 Goldbeter A. Pourquie O. (2008). Modeling the segmentation clock as a network of coupled oscillations in the Notch, Wnt and FGF signaling pathways. J. Theor. Biol. 252, 574–585.10.1016/j.jtbi.2008.01.00618308339 Gurdon J. B. Bourillot P. Y. (2001). Morphogen gradient interpretation. Nature 413, 797–803.10.1038/3510150011677596 Hailfinger S. Jaworski M. Braeuning A. Buchmann A. Schwarz M. (2006). Zonal gene expression in murine liver: lessons from tumors. Hepatology 43, 407–414.10.1002/hep.2108216496347 Hao H. X. Xie Y. Zhang Y. Charlat O. Oster E. Avello M. (2012). ZNRF3 promotes Wnt receptor turnover in an R-spondin-sensitive manner. Nature 485, 195–200.10.1038/nature1101922575959 Hecht A. Kemler R. (2000). Curbing the nuclear activities of beta-catenin. Control over Wnt target gene expression. EMBO Rep. 1, 24–28.10.1093/embo-reports/kvd01211256619 Hovanes K. Li T. W. Munguia J. E. Truong T. Milovanovic T. Lawrence Marsh J. (2001). Beta-catenin-sensitive isoforms of lymphoid enhancer factor-1 are selectively expressed in colon cancer. Nat. Genet. 28, 53–57.10.1038/ng0501-5311326276 Huber O. Korn R. McLaughlin J. Ohsugi M. Herrmann B. G. Kemler R. (1996). Nuclear localization of beta-catenin by interaction with transcription factor LEF-1. Mech. Dev. 59, 3–10.10.1016/0925-4773(96)00597-78892228 Hurvich C. M. Tsai C. L. (1995). Model selection for extended quasi-likelihood models in small samples. Biometrics 51, 1077–1084.10.2307/25330067548692 Jaeger J. Irons D. Monk N. (2008). Regulative feedback in pattern formation: towards a general relativistic theory of positional information. Development 135, 3175–3183.10.1242/dev.01869718776142 Jamora C. Dasgupta R. Kocieniewski P. Fuchs E. (2003). Links between signal transduction, transcription and adhesion in epithelial bud development. Nature 422, 317–322.10.1038/nature0145812646922 Jensen P. B. Pedersen L. Krishna S. Jensen M. H. (2010). A Wnt oscillator model for somitogenesis. Biophys. J. 98, 943–950.10.1016/j.bpj.2009.12.156720303851 Jho E. H. Zhang T. Domon C. Joo C. K. Freund J. N. Costantini F. (2002). Wnt/beta-catenin/Tcf signaling induces the transcription of Axin2, a negative regulator of the signaling pathway. Mol. Cell. Biol. 22, 1172–1183.10.1128/MCB.22.4.1172-1183.200211809808 Jungermann K. Kietzmann T. (1996). Zonation of parenchymal and nonparenchymal metabolism in liver. Annu. Rev. Nutr. 16, 179–203.10.1146/annurev.nu.16.070196.0011438839925 Kimelman D. Xu W. (2006). Beta-catenin destruction complex: insights and questions from a structural perspective. Oncogene 25, 7482–7491.10.1038/sj.onc.121005517143292 Klaus A. Birchmeier W. (2008). Wnt signalling and its impact on development and cancer. Nat. Rev. Cancer 8, 387–398.10.1038/nrc238918432252 Kofahl B. Wolf J. (2010). Mathematical modelling of Wnt/beta-catenin signalling. Biochem. Soc. Trans. 38, 1281–1285.10.1042/BST038128120863299 Lee E. Salic A. Kruger R. Heinrich R. Kirschner M. W. (2003). The roles of APC and Axin derived from experimental and theoretical analysis of the Wnt pathway. PLoS Biol. 1:e10.10.1371/journal.pbio.0000010 Leung J. Y. Kolligs F. T. Wu R. Zhai Y. Kuick R. Hanash S. (2002). Activation of AXIN2 expression by beta-catenin-T cell factor. J. Biol. Chem. 277, 21657–21665.10.1074/jbc.M20013920011940574 Logan C. Y. Nusse R. (2004). The Wnt signaling pathway in development and disease. Annu. Rev. Cell Dev. Biol. 20, 781–810.10.1146/annurev.cellbio.20.010403.11312615473860 Lukacs P. M. Burnham K. P. Anderson D. R. (2010). Model selection bias and Freedman’s paradox. Ann. Inst. Stat. Math. 62, 117–125.10.1007/s10463-009-0234-4 Lustig B. Jerchow B. Sachs M. Weiler S. Pietsch T. Karsten U. (2002). Negative feedback loop of Wnt signaling through upregulation of conductin/Axin2 in colorectal and liver tumors. Mol. Cell. Biol. 22, 1184–1193.10.1128/MCB.22.4.1184-1193.200211809809 MacDonald B. T. Tamai K. He X. (2009). Wnt/beta-catenin signaling: components, mechanisms, and diseases. Dev. Cell 17, 9–26.10.1016/j.devcel.2009.06.01619619488 Mao B. Wu W. Li Y. Hoppe D. Stannek P. Glinka A. (2001). LDL-receptor-related protein 6 is a receptor for Dickkopf proteins. Nature 411, 321–325.10.1038/3507710811357136 Meinhardt H. (2009). Models for the generation and interpretation of gradients. Cold Spring Harb. Perspect. Biol. 1, a001362.10.1101/cshperspect.a00136220066097 Melen G. J. Levy S. Barkai N. Shilo B. Z. (2005). Threshold responses to morphogen gradients by zero-order ultrasensitivity. Mol. Syst. Biol. 1, 2005.0028.10.1038/msb410003616729063 Mirams G. R. Byrne H. M. King J. R. (2010). A multiple timescale analysis of a mathematical model of the Wnt/beta-catenin signalling pathway. J. Math. Biol. 60, 131–160.10.1007/s00285-009-0262-y19288106 Murray P. J. Kang J. W. Mirams G. R. Shin S. Y. Byrne H. M. Maini P. K. (2010). Modelling spatially regulated beta-catenin dynamics and invasion in intestinal crypts. Biophys. J. 99, 716–725.10.1016/j.bpj.2010.05.01620682248 Nahmad M. Lander A. D. (2011). Spatiotemporal mechanisms of morphogen gradient interpretation. Curr. Opin. Genet. Dev. 21, 726–731.10.1016/j.gde.2011.10.00222033220 Niida A. Hiroko T. Kasai M. Furukawa Y. Nakamura Y. Suzuki Y. (2004). DKK1, a negative regulator of Wnt signaling, is a target of the [beta]-catenin//TCF pathway. Oncogene 23, 8520–8526.10.1038/sj.onc.120789215378020 Ovejero C. Cavard C. Perianin A. Hakvoort T. Vermeulen J. Godard C. (2004). Identification of the leukocyte cell-derived chemotaxin 2 as a direct target gene of beta-catenin in the liver. Hepatology 40, 167–176.10.1016/S0168-8278(04)90571-515239100 Roose J. Huls G. Van Beest M. Moerer P. Van Der Horn K. Goldschmeding R. (1999). Synergy between tumor suppressor APC and the beta-catenin-Tcf4 target Tcf1. Science 285, 1923–1926.10.1126/science.285.5435.192310489374 Saegusa M. Hashimura M. Kuwata T. Hamano M. Okayasu I. (2005). Upregulation of TCF4 expression as a transcriptional target of beta-catenin/p300 complexes during trans-differentiation of endometrial carcinoma cells. Lab. Invest. 85, 768–779.10.1038/labinvest.370027315806138 Schütze J. Wolf J. (2010). Spatio-temporal dynamics of glycolysis in cell layers. A mathematical model. Biosystems 99, 104–108.10.1016/j.biosystems.2009.10.00219837130 Schwanhausser B. Busse D. Li N. Dittmar G. Schuchhardt J. Wolf J. (2011). Global quantification of mammalian gene expression control. Nature 473, 337–342.10.1038/nature1009821593866 Takigawa Y. Brown A. M. (2008). Wnt signaling in liver cancer. Curr. Drug Targets 9, 1013–1024.10.2174/13894500878678612718991612 Tan C. W. Gardiner B. S. Hirokawa Y. Layton M. J. Smith D. W. Burgess A. W. (2012). Wnt signalling pathway parameters for mammalian cells. PLoS ONE 7:e31882.10.1371/journal.pone.0031882 Torre C. Perret C. Colnot S. (2010). Molecular determinants of liver zonation. Prog. Mol. Biol. Transl. Sci. 97, 127–150.10.1016/B978-0-12-385233-5.00005-221074732 Tyson J. J. Chen K. C. Novak B. (2003). Sniffers, buzzers, toggles and blinkers: dynamics of regulatory and signaling pathways in the cell. Curr. Opin. Cell Biol. 15, 221–231.10.1016/S0955-0674(03)00017-612648679 Vadlamudi U. Espinoza H. M. Ganga M. Martin D. M. Liu X. Engelhardt J. F. (2005). PITX2, beta-catenin and LEF-1 interact to synergistically regulate the LEF-1 promoter. J. Cell. Sci. 118, 1129–1137.10.1242/jcs.0170615728254 van Leeuwen I. M. M. Byrne H. M. Jensen O. E. King J. R. (2007). Elucidating the interactions between the adhesive and transcriptional functions of [beta]-catenin in normal and cancerous cells. J. Theor. Biol. 247, 77–102.10.1016/j.jtbi.2007.01.01917382967 Wawra C. Kuhl M. Kestler H. A. (2007). Extended analyses of the Wnt/beta-catenin pathway: robustness and oscillatory behaviour. FEBS Lett. 581, 4043–4048.10.1016/j.febslet.2007.07.04317678900 Whangbo J. Kenyon C. (1999). A Wnt signaling system that specifies two patterns of cell migration in C. elegans. Mol. Cell 4, 851–858.10.1016/S1097-2765(00)80394-9 Wolpert L. (2011). Positional information and patterning revisited. J. Theor. Biol. 269, 359–365.10.1016/j.jtbi.2010.10.03421044633 Zecca M. Basler K. Struhl G. (1996). Direct and long-range action of a wingless morphogen gradient. Cell 87, 833–844.10.1016/S0092-8674(00)81991-18945511 Zucman-Rossi J. Benhamouche S. Godard C. Boyault S. Grimber G. Balabaud C. (2007). Differential effects of inactivated Axin1 and activated beta-catenin mutations in human hepatocellular carcinomas. Oncogene 26, 774–780.10.1038/sj.onc.120982416964294 Research Topic:

From structural to molecular systems biology: experimental and computational approaches to unravel mechanisms of kinase activity regulation in cancer and neurodegeneration

 Topic Editor:

Matteo Barberis, Humboldt University Berlin, Germany; Max Planck Institute for Molecular Genetics, Berlin, Germany