Cellular responding DNA damage: An improved modeling of P53 gene regulatory networks under ion radiation (IR)
Introduction
Under genome stresses such as DNA damage, hypoxia, and aberrant oncogene signals, cells can trigger their internal defensive mechanisms in fighting against these perturbations from outsides [1], [2], [3]. One critical response is the activation of the tumor suppressor P53, which can control the transcription and translation of series downstream genes. By regulating the complicated signal pathways, P53 induce or repress the expression level of the vital genes within the cell, leading to cell cycle arrest and cell apoptosis in response to genome stresses induced by abnormal signals from outsides [5]. Induction of apoptosis indirectly triggered by P53 is pivotal for eliminating abnormal cells with genome damage by irradiation or chemotherapeutic drugs [4], [5]. Whereas, abnormalities in the P53 tumor suppressor have been identified in over 60% of human cancers and the status of P53 within tumor cells has been proposed to be one of the determinant response to anticancer therapies [3]. Furthermore, controlled radiotherapy studies show the existence of a strong biologic basis for considering P53 status as a radiation predictor [6], [7], therefore, P53 status in tumor cell can be considered as a predictor for long-term biochemical control during and after radiotherapy [7], [8].
The combined approaches of systems analysis, control theory, and computer science can stimulate new approaches to simulate the investigation of the complicated mechanisms of cellular responding genome stresses [9], [10]. These methods provide a good link between the diverging areas of biomedicine and mathematics [11]. Using differential equations and graphic approaches to study various dynamical and kinetic processes of biological systems can provide useful insights, as indicated by many previous studies on a series of important biological topics, such as enzyme-catalyzed reactions, etc. [11], [12], [13].
Recently, several models have been proposed to explain the damped oscillations of P53 in cell populations [14], [15], [16], [17], [18]. However, under continuous effect of acute IR, the complicated mechanisms of cellular response in fighting against DNA damage need to be further addressed at single cell level. To simulate the kinetics of cellular self-defense in response to different IR dose, an improved model of P53 gene regulatory networks is proposed based on the previous models [14], [15], [16], [17], [18] and latest biomedical studies [19], [20]. In our model, the dynamics of double-strand breaks (DSBs) generating and their repairing, ATM and ARF activation, as well as the regulatory kinetics of P53–MDM2 feedback loop are implemented by using a set of differential equations. Meanwhile, the plausible outcomes of cellular responding DNA damage are presented versus continuous radiation time within different IR dose domains.
Section snippets
Model implementation
The scheme of the integrated model is given in Fig. 1. Compared with previous models [14], [15], [16], [17], [18], our model contains more vital components which is most related with cellular response kinetics, such as the over-expression of oncogenes, ARF activation and some signal pathways [20], [21]. As continuous IR is applied into a cell, the resulting DSBs occur stochastically, and then form the DSB-protein complexes (DSBCs) after interacting with the DNA repair proteins around damage
Simulation results
To ensure the rightness of the simulation results, we consider the fact that the valid parameter sets should obey the following rules [10], [14], [15], [16], [17]: (1) the model must contain oscillations. This is important as there has been experimental evidence that oscillations occur between P53 and MDM2 after cell stress; (2) the mechanism used to mathematically describe the degradation of P53 by MDM2 is accurate only for low concentrations of P53; (3) the concentration of P53∗ is much
Conclusion
By using a set of differential equations, we proposed an improved model of P53 gene regulatory network under different IR dose domains. The modules for DSBs generation and repair, ATM and ARF activation, as well as P53–MDM2 feedback loop are implemented. In our simulations, ATM and ARF exhibits a strong sensitivity behavior in response to acute IR, fully consistent with the necessary outcomes to efficiently transfer the stress signal to the P53–MDM2 feedback loop and trigger cellular responding
Acknowledgements
This work is supported in part by Doctorial Foundation from Education Committee (20060255006), China, and National Natural Science Foundation of China (60661003), China. Especially, we thank so much for the helps of Prof. G.C. Chou in Gordon Life Science Institute, USA, and helpful discussions with Prof. G.P. Zhou in the Research Center of Structural Biology Harvard Medical School.
References (26)
- et al.
Molecular interaction map of the P53 and Mdm2 logic elements, which control the Off–On switch of P53 in response to DNA damage
Biochem. Biophys. Res. Commun.
(2005) - et al.
Cellular responses to ionizing radiation damage
Int. J. Radiat. Oncol. Biol. Phys.
(2001) - et al.
The role of P53 in radiation therapy outcomes for favorable-to-intermediate-risk prostate cancer T
Radiat. Oncol. Biol. Phys.
(2002) - et al.
A mathematical model of P53 gene regulatory networks under radiotherapy
BioSystems
(2007) - et al.
ARF and ATM/ATR cooperate in P53-mediated apoptosis upon oncogenic stress
Biochem. Biophys. Res. Commun.
(2005) - et al.
DNA double strand break repair signalling: the case of RAD51 post-translational regulation
Cell. Signal.
(2002) Decision making by P53: life, death and cancer
Cell Death Differ.
(2003)- et al.
Surfing the P53 network
Nature
(2000) - et al.
Principles and Practice of Radiation Oncology
(1998) Predicting response to cancer chemotherapy: the role of P53
Cell Tissue Res.
(1998)
The P53 protein family and radiation sensitivity: yes or no?
Cancer Metastasis Rev.
Cancer: a challenge for control theory and computer modeling
Eur. J. Cancer
Modeling and simulation of genetic regulatory systems: a literature review
Comput. Biol.
Cited by (15)
A mathematical study of the robustness of G2/M regulatory network in response to DNA damage with parameters sensitivity
2014, Applied Mathematics and ComputationCitation Excerpt :With incorporation of ATM, an oscillation behavior occurred in this feedback loop after a sufficiently strong DNA damage signal [7]. The primary phosphorylation of p53 on Ser-15 is induced by Wip1, causing an arrest of the cell cycle through transcriptional activation of p21 [8]. Phosphorylated by Chk1, Cdc25C is exported to the cytoplasm by 14-3-3.
Mathematical modeling of p53 pulses in G2 phase with DNA damage
2014, Applied Mathematics and ComputationCitation Excerpt :Fig. 3 also shows the details of the G2 signaling pathway. To ensure the accuracy of the simulated results and the robustness of the presented model, parameters are estimated according to two rules [19]: (i) the model must contain oscillations; and (ii) the mechanism used to mathematically describe the degradation of P53 by MDM2 is accurate only for low concentration of P53. Existing parameters described in previous reports of Zhang et al. [11] and Batchelor et al. [13] are integrated into our model.
A novel method for modelling cellular response genome stress by combined geneenvironment network (GEN) and kinetic theory framework
2014, IFAC Proceedings Volumes (IFAC-PapersOnline)A mathematical analysis of DNA damage induced G2 phase transition
2013, Applied Mathematics and ComputationCitation Excerpt :This leads to p21 binds to MPF to form p21/MPF complex, and retards the progression of cell cycle [7,21]. The p53-Mdm2 feedback loop formed by p53 and Mdm2 plays an essential role in responding to DNA damage signal [22,23]. With incorporation of ATM, the p53-Mdm2 negative feedback loop leads to an oscillations behavior after a sufficiently strong DNA damage signal.
Stability analysis of standard genetic regulatory networks with time-varying delays and stochastic perturbations
2011, NeurocomputingCitation Excerpt :The oscillation of GRNs is the main reason causing immunity decline that leads to illness. Therefore, the study on the stability of GRNs could further reveal life activity rule from the viewpoint of system biology, and also contribute to both the drug design and disease diagnosis [5–7]. In biological systems especially GRNs, time delays are unavoidable due primarily to the slow process of transcription, translation, translocation as well as the finite switching speed of amplifiers.
Dynamic modeling of cellular response to DNA damage based on p53 stress response networks
2009, Progress in Natural Science