|Home | About | Journals | Submit | Contact Us | Français|
E3 ubiquitin ligases are important cellular components that determine the specificity of proteolysis in the ubiquitin-proteasome system. However, an increasing number of studies have indicated that E3 ubiquitin ligases also participate in transcription. Intrigued by the apparently paradoxical functions of E3 ubiquitin ligases in both proteolysis and transcriptional activation, we investigated the underlying design principles using mathematical modeling. We found that the antagonistic functions integrated in E3 ubiquitin ligases can prevent any undesirable sustained activation of downstream genes when E3 ubiquitin ligases are destabilized by unexpected perturbations. Interestingly, this design principle of the system is similar to the operational principle of a safety interlock device in engineering systems, which prevents a system from abnormal operation unless stability is guaranteed.
The ubiquitin-proteasome system (UPS) is primarily known as an important part of the protein degradation (proteolysis) machinery. The proper and timely degradation of proteins is critical to maintaining the homeostasis of biological systems, and a failure of the UPS can cause abnormal development1 and devastating diseases such as cancer2,3,4, neurodegeneration5, cardiac diseases6,7, autoimmunity, and inflammatory diseases8,9.
UPS-dependent proteolysis requires the attachment of multiple ubiquitin molecules to a target protein, which results in the subsequent degradation of the polyubiquitinated target protein through the 26S proteasome complex (Fig. 1a)3,10,11,12,13. The ubiquitination step involves at least three enzymes: a ubiquitin-activating enzyme (E1), a ubiquitin-conjugating enzyme (E2), and a ubiquitin-protein ligase (E3). E1 activates a ubiquitin molecule via an ATP-dependent reaction, forming a high-energy thioester linkage between E1 and a ubiquitin. The activated ubiquitin is transferred to E2; subsequently, E3, supported by E2, attaches the ubiquitin to a lysine residue of the substrate protein or a lysine residue of the previously attached ubiquitin. E3, in particular, plays a key role in the recognition of the amino acid sequence of a target protein, which is called a “degron.” Thus, E3 determines the specificity for the “death” of the target protein molecule. Finally, the polyubiquitinated protein is recognized by a specific subunit of the 19S regulatory particle of the proteasome and is degraded by the 20S core subunit in an ATP-dependent manner.
Several recent studies have demonstrated that the UPS plays a role in transcriptional activation (Fig. 1b)13. For example, the UPS was found to activate a transcription factor on the promoter of a target gene via ubiquitination14,15,16,17. Other studies have indicated that E3 ubiquitin ligases can act as transcription cofactors in the absence of ubiquitin ligase activity18,19. In addition, it was found that a certain type of polyubiquitination does not promote the degradation of a transcription factor but rather stabilizes it20,21. Figure 2a presents examples in which E3 ubiquitin ligases facilitate both proteolysis and transcription in various signaling pathways. (1) Skp2 not only promotes the degradation of Myc but also enhances the transcription of a subset of Myc target genes with or without ubiquitin ligase activity14,15,19. (2) The regulatory relationship of β-catenin and β-TrCP1 (also known as Fbw1a) in the canonical Wnt pathway22 is similar to that of Myc and Skp218; the major difference is that the Wnt/β-catenin pathway involves phosphorylation for the recognition of β-catenin by β-TrCP123. (3) β-TrCP1 also facilitates the proteolysis of Smad4 and the Smad3/Smad4 complex24,25,26, whereas it participates in the transcriptional activation of Smad3 without ubiquitin ligase activity18. (4) For Smad2 and Wwp2 (also known as AIP2), Wwp2 oppositely regulates TGF-β signal transduction to its target genes by participating in both the proteolysis of Smad227,28 and the transcriptional activity of Goosecoid (Gsc)17 or Sox929 (see Supplementary Methods (section S1) for details).
We were intrigued as to why a single protein, i.e., the E3 ubiquitin ligase, participates in both proteolysis and transcriptional activation (i.e., the downregulation and upregulation, respectively, of a signal transduction pathway), which are apparently opposing reactions. In particular, we wondered what would be the advantage of simultaneously promoting both processes via this single component of the UPS. To decipher the design principles underlying these paradoxical roles of E3 ubiquitin ligases, we developed a generic mathematical model that represents the essential dynamics of the signaling pathways exemplified in Figure 2a, which we call integrated transcription and UPS-dependent degradation (ITUD)18 (Fig. 2b). For comparison, we considered a simple negative feedback loop (SNFL) and dissociated transcription and UPS-dependent degradation (DTUD) as alternative models of the ITUD (Fig. 2b). Based on a mathematically controlled comparison of these models, we found that only the ITUD system can suppress abnormal responses under unexpected downregulation of the E3 ubiquitin ligase by which the cognate transcription factor might be excessively stabilized. We further confirmed these characteristics of ITUD at the cell population level using in silico cell proliferation and migration experiments. Interestingly, the design principle that underlies the paradoxical roles of the E3 ubiquitin ligase in the ITUD model is similar to the operational principle of a safety interlock device in engineering systems in which any unexpected perturbation of the safety device results in the system shutdown to prevent abnormal operation.
ITUD is the core structure of the signaling networks presented in Figure 2a, which includes only the essential components for investigation of the underlying design principles of biological systems30,31,32,33,34,35,36 (see Methods and Supplementary Table S1). The ITUD system consists of five components: the signal (S), the transcription factor (T), polyubiquitinated T (ubT), E3 ubiquitin ligase (E3), and the target gene product (P) (Fig. 2b and 2c). The signal, S, is an extracellular ligand or a subcellular biomolecule, which increases the amount of T. T is a transcription factor that promotes the expression of the target gene P. It is assumed that T can undergo UPS-dependent proteolysis promoted by E3 in both cytoplasm and nucleus. ubT is the polyubiquitinated form of T, which undergoes proteasomal degradation. E3 is an E3 ubiquitin ligase that not only facilitates the ubiquitin-proteasomal degradation of T but also ‘paradoxically' aids in the transcriptional function of T. It is assumed that the stability of T is not affected by either the proteolytic or the nonproteolytic function of E3 in transcription37,38 (see Supplementary Methods (section S1 and S2)). E3 is also transactivated by T, and therefore, E3 and T form a negative feedback loop by which the appropriate level of T is maintained. For instance, Skp2 is a direct target gene of c-Myc in leukemia39 and melanoma cells40. P is the final product, or the output of the system, whose expression is activated by both T and E3.
We employed the method of mathematically controlled comparison (MCC), which identifies the essential characteristics of a mathematical model by comparison with alternative models according to a set of mathematical criteria and statistical methods41,42,43,44. To understand the essential characteristics of ITUD, we compared it with two alternative systems having different configurations: SNFL and DTUD (Fig. 2b). SNFL represents a negative feedback loop where E3 does not participate in the transcriptional activation of P. In SNFL, T transactivates E3 and the enhanced E3 only regulates T through UPS-dependent proteolysis as a feedback control. In DTUD, the proteolysis of T and the transcriptional activation of P by T are processed by different components: E3 only participates in the degradation of T, and another component, I (i.e., the intermediate node), assumes the role of E3 for the transcriptional activation of P. We developed mathematical models for SNFL, DTUD, and ITUD such that the mathematical terms of the biochemical reactions and their common parameters were established to be comparable on the basis of MCC (Supplementary Table S1, Table S2, and Supplementary Methods (section S2)).
We wondered what distinguishes ITUD from the alternative systems. To answer this question, we explored the critical determinants that characterize ITUD. Specifically, we analyzed how these three systems differently respond to an identical small perturbation of each kinetic parameter (i.e., local sensitivity analysis, LSA), in which the system output is the steady-state level of P (Pss). Among the 20 common kinetic parameters, five parameters (ID: 11–15) exhibited a noticeable difference (i.e., <20%) in the ratio of the sensitivities between ITUD and the other two systems (Table 1 and Fig. 3). The ratio of sensitivities denotes the ratio of how the steady-state level of P changes between the two systems in response to a small perturbation on the same parameter. So, a small ratio of the sensitivities means that the two systems have a large difference in their responses to the change of the parameter value. To overcome the limitations of using specific values for a parameter set, we further performed LSA based on MCC over a wide range of parameter values (see Methods). The five parameters (ID: 11–15) also exhibited different patterns as compared to the remaining parameters (Supplementary Fig. S1). These five parameters are involved in the regulation of the level of E3 and, thus, are important for the stability of E3; i.e., ITUD and the other two systems differently respond to a change in the level of E3. Specifically, βE3 (ID: 11) and αE3 (ID: 12) represent the basal production and degradation rates, respectively, which determine the basal level of E3. βTE3 (ID: 13), KTE3 (ID: 14), and nTE3 (ID: 15) represent the maximum transactivation rate, half-maximal occupancy, and Hill-coefficient, respectively, in the regulation of E3 by T. Among the five parameters, βTE3 (ID: 13) and KTE3 (ID: 14) are the critical determinants because βE3 (ID: 11) and αE3 (ID: 12) do not affect the signal transduction from S to E3, and the Hill coefficient represents the cooperativity of ligand binding45, which is not within the primary scope of this study. Therefore, the regulation (represented by the two kinetic parameters, βTE3 and KTE3) of E3 by T distinguishes ITUD from the two alternative systems, and we can conclude that the two kinetic parameters are critical determinants.
We further investigated how the critical determinants (i.e., βTE3 and KTE3) characterize ITUD in comparison to the two alternative systems. In contrast to the small discrete perturbation in LSA, we perturbed E3 by continuously varying βTE3 or KTE3 such that E3 was upregulated initially and downregulated subsequently, and measured the steady-state levels of the system components. Pss rose gradually until βTE3 was decreased to the nominal value of 1.0 in all systems (Fig. 4a). However, below the nominal value, Pss was dramatically decreased in ITUD, whereas it is finally saturated in the other two systems. For KTE3, the increasing phase of Pss in ITUD was observed within a small range, and the decreasing phase began near the nominal value of 0.5 (Fig. 4b). This result suggests that ITUD exhibits a biphasic response with respect to the factors that alter the stability of E3 and that the output level is lowered only in ITUD when the amount of E3 is reduced by perturbations.
To explore the biphasic characteristics of ITUD with respect to the level of E3, we applied the steady-state analysis solely on the ITUD system through direct manipulation of the steady-state E3 level (E3ss). The regulation of E3 by T was deliberately removed to modulate the E3ss irrespective of the signal S (Fig. 4c). At the higher E3ss, the outputs of all three systems decreased due to the enhanced proteolysis facilitated by an excess of E3. However, at the lower E3ss, only ITUD demonstrated a decreasing phase, whereas the other two systems exhibited saturated responses (Fig. 4d and 4e). Thus, ITUD demonstrates a differential behavior as compared to the alternative systems at a lower E3 level.
We also examined how binding affinities (or reaction intensities) of E3 for T in proteolysis (Km_E3T) and transcription (KE3P) affect the shapes of the response curves. In all three systems, increasing Km_E3T (i.e., decreasing proteolysis) reduced the sharpness of the curves, whereas KE3P primarily altered the amplitudes of the curves. In addition, the effect of Km_E3T on the response curves was large at the higher E3ss, whereas that of KE3P was prominent at the lower E3ss. Similar patterns of the biphasic response were also observed in the intact ITUD (Supplementary Fig. S2). These results suggest that the biphasic characteristics of the ITUD system are tunable by modulating the regulatory mechanisms involved in proteolysis and transcription such as posttranslational modifications or subcellular localization.
Having identified the critical determinants that characterize the biphasic responses of ITUD at the steady state, we sought to determine how the critical determinants affect the temporal dynamics. Hence, we observed the temporal profile of the three systems under the perturbation of the critical determinants. As the signal S increased, T was activated, and subsequently the levels of E3, ubT, and P were also increased in the control condition (no perturbation) (Fig. 5a). The three systems exhibited the identical temporal dynamics for T, ubT, and E3 because they share the negative feedback loop comprising T, ubT, and E3. However, the output responses of DTUD and ITUD were slower than that of SNFL, in which the response time was defined as the time required for a system output to exceed the activation threshold (i.e., PThreshold = 0.5; approximately half of the maximum level of P at the control condition) after the input signal was applied. This slower response is attributable to a “sign-sensitive delay” of the AND-gated coherent feed-forward loop46,47 embedded in DTUD and ITUD. The presence of a sign-sensitive delay means that the output response is delayed to an extent almost identical to the duration for the expression of intermediate components such as E3 and I if the input signal is increased (plus), whereas it is not delayed when the input is decreased (minus). By contrast, when S was diminished T was rapidly degraded by E3 and exhibited undershooting dynamics. The decreasing dynamics of P was similar in all three systems. Notably, the three systems also exhibited similar overall dynamics over a wide range of kinetic parameter values with or without noise in the signal (Supplementary Fig. S3, control).
In the perturbation condition, we perturbed the critical determinants (i.e., βTE3 and KTE3) toward the downregulation of E3. The perturbation for the upregulation of E3 was not considered because it did not noticeably distinguish ITUD at the higher E3ss in the steady-state analysis (Fig. 4). When βTE3 was decreased by 80% (0.2-fold), the outputs of SNFL and DTUD nearly reached the activation threshold even at a low signal S level (Fig. 5b). By contrast, ITUD was irresponsive to the signal. When KTE3 was increased by 150% (2.5-fold), SNFL and DTUD exhibited similar overactivated dynamics as in the βTE3 perturbation case (Fig. 5b), in which the outputs exceeded the activation threshold even at the low S level (Fig. 5c). However, ITUD exhibited an attenuated response (a response time of 2.09 h in comparison to 1.33 h in the control).
The temporal dynamics under the perturbation of the critical determinants were also investigated using a noisy signal because biological environments fluctuate and are noisy48,49. The dynamics of SNFL and DTUD under the perturbation of βTE3 or KTE3 were similar to those of the control condition in Figure 5a (Fig. 5d and 5e). However, ITUD under the perturbation of βTE3 with a noisy signal was irresponsive, as shown in Figure 5b (Fig. 5d). For KTE3, the output of ITUD failed to exceed the activation threshold (Fig. 5e), suggesting that ITUD lost its functionality under this condition. The difference between ITUD and the other two systems was more evident upon variation of the kinetic parameters with or without noise in the signal (Supplementary Fig. S3). Therefore, the temporal dynamics results suggest that the speed and amplitude of the system response is attenuated only in ITUD, whereas the other two systems are overactivated when the E3 ubiquitin ligase is destabilized and thus fail to regulate the transcription factor.
The analyses of the steady-state response and the temporal dynamics suggest that ITUD is characterized by the critical determinants that dominate the stability of E3 under the effect of signal S at the single-cell level. To understand the characteristics of ITUD at the cell population level, we performed in silico analysis of the cell population dynamics50,51,52,53,54. Because the downstream genes or cellular functions activated by the signaling pathways in Figure 2a are primarily related to cell proliferation55,56, migration, or metastasis19,57, we designed the simulations for two in silico experiments: 1) cell proliferation and 2) cell migration. The SNFL, DTUD, or ITUD system is implemented in a cell model (see Supplementary Methods (section S3)) as a subcellular pathway (Fig. 6a). In a two-dimensional simulation space, the cell population consists of five types of cells: 1) normal cells, 2) active normal cells, 3) perturbed cells, 4) active perturbed cells, and 5) stromal cells (Fig. 6b and Supplementary Methods (section S3)). Stromal cells secrete the signal molecule S but do not possess any subcellular pathway. The other cells receive the signal to be activated, and only the active cells are capable of cell division or chemotaxis-driven migration. The activation of each cell is determined by measuring the output of its subcellular pathway (i.e., the level of P) (Fig. 6a). The cell population was configured to receive a stochastic signal and perturbation of the critical determinants toward the downregulation of E3. Thus, the subcellular pathway in each cell exhibited similar temporal dynamics, as shown in Figure 5d and 5e.
In the absence of perturbation, the cell populations of all three systems grew rapidly but gradually declined after reaching the maximum population due to cellular senescence and apoptosis (Fig. 6c, left and Supplementary Fig. S4a). Conversely, in the perturbation condition, the cell populations of SNFL and DTUD continuously increased despite cellular senescence and apoptosis. These increases in the SNFL and DTUD populations are akin to tumorigenic cell proliferation. However, the ITUD population exhibited dynamics with a smaller peak in comparison to the control condition (Fig. 6c, right and Supplementary Fig. S4b). Notably, when the perturbation effect of DTUD was reduced to half by introducing a ‘buffering effect' (DTUD-BE) through the intermediate node (I), the population dynamics was similar to that of ITUD in the perturbation condition with a larger peak (Fig. 6c, right). This result demonstrates that a ‘semi-integration' of proteolysis and transcription in DTUD-BE under the perturbation (i.e., the perturbation of both I and E3) can prevent uncontrolled proliferation that was observed in the cell population of the intact DTUD. The results of the in silico cell proliferation dynamics are also provided as snapshots and movies (Fig. 6d and Supplementary Movies S1–S2).
For in silico cell migration experiments, cells were laid in a 200 μm × 80 μm space and were attracted by signal S that was secreted by stromal cells (Fig. 6e). We varied the probability of the perturbation occurrence (i.e., P1, P2, and P3; see Supplementary Methods (section S3)) and observed the cell migration. When the perturbation occurred more frequently, the mean displacement of cells possessing SNFL or DTUD was slightly increased. Conversely, the cell population of ITUD exhibited a significantly reduced migration (Fig. 6f and Supplementary Movies S3–S4). We also tracked the migration of a single cell for the most frequent perturbation (P3). In the 25-h simulation, SNFL and DTUD cells moved toward stromal cells, whereas the ITUD cell wandered around its initial position (Fig. 6g). This result suggests that cell migration can also be attenuated by ITUD if the E3 ubiquitin ligase is destabilized and thereby fails to regulate the transcription factor that drives cell migration.
To experimentally validate our results obtained from the in silico analysis at the population level, we performed in vitro cell proliferation and migration experiments in HEK293T cells. These experiments were confined to ITUD only since there is a limitation in preparing experimental systems that exactly correspond to the three mathematical models (see Supplementary Note 1). The perturbation of the critical determinants was provided through the knockdown of E3 ubiquitin ligases (e.g., Wwp2 and β-TrCP1) or a protein that regulates the production of E3 proteins (e.g., CRD-BP) (Fig. 2a). Subsequently, cell proliferation or migration was observed, as shown in Figure 6.
The proliferation of HEK293T cells was significantly reduced when we transiently knocked down Wwp2, β-TrCP1, or CRD-BP (after 96-h culture: 65.2%, 85.9%, and 62.5% of the control, respectively) (Fig. 7a). For cell migration, we considered Wwp2 in the Smad2/Gsc pathway (Fig. 2a, fourth panel). Gsc has been reported to promote EMT (epithelial-mesenchymal transition), migration, and metastasis of breast cancer cells57. In addition, Wwp2 transactivates Gsc17. Thus, we knocked down Wwp2 in HEK293T cells, expecting that the transcriptional effects of Gsc on cell migration would be attenuated through the downregulation of Wwp2. The migration of HEK293T cells was found to be significantly reduced through the knockdown of Wwp2 (Fig. 7b). The efficiency of RNAi-mediated knockdown of Wwp2, β-TrCP1, and CRD-BP transcripts is shown in Supplementary Figure S5a–c.
We also observed the expression of target genes that were expected to be affected by the knockdown of Wwp2, β-TrCP1, or CRD-BP in the cell proliferation and migration assays. The known pro-proliferative gene c-Jun58 was selected as a target gene (P) of Gsc(T)/Wwp2(E3) and Smad3(T)/β-TrCP1(E3)59. The knockdown of Wwp2 or β-TrCP1 resulted in the reduced expression of c-Jun (Supplementary Fig. S5a and S5b). The knockdown of CRD-BP, which is in the Wnt/β-catenin pathway (Fig. 2a, second panel), efficiently decreased β-TrCP1 and thereby the expression of CyclinD1, a pro-proliferative target gene of β-catenin (Supplementary Fig. S5c). For the cell migration assay, Wwp2 knockdown upregulated FGF2 which was reported to inhibit cell migration in breast cancer cells60,61 and MEF cells62 (Supplementary Fig. S5a). These experimental results are consistent with our simulation results indicating that cell proliferation and migration are significantly reduced through the downregulation of E3 in ITUD. However, the results are not entirely attributable to the effects of ITUD because E3 proteins possess many substrates, and, thus, might have pleiotropic effects on cellular functions3,4,63,64,65.
We further investigated whether the decreasing phase of ITUD caused by the enhanced proteolysis of T by an excess of E3 (Fig. 4d and 4e) can also be observed in the target gene expression and the actual cell proliferation. For this purpose, we overexpressed CRD-BP to stabilize β-TrCP1 mRNA in HEK293T cells, which corresponds to an increase of βTE3 in ITUD. Consequently, the overexpression of CRD-BP reduced both the proliferation rate of HEK293T cells (Supplementary Fig. S5d) and CyclinD1 expression (Supplementary Fig. S5e), which is consistent with our simulation results.
Perturbation of the critical determinants for both downregulation and upregulation of E3 (corresponding to the increase and decrease of T, respectively) reduced the output of the ITUD system (Fig. 8a). The maximum output of ITUD is achieved by an intricate balance between T and E3. Thus, ITUD exhibits biphasicity with respect to the stability of the E3 ubiquitin ligase. We suggest two molecular mechanisms underlying the biphasicity of ITUD. The first mechanism involves complex formation or dimerization of the E3 ubiquitin ligase (Fig. 8b, top). Some E3 ubiquitin ligases must form a complex to achieve their proteolytic function. β-TrCP1 or Skp2, for example, interacts with Skp1 to form an SCF complex for proteolysis3, whereas these proteins can also participate in transcription without complex formation18,19. In addition, dimerization is necessary for some E3 proteins to achieve ubiquitin ligase activity (Table 2). This evidence indicating the necessity of the dimerization of and complex formation by E3 proteins for their proteolytic function suggests that the maximum, or the optimal response of the ITUD system, may be established by a subtle balance between the monomeric form for transcription and the oligomeric form for proteolysis. The second mechanism is based on subcellular localization (Fig. 8b, bottom). Several isoforms of β-TrCP1 and β-TrCP2 have been reported to predominantly localize in the nucleus66,67. Assuming that proteolysis is facilitated primarily in the cytoplasm, transcriptional activation may be initially promoted by a small amount of E3 proteins that preferentially localize in the nucleus. Eventually, the accumulation of the E3 ubiquitin ligase in the cytoplasm and nucleus would enhance proteolysis, which prevents the transcription factor from entering the nucleus to participate in transcription, thereby suppressing transcriptional activation (these two mechanisms can be classified as Type I and Type IV biphasic regulation according to Levchenko et al.68).
We suppose that the characteristics of the biphasic response in ITUD might be modulated by regulating the extent of phosphorylation/dephosphorylation or nuclear-cytoplasmic shuttling. For example, the proteolysis of β-catenin by β-TrCP1 requires the phosphorylation of β-catenin, and the phosphorylation can be reversed by a phosphatase such as PP2A. Therefore, the intensity of phosphorylation/dephosphorylation of the transcription factor by an upstream enzyme determines both the timing and the extent of E3-mediated proteolysis, which consequently alters the biphasic response (Fig. 8c, left). This effect was observed by varying Km_E3T in ITUD (Fig. 4d). Alternatively, nuclear-cytoplasmic shuttling is another possible mechanism that exemplifies the effects of varying KE3P in ITUD (Fig. 4e). For example, the nuclear-cytoplasmic shuttling of Skp2 is regulated by phosphorylation of Akt69 or the acetylation of p30070. The upstream signals of Akt or p300 may modulate the biphasic response by regulating the nuclear-cytoplasmic shuttling of Skp2 for its transcriptional activation (Fig. 8c, right).
Biological circuits employ elaborate structures and dynamic properties to achieve specialized functions31. Interestingly, such circuits also possess mechanisms to protect the biological systems from intrinsic or extrinsic perturbations. A representative example is the DNA repair system, which ensures the integrity of genetic information71. Another intriguing example is the responsive backup circuit of genetic redundancy, in which any problem that reduces the stability of the activator is buffered by a functionally redundant partner to guarantee a robust output72,73,74. This principle is comparable to the redundancy of critical components for system reliability in engineering. We suggest that the apparently paradoxical roles of E3 ubiquitin ligases in ITUD may represent a safety interlock device that prevents the overactivation of a molecular or cellular function due to the destabilization of E3. A safety interlock device does not allow the entire system to operate unless the safety condition is satisfied. A familiar example of a safety interlock device is the door of a washing machine (Fig. 8d, left). In the washing machine, the door is interlocked with the motor so that the motor can operate only if the door safety is guaranteed. E3 and the transcription factor in ITUD correspond to the door and motor of the washing machine, respectively (Fig. 8d, right). So, E3 ubiquitin ligase functions as a safety interlock device which prevents a devastating overactivation as observed in SNFL and DTUD (Fig. 5 and and6),6), by controlling an appropriate level of the transcription factor through UPS-dependent proteolysis. If E3 is not functional due to a certain perturbation in ITUD, transcription should be attenuated in this case since ‘proteolysis for safety' is not guaranteed. Together, we find that both washing machine and ITUD can cope with the breakdown of safety device by shutting down the system to avoid its abnormal operation.
The design principles of ITUD may also represent an important mechanism used to prevent potentially dangerous functions in multicellular organisms (Fig. 8e). Proliferation or migration of cells in a tissue should be tightly regulated because the overactivation of such cellular functions can disrupt tissue homeostasis and even result in death. A typical example is the manifestation of abnormal cellular physiology in cancer development75. The strategy of ITUD in support of tissue homeostasis is that of “if impossible to control potentially dangerous cellular functions such as cell division or migration, it would be better not to activate it.” If the stability of E3 ubiquitin ligase in a cell collapses owing to any perturbation, ITUD suppresses transcriptional activation (Fig. 5b–e) and, therefore, the corresponding cellular functions (Fig. 6c and 6f). This preventive strategy of ITUD may represent one of the mechanisms for cellular dormancy under harsh environments76. Moreover, the paradoxical bifunctionality of E3 ubiquitin ligases presents an interesting example of the utility of paradoxical components77 to ensure the prevention of disastrous consequences in the cell population.
The consideration of E3 ubiquitin ligases as targets for cancer treatment relates to the successful application of the proteasome inhibitor bortezomib (commercialized name, Velcade®) for multiple myeloma78. The rationale for the development of therapeutics targeting E3 proteins is based on enhancing the specificity of the drug to reduce the side effects due to proteasome inhibition. However, targeting an E3 protein, particularly the SCF E3 ubiquitin ligase, can cause unpredictable effects because a single E3 is able to promote the proteolysis of both oncogenes and tumor-suppressor genes3. Therefore, we suggest a combinatorial therapy to downregulate both oncogenes and cognate E3 proteins by exploiting the E3 function as the safety interlock device in the ITUD system. Targeting of a single oncogene or a specific signal transduction pathway based on the rationale of “oncogene addiction” is one of the basic approaches in cancer treatment79. However, this approach typically includes drawbacks such as dose toxicity80 and drug resistance81. To mitigate such problems in a given ITUD signaling pathway, we can distribute the perturbation strength from an addictive oncogene to the E3 ubiquitin ligases to reduce toxicity and to diversify targets to block as many escape routes as possible. A toy example of the steady-state analysis in ITUD system indicated that 114% (i.e., 57% for each T and E3) of the single perturbation strength on T is required to inhibit both T and E3 to attain an identical result (Fig. 8f and see Methods for further details). If the dose required to cause this perturbation in clinical treatment is known, our approach can represent a promising strategy for the treatment of cancer or other obstinate diseases.
In this study, the role of E3 in transcription was limited to the transcriptional activation. However, E3 ubiquitin ligases can also participate in transcriptional repression independent of proteolysis82,83,84. In yeast, Met30, an SCF E3 ubiquitin ligase, not only promotes the UPS-dependent proteolysis of Met4, a transcriptional activator85,86,87,88,89, but also represses the transcriptional activity of Met4 through non-proteolytic ubiquitination84. In this case, the roles of E3 in transcription and proteolysis are not paradoxical, but consistent in inhibiting the function of the cognate transcription factor. On the other hand, BRCA1, a tumor suppressor E3 in breast cancers, can activate or repress transcription pre-initiation complex depending on the context82,90. Therefore, the effects of BRCA1 should be carefully interpreted in consideration of the cellular environment. Although the design principle of ITUD that we unraveled in this study cannot be generalized to expound all the complicated roles of E3 in transcription and proteolysis, it is still useful to get a new insight into the apparently paradoxical roles of E3. We need to further investigate the design principles underlying the other roles of E3 as in the case of Met30.
Numerical solutions and random numbers were obtained using MATLAB® 2012a (64 bit). The fsolve and ode15s/ode45 functions were used for steady-state analysis and temporal dynamics, respectively. For the noise signal in temporal dynamics (Fig. 5d and 5e), we generated a random number that is uniformly distributed on (0, 1) as an amplitude of the signal for each time point using rand function.
Mathematical models of SNFL, DTUD, and ITUD are systems of ordinary differential equations as follows (the ordinary differential equations for T, ubT, E3 of SNFL and DTUD are the same as those of ITUD):
The individual biological processes are described by typical kinetic terms such as Michaelis-Menten type kinetics and Hill equation. Refer to Supplementary Methods (section S2) for model derivation and major assumptions.
In this study, local sensitivity was defined as the relative change in the system output (i.e., the steady-state level of P) when a small perturbation (±1%) was imposed on a single kinetic parameter92. The relative local sensitivity coefficient was approximated by the following formula:
where s is the system output, and pi is a parameter of ID number i. The steady-state levels of the system components were numerically solved, and thus, did not rely on the use of power-law formalism and an analytic solution93.
We sampled 5000 sets of parameter values for the steady-state analysis over a wide range of parameter values (Supplementary Fig. S1). The parameter values were generated using a uniform random variable on [0.01, 100]. Subsequently, we calculated M and R, which are defined as
where M is the relative sensitivity coefficient, as previously defined, and R is the ratio of two Ms between ITUD and an alternative system such as SNFL and DTUD. Therefore, R represents the ratio of the rate of change in the system outputs as follows:
The external equivalence of MCC was achieved for each parameter set (see Supplementary Methods (section S2)). The 5000 data entities of a pair (M, R) were sorted according to MITUD, and the moving quantile method with a window size of 1000 was applied as in the method of Alves and Savageau42, in which the median of R (i.e., quantile 0.5) and its corresponding MITUD were obtained from a subset of the data extracted by the sliding window. Thus, we finally obtained 4001 data entities for <R> and <M> (the total number of windows is “sample size − window size + 1”; i.e., 5000 − 1000 + 1). To exclude extreme cases, 10% of each end in the total windows was trimmed (i.e., 80% of 4001 windows 3200 windows remaining), and density plots were generated (Supplementary Fig. S1).
We first applied the perturbation on T by disrupting the signal flow from S to T (single perturbation), which was implemented by increasing the parameter value of Km_ST (ID: 4). Next, both values of Km_ST (ID: 4) and KTE3 (ID: 14) were increased for the inhibition of both T and E3 (dual perturbation). We defined a perturbation strength (PS) as the percentage change of a kinetic parameter value to the perturbation of a system component (e.g., T or E3) as follows:
where perturbed(pi) and nominal(pi) are perturbed and nominal value of parameter pi, respectively. Therefore, the single perturbation strength (PSs) of Km_ST, when it is perturbed from the nominal value 0.5 (Supplementary Table S2) to 1.5, is calculated as follows:
In this case, Pss under single perturbation is 0.0721, which is decreased by 91.8% (0.08 fold) from Pss without any perturbation (i.e., 0.8781). The relative perturbation strength (RPS) was obtained based on the ratio between the single and dual perturbation strengths as follows:
where PSd(pi) represents the dual perturbation strength for a given parameter pi with the inhibition of both T and E3. For example, if we change Km_ST from 0.5 to 1.5 as the single perturbation and both Km_ST and KTE3 from 0.5 to 0.75 as the dual perturbation, the relative perturbation strength in Figure 8f is calculated as follows:
For simplicity, the perturbation strengths of Km_ST and KTE3 were set to be identical in the dual perturbation. So, RPS was calculated as follows:
where PSd = PSd(Km_ST) = PSd(KTE3). The relative output in Figure 8f represents the ratio of Pss between the single and dual perturbation. Pss under single perturbation is 0.0721, so the relative output 0.5 means that Pss under dual perturbation is 0.5 × 0.0721 = 0.03605.
Refer to Supplementary Methods (section S4) for cell culture and transfection, RNA isolation and qRT-PCR, cell proliferation assay, and cell migration assay.
This work was supported by grants from the National Research Foundation of Korea (NRF), which is funded by the Korean Government, the Ministry of Science, ICT & Future Planning (MSIP) (2010-0017662 and 2013M3A9A7046303). This work was also supported by the “Systems biology infrastructure establishment grant” provided by Gwangju Institute of Science and Technology in 2014. We thank Tae-Hwan Kim and Taejun Park for their critical reading and valuable comments.
The authors declare no competing financial interests.
Author Contributions K.-H.C. designed the project and supervised the research; D.L. and K.-H.C. performed the mathematical modeling and analysis; M.K. performed the experiments; and D.L. and K.-H.C. wrote the manuscript.