plasticity of TGF-β signaling
© Cellière et al; licensee BioMed Central Ltd. 2011
Received: 20 June 2011
Accepted: 3 November 2011
Published: 3 November 2011
Skip to main content
© Cellière et al; licensee BioMed Central Ltd. 2011
Received: 20 June 2011
Accepted: 3 November 2011
Published: 3 November 2011
The family of TGF-β ligands is large and its members are involved in many different signaling processes. These signaling processes strongly differ in type with TGF-β ligands eliciting both sustained or transient responses. Members of the TGF-β family can also act as morphogen and cellular responses would then be expected to provide a direct read-out of the extracellular ligand concentration. A number of different models have been proposed to reconcile these different behaviours. We were interested to define the set of minimal modifications that are required to change the type of signal processing in the TGF-β signaling network.
To define the key aspects for signaling plasticity we focused on the core of the TGF-β signaling network. With the help of a parameter screen we identified ranges of kinetic parameters and protein concentrations that give rise to transient, sustained, or oscillatory responses to constant stimuli, as well as those parameter ranges that enable a proportional response to time-varying ligand concentrations (as expected in the read-out of morphogens). A combination of a strong negative feedback and fast shuttling to the nucleus biases signaling to a transient rather than a sustained response, while oscillations were obtained if ligand binding to the receptor is weak and the turn-over of the I-Smad is fast. A proportional read-out required inefficient receptor activation in addition to a low affinity of receptor-ligand binding. We find that targeted modification of single parameters suffices to alter the response type. The intensity of a constant signal (i.e. the ligand concentration), on the other hand, affected only the strength but not the type of the response.
The architecture of the TGF-β pathway enables the observed signaling plasticity. The observed range of signaling outputs to TGF-β ligand in different cell types and under different conditions can be explained with differences in cellular protein concentrations and with changes in effective rate constants due to cross-talk with other signaling pathways. It will be interesting to uncover the exact cellular differences as well as the details of the cross-talks in future work.
Transforming growth factor beta (TGF-β) signaling has been implicated as an important regulator of almost all major cell behaviors, including proliferation, differentiation, cell death, and motility . Which response is induced or repressed depends on the cell type and context in which the signal is received.
The complexity of the biological outcomes elicited by TGF-β stands in stark contrast to the apparent simplicity of the signaling cascade. In response to TGF-β, type 1 (ALKs 1-7 in humans) and type 2 receptors (ActR-IIA, ActR-IIB, BMPR-II, AMHR-II and TbR-II, in humans) form complexes and the constitutively active type 2 serine/threonine kinase phosphorylates the type 1 receptor. The activated type 1 receptor transduces the signal into the cell by phosphorylating the regulatory Smads (R-Smad: Smad 2 and 3 in case of the TGF-β subfamily, and Smad 1, 5 and 8 for the BMP subfamily). Once activated R-Smads form homomeric complexes and heteromeric complexes with the common Smad, Co-Smad (Smad 4) . Smads continuously shuttle between nucleus and cytoplasm . TGF-β signaling biases Smad localisation to the nucleus  where Smad complexes associate with chromatin and regulate the transcription of hundreds of genes . Signal termination is achieved through continuous dephosphorylation of the R-Smad (mainly in the nucleus ) and induction of inhibitory Smads (I-Smad: Smad 6 for the BMP subfamily, and Smad7 for the TGF-β subfamily). I-Smads act through diverse mechanisms: by targeting active receptor for proteasomal degradation [6, 7], inducing receptor dephosphorylation  and competing with R-Smad for the receptor binding site . Rapid shuttling and inactivation enables a continuous sensing of the extracellular ligand concentrations . This is likely to be particular important when members of the TGF-β ligand family acts as morphogen and determine cell-fate in a concentration-dependent manner.
Beyond the core components of this signaling pathway many other factors modulate the signal and thereby contribute to the versality of the response. At the membrane level, the access to receptor is controlled by soluble proteins that sequester TGF-β ligand (i.e. decorin) , and by membrane-bound co-receptors that promote binding (i.e. betaglycan) . The receptor activity is further regulated by several receptor internalization routes , and by receptor turnover. Intracellularly, many processes require auxiliary proteins (i.e. SARA for the binding of R-Smad to the receptor and Schnurri for the binding of the R-Smad/Co-Smad complex to the DNA binding element) [2, 13]. The restriction of those auxiliary factors to specific cell-types will make the response cell context dependent . Diversity can also be generated by the huge number of different possible combinations of type 1 and type 2 receptors  and the multiple crosstalks of the TGF-β signaling cascade with other pathways. One example of regulation by cross-talk is the phosphorylation of R-Smads in the linker region by Ras-activated MAPK , calcium-calmodulin-dependent protein kinase II  or CDKs . Phosphorylation reduces the transcriptional activity of the R-Smad .
Several mathematical models have been developed to gain further insights into the complex TGF-β-dependent signaling network . An early model by Clarke and co-workers (2006)  focused on the nuclear accumulation of Smad complexes. Their conclusion on the central role of the imbalance between R-Smad phosphorylation and dephosphorylation rates were confirmed by a more detailed model by Schmierer et al. (2008) . Experiments suggest that the duration of the response to a ligand stimulation strongly impacts on the cellular response. Thus epithelial cells that elicit sustained nuclear Smad complex accumulation respond to TGF-β with cell growth arrest, whereas pancreatic tumor cells that elicit a transient response continue proliferating (while keeping other TGF-β induced behaviors) . Much theoretical work therefore focused on how sustained, transient, or switch-like responses could be obtained by adjusting the receptor dynamics, ligand depletion, and the I-Smad dependent negative feedback. Melke et al. (2006)  focused on the potential role of I-Smads in generating transient responses while Vilar et al. (2006) focused on the receptor dynamics to explain the occurrence of both transient and sustained responses. Zi et al. (2007)  included a simple model of the Smad dynamics and highlighted the importance of the balance between clathrin-dependent endocytosis and non-clathrin mediated endocytosis. All pathway elements were finally brought together by Chung et al. (2009)  in a more comprehensive model, used to examine the contradictory roles of TGF-β in cancer progression. Lately Zi et al. (2011)  published a study that highlights the potential of TGF-β ligand depletion in converting short-term graded signaling responses into long-term switch-like responses. Unlike for other pathways oscillations have not yet reported for the TGF-β signaling pathway [26, 27]. TGF-β type ligands are also acting as morphogens, and the response to these appears to be proportional. Recently, Paulsen and co-workers published a study on the impact of synexpression of the feedback inhibitors BAMBI, Smad6, and Smad7 on the read-out of morphogen gradients during embryogenesis .
While the many published studies explain the different behaviours for the different situations for which they are observed and highlight the many mechanisms that enable the different response types it remains largely unclear how easily the response type can be changed. We wondered how the TGF-β signaling pathway accomplishes the flexibility in its responses and which and how many parameters have to be altered for cells to respond differently. To efficiently explore the canonical response we focused on the core signaling architecture, and did not consider the detailed receptor dynamics and cross-talks in the model; they are included indirectly through the parameters that they modulate. We explored the response types and in particular changes in the response type as we explored the parameter values within biologically meaningful ranges. We find that relatively small changes in single parameters can alter the response. Cellular protein concentrations are a particular powerful point of control and this explains how different cell types can show different responses. Importantly we also identify key parameters that affect the response and we can relate these to observed points of cross-talk between signaling pathways. The particular architecture of the TGF-β network thus allows for the great flexibility in the response.
Several models for the TGF-β signaling network have been developed that focus on different aspects of the TGF-β signaling network, i.e. the receptor dynamics [23, 29], the shuttling between the cytoplasm and the nucleus , and the negative feedback via the I-Smad (Smad7/Dad) . These different aspects have lately been combined in a model that addresses differences in TGF-β signaling between normal and cancerous cells . The models of the TGF-β signaling pathway showed that stimulation could result in either transient and sustained responses dependent on the choice of parameters [3, 22–25, 29]. Transient responses could be obtained through complex receptor dynamic , the I-Smad-mediated negative feedback [3, 22], or ligand depletion . Negative feedbacks can in principle also give rise to oscillatory behaviour. We wondered whether all three qualitative behaviours (sustained, transient, or oscillatory response) could be obtained already with the most simple intracellular feedback mechanism, and how these behaviours would depend on the parameters. Since the more complex interactions (that we ignore) effectively modulate the parameter values in our model an in-depth understanding of the parameter dependencies in the simple model should also enable a better understanding of the complex network interactions that are found in the cell. The different response types can also (trivially) be obtained by modulating the protein concentrations accordingly. We, however, keep the concentrations of receptors, ligand, R-Smad and Co-Smad constant and thus include these effects only indirectly as changes in the effective binding rates.
We translated those interactions into sets of ODEs using the law of mass action where appropriate. To reduce the complexity of the model we also employed Hill functions to describe the regulation by cooperative interactions. To efficiently investigate the impact of changes in total concentration of receptors, R-Smad, and Co-Smad we used a total concentration rather than production and degradation rates for these species.
To respond to TGF-β cells must be able to detect changes in the ligand concentration and convert the differences into different transcriptional responses. Transcriptional activity is determined by the concentration of transcription factors in the nucleus. We therefore monitor the nuclear concentration of R-Smad/Co-Smad complexes as a measure of transcriptional activity, in response to a change in the extracellular TGF-β concentration.
We are interested in the signaling capacity of the TGF-β pathway within its physiological limits. These physiological limits are set by the plausible range that the parameter values can take. We established a likely range for each parameter value based on available data and estimates (Additional file 2, Table S1 and Table S2). While previous measurements and estimates are necessarily of limited accuracy and differences are likely to exist between different cells and different cell types [3, 22–24, 29] we expect that basing ourselves on the available data will not too much distort the ranges that we screen. Most parameters were varied over 3 or 4 orders of magnitude, centered around the mean of values found in the literature. Since there are no good estimates for the I-Smad expression rates k14 and k15 were varied over 5 orders of magnitude. The rates of phosphorylation and dephosphorylation (k7 and k13) were varied only over two orders of magnitude because a large fraction of the simulations failed when these rate constants were varied over a wider range. To avoid a bias to the few parameter sets that do not lead to extreme dynamics we had to constrain these two parameters to only vary over two orders of magnitude. To determine the possible range of pathway responses to a defined stimulus, we carried out 106 independent simulations with parameter values randomly picked from a uniform logarithmic distribution of parameter values within the set ranges (as discussed in Geier et al. ) and compared the predicted nuclear concentration of R-Smad/Co-Smad complexes in response to the ligand stimulus. In a first step, we let the system equilibrate for 1 hour with almost no ligand (concentration of 10-6 pM to avoid failure of the solver) and initial cellular concentrations TGFβR = 1 nM, Smad = 60 nM and Co-Smad = 100 nM. We then used the steady-state value of the first step and solved the simulations for 10 hours with a constant ligand concentration of 200 pM. Using MATLAB's ode15 s routine the 106 simulations took in total approximately 140 hours of CPU time.
After the initial peak the response must retain at least 90% of its maximal value (called Opeak). To exclude slowly increasing responses we further require that 90% of the peak value Opeak is reached within less than 7200 s (2 hours).
After the initial peak the response must drop to levels lower than 10% of its peak value Opeak within less than two hours and the final value (after 10 hours) (called Oend) must be lower than 0.1 × θ.
After the initial peak the amplitude (difference between the local maximum and local minimum) must exceed 0.1 × θ at least 4 times.
The fifth amplitude must be less than half the second amplitude.
The fifth amplitude must be higher than half the second amplitude.
We characterized the long-term behaviour of oscillations based on the relative amplitudes of the second and fifth peak because the first peak can be particularly high (Figure 2D), and most dampened simulations have no more than five peaks.
Quantitative data on the physiological concentrations of the cellular proteins and the transcription factor complex (nuclear Smad/CoSmad complex) do not exist, and we therefore had to set our detection threshold arbitrarily to θ = 10 pM when analysing a unique constant stimulus with 200 pM TGF-β ligand. When the response to several ligand concentrations or with several protein concentrations was studied we used the maximal response value as θ. Simulations were run for 10 hours. In case if oscillations, if the amplitude of oscillations was still larger than 0.1 × θ after 10 hours, the simulation was continued until the oscillations vanished, but for a maximum 100 hours. In this way we avoid any impact of period length on the classification of oscillations, and the length of the period indeed does not bias our characterisation of oscillations to dampened or sustained oscillatory behaviour (Additional file 3, Fig. S5C). The time thresholds 2 hours and 10 hours were chosen based on experimental data .
In response to a sustained stimulus (200 pM ligand) our simple model for TGF-β signaling can give rise to sustained (Figure 2A), transient (Figure 2B), or oscillatory (Figure 2C, D) responses. The sustained/transient distinction is particularly relevant, as it has been shown that those two qualitative behaviors are related to the growth inhibitory effect of TGF-β . To better understand the conditions for these different behaviours we sought to identify parameter families that would give rise to a certain response type. We hoped that a comparison of those families would reveal the critical parameters that determine the response type. To that end we screened a large number of parameter sets and classified them according to their responses as described in detail in the Materials and Methods section.
Our criteria in Figure 2 are very strict (i.e. the speed of responses is an important criterium) and there is a wide undefined range between sustained and transient responses. As a consequence most parameter sets (57.5%) do not fall into any of the defined categories (Additional file 3, Fig. S5A, black sets). Of those that can be classified most (25.5% of the parameter sets tested) led to no response (Additional file 3, Fig. S5A, marked in grey). Among the "responsive" parameter sets most lead to sustained responses (14.8% of the parameter sets tested) (Additional file 3, Fig. S5A, marked in red) while transient responses are observed less frequently (2.2%of the parameter sets tested) (Additional file 3, Fig. S5A, marked in green). All three behaviors have previously been observed in various models of TGF-β signaling. We find that in addition in a minority of cases (306 simulations, i.e 0.046%) also oscillatory responses can be produced (Additional file 3, Fig. S5B). Even though the number of sets that give rise to oscillations in the concentration of nuclear transcription factor complexes is small, these may occupy a sufficiently dense subspace in the parameter space to be physiologically relevant. The oscillations can either be sustained or dampened, depending on how fast their amplitude decays (Figure 2C, D). As expected sustained oscillations have a larger number of peaks (Additional file 3, Fig. S5B). While the period of the oscillations is not biased to dampened or sustained oscillatory behaviour (Additional file 3, Fig. S5C), the duration (the time until oscillations vanish) depends on both the number of peaks and the duration that tends to be higher for sustained oscillations (Additional file 3, Fig. S5D).
Oscillatory behavior has been reported for a number of other signaling pathways (i.e. the ERK cascade ), but so far no experimental evidence exists for oscillations in the TGF-β pathway. However, standard biochemical experiments average over a large number of non-synchronized cells. If the nuclear concentration of transcription factor indeed oscillated, only sophisticated single-cell assays would reveal these.
Transient and oscillatory responses are similar in that the response must decay quickly in spite of the continuous presence of ligand. A similar comparison of the parameter ranges that permit transient (green) or oscillatory (blue and magenta) responses (Figure 3B) indeed reveals that similar restrictions apply (i.e. large shuttling rate between cytoplasm and nucleus, k8, and strong negative feedback, k14, k15, k17, k18, and k19). However, in case of oscillations the response restarts and in addition we indeed notice a strong restriction of the rate of ligand-receptor binding k2 in case of oscillatory responses. Figure 4B shows the clear segregation of the parameter sets that give rise to "transient" (green), dampened (blue) or sustained (purple) oscillatory responses in a plane spanned by the receptor-ligand binding rate k2 and the speed of I-Smad turn-over (k16 × k17 × k19). Oscillations are observed only when k2 is small such that ligand binds slowly to its receptor and the pool of free receptor is depleted gradually (Additional file 3, Fig. S6A). As a consequence free receptor is still available when I-Smad has downregulated the response and ligand can still trigger a further response. Conversely, if k2 is large, receptors are rapidly bound to the ligand (whose concentration is constant in this model), and once the response has terminated, there is no free receptor available to induce a new response (Additional file 3, Fig. S6B). k2, and thus the speed with which the free receptor concentration decreases, critically determine the dampening of oscillations. A 10-fold change in the value of k2 can transform sustained oscillations in highly dampened ones (compare panels A and B in Additional file 3, Fig. S6). Rapid degradation of I-Smads is important for sustained oscillations because otherwise all receptors become rapidly sequestered by I-Smad and the response is terminated. Accordingly inclusion of receptor endocytosis and recycling to the membrane (cycling) combined with the removal of the I-Smad would allow further oscillatory cycles.
Each parameter in our simple model integrates the effects of many further interactions as may also arise from cross-talk. Thus it has been shown that the phosphorylation of Smad in its linker region by Ras-activated MAPK induces a cytoplasmic retention of R-Smads , which in our system would be represented by a lower shuttling rate into the nucleus (k8). Interestingly k8 indeed strongly influenced the response type. Another parameter that appears to be important in determining the response characteristics is the binding rate of TGF-β to the receptor (k2). Our description of the processes at the cell membrane is very simple and thus k2 has also to take in account the regulation of TGF-β outside the cell by soluble sequestering factors and membrane-bound co-receptors as well as processes that affect the receptor density on the cell membrane. Those auxiliary factors play therefore a crucial role in the TGF-β pathway flexibility. We should stress that all parts of the parameter space should be readily reachable for the cell and small adjustments in the parameter values should thus be sufficient to alter the response type.
Finally we wondered how different ligand concentrations would affect the cellular response. The impact of different TGF-β concentrations have already been studied by Clarke et al.  and Zi et al. , but there the results were strongly affected by ligand depletion since the TGF-β concentrations were allowed to go down over time because of internalization and degradation. We were interested how different, but constant, stimuli would affect the response - the effect of ligand depletion can then be deduced as response to decreasing ligand concentrations. As we varied the ligand concentration between 0.2 pM and 20 nM we noticed that only for a very small fraction of parameter sets (less than 0.5%) the response type changed qualitatively as the concentration varied. The parameter sets were not clustered and a further increase by restricting parameter ranges (as in case of the cellular protein concentrations) could not be achieved. Even though changes in the TGF-β concentration cannot switch the response type in our simulations, the duration of the response increases with increasing TGF-β concentrations as previously observed by Zi et al. . This increase was, however, insufficient to alter the response type according to our definitions.
For parameter sets that give rise to oscillatory responses, changing the input strength and shape does not influence the period of oscillation but modulates the evolution of the oscillations amplitudes (data not shown). When exposed to sustained, high TGF-β concentrations the amplitude of oscillations starts to decay from the beginning. When the TGF-β concentration raises progressively, the amplitude of oscillation first raises and then decays, reflecting two competing phenomena: the amplitude of oscillations tends to be proportional to the input, but at the same time the sequestration of the receptor by the inhibitor leads to a dampening of the amplitude.
We next investigated in how far the kinetic parameters can influence the saturation concentration (Figure 6B, C) and the maximal output value at saturation (Additional file 3, Fig. S8 and Additional file 3, Fig. S9). For transient responses it is mainly the rate of ligand-receptor binding, k2, that determines the saturation concentration (Figure 6D and Additional file 3, Fig. S10). In case of slow binding higher concentrations of ligand are required to saturate the receptors. The saturation concentration for sustained responses are determined both by the receptor-ligand binding rate, k2, and by the cytoplasm-nucleus shuttling rate, k8 (Figure 6E and Additional file 3, Fig. S10). Fast shuttling enables more rapid deactivation of Smads as based on observations by Hill and coworkers  dephosphorylation is restricted to the nucleus in our model. As discussed above k2 and k8 have both been reported to be modulated by other processes. The saturation concentration can therefore also be adjusted by cross-talk.
The different saturation concentrations are likely important for the TGF-β response as different genes can be activated or repressed depending on the nuclear Smad complex concentration. While the mechanism by which different concentrations of the nuclear transcription factor complex translate into different transcriptional responses has not been resolved, likely mechanisms include promotor selection based on differences in the promoter binding-site affinities, cross-repression, and the establishment of a reciprocal of repressor gene expression [36, 37].
A response is faithful if the output is proportional to the input over time, i.e. y output (t) = α × y input (t), where α is the proportionality coefficient. This requires (i) that the output adapts rapidly to changes in the input, and (ii) that the response does not saturate, i.e. max(y output (t)) < max(O peak ), which is the case if the proportionality coefficient α is low and/or the maximal response value max(O peak ) is high. Those requirements are reflected in the constraints on the parameter values (Figure 7B-C) for faithful responses, i.e. a low binding rate of TGF-β to its receptor and a low phosphorylation rate prevent early saturation of the output, while a relative weak feedback and a low binding rate of the I-Smad to the receptor prevent a premature termination of the response. We have previously discussed the regulation of the binding rate of TGF-β to its receptor, k2 and thus now focus on the feedback. The different I-Smads have been shown to vary in their effects. Thus Dad, the Drosophila I-Smad, appears to interfere mainly with the BMP-like pathways (Tkv and Sax receptor dependent pathways) but not the Activin-like Babo-dependent pathway . Inhibition by vertebrate Smad6 and Smad7 can be achieved by sequestration, enhanced degradation, or an impact on phosphorylation. The different processes likely have different efficiencies and this will determine the efficiency of the negative feedback.
Our results indicate that under certain parameter restrictions the extracellular concentration is directly reflected in the output concentration. In that case, TGF-β can act as a morphogen, conveying positional information and determining cell-fate, subjected to the set of activated and repressed genes.
The duration of the signaling response is thought to be an important factor influencing the cell's phenotypic response to TGF-β. We have employed a very simple model of the TGF-β network to better understand the mechanistic basis of the observed signaling plasticity. We find that the qualitative response (transient, sustained, oscillations, proportional responses) to a constant ligand exposure can indeed be changed by altering the value of a single parameter value. Since we consider a simple model each parameter value represents a wider range of processes and our observation thus implies that both changes in protein concentration as well as cross-talk between signaling pathways can alter the qualitative response to a TGF-β stimulus. Many more complicated models for TGF-β signaling as well as for other signaling networks have been proposed already. To better understand the regulatory impact of cross-talk it will be important to connect experimentally validated models for the TGF-β network also to those for other pathway models. While many kinetic parameters have been measured an important parameter that remains often unmeasured is the protein concentrations. To better predict the responses in different cell types it will be important to obtain quantitative information on protein abundance in different cell types - and eventually in individual cells.
This research was supported by an iPhD SystemsX grant to Georgios Fengos.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.