Search tips
Search criteria 


Logo of ploscompComputational BiologyView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS Comput Biol. 2006 January; 2(1): e3.
Published online 2006 January 27. Prepublished online 2005 December 9. doi:  10.1371/journal.pcbi.0020003
PMCID: PMC1356091

Signal Processing in the TGF-β Superfamily Ligand-Receptor Network

Adam Arkin, Editor


The TGF-β pathway plays a central role in tissue homeostasis and morphogenesis. It transduces a variety of extracellular signals into intracellular transcriptional responses that control a plethora of cellular processes, including cell growth, apoptosis, and differentiation. We use computational modeling to show that coupling of signaling with receptor trafficking results in a highly versatile signal-processing unit, able to sense by itself absolute levels of ligand, temporal changes in ligand concentration, and ratios of multiple ligands. This coupling controls whether the response of the receptor module is transient or permanent and whether or not different signaling channels behave independently of each other. Our computational approach unifies seemingly disparate experimental observations and suggests specific changes in receptor trafficking patterns that can lead to phenotypes that favor tumor progression.


Sensing of the environment by cells relies extensively on receptors that bind extracellular molecules and trigger intracellular responses. The TGF-β pathway transduces a broad range of extracellular signals into transcriptional responses that affect many cellular processes, including cell growth, apoptosis, differentiation, homeostasis, and morphogenesis. It is used, for instance, to control the precise patterns and forms that arise during development, and its malfunction contributes to a wide variety of diseases and developmental disorders. Here the authors develop a concise computational model of the TGF-β pathway and show that the first layer of communication with the environment, the ligand-receptor network, is not merely a passive transducer of signals but rather embeds properties that makes it a signal processing unit. Receptors traffic between different cellular compartments from which they signal distinctly, leading to an unexpected richness of types of behavior that is not apparent from the simplicity of the typical cartoon representations of this pathway. At the receptor level, the system can select among different functioning modes to sense absolute levels of ligand, temporal changes in ligand concentration, and ratios of multiple ligands. This extra level of regulation can explain a wide variety of phenomena and leads to a unified interpretation of seemingly disparate experimental observations.


The TGF-β signal transduction pathway follows an apparently straightforward downstream cascade, progressing sequentially from the interaction of ligands with transmembrane receptors, through phosphorylation of mediator Smad proteins, to transcriptional responses (Figure 1). The simple logic of this signal transduction cascade strongly contrasts with the molecular complexity of the cellular processes involved and the wide diversity of responses triggered.

Figure 1
Formation of Receptor Hetero-Tetramers

At the molecular level, there is an intricate signal transduction machinery that integrates signals from the 42 known ligands of the TGF-β superfamily, funnels them through the two principal regulatory Smad (R-Smad) channels (Smad1/5/8 or Smad2/3), and subsequently leads to the widespread transcriptional control of more than 300 target genes in a cell-context dependent manner [1] (see Figure 2). The components of this machinery include the members of the two main receptor families (type I and type II receptors), a myriad of adaptor proteins, and the trafficking apparatus of the cell, which shuttles proteins between different subcellular compartments. Each ligand induces the formation of a receptor complex with type I and type II receptors, which then signal through one of the two Smad channels [2,3]. The ability of most ligands to bind several type I and type II receptors results in a complex ligand-receptor interaction network (Figure 2).

Figure 2
Interactions among the Ligands of the TGF-β Superfamily and Their Receptors

At the phenotypic level, the responses are extremely diverse. The members of the TGF-β superfamily act prototypically as potent negative growth regulators, but, depending on the cell type and context, they can also induce differentiation, apoptosis, cell migration, adhesion, and extracellular matrix deposition. TGF-β itself is of particular interest in cancer research. In epithelial cells, it suppresses cellular growth by inducing G1 arrest (mediated by transcriptional activation of p15 and p21) [4], and its inactivation contributes to tumorigenesis. The versatility of the pathway in eliciting different types of behavior is perhaps best epitomized by the pervasive, rather paradoxical ability of TGF-β to change its function from suppressor to promoter of growth in epithelial cells during tumor progression [4,5].

Current theories for explaining the variety of responses to members of the TGF-β superfamily of ligands focus mainly on the downstream transcriptional regulatory networks they activate: transcriptional cofactors of the R-Smads are expressed at different levels in a cell-specific manner, thereby modifying downstream responses. In fact, the role reversal of TGF-β from negative to positive growth regulator has been found to be associated with a phenotypic change known as epithelial-to-mesenchymal transition, in which cells change the cofactors recruited by the R-Smads and acquire motile phenotypes [57].

It is striking, however, that such a variety of complex responses and intricate molecular components are connected through just two Smad channels by such a simple downstream signal transduction cascade. There is a richness of experimental observations that are difficult to reconcile with this observation. In particular, whether TGF-β acts as a growth suppressor or promoter can depend on whether the tumor cells were grown in vitro or in vivo [8]. In these two different situations, the extracellular context determines the way in which cells respond to TGF-β. It has been suggested that TGF-β can suppress the growth of cells around the tumor, that it can shut down locally the immune system, and that it can promote angiogenesis. All these paracrine effects would help the growth of the tumor in vivo, where it has to compete with neighboring cells. So far, although appealing, none of these mechanisms has been identified as an alternative cause of the TGF-β role reversal.

The most direct way in which the extracellular context can affect the functioning of the TGF-β pathway is through signaling of other ligands of the TGF-β superfamily. As we have mentioned, ligands and receptors form a complex interaction network, where multiple ligands share receptors, potentially coupling their signaling. All these interactions are in turn coupled to receptor trafficking, which is known to be a mechanism that regulates signal transduction [9,10]. Trafficking has been investigated in detail in many signal transduction pathways, such as the epidermal growth factor receptor (EGFR) and G protein-coupled receptor (GPCR) pathways [1113]. The typical way in which trafficking and signaling are coupled is by the induction of receptor internalization upon ligand binding and receptor activation, as for instance in the EGFR and GPCR pathways. After internalization, receptors can activate other signaling pathways, be modified in specific ways, and be targeted for degradation or recycling back to the plasma membrane.

A peculiarity of the TGF-β pathway is that receptors are constitutively internalized, even in the absence of ligand [14,15]. The trafficking route that the receptors follow, however, depends on whether or not they are in a signaling complex (Figure 3). Different routes will trigger different signaling outcomes and affect how receptors are degraded. Therefore, although the explicit implementation of the coupling is different than in the EGFR and GPCR pathways, receptor trafficking and signaling are also tightly coupled in the TGF-β pathway.

Figure 3
Signaling and Trafficking in the TGF-β Pathway

Taking all the current experimental data together, it is clear that many details of the underlying processes remain largely unknown. Recent experiments [14,15], however, provide key information that allows computational modeling to bridge the gap between potential molecular mechanisms and experimentally observable behavior. The TGF-β pathway is at a discovery stage where modeling can help to gain “functional” intuition.

Here we characterize computationally the diverse potential types of behavior that the pathway structure itself can confer on the system. The types of behavior include responses to persistent changes in ligand concentration that can be transient or sustained and simultaneous responses to multiple ligands that can be passed downstream independently of or dependently on each other.

A sustained response implies that the steady-state signaling activity is a function of the ligand concentration. In this case, the higher the ligand concentration, the higher the activity of the pathway. For a transient response that precisely returns to the prestimulus level, in contrast, the steady-state activity is always the same and the pathway can only detect changes in ligand concentration. When multiple ligands signal in a dependent fashion, the extent of the coupling can be such that one ligand can suppress the effects of another one. In this regime, the pathway does not detect ligand concentrations but ratios of concentrations.

As we show here, all these types of behavior can be present in the TGF-β pathway. Which specific one is selected is determined by the interplay between trafficking and signaling. Thus the pathway can be set to detect, at the receptor level, absolute levels of ligand, temporal changes in ligand concentration, and ratios of multiple ligands. Such flexibility in the pathway behavior can lead to diverse physiological outcomes that have been associated with facilitated tumor progression.


A Concise Computational Model

In order to study the signal processing potential of the ligand-receptor network coupled with receptor trafficking, we assemble all the essential elements into a concise mathematical model that captures the logic of the underlying processes. The main goal is to represent as much complexity as possible through a small number of quantities that have direct experimental interpretation.

The essential elements we consider are (Figure 4): (i) Ligands induce the formation of receptor complexes with type II and type I receptors. (ii) Receptors and ligand-receptor complexes can be present in two spatially distinct compartments: plasma membrane and internalized endosomes. (iii) The signaling activity is proportional to the number of ligand-receptor complexes that are present in the internalized endosomes [16]. (iv) Receptors and ligand-receptor complexes are continuously internalized into endosomes and recycled back to the plasma membrane [14, 15]. (v) Receptor degradation has a constitutive contribution, which is the same for free receptors and ligand-receptor complexes [14,15]. (vi) Receptor degradation has a ligand-induced contribution, which affects only receptors that have been complexed with ligands [14,15].

Figure 4
Two-Compartment Model of Receptor Trafficking and Signaling

We use these elements to develop a mathematical model based on rate equations that describe the dynamics of both how different molecular species transform into each other and how they traffic between the different cellular compartments. We assume that internalization, recycling, and degradation rates are proportional to the number of receptors or ligand-receptor complexes; and that the formation of ligand-receptor complexes is proportional to the ligand concentration and to the type I and type II receptor concentrations.

In a first step toward characterizing the effects of the coupling of signaling with receptor trafficking, we consider that only a single type of ligand is present. Explicitly, we study how the components of the canonical TGF-β pathway—one ligand (TGF-β) and two receptors (Alk5 and TGFβRII), as emphasized in Figure 2—respond to changes in ligand concentration. The mathematical equations are shown in Figure 4. For typical trafficking rates (see Materials and Methods), this model closely reproduces the typical time courses of Smad phosphorylation upon addition of ligand (Figure 5A and and55B).

Figure 5
Typical Time Courses of the Number of Active Receptor Complexes upon Addition of TGF-β

The computational model can be used also to analyze how different parameters affect the behavior of the system. For instance, the time during which the signaling activity rises is related to the time required for internalization and recycling of the receptors. Thus, the signal will peak, or stop rising, at about 30–60 min after ligand addition. If the internalization and recycling rates are changed, the position of the peak changes accordingly (Figure 5C and and55D).

Other transmembrane receptor pathways, such as the EGFR pathway, have much faster kinetics; the EGFR pathway reaches peak activity as fast as 5 min after stimulation [17]. The main reason for these differences is that most of the EGF receptors are present in the plasma membrane and they are ready to signal upon the addition of ligand. A similar kinetics is also observed for many GPCRs. In the TGF-β pathway, internalization occurs continuously, and only about 5%–10% of the receptors are present in the plasma membrane at a given time [15]. The remaining 90%–95% of the receptors are internalized in endosomes. Receptors need to be recycled from the endosomes back to the plasma membrane in order to be able to interact with the ligand; and this process takes about 30 min on average.

We can use the computational model to study the effects of different mechanisms, such as different forms of receptor degradation, on the behavior of the system. It has been observed that the addition of ligand can stimulate the degradation of the receptors in two different ways. On the one hand, receptors complexed with ligand in the plasma membrane can be internalized through a lipid raft–caveolar degradation pathway without becoming active signalers [15]. The behavior obtained in this case is the one we have already described (Figure 5A). On the other hand, ligand-receptor complexes can follow the standard internalization clathrin pathway, signal, and then be targeted for degradation upon returning to the plasma membrane [14]. Interestingly, when this mechanism is explicitly modeled, the behavior obtained (Figure 5E) is qualitatively the same as the previous one (Figure 5A). Likewise, when both degradation mechanisms are considered together, a similar type of behavior is also obtained (Figure 5F). Thus, at this level of detail and for this range of parameter values, different mechanisms that implement ligand-induced degradation can lead to similar behavior.

Control of the Signal: Transient versus Permanent Responses

How is it possible to modify the form in which the system responds to changes in TGF-β concentration? A mathematical analysis of the model (see below and Materials and Methods for details) indicates that the key quantity that determines the qualitative behavior of the pathway is the ratio of the constitutive to the ligand-induced rate of degradation, referred to, in short, as the constitutive-to-induced degradation ratio (CIR). This quantity compares the rates of two degradation processes and, in general, does not have a simple expression in terms of rate constants.

Depending on the CIR, a permanent change in ligand concentration can elicit responses between two extremes (Figure 6). For low CIR, the ligand-induced degradation process dominates and there is a transient increase in signaling activity that returns to pre-stimulus levels (Figure 6A and and6B).6B). For high CIR, the constitutive degradation process dominates and there is a permanently elevated level of signaling activity (Figure 6E and and66F).

Figure 6
Control of the Kinetic Signaling Behavior

For intermediate CIR, the behavior of the system is a mixture of both limiting types of behavior, with transient and permanent components (Figure 6C and and6D).6D). The precise parameter values influence the amplitude and characteristic time of the response (see for instance Figure 5C and and5D),5D), but its qualitative shape, that is, whether the response is transient (Figure 6A and and6B)6B) or permanent (Figure 6E and and6F),6F), depends only on the CIR.

The intuitive explanation of such types of behavior is as follows (for a detailed mathematical analysis, see Materials and Methods). The probability for a receptor to bind the ligand, and therefore to become active, increases with the ligand concentration. If the ligand does not induce the degradation of receptors, the number of receptors remains constant and the total activity increases when the ligand concentration increases. If the ligand induces the degradation of the receptors, the number of receptors starts to decrease after ligand addition, which will eventually attenuate the signal. At steady state, the production and degradation of receptors equal each other. In the limit of the CIR going to 0, the signal adapts completely because degradation is proportional to the activation of receptors, and therefore activation is also proportional to the production of receptors. Thus, it is the receptor production rate, not the ligand concentration, that determines the steady-state signaling activity.

There are clear examples in other signal transduction pathways that show that these two limiting types of behavior can potentially lead to different physiological outcomes. For instance, transient activation of the MAPK cascade by EGF leads to cell proliferation. In contrast, permanent activation of the MAPK cascade by NGF leads to cell differentiation. In both cases, activation of the MAPK cascade induces the expression of a negative regulator that shuts down the activity of this cascade. The differences between EGF and NGF have been attributed to additional pathways activated by NGF that can prevent the inactivation of the MAPK cascade [18]. Our model shows that such transient and permanent types of behavior can also be achieved by just changing the trafficking patterns, in particular by adjusting the CIR, without the need for explicitly expressing a negative regulator to shut down the cascade after signaling.

Remarkably, the duration of the signaling activity also seems to affect the physiological outcomes triggered by TGF-β [19]. Epithelial cells that are sensitive to the antiproliferative effects of TGF-β (HaCaT and Colo-357) have sustained activity of more than 6 h. In contrast, pancreatic tumor cell lines (PT45 and Panc-1) show short transient activity of about 1–2 h. Such a short transient confers resistance to the antiproliferative effects of TGF-β but maintains other responses to TGF-β that can lead to increased malignancy and invasiveness [19]. In our model, those differentiated types of behavior arise naturally for different trafficking patterns. In particular, short transients and sustained responses imply a low and a high CIR, respectively.

Control of the Signals: Coupled versus Uncoupled Channels

In vivo conditions, in contrast to those typical of in vitro experiments, expose cells to complex environments with many different growth factors. When multiple ligands of the TGF-β superfamily are present at the same time, they are likely to affect each other's signaling (Figure 2). To study how multiple simultaneous input signals are integrated into coordinated transcriptional responses, we extend our computational model to consider two ligands that signal through two different type II receptors and a shared common type I receptor (see Materials and Methods for the mathematical equations). This is the simplest case of signal integration.

Intuitively, one should expect signals to be coupled when the shared receptor is saturated with ligands and uncoupled when ligand concentrations are low. At saturation, increasing the concentration of one ligand, and thus the concentration of the corresponding ligand-receptor complex, will take the shared receptor away from the complex formed by the other ligand, thus decreasing its signaling.

A mathematical analysis of the model (see Materials and Methods) indicates that even when the receptors are far from ligand-saturating conditions, it is possible for signals to affect each other. The key element is again receptor trafficking. In essence, the coupling arises because the induction of degradation of the common receptor by one ligand attenuates the effects of the other ligand, which also requires the common receptor to transduce the signal.

For pathways working away from receptor saturation, the interplay between trafficking and signaling determines how multiple simultaneous signals are passed downstream. As in the single ligand case, there are two extreme types of behavior (Figure 7): For low CIR, the ligand-induced degradation process dominates and signals are completely coupled (Figure 7A and and7B).7B). For high CIR, the constitutive degradation process dominates and signaling is uncoupled (Figure 7E and and77F).

Figure 7
Control of Signal Integration

When the signals are completely coupled, the steady-state number of all ligand-receptor complexes remains constant and is independent of the ligand concentration. In this case, increasing one signal will decrease the other one by the same amount. When signals are uncoupled, the numbers of each species of ligand-receptor complexes change independently of each other. In general, for intermediate CIR, signals will show some degree of coupling (Figure 7C and and7D).7D). These results demonstrate that changes in trafficking patterns, and the corresponding degradation, can alter the way in which the pathway integrates multiple, simultaneous signals.

The completely coupled case is especially interesting because it indicates that one ligand can potentially inhibit the effects of another one. Ligand-induced degradation is thus not only a mechanism for achieving transient responses, but also for coupling multiple signals. The fact that TGF-β can signal not only via Alk5 but also via Alk1 and Alk2 (Figure 2) potentially couples TGF-β signals to those of Activin A, BMP 6, BMP 7, and MIS. Thus, if TGF-β loses its growth suppressor properties, it could promote growth by inhibiting other growth suppressor pathways. For instance, there are dominant negative TGFβRII mutants that when overexpressed attenuate the response to TGF-β [20]. The presence of any of these mutant receptors and TGF-β results in the formation of futile receptor complexes that can target receptors for degradation, or take receptors away, which otherwise would be available to transduce the signals of other members of the TGF-β superfamily.

Context-Dependent Response to TGF-β

The role reversal of TGF-β from negative to positive growth regulator is a widespread feature of tumor progression and is often associated with endogenous overexpression of TGF-β. As we have mentioned in the introduction, it is associated in some situations with the epithelial-to-mesenchymal transition [6,7]. Under these conditions, the transcriptional program of tumor cells changes so that the Smad-activated genes promote rather than repress growth.

In other situations such a transition does not seem to be present. It has been observed in breast, prostate, and colon cancer cell lines that the action of TGF-β as growth promoter or suppressor depends on whether the cells were grown in an in vitro environment or in vivo in mouse xenografts [8,21,22]. The reasons for such a change remain largely unknown. It has been speculated that it could be a consequence of the effects of TGF-β on the in vivo microenvironment of the tumor cells. Another possibility is that other growth factors, such as EGF, affect how TGF-β is ultimately coupled to the cell cycle. Our model explicitly shows that the role reversal is a potentially intrinsic consequence of the design of the ligand-receptor interaction network and trafficking machineries and that it could be the result of TGF-β attenuating the effects of growth-suppressing signals of other members of the TGF-β superfamily that might be present in the in vivo cell environment.

Simultaneous Perfect Adaptation and Coupled Signaling

Our model also indicates that the conditions that give rise to completely coupled integration of multiple signals are the same that, in a single-ligand system, cause the signaling activity to completely adapt to its prestimulus level. Remarkably, this concurs with observations in prostate cancer cell lines, which show that the in vivo context can not only make TGF-β a growth promoter but also that the in vitro response to TGF-β is transient [8]. One should expect the extracellular environment of growth factors to be more complex in vivo than in vitro. This relationship between in vivo and in vitro behavior and its connection to the coupling between receptor trafficking and signaling underscores the importance of understanding how signal transduction pathways are embedded within the cellular microenvironment under physiologically relevant conditions. Not only mutations in the canonical pathway but also changes in trafficking patterns can move the pathway to a different functioning point.

The qualitative results of our model, such as the regimes leading to transient and permanent responses as well as to completely coupled and uncoupled modes of signal integration, do not depend on the details of the model but on general properties (see Materials and Methods). Thus the main ideas are also relevant to other signal transduction pathways that are coupled to receptor trafficking. In particular, revisiting the experimental data, one can see that the interplay between adaptation and signal integration (Materials and Methods) is also present in the EGFR pathway (see Figures 4 and and55 of reference [23]), in which down-regulation of erbB-2 by EGF concurs with adaptation of the signal transmitted by EGFR.


Cellular functions are controlled by networks of interacting molecules that operate at different levels of organization [2426]. Here, we have developed a concise computational model of the TGF-β pathway that shows that the receptors for the TGF-β superfamily of ligands are not just passive signal transducers. They are organized in a network that is able to process the signals before passing them downstream. Changes in receptor trafficking patterns can modify the type of behavior of the pathway in response to single and multiple ligand inputs. Already at the receptor level, the pathway can detect absolute levels of ligands, temporal changes in ligand concentration, and ratios of multiple ligands. This extra level of regulation can explain a wide variety of phenomena, such as the counterintuitive role reversal of TGF-β from suppressor to promoter of growth, and leads to a unified interpretation of seemingly disparate experimental observations. A key quantity that determines the qualitative behavior of the pathway is the CIR of the receptors. For low CIR, the pathway responds transiently to sustained changes in ligand concentration, and the signaling activities of multiple simultaneous ligands become dependent on each other. Ligand-induced degradation is thus not only a mechanism for achieving transient responses and perfect adaptation, but also for coupling multiple signals.

Various experiments can be designed to test the predicted types of behavior. The most direct evidence would come from biochemical measurements of ligand-induced and constitutive degradation rates [14,15]. The observed degradation rates can then be related to measurements of the levels of signaling activity. For instance, decreasing ligand-induced degradation by blocking the lipid raft–caveolar pathway with nystatin, as in reference [15], should bias the system behavior toward permanent responses to step changes in ligand concentration. Molecular interventions, such as RNAi against the mRNA of proteins involved in the trafficking processes, and their effects on signaling activity would provide indirect evidence that can be related to the details of the model. Because of their well-established role in intracellular trafficking, small GTPase rab proteins [27] are prime candidates for this type of approach. Certain experimental observations [8,19] in conjunction with our model suggest that cells with phenotypes that favor tumor progression have low CIR. Further experiments could test the extent of this correlation by comparing trafficking rates and signaling activity between different cancerous and non-cancerous cell lines.

Materials and Methods

Internalization rate.

It has been reported in Table 1 of reference [15] that after 15 min of labeling the receptors at the plasma membrane, only 8%, 6%, 4%, and 2% of the labeled receptors remain at the plasma membrane (the different percentage values correspond to different experimental conditions). The remaining labeled receptors have been internalized in either caveolin-1 positive or caveolin-1 negative vesicles. By using the formula An external file that holds a picture, illustration, etc.
Object name is pcbi.0020003.ex001.jpg , where ki is the internalization rate and ft the fraction of labeled receptors that remain at the plasma membrane after a time t, we obtain internalization rates of 1/5.9, 1/5.3, 1/4.7, and 1/3.8 min−1, respectively. We have chosen ki = 1/(3 min) for comparison with experimental data in Figure 5B.

Figure 4 of reference [14] shows that 1.7% of the total number of receptors is internalized per minute. When this value is rescaled by the fraction of receptors in the plasma membrane, it translates into 18% of surface receptors internalized per minute. This implies that this internalization rate is ki = 1/(5.3 min), which is similar to the results obtained from reference [15]. The details of this rescaling are as follows. Mathematically, the internalization rate constant is defined as An external file that holds a picture, illustration, etc.
Object name is pcbi.0020003.ex002.jpg , where Nsur is the number of receptors at the plasma membrane. Mitchell et al. [14] measured An external file that holds a picture, illustration, etc.
Object name is pcbi.0020003.ex003.jpg , where Ntot is the total number of receptors. The fact that most of the receptors are internalized, so that Nsur≈0.1Ntot, leads to ki≈10kMitchell The factor 0.1 results from the fact that the internalization rate is about 10× higher than the recycling rate. Therefore, under stationary conditions, the number of internalized receptors is 10× higher than the number of receptors at the plasma membrane.

Ligand-induced degradation rate.

Active ligand-receptor complexes in lipid raft–caveolar compartments can recruit Smad7-Smurf2 [28], which then targets them for degradation [15]. Reference [15] shows that receptors are internalized through the clathrin pathway and lipid raft–caveolar compartments with similar rates. We have chosen klid = 1/(4min) for comparison with experimental data in Figure 5B.

Constitutive degradation rate.

Figure 3 of reference [15] shows that when the lipid raft–caveolar pathway is blocked with nystatin, only ~30% of the initially labeled receptors remain in the cell after 8 h. This gives a characteristic degradation time of ~400 min with respect to the total number of receptors. Rescaling this number to the plasma membrane receptors, we obtain An external file that holds a picture, illustration, etc.
Object name is pcbi.0020003.ex004.jpg . We have chosen kcd = 1/(36 min) for comparison with experimental data in Figure 5B.

Recycling rate.

Figure 3 of reference [14] shows that after about 30 min cells stop secreting internally labeled TGF-β receptors. This recycling rate is similar to that for the EGF receptor. We have chosen kr = 30 min for comparison with experimental data in Figure 5B.

Steady and quasi-steady state analysis.

Here we study mathematically the properties of the steady state of the system with a single ligand. By equating to 0 the derivatives in the model equations of Figure 4, we obtain that the steady-state number of internalized ligand-receptor complexes is

equation image

where the steady-state number of type I and type II receptors at the plasma membrane are obtained by solving the equations

equation image


equation image

The solution of these equations is

equation image

which leads to

equation image

For low values of c, Equation 5 reduces to

equation image

which indicates that the steady-state number of internalized ligand-receptor complexes is proportional to the ligand concentration and the production of each receptor type.

For high values of c, in contrast, we obtain

equation image


equation image

Therefore, for high ligand concentration or low constitutive degradation, the steady-state number of internalized ligand-receptor complexes is controlled by the receptor with the smallest production rate and this number does not depend on the ligand concentration.

The case of high c and low constitutive degradation is especially interesting because the steady-state signal does not depend on the ligand concentration, even when the ligand is present in small quantities. An important question to address now is: can the system detect changes in concentration in this regime? When the recycling rate is much lower than the internalization rate, the number of ligand-receptor complexes in the plasma membrane equilibrates faster than all the other variables. Therefore assuming quasi-equilibrium in this variable, An external file that holds a picture, illustration, etc.
Object name is pcbi.0020003.ex005.jpg , we obtain that, upon changes in the ligand concentration (Δ[l]), the changes in the number of ligand-receptor complexes in the plasma membrane (Δ[lRIRII]) follow the equation

equation image

Note that we have assumed that the number of receptors in the plasma membrane is conserved at these time scales.

For small changes in ligand concentration, we obtain

equation image

This expression indicates that for high c, low constitutive degradation, and slow recycling (compared to internalization), the system can detect changes in ligand concentration while keeping a steady-state signal that does not depend on ligand concentration.

Two-compartment model of receptor trafficking for two ligands.

The equations for a system with two ligands with concentrations [l1] and [l2] and are:

equation image
equation image
equation image
equation image
equation image
equation image
equation image
equation image
equation image
equation image

The variables [RI], [RII, 1], and [RII, 2], are the numbers of type I and type II receptors in the plasma membrane; and [l 1 RIRII, 1] and [l 2 RIRII, 2] refer to the corresponding ligand-receptor complexes. The overline indicates internalized receptors and ligand-receptor complexes. The signaling activity triggered by each ligand is assumed to be proportional to the corresponding number of internalized ligand-receptor complexes. kα is the rate constant of ligand-receptor complex formation; pRI, pRII, 1, and pRII, 2 are the rates of receptor production; ki, kr, kcd, and klid, are the internalization, recycling, constitutive degradation, and ligand-induced degradation rate constants; and α is the fraction of active receptors that are recycled back to the plasma membrane and can interact again with the ligand.

Coupled signaling and perfect adaptation.

Computational modeling offers precise insights into the functioning of the TGF-β pathway. It is possible to go a step further and generalize the conditions that give rise to different qualitative types of behavior.

Let us consider two ligands (l 1 and l 2), one type I receptor (RI), and two type II receptors (RII, 1 and RII, 2). The type I receptor is shared among the two ligand-receptor complexes [l 1 RI RII, 1] and [l 2 RI RII, 2]. The following conservation equations refer to the common type I receptor at steady-state.

Under stationary conditions, the number of receptors produced (by gene expression) is equal to the number of receptors degraded:

equation image

where p is the receptor production rate; and dconst and dlid are the constitutive and ligand-induced degradation rates, respectively. Assuming that a fraction δ of the activated receptors is degraded through a ligand-induced degradation process, we can express dlid as

equation image

where ia 1 and ia 2 are the rates of formation of the ligand-receptor complexes ([l 1 RI RII, 1] and [l 2 RI RII, 2], respectively). Therefore,

equation image

We explicitly consider two limiting cases:

In the first case, there is no ligand-induced degradation (dconst > 0, δ = 0). Therefore, we have dlid = 0, which leads to p = dconst. Because dconst [equivalent] dconst(RT) is a function of the total number of receptors RT, the previous condition indicates that the number of receptors remains constant RT = dconst −1(p), where dconst −1 is the inverse function of dconst. For instance, if the constitutive receptor degradation follows first order kinetics, dconst = γRT, then RT = p/γ. Under these conditions, if the rate of formation of complexes (ia 1 + ia 2) is small (for instance, for low ligand concentrations) compared with constitutive internalization and degradation, there is no coupling between signaling channels.

In the second case, there is only ligand-induced degradation (dconst = 0, δ > 0). Therefore, we have p = δ (ia 1 + ia 2), which implies that the formation of one ligand-receptor complex excludes the formation of the other one. In this case, the number of receptors in the plasma membrane does not remain constant, but is adjusted so that for a given ligand concentration the rate of formation of complexes (ia1 + ia2) remains equal to p/δ. As an explicit example, the kinetics ia 1 = l 1 RI RII, 1 and ia 2 = l 1 RIRII, 2 implies RI = p/(l 1 RII, 1 + l 2 RII, 2). Under these conditions, the completely coupled mode of signal integration arises even for low ligand concentrations.

The conditions that lead to the completely coupled mode of signal integration also lead to perfect adaptation for a single ligand. Consider, for example, ia 2 = 0. When there is only ligand-induced degradation, the fact that the rate of formation of complexes remains constant implies p = δia 1. The steady state of the system is fixed irrespective of the ligand concentration. Consequently, changes in ligand concentration can only elicit transient responses that completely adapt to the prestimulus level and the system exhibits perfect adaptation [29].


We are indebted to Joan Massagué and Van Le for invaluable help during the early stages of this work. We are also very grateful to Gary Bader, Debora Marks, Leonor Saiz, Nikolaus Schultz, Wenying Shou, and Stas Shvartsman for numerous discussions, comments, and suggestions. CS was funded in part by the Alfred W. Bressler Scholars Endowment Fund.


constitutive-to-induced degradation ratio
epidermal growth factor receptor
G protein-coupled receptor
regulatory Smad


Author contributions. JMGV conceived and designed the experiments, performed the experiments, and analyzed the data. RJ and CS contributed reagents/materials/analysis tools. JMGV and RJ wrote the paper. CS critically edited the manuscript.

Competing interests. The authors have declared that no competing interests exist.

A previous version of this article appeared as an Early Online Release on December 9, 2005 (DOI: 10.1371/journal.pcbi.0020003.eor).


  • Kang Y, Siegel PM, Shu W, Drobnjak M, Kakonen SM, et al. A multigenic program mediating breast cancer metastasis to bone. Cancer Cell. 2003;3:537–549. [PubMed]
  • Yamashita H, ten Dijke P, Franzen P, Miyazono K, Heldin CH. Formation of hetero-oligomeric complexes of type I and type II receptors for transforming growth factor-beta. J Biol Chem. 1994;269:20172–20178. [PubMed]
  • Goumans MJ, Valdimarsdottir G, Itoh S, Lebrin F, Larsson J, et al. Activin receptor-like kinase (ALK)1 is an antagonistic mediator of lateral TGFbeta/ALK5 signaling. Mol Cell. 2003;12:817–828. [PubMed]
  • Siegel PM, Massague J. Cytostatic and apoptotic actions of TGF-beta in homeostasis and cancer. Nat Rev Cancer. 2003;3:807–821. [PubMed]
  • Cui W, Fowlis DJ, Bryson S, Duffie E, Ireland H, et al. TGFbeta1 inhibits the formation of benign skin tumors, but enhances progression to invasive spindle carcinomas in transgenic mice. Cell. 1996;86:531–542. [PubMed]
  • Gotzmann J, Mikula M, Eger A, Schulte-Hermann R, Foisner R, et al. Molecular aspects of epithelial cell plasticity: Implications for local tumor invasion and metastasis. Mutat Res. 2004;566:9–20. [PubMed]
  • Thiery JP. Epithelial-mesenchymal transitions in development and pathologies. Curr Opin Cell Biol. 2003;15:740–746. [PubMed]
  • Steiner MS, Barrack ER. Transforming growth factor-beta 1 overproduction in prostate cancer: Effects on growth in vivo and in vitro. Mol Endocrinol. 1992;6:15–25. [PubMed]
  • Gonzalez-Gaitan M. Signal dispersal and transduction through the endocytic pathway. Nat Rev Mol Cell Biol. 2003;4:213–224. [PubMed]
  • Lauffenburger DA, Linderman JJ. Receptors: Models for binding, trafficking, and signalling. New York: Oxford University Press; 1993. 365. p.
  • Wiley HS, Shvartsman SY, Lauffenburger DA. Computational modeling of the EGF-receptor system: A paradigm for systems biology. Trends Cell Biol. 2003;13:43–50. [PubMed]
  • Tan CM, Brady AE, Nickols HH, Wang Q, Limbird LE. Membrane trafficking of G protein-coupled receptors. Annu Rev Pharmacol Toxicol. 2004;44:559–609. [PubMed]
  • Bohm SK, Grady EF, Bunnett NW. Regulatory mechanisms that modulate signalling by G-protein-coupled receptors. Biochem J. 1997;322((Pt 1)):1–18. [PubMed]
  • Mitchell H, Choudhury A, Pagano RE, Leof EB. Ligand-dependent and -independent transforming growth factor-beta receptor recycling regulated by clathrin-mediated endocytosis and Rab11. Mol Biol Cell. 2004;15:4166–4178. [PMC free article] [PubMed]
  • Di Guglielmo GM, Le Roy C, Goodfellow AF, Wrana JL. Distinct endocytic pathways regulate TGF-beta receptor signalling and turnover. Nat Cell Biol. 2003;5:410–421. [PubMed]
  • Inman GJ, Nicolas FJ, Hill CS. Nucleocytoplasmic shuttling of Smads 2, 3, and 4 permits sensing of TGF-beta receptor activity. Mol Cell. 2002;10:283–294. [PubMed]
  • Yarden Y, Sliwkowski MX. Untangling the ErbB signalling network. Nat Rev Mol Cell Biol. 2001;2:127–137. [PubMed]
  • Vaudry D, Stork PJ, Lazarovici P, Eiden LE. Signaling pathways for PC12 cell differentiation: Making the right connections. Science. 2002;296:1648–1649. [PubMed]
  • Nicolas FJ, Hill CS. Attenuation of the TGF-beta-Smad signaling pathway in pancreatic tumor cells confers resistance to TGF-beta-induced growth arrest. Oncogene. 2003;22:3698–3711. [PubMed]
  • Dumont N, Arteaga CL. A kinase-inactive type II TGFbeta receptor impairs BMP signaling in human breast cancer cells. Biochem Biophys Res Commun. 2003;301:108–112. [PubMed]
  • Ye SC, Foster JM, Li W, Liang J, Zborowska E, et al. Contextual effects of transforming growth factor beta on the tumorigenicity of human colon carcinoma cells. Cancer Res. 1999;59:4725–4731. [PubMed]
  • Tobin SW, Douville K, Benbow U, Brinckerhoff CE, Memoli VA, et al. Consequences of altered TGF-beta expression and responsiveness in breast cancer: Evidence for autocrine and paracrine effects. Oncogene. 2002;21:108–118. [PubMed]
  • Worthylake R, Opresko LK, Wiley HS. ErbB-2 amplification inhibits down-regulation and induces constitutive activation of both ErbB-2 and epidermal growth factor receptors. J Biol Chem. 1999;274:8865–8874. [PubMed]
  • Hartwell LH, Hopfield JJ, Leibler S, Murray AW. From molecular to modular cell biology. Nature. 1999;402:C47–C52. [PubMed]
  • Vilar JM, Guet CC, Leibler S. Modeling network dynamics: The lac operon, a case study. J Cell Biol. 2003;161:471–476. [PMC free article] [PubMed]
  • Stelling J, Sauer U, Szallasi Z, Doyle FJ, III, Doyle J. Robustness of cellular functions. Cell. 2004;118:675–685. [PubMed]
  • Zerial M, McBride H. Rab proteins as membrane organizers. Nat Rev Mol Cell Biol. 2001;2:107–117. [PubMed]
  • Monteleone G, Pallone F, MacDonald TT. Smad7 in TGF-beta-mediated negative regulation of gut inflammation. Trends Immunol. 2004;25:513–517. [PubMed]
  • Csete ME, Doyle JC. Reverse engineering of biological complexity. Science. 2002;295:1664–1669. [PubMed]
  • Massague J. TGF-beta signal transduction. Annu Rev Biochem. 1998;67:753–791. [PubMed]
  • Sun PD, Davies DR. The cystine-knot growth-factor superfamily. Annu Rev Biophys Biomol Struct. 1995;24:269–291. [PubMed]
  • de Caestecker M. The transforming growth factor-beta superfamily of receptors. Cytokine Growth Factor Rev. 2004;15:1–11. [PubMed]
  • Shi Y, Massague J. Mechanisms of TGF-beta signaling from cell membrane to the nucleus. Cell. 2003;113:685–700. [PubMed]

Articles from PLoS Computational Biology are provided here courtesy of Public Library of Science