- Research article
- Open Access

# Post-translational regulation enables robust p53 regulation

- Yong-Jun Shin
^{1}, - Kai-Yuan Chen
^{1}, - Ali H Sayed
^{2}, - Brandon Hencey
^{3}and - Xiling Shen
^{1}Email author

**7**:83

https://doi.org/10.1186/1752-0509-7-83

© Shin et al.; licensee BioMed Central Ltd. 2013

**Received:**24 August 2012**Accepted:**23 August 2013**Published:**30 August 2013

## Abstract

### Background

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.

### Results

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.

### Conclusions

Our analyses suggest that the p53-Mdm2 feedback favors regulatory mechanisms that provide the optimal trade-offs for dynamic control.

## Keywords

- Feedback control theory
- p53-Mdm2 feedback loop
- Robustness
- Disturbance rejection

## Background

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 [3]. Intrinsic noise arises from the inherent randomness during transcription and translation, key processes for gene expression [4]. 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 [3].

Nevertheless, gene networks are usually able to perform their regulatory functions under the influence of such disturbances [5], 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 [21] (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 [16] (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 [7]. 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 [21] (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 [22].

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 [23], 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 [23]. 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.

*X*(

*s*),

*Y*(

*s*), and

*Z*(

*s*) denote the Laplace transforms of

*x*(

*t*),

*y*(

*t*), and

*z*(

*t*) respectively.

*E*(

*s*) is the difference, or error, between the input,

*p*

_{ xy }

*X*(

*s*), and the output,

*p*

_{ zy }

*Z(s*).

*D*(

*s*) represents any disturbance exerted on p53.

*G*

_{ Y }(

*s*) and

*G*

_{ Z }(

*s*) are the transfer functions that represent the Mdm2 and p53 systems as we showed previously [24]:

*E*(

*s*) can be expressed as [16]:

*Z*(

*s*) can be expressed as:

*Z*(

*s*) into Eq. 5, we obtain:

*E(s)*from

*D*(

*s*), the Laplace transform of the disturbance signal. We denote this term as

*E*

_{ D }(

*s*) and its corresponding time domain function as

*e*

_{ D }(

*t*):

*D*(s) = 1/s), we can determine the steady-state error due to the disturbance as follows:

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

*ζ*), we derive a second-order transfer function representative of the block diagram in Figure 2A to determine its damping ratio (

*ζ*). Substituting Eq. 5 for

*E*(

*s*) in Eq. 6, the transfer function

*G*(

*s*), which directly relates the input

*X*(

*s*) to the output

*Z*(

*s*), can be expressed as:

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.

*p*

_{ z }affect the overshoot? If we take the partial derivative of the damping ratio in regard to

*p*

_{ z }, we get:

The second-order partial derivative of Eq. 14 shows that *ζ* reaches the minimum value at ${p}_{z}={p}_{y}-\frac{2{p}_{\mathit{yz}}{p}_{\mathit{zy}}}{{p}_{y}}({p}_{y},{p}_{z}>0\to \underset{y}{\overset{2}{p}}>2{p}_{\text{zy}}{p}_{\text{zy}})$. When ${p}_{z}<{p}_{y}-\frac{2{p}_{\mathit{yz}}{p}_{\mathit{zy}}}{{p}_{y}}$, *ζ* decreases with *p*_{
z
}; when ${p}_{z}>{p}_{y}-\frac{2{p}_{\mathit{yz}}{p}_{\mathit{zy}}}{{p}_{y}}$, *ζ* increases with p_{z}. This suggests that when ${p}_{z}>{p}_{y}-\frac{2{p}_{\mathit{yz}}{p}_{\mathit{zy}}}{{p}_{y}}$, increasing p_{z} 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
}.

*p*

_{ yz }= 0.8 h

^{-1}[23],

*p*

_{ zy }= 0.8 h

^{-1}[23],

*p*

_{ y }= 2.0 h

^{-1}(half-life ≈ 20 min) [25],

*p*

_{ z }= 2.0 h

^{-1}(half-life ≈ 20 min) [26]). As

*p*

_{ yz }increases from 0.5 to 3.0, the steady-state error decreases while the %OS increases (the damping ratio decreases) (Figure 3A, left panel). Hence there is a trade-off between steady-state error and transient overshoot (Figure 3A, right panel), consistent with the previous analysis. Any further decrease of the steady-state error from the system operating point (

*p*

_{ yz }= 0.8) will have to pay a hefty penalty for transient overshoot. On the other hand, as

*p*

_{ z }is increased from 0.5 to 3.0, the %OS initially increases but then decreases (Figure 3B, left panel), exactly as Eqs. 14 and 15 have indicated. The system operates at a point (

*p*

_{ y }= 2.0,

*p*

_{ z }= 2.0) where ${p}_{z}\ge {p}_{y}-\frac{2{p}_{\mathit{yz}}{p}_{\mathit{zy}}}{{p}_{y}}$, which allows the system to reduce both the steady-state error and %OS through

*p*

_{ z }(Figure 3B, right panel). The combined effect of

*p*

_{ yz }and

*p*

_{ z }is shown in Figure 3C and D. The steady-state error can be decreased by either increasing

*p*

_{ yz }or

*p*

_{ z }(Figure 3C), but increasing

*p*

_{ yz }can increase %OS (Figure 3D), as shown by the arrows in the figures. Therefore, by modulating both

*p*

_{ yz }and

*p*

_{ z }, post-translational suppression can achieve better steady-state error and %OS than transcriptional suppression, which only modulates

*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

*p*

_{ yz }and

*p*

_{ z }in Figure 3 demonstrated the robustness trade-off between steady-state and transient disturbance rejection, but it will be useful to examine the trade-off within the range of parameter fluctuation that a realistic cell has to experience. To estimate the time-varying parameter values from time-series data of the p53-Mdm2 network [21], we constructed a corresponding discrete-time model, which can be described by Eqs. 16–18:

*i*), y(

*i*), and z(

*i*) represent ATM, Mdm2, and p53 levels respectively. As before, the ATM level is assumed to be low and constant under normal conditions.

*w*

_{ xy }represents the suppression of Mdm2 by ATM

*, w*

_{ yz }represents the suppression of p53 by Mdm2,

*w*

_{ zy }represents the transcriptional activation of Mdm2 by p53,

*w*

_{ y }represents the stability of Mdm2 (1-

*w*

_{ y }represents Mdm2 degradation), and

*w*

_{ z }represents the stability of p53 (1-

*w*

_{ z }represents p53 degradation). Note that

*w*

_{ y }and

*w*

_{ z }represent stability rather than degradation, opposite to

*p*

_{ y }and

*p*

_{ z }. Block diagram representation of the discrete-time feedback model is shown in Figure 2D. Consistent with the continuous model, the discrete-time model confirms that the steady-state error can be decreased by increasing

*w*

_{ yz }or decreasing

*w*

_{ z }(Eq. 19, also see Eq. S11 in Additional file 1). Note that decreasing

*w*

_{ z }means increasing p53 degradation and corresponds to increasing

*p*

_{ yz }.

*w*

_{ yz }and

*w*

_{ z }affect %OS.

*G*(

*z*), which represents the feedback system shown as the block diagram in Figure 2D (see Additional file 1):

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 [30]. 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

*y*(t) has a constant value of

*Y*. Eq. 21 suggests that

*p*

_{ z }but not

*p*

_{ yz }determines the response time (the time needed for

*z*(

*t*) to reach the half steady-state value) (Figure 5A and B). However, the increase in speed comes at a cost; it decreases the response amplitude as shown in Figure 5B. Interestingly, increasing both

*p*

_{ yz }and

*p*

_{ z }allows the system to achieve a faster response time with constant steady-state amplitude (Figure 5C), because the respective increases of

*p*

_{ yz }and

*p*

_{ z }offset each other at the steady level (Eq. 21). Therefore, by modulating both

*p*

_{ yz }and

*p*

_{ z }, post-translational suppression can generate more prompt and consistent p53 response to stress conditions. Indeed, cells respond to DNA damage by modulating the post-translational degradation of p53.

### 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. [23]. Even though this model matches experimental measurements [23], 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 *K*_{1} is the dissociation constant. Mdm2-mediated p53 ubiquitination and degradation are modeled by a Michaelis-Menten (MM) function as described by Xu et al. [31]. *γ* is the reaction rate of p53 ubiquitination by Mdm2, and *K*_{2} is the saturation constant.

_{old}is the p53 steady-state level before perturbation, and SS

_{new}is the p53 steady-state level after perturbation (Figure 6A).

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 p_{z} 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).

### Stochastic simulation

*γ*, the degradation rate of Mdm2-ubiquitinated p53, while keeping the p53 level constant, to investigate how the negative feedback affects p53 variance (Figure 7C). The stochastic noise was then calculated as the variance normalized by the mean (Figure 7D). These stochastic simulations reveal that stochastic noise on p53 increases when

*γ*increases. This suggests that while Mdm2-mediated p53 degradation improves disturbance rejection, it pays the price of amplifying stochastic noise. Hence, the p53-Mdm2 negative feedback loop has to balance the trade-off between disturbance rejection and stochastic noise.

## Conclusions

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 [36], 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.

## Methods

### Computational methods

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 [21] 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 [21] 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.

## Declarations

### Acknowledgements

This work was supported by National Institutes of Health grant R01GM95990 and National Science Foundation grants NSF 1137269, CCF-1011918 and CCF-0942936.

## Authors’ Affiliations

## References

- 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

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.