The role of time delays in P53 gene regulatory network stimulated by growth factor

Interplay of p53 and XIAP protein dynamics orchestrates cell fate in response to chemotherapy

Chemotherapeutic drugs are used to treat almost all types of cancer, but the intended response, i.e., elimination, is often incomplete, with a subset of cancer cells resisting treatment. Two critical factors play a role in chemoresistance: the p53 tumour suppressor gene and the X-linked inhibitor of apoptosis (XIAP). These proteins have been shown to act synergistically to elicit cellular responses upon DNA damage induced by chemotherapy, yet, the mechanism is poorly understood. This study introduces a mathematical model characterising the apoptosis pathway activation by p53 before and after mitochondrial outer membrane permeabilisation upon treatment with the chemotherapy Doxorubicin (Dox). "In-silico" simulations show that the p53 dynamics change dose-dependently. Under medium to high doses of Dox, p53 concentration ultimately stabilises to a high level regardless of XIAP concentrations. However, caspase-3 activation may be triggered or not depending on the XIAP induction rate, ultimately determining whether the cell will perish or resist. Consequently, the model predicts that failure to activate apoptosis in some cancer cells expressing wild-type p53 might be due to heterogeneity between cells in upregulating the XIAP protein, rather than due to the p53 protein concentration. Our model suggests that the interplay of the p53 dynamics and the XIAP induction rate is critical to determine the cancer cells' therapeutic response.

Mathematical Modelling of p53 Signalling during DNA Damage Response: A Survey

No gene has garnered more interest than p53 since its discovery over 40 years ago. In the last two decades, thanks to seminal work from Uri Alon and Ghalit Lahav, p53 has defined a truly synergistic topic in the field of mathematical biology, with a rich body of research connecting mathematic endeavour with experimental design and data. In this review we survey and distill the extensive literature of mathematical models of p53. Specifically, we focus on models which seek to reproduce the oscillatory dynamics of p53 in response to DNA damage. We review the standard modelling approaches used in the field categorising them into three types: time delay models, spatial models and coupled negative-positive feedback models, providing sample model equations and simulation results which show clear oscillatory dynamics. We discuss the interplay between mathematics and biology and show how one informs the other; the deep connections between the two disciplines has helped to develop our understanding of this complex gene and paint a picture of its dynamical response. Although yet more is to be elucidated, we offer the current state-of-the-art understanding of p53 response to DNA damage.