Post-translational regulation enables robust p53 regulation
© Shin et al.; licensee BioMed Central Ltd. 2013
Received: 24 August 2012
Accepted: 23 August 2013
Published: 30 August 2013
The tumor suppressor protein p53 plays important roles in DNA damage repair, cell cycle arrest and apoptosis. Due to its critical functions, the level of p53 is tightly regulated by a negative feedback mechanism to increase its tolerance towards fluctuations and disturbances. Interestingly, the p53 level is controlled by post-translational regulation rather than transcriptional regulation in this feedback mechanism.
We analyzed the dynamics of this feedback to understand whether post-translational regulation provides any advantages over transcriptional regulation in regard to disturbance rejection. When a disturbance happens, even though negative feedback reduces the steady-state error, it can cause a system to become less stable and transiently overshoots, which may erroneously trigger downstream reactions. Therefore, the system needs to balance the trade-off between steady-state and transient errors. Feedback control and adaptive estimation theories revealed that post-translational regulation achieves a better trade-off than transcriptional regulation, contributing to a more steady level of p53 under the influence of noise and disturbances. Furthermore, post-translational regulation enables cells to respond more promptly to stress conditions with consistent amplitude. However, for better disturbance rejection, the p53- Mdm2 negative feedback has to pay a price of higher stochastic noise.
Our analyses suggest that the p53-Mdm2 feedback favors regulatory mechanisms that provide the optimal trade-offs for dynamic control.
KeywordsFeedback control theory p53-Mdm2 feedback loop Robustness Disturbance rejection
Gene networks are constantly subject to noise or fluctuations, which originate from variations in transcription, translation, and environmental conditions. The stochastic nature of gene networks has been the focus of many studies (reviewed in [1, 2]). There are at least three types of fluctuations that affect gene network dynamics: 1) intrinsic noise or fast fluctuations, 2) extrinsic noise or slow fluctuations, and 3) periodic DNA replication-dependent oscillations . Intrinsic noise arises from the inherent randomness during transcription and translation, key processes for gene expression . Extrinsic noise arises from the factors that universally affect the expression of all genes in a given cell, such as variations in the number of RNA polymerase, ribosome, etc. [1, 2]. The third type of fluctuation is due to periodic DNA replication in growing and dividing cells .
Nevertheless, gene networks are usually able to perform their regulatory functions under the influence of such disturbances , which provokes the question: How do they manage to achieve this remarkable robustness? In control theory, it is known that feedback, a situation in which two (or more) dynamical sub-systems are connected in a way that their dynamics are coupled, can make a system resilient towards disturbances [6, 7]. A well-known example of feedback in the context of gene networks is negative auto-regulation, in which a transcription factor represses the transcription of its own gene and reduces the effects of noise exerted on the transcription process [8–10].
The feedback is suppressed to stop disturbance rejection during DNA-damage. Within minutes of exposure to DNA-damaging agents (UV, X-rays, etc.), the p53 protein level increases rapidly without any significant change in the p53 mRNA level, for p53 is stabilized and no longer degraded by Mdm2 at the normal rate [17, 18]. The suppression of the feedback and the stoppage of disturbance rejection allow external factors such as ATM to modulate p53 in response to stress [19, 20], which can result in pulses or sustained oscillation of p53  (Figure 1B). Therefore, the feedback mechanism is adaptive – under normal conditions the feedback rejects disturbance to maintain a low steady level of p53 whereas upon DNA damage, the feedback is inactivated to allow pulses or oscillation  (Figure 1C).
Interestingly, even though the negative feedback is stronger in normal conditions to reject disturbances, it is well known that feedback can also cause instability, which leads to transient fluctuation (overshoot) and oscillation . Both the strength and the delay of a negative feedback can contribute to instability. As previous measurements showed, the transcriptional regulation of Mdm2 by p53 has a noticeable delay  (Figure 1B). Therefore, to maintain a steady p53 level in normal conditions, the p53-Mdm2 feedback faces a dilemma: the stronger the feedback is to reject disturbances at the steady state, the more likely the feedback will become unstable and cause transient or sustained fluctuations. The feedback has to carefully balance the trade-off between steady-state and transient errors for disturbance rejection .
Does the p53-Mdm2 feedback adopt any strategy to optimize this trade-off? It is intriguing that in this negative feedback loop Mdm2 downregulates p53 through a post-translational mechanism (protein degradation), which is not as energy efficient as transcriptional repression because p53 is being produced and actively degraded simultaneously. However, post-translational regulation has distinct dynamic properties and is relatively faster than transcriptional regulation, so it is conceivable that post-translational regulations provides an advantage over transcriptional regulation in terms of the robustness-stability trade-off despite its less energy efficiency. To test this hypothesis, we used techniques from feedback control and adaptive estimation theories to analyze the p53-Mdm2 feedback loop.
Results and Discussion
p53-Mdm2 feedback model
where x(t), y(t), and z(t) represent ATM, Mdm2, and p53 levels, respectively. In Eq. 1, since ATM is inactive under normal conditions, the ATM levels are assumed to be low and constant and the rate of change is zero. Eq. 2 is composed of the p53-dependent production of Mdm2 (first term), ATM-dependent suppression of Mdm2 (second term), and Mdm2 degradation (third term). Eq. 3 is composed of the Mdm2-dependent suppression of p53 (first term) and p53 degradation (second term). Note that all the parameter values are supposed to be positive. Following the practice in the published model , Eq. 3 does not include the constant basal production rate because it has no effect on the frequency domain and disturbance rejection analysis we will perform next . A discrete-time model for the p53-Mdm2 feedback system was also built for parameter estimation, which will be described in later sections of the paper.
Feedback reduces steady-state error
Eq. 3 is rather generic and can represent alternative p53 suppression mechanisms by Mdm2. A transcriptional mechanism would only affect p yz in the first term because the strength of transcriptional suppression only depends on the number of Mdm2 molecules (suppressor) but not on the number of p53 protein molecules. In contrast, a post-translational mechanism would affect both p yz (first term) and p z (second term), because the degradation rate of p53 depends on both the number of p53 molecules and the number of Mdm2 molecules. This generic equation enables us to derive a common set of equations to compare the loop dynamics between transcriptional and translational regulation.
According to Eq. 9, increasing either p yz or p z will reduce the steady-state error, because both terms are only in the denominator. Therefore, either a stronger negative feedback through the Mdm2 suppression of p53 (p yz ) or less stable p53 (p z ) can reduce the steady-state error. Interestingly, as stated previously, post-translational degradation of p53 by Mdm2 increases both p yz and p z (degradation rate), while a hypothetical transcriptional suppression would only increase p yz . Hence post-translational suppression of p53 may be more efficient at reducing steady-state error than transcriptional suppression, even though it is less energy efficient.
Trade-off between steady-state and transient errors
The typical time responses (unit step responses) for a second-order system with different damping ratios are shown in Figure 2B (p z = p y = 0.01 min-1). The %OS, or overshoot of the p53 level, decreases when the damping ratio increases (Figure 2C).
According to Eq. 12, increasing p yz monotonically reduces the damping ratio and increases %OS, for p yz only appears in the denominator. Since, transcriptional suppression can only influence p yz , a transcriptional negative feedback is limited by the trade-off between steady-state error and transient overshoot – increasing the strength of the negative feedback reduces steady-state error at the expense of increasing transient overshoot.
The second-order partial derivative of Eq. 14 shows that ζ reaches the minimum value at . When , ζ decreases with p z ; when , ζ increases with pz. This suggests that when , increasing pz can decrease both the steady-state error and overshoot (which is inverse to the damping ratio). This insight potentially explains why cells choose to spend extra energy producing and then actively degrading p53 under normal conditions – it rejects disturbances by reducing both the steady-state error and the transient overshoot. Furthermore, post-translational suppression influences both p yz and p z , because the Mdm2-mediated degradation rate of p53 depends on both Mdm2 and p53 concentrations. This suggests that post-translational suppression can achieve a better trade-off than transcriptional suppression, which can only influence p yz .
Interestingly, the fact that p z can improve both steady-state error and %OS seems to suggest that cells should keep increasing p z to achieve ever better robustness. However, it is worth noting that there is another trade-off factor that ultimately comes into the picture, which is the energy cost. Higher p z means that the cells have to spend more resources to rapidly produce and then degrade p53 in a “futile” cycle, so the degradation rate cannot be increased indefinitely. The cell has to balance between robustness and energy spending to achieve the optimal operating point.
Altogether, our analytical and quantitative analyses suggest that cells spend extra energy to produce and degrade p53 simultaneously under normal conditions in order to maintain a more robust p53 level. By modulating both p yz and p z , post-translational suppression of p53 enables the p53-Mdm2 feedback to achieve a better trade-off by reducing both steady-state errors and transient overshoots. On the contrary, transcriptional suppression only modulates p yz , which does not improve the overall trade-off.
Robustness analyses with pole plots
To investigate the transient behavior and stability of G(z), we plotted its poles using the MATLAB Robust Control Toolbox (Figure 4B-D), which can be reproduced using the supplementary MATLAB file (Additional file 6). For each map, 100 points were calculated based on Monte Carlo sampling of the estimated parameter range in Figure 4A. We first plotted the poles for the DNA damage condition, under which the Mdm2 suppression of p53 (the negative feedback) is weakened to stop disturbance rejection (Figure 4B). Consistent with the weakening of the feedback, the plot shows that the poles are mostly real and the system is stable, with damping ratios ranging from 0.8 to 1.0.
Western blot measurements of the total and ubiquitinated p53 levels showed that the suppression of p53 is 3.67 fold higher under the normal condition compared to the DNA damage condition . We first increased w yz by 3.67 fold (from 0.1656±0.1169 to 0.6708±0.1169) to evaluate how transcriptional suppression will affect transient overshoot and stability. The increased w yz values shifted poles higher on the map compared to the DNA damage condition (Figure 4C). The bigger imaginary components of the poles indicate that the damping ratios are lower (0.4 – 0.8) and the corresponding %OS values are higher for the system, therefore confirming our previous conclusion that transcriptional suppression reduces the steady-state error at the expense of transient overshoot and stability.
However, if the increase of w yz is accompanied by a decrease of w z (stability of p53) by 3.67 fold (0.3627±0.0825 to 0.0998±0.0825), which approximates post-translational suppression, the poles are located closer to the real axis. The smaller imaginary components of the poles indicate that the damping ratios are greater (the %OS values are smaller). Therefore, post-translational suppression can reduce the steady-state error without a hefty penalty of transient overshoot (%OS) and stability, unlike transcriptional suppression.
Post-translational regulation enables faster responses with predictable level shifts
Validation by a non-linear, mechanistic model
So far, our analyses have been performed using the linear p53-MDM2 feedback model (Eq. 1–3) modified from Geva-Zatorsky et al. . Even though this model matches experimental measurements , it does not capture the non-linear aspects of the negative feedback loop. Therefore, the linear model, and hereby its analyses, can only be applied to the first order approximation.
where, [mdm2 mRNA ], [Mdm2], [p53 mRNA ], and [P53] denote the mRNA and protein levels of p53 and Mdm2. β1 and β3 denotes transcription rates, and β2 and β4 denote translation rates of Mdm2 and p53. α s (s = 1,2,3,4) denote mRNA and protein degradation rates. Transcriptional activation of Mdm2 by p53 is modeled by a Hill function, where n is the hill coefficient and K1 is the dissociation constant. Mdm2-mediated p53 ubiquitination and degradation are modeled by a Michaelis-Menten (MM) function as described by Xu et al. . γ is the reaction rate of p53 ubiquitination by Mdm2, and K2 is the saturation constant.
When γ, the degradation rate of p53 by Mdm2-mediated ubiquitination, increases, the %OS initially increases but then decreases, which allows the system to reduce both the steady-state error and %OS through gamma (Figure 6B). This plot is similar to the previous analysis with the linear model when pz is increased (Figure 3B). Therefore, the mechanistic model confirmed the insight from the linear model that post-translational suppression enables the p53-Mdm2 negative feedback to reduce both steady-state error and %OS.
We explored the parameter space to examine the different operating regions of the non-linear model, and found that the above tradeoff trend is generally preserved even though the absolute values of steady-state error and %OS vary (Figure 6C). This is probably due to the fact that the model is roughly piece-wise linear at each operating point when given a modest disturbance, so that higher-order effects do not dominate the first-order behavior predicted by the linear model.
where τ denotes the time delay between transcription and translation.
Simulations with varying time delay reveal that delay does not affect steady-state error but increases %OS, which is consistent with the proposition that time delay decreases stability of negative feedback loops (Figure 6D). A longer delay also increases the likelihood of sustained oscillation (Figure 6E).
In summary, our analyses indicate that cells maintain robust p53 levels and reject disturbances by simultaneously producing and degrading p53, even though this process is more energy intensive. Within the p53-Mdm2 negative feedback loop, post-translational suppression of p53 by Mdm2 achieves a better trade-off between steady-state and transient errors than transcriptional suppression, which potentially explains why the former has been experimentally observed in the cell. Furthermore, post-translational suppression enables p53 to respond faster to stress conditions with a more predictable level shift. Understanding these nuances allows us to appreciate the complexity of regulatory networks, which will potentially lead to better therapeutics.
Besides ubiquitination, p53 activity is also heavily regulated by phosphorylation and nucleocytoplasmic shuttling [32–35]. Interestingly, these post-translational processes provide similar beneficial tradeoff between steady-state errors and overshoot, while being faster and less energy intensive than protein degradation. It is probably not a coincidence then that p53 has multiple phosphorylation sites and is shuttled alongside Mdm2 between nucleus and cytoplasm. Therefore, protein degradation is only one of several post-translational mechanisms that enhance the robustness of the system.
Increasing degradation rather than reducing production might be a common strategy evolved by biological systems for robustness. For example, hematopoietic stem cells continuously go through apoptosis , which seems energy-inefficient and futile as the protein degradation we discussed. However, our analysis would suggest that regulating apoptosis rather than cell division may enable the stem cell population to become more robust to disturbances and respond faster to changes. The fact that biological systems employ mechanisms for robustness at many different levels raises an interesting question - how does robustness at each level contribute to the overall robustness of the whole system? Undoubtedly challenging, attempts to answer this question will help unravel the underlying design principles of complex biological systems.
Ordinary differential and difference equations were used for physics-based modeling of the p53-Mdm2 feedback loop (see Additional file 1 for the derivation of the p53-Mdm2 discrete-time model from underlying physics). Parameter ranges of the discrete-time model were estimated using the Normalized Least Mean Squares (NLMS) method detailed in Additional file 1. The experimental data  and MATLAB (Mathworks, USA) codes are provided as (Additional files 1, 2, 3, 4 and 5). Steady-state and transient error analysis was performed using Laplace and Z-domain analysis (see Additional file 1 for steady-state error analysis using the discrete-time model). The Monte Carlo method and MATLAB Robust Control Toolbox (Mathworks, USA) were used for the estimation-based robustness analysis (see Additional file 6).
Image extraction and fluorescence quantification
285 Image frames were extracted from the video file  and the fluorescence quantification of p53 and Mdm2 was carried out using the LabVIEW Vision Assistant 2010 (National Instruments, USA). We manually marked the location of each cell nucleus in each frame and 285 data points were obtained for each protein.
This work was supported by National Institutes of Health grant R01GM95990 and National Science Foundation grants NSF 1137269, CCF-1011918 and CCF-0942936.
- Kaern M, Elston TC, Blake WJ, Collins JJ: Stochasticity in gene expression: From theories to phenotypes. Nat Rev Genet. 2005, 6: 451-464. 10.1038/nrg1615.PubMedView ArticleGoogle Scholar
- Raj A, van Oudenaarden A: Nature, nurture, or chance: Stochastic gene expression and its consequences. Cell. 2008, 135: 216-226. 10.1016/j.cell.2008.09.050.PubMedPubMed CentralView ArticleGoogle Scholar
- Rosenfeld N, Young JW, Alon U, Swain PS, Elowitz MB: Gene regulation at the single-cell level. Science. 2005, 307: 1962-1965. 10.1126/science.1106914.PubMedView ArticleGoogle Scholar
- Dunlop MJ, Cox RS, Levine JH, Murray RM, Elowitz MB: Regulatory activity revealed by dynamic correlations in gene expression noise. Nat Genet. 2008, 40: 1493-1498. 10.1038/ng.281.PubMedPubMed CentralView ArticleGoogle Scholar
- Jen E: Robust design : Repertoire of biological, ecological, and engineering case studies. 2005, New York: Oxford University PressGoogle Scholar
- Cosentino C, Bates D: Feedback control in systems biology. 2012, Boca Raton, Florida: CRC PressGoogle Scholar
- Åström KJ, Murray RM: Feedback systems: An introduction for scientists and engineers. 2008, Princeton: Princeton University PressGoogle Scholar
- Becskei A, Serrano L: Engineering stability in gene networks by autoregulation. Nature. 2000, 405: 590-593. 10.1038/35014651.PubMedView ArticleGoogle Scholar
- Simpson ML, Cox CD, Sayler GS: Frequency domain analysis of noise in autoregulated gene circuits. Proc Natl Acad Sci U S A. 2003, 100: 4551-4556. 10.1073/pnas.0736140100.PubMedPubMed CentralView ArticleGoogle Scholar
- Dublanche Y, Michalodimitrakis K, Kummerer N, Foglierini M, Serrano L: Noise in transcription negative feedback loops: Simulation and experimental analysis. Mol Syst Biol. 2006, 2: 41-PubMedPubMed CentralView ArticleGoogle Scholar
- Piette J, Neel H, Marechal V: Mdm2: Keeping p53 under control. Oncogene. 1997, 15: 1001-1010. 10.1038/sj.onc.1201432.PubMedView ArticleGoogle Scholar
- Vogelstein B, Lane D, Levine AJ: Surfing the p53 network. Nature. 2000, 408: 307-310. 10.1038/35042675.PubMedView ArticleGoogle Scholar
- Levine AJ, Oren M: The first 30 years of p53: Growing ever more complex. Nat Rev Cancer. 2009, 9: 749-758.PubMedPubMed CentralView ArticleGoogle Scholar
- Kaina B: DNA damage-triggered apoptosis: Critical role of DNA repair, double-strand breaks, cell proliferation and signaling. Biochem Pharmacol. 2003, 66: 1547-1554. 10.1016/S0006-2952(03)00510-0.PubMedView ArticleGoogle Scholar
- Roos WP, Kaina B: DNA damage-induced cell death by apoptosis. Trends Mol Med. 2006, 12: 440-450. 10.1016/j.molmed.2006.07.007.PubMedView ArticleGoogle Scholar
- Shin YJ, Hencey B, Lipkin SM, Shen X: Frequency domain analysis reveals external periodic fluctuations can generate sustained p53 oscillation. PLoS One. 2011, 6: e22852-10.1371/journal.pone.0022852.PubMedPubMed CentralView ArticleGoogle Scholar
- Maltzman W, Czyzyk L: UV irradiation stimulates levels of p53 cellular tumor antigen in nontransformed mouse cells. Mol Cell Biol. 1984, 4: 1689-94.PubMedPubMed CentralView ArticleGoogle Scholar
- Price BD, Calderwood SK: Increased sequence-specific p53-DNA binding activity after DNA damage is attenuated by phorbol esters. 1993, BASINGSTOKE: ONCOGENE, 8: 3055Google Scholar
- Kastan MB, Zhan Q, El-Deiry WS, Carrier F, Jacks T, et al.: A mammalian cell cycle checkpoint pathway utilizing p53 and GADD45 is defective in ataxia-telangiectasia. Cell. 1992, 71: 587-597. 10.1016/0092-8674(92)90593-2.PubMedView ArticleGoogle Scholar
- Gannon HS, Woda BA, Jones SN: ATM phosphorylation of Mdm2 Ser394 regulates the amplitude and duration of the DNA damage response in mice. Cancer Cell. 2012, 21: 668-679. 10.1016/j.ccr.2012.04.011.PubMedPubMed CentralView ArticleGoogle Scholar
- Geva-Zatorsky N, Rosenfeld N, Itzkovitz S, Milo R, Sigal A, et al.: Oscillations and variability in the p53 system. Mol Syst Biol. 2006, 2: 2006.0033-PubMedPubMed CentralView ArticleGoogle Scholar
- Chandra FA, Buzi G, Doyle JC: Glycolytic oscillations and limits on robust efficiency. Science Science. 2011, 333: 187-192.PubMedGoogle Scholar
- Geva-Zatorsky N, Dekel E, Batchelor E, Lahav G, Alon U: Fourier analysis and systems identification of the p53 feedback loop. Proc Natl Acad Sci U S A. 2010, 107: 13550-13555. 10.1073/pnas.1001107107.PubMedPubMed CentralView ArticleGoogle Scholar
- Shin YJ, Bleris L: Linear control theory for gene network modeling. PLoS One. 2010, 5: e12785-10.1371/journal.pone.0012785.PubMedPubMed CentralView ArticleGoogle Scholar
- Stommel JM, Wahl GM: Accelerated MDM2 auto-degradation induced by DNA-damage kinases is required for p53 activation. EMBO J. 2004, 23: 1547-56. 10.1038/sj.emboj.7600145.PubMedPubMed CentralView ArticleGoogle Scholar
- Ju J, Schmitz JC, Song B, Kudo K, Chu E: Regulation of p53 expression in response to 5-fluorouracil in human cancer RKO cells. Clin Cancer Res. 2007, 13: 4245-4251. 10.1158/1078-0432.CCR-06-2890.PubMedView ArticleGoogle Scholar
- Shin YJ, Sayed AH, Shen X: Adaptive models for gene networks. PloS one. 2012, 7: e31657-10.1371/journal.pone.0031657.PubMedPubMed CentralView ArticleGoogle Scholar
- Sayed AH: Adaptive Filters. 2008, New Jersey: WilleyView ArticleGoogle Scholar
- Franklin GF, Powell JD: Digital control of dynamic systems. 1980, Reading, Mass.: Addison-Wesley Pub. Co.Google Scholar
- Cheng Q, Cross B, Li B, Chen L, Li Z, et al.: Regulation of MDM2 E3 ligase activity by phosphorylation after DNA damage. Mol Cell Biol. 2011, 31: 4951-4963. 10.1128/MCB.05553-11.PubMedPubMed CentralView ArticleGoogle Scholar
- Xu L, Qu Z: Roles of protein ubiquitination and degradation kinetics in biological oscillations. PloS one. 2012, 7 (4): e34616-10.1371/journal.pone.0034616.PubMedPubMed CentralView ArticleGoogle Scholar
- O'Keefe K, Li H, Zhang Y: Nucleocytoplasmic shuttling of p53 is essential for MDM2-mediated cytoplasmic degradation but not ubiquitination. Mol Cell Biol. 2003, 23: 6396-405. 10.1128/MCB.23.18.6396-6405.2003.PubMedPubMed CentralView ArticleGoogle Scholar
- Tao W, Levine AJ: Nucleocytoplasmic shuttling of oncoprotein Hdm2 is required for Hdm2-mediated degradation of p53. Proc Natl Acad Sci. 1999, 96: 3077-3080. 10.1073/pnas.96.6.3077.PubMedPubMed CentralView ArticleGoogle Scholar
- Ashcroft M, Kubbutat MH, Vousden KH: Regulation of p53 function and stability by phosphorylation. Mol Cell Biol. 1999, 19: 1751-8.PubMedPubMed CentralView ArticleGoogle Scholar
- Kazuyasu S, Herrera JE, Shin’ichi S, Toru M, Michael B, Alex V, Anderson CW, Anderson CW, Ettore A: DNA damage activates p53 through a phosphorylation–acetylation cascade. Genes Dev. 1998, 12 (18): 2831-2841. 10.1101/gad.12.18.2831.View ArticleGoogle Scholar
- Alenzi FQ, Ahmad SY, Alenazi BQ, Salem ML, Al-Jabri AA, et al.: The haemopoietic stem cell: Between apoptosis and self renewal. Yale J Biol Med. 2009, 82: 7-18.PubMedPubMed CentralGoogle Scholar