Recurrent initiation: a mechanism for triggering p53 pulses in response to DNA damage

Chromosome Instability in Fanconi Anemia: From Breaks to Phenotypic Consequences

Fanconi anemia (FA), a chromosomal instability syndrome, is caused by inherited pathogenic variants in any of 22 FANC genes, which cooperate in the FA/BRCA pathway. This pathway regulates the repair of DNA interstrand crosslinks (ICLs) through homologous recombination. In FA proper repair of ICLs is impaired and accumulation of toxic DNA double strand breaks occurs. To repair this type of DNA damage, FA cells activate alternative error-prone DNA repair pathways, which may lead to the formation of gross structural chromosome aberrations of which radial figures are the hallmark of FA, and their segregation during cell division are the origin of subsequent aberrations such as translocations, dicentrics and acentric fragments. The deficiency in DNA repair has pleiotropic consequences in the phenotype of patients with FA, including developmental alterations, bone marrow failure and an extreme risk to develop cancer. The mechanisms leading to the physical abnormalities during embryonic development have not been clearly elucidated, however FA has features of premature aging with chronic inflammation mediated by pro-inflammatory cytokines, which results in tissue attrition, selection of malignant clones and cancer onset. Moreover, chromosomal instability and cell death are not exclusive of the somatic compartment, they also affect germinal cells, as evidenced by the infertility observed in patients with FA.

Wip1 phosphatase deficiency impairs spatial learning and memory

Spatial learning and memory are typically assessed to evaluate hippocampus-dependent cognitive and memory functions in vivo. Protein phosphorylation and dephosphorylation by kinases and phosphatases play critical roles in spatial learning and memory. Here we report that the Wip1 phosphatase is essential for spatial learning, with knockout mice lacking Wip1 phosphatase exhibiting dysfunctional spatial cognition. Aberrant phosphorylation of the Wip1 substrates p38, ATM, and p53 were observed in the hippocampi of Wip1^-/- mice, but only p38 inhibition reversed impairments in long-term potentiation in Wip1-knockout mice. p38 inhibition consistently ameliorated the spatial learning dysfunction caused by Wip1 deficiency. Our results demonstrate that deletion of Wip1 phosphatase impairs hippocampus-dependent spatial learning and memory, with aberrant downstream p38 phosphorylation involved in this process and providing a potential therapeutic target.

Rucaparib Treatment Alters p53 Oscillations in Single Cells to Enhance DNA-Double-Strand-Break-Induced Cell Cycle Arrest

DNA double strand breaks induce oscillatory expression of the transcription factor p53 that is dependent on ataxia telangiectasia mutated (ATM) activity and the rate of double strand break resolution. Although p53 dynamics are known to play a role in the regulation of cell fate determination, the consequences of the variability in dynamics associated with differences in repair rates and utilized repair pathways are unknown. Using single-cell time-lapse microscopy, we found that disruption of specific repair pathways has distinct impacts on p53 dynamics. The small-molecule rucaparib, an inhibitor of the alternative end-joining-associated protein poly (ADP-ribose) polymerase (PARP), increased p53 pulse duration, altering the temporal expression of multiple p53 target genes. As a result, combination treatments of the radiomimetic drug neocarzinostatin with rucaparib drove prolonged growth arrest beyond that of DNA damage alone. This study highlights how pharmacological manipulation of DNA repair pathways may be used to alter p53 dynamics to enhance therapeutic regimens.

Towards DNA-damage induced autophagy: A Boolean model of p53-induced cell fate mechanisms

How a cell determines a given phenotype upon damaged DNA is an open problem. Cell fate decisions happen at cell cycle checkpoints and it is becoming clearer that the p53 pathway is a major regulator of cell fate decisions involving apoptosis or senescence upon DNA damage, especially at G1/S. However, recent results suggest that this pathway is also involved in autophagy induction upon DNA damage. To our knowledge, in this work we propose the first model of the DNA damage-induced G1/S checkpoint contemplating the decision between three phenotypes: apoptosis, senescence, and autophagy. The Boolean model is proposed based on experiments with U87 glioblastoma cells using the transfection of miR-16 that can induce a DNA damage response. The wild-type case of the model shows that DNA damage induces the checkpoint and the coexistence of the three phenotypes (tristable dynamics), each with a different probability. We also predict that the positive feedback involving ATM, miR-16, and Wip1 has an influence on the tristable state. The model predictions were compared to experiments of gain and loss of function in other three different cell lines (MCF-7, A549, and U2OS) presenting agreement. For p53-deficient cell lines such as HeLa, H1299, and PC-3, our model contemplates the experimental observation that the alternative AMPK pathway can compensate this deficiency. We conclude that at the G1/S checkpoint the p53 pathway (or, in its absence, the AMPK pathway) can regulate the induction of different phenotypes in a stochastic manner in the U87 cell line and others.

Oscillatory Dynamics of p53-Mdm2 Circuit in Response to DNA Damage Caused by Ionizing Radiation

Although the dynamical behavior of the p53-Mdm2 loop has been extensively studied, the understanding of the mechanism underlying the regulation of this pathway still remains limited. Herein, we developed an integrated model with five basic components and three ubiquitous time delays for the p53-Mdm2 interaction in response to DNA damage following ionizing radiation (IR). We showed that a sufficient amount of activated ATM level can initiate the p53 oscillations with nearly the same amplitude over a wide range of the ATM level; a proper range of p53 level is also required for generating the oscillations, for too high or too low levels it would fail to generate the oscillations; and increased Mdm2 level leads to decreased amplitude of the p53 oscillation and reduced expression of the p53 activity. Moreover, we found that the negative feedback loop formed between p53 and nuclear Mdm2 plays a dominant role in determining the p53 dynamics, whereas when interaction strength of the negative feedback loop becomes weaker, the positive feedback loop formed between p53 and cytoplasmatic Mdm2 can induce different types of dynamics. Furthermore, we demonstrated that the total time delay required for protein production and nuclear translocation of Mdm2 can induce p53 oscillations even when the p53 level is at a certain stable high steady state or at a certain stable low steady state. In addition, the two important features of the oscillatory dynamics-amplitude and period-can be controlled by such time delay. These results are in agreement with multiple experimental observations and may enrich our understanding of the dynamics of the p53 network.

A Switch in p53 Dynamics Marks Cells That Escape from DSB-Induced Cell Cycle Arrest

Cellular responses to stimuli can evolve over time, resulting in distinct early and late phases in response to a single signal. DNA damage induces a complex response that is largely orchestrated by the transcription factor p53, whose dynamics influence whether a damaged cell will arrest and repair the damage or will initiate cell death. How p53 responses and cellular outcomes evolve in the presence of continuous DNA damage remains unknown. Here, we have found that a subset of cells switches from oscillating to sustained p53 dynamics several days after undergoing damage. The switch results from cell cycle progression in the presence of damaged DNA, which activates the caspase-2-PIDDosome, a complex that stabilizes p53 by inactivating its negative regulator MDM2. This work defines a molecular pathway that is activated if the canonical checkpoints fail to halt mitosis in the presence of damaged DNA.

Dynamical gene regulatory networks are tuned by transcriptional autoregulation with microRNA feedback

Concepts from dynamical systems theory, including multi-stability, oscillations, robustness and stochasticity, are critical for understanding gene regulation during cell fate decisions, inflammation and stem cell heterogeneity. However, the prevalence of the structures within gene networks that drive these dynamical behaviours, such as autoregulation or feedback by microRNAs, is unknown. We integrate transcription factor binding site (TFBS) and microRNA target data to generate a gene interaction network across 28 human tissues. This network was analysed for motifs capable of driving dynamical gene expression, including oscillations. Identified autoregulatory motifs involve 56% of transcription factors (TFs) studied. TFs that autoregulate have more interactions with microRNAs than non-autoregulatory genes and 89% of autoregulatory TFs were found in dual feedback motifs with a microRNA. Both autoregulatory and dual feedback motifs were enriched in the network. TFs that autoregulate were highly conserved between tissues. Dual feedback motifs with microRNAs were also conserved between tissues, but less so, and TFs regulate different combinations of microRNAs in a tissue-dependent manner. The study of these motifs highlights ever more genes that have complex regulatory dynamics. These data provide a resource for the identification of TFs which regulate the dynamical properties of human gene expression.

Model-based optimization of combination protocols for irradiation-insensitive cancers

Alternations in the p53 regulatory network may render cancer cells resistant to the radiation-induced apoptosis. In this theoretical study we search for the best protocols combining targeted therapy with radiation to treat cancers with wild-type p53, but having downregulated expression of PTEN or overexpression of Wip1 resulting in resistance to radiation monotherapy. Instead of using the maximum tolerated dose paradigm, we exploit stochastic computational model of the p53 regulatory network to calculate apoptotic fractions for both normal and cancer cells. We consider combination protocols, with irradiations repeated every 12, 18, 24, or 36 h to find that timing between Mdm2 inhibitor delivery and irradiation significantly influences the apoptotic cell fractions. We assume that uptake of the inhibitor is higher by cancer than by normal cells and that cancer cells receive higher irradiation doses from intersecting beams. These two assumptions were found necessary for the existence of protocols inducing massive apoptosis in cancer cells without killing large fraction of normal cells neighboring tumor. The best found protocols have irradiations repeated every 24 or 36 h with two inhibitor doses per irradiation cycle, and allow to induce apoptosis in more than 95% of cancer cells, killing less than 10% of normal cells.

Wild-type p53-induced phosphatase 1 promotes vascular smooth muscle cell proliferation and neointima hyperplasia after vascular injury via p-adenosine 5'-monophosphate-activated protein kinase/mammalian target of rapamycin complex 1 pathway

OBJECTIVES

Vascular smooth muscle cell (VSMC) proliferation is a crucial cause of vascular neointima hyperplasia and restenosis, thus limiting the long-term efficacy of percutaneous vascular intervention. We explored the role of wild-type p53-induced phosphatase 1 (Wip1), a potent regulator of tumorigenesis and atherosclerosis, in VSMC proliferation and neointima hyperplasia.

METHODS AND RESULTS

Animal model of vascular restenosis was established in wild type C57BL/6J and VSMC-specific Tuberous Sclerosis 1 (TSC1)-knockdown mice by wire injury. We observed increased protein levels of Wip1, phospho (p)-S6 Ribosomal Protein (S6), p-4EBP1 but decreased p-adenosine 5'-monophosphate-activated protein kinase (AMPK)α both in carotid artery at day 28 after injury and in VSMCs after 48 h of platelet derived growth factor-BB (PDGF-BB) treatment. By using hematoxylin-eosin staining, Ki-67 immunohistochemical staining, cell counting kit-8 assay and Ki-67 immunofluorescence staining, we found Wip1 antagonist GSK2830371 (GSK) or mammalian target of rapamycin complex 1 (mTORC1) inhibitor rapamycin both obviously reversed the neointima formation and VSMC proliferation induced by wire injury and PDGF-BB, respectively. GSK also reversed the increase in mRNA level of Collagen I after wire injury. However, GSK had no obvious effects on VSMC migration induced by PDGF-BB. Simultaneously, TSC1 knockdown as well as AMPK inhibition by Compound C abolished the vascular protective and anti-proliferative effects of Wip1 inhibition. Additionally, suppression of AMPK also reversed the declined mTORC1 activity by GSK.

CONCLUSION

Wip1 promotes VSMC proliferation and neointima hyperplasia after wire injury via affecting AMPK/mTORC1 pathway.

Oscillatory behavior of p53-Mdm2 system driven by transcriptional and translational time delays

Biological experiments clarify that p53-Mdm2 module is the core of tumor network and p53 oscillation plays an important role in determining the tumor cell fate. In this paper, we investigate the effect of time delay on the oscillatory behavior induced by Hopf bifurcation in p53-Mdm2 system. First, the stability of the unique positive equilibrium point and the existence of Hopf bifurcation are investigated by using the time delay as the bifurcation parameter and by applying the bifurcation theory. Second, the explicit criteria determining the direction of Hopf bifurcation and the stability of bifurcating periodic solutions are developed based on the normal form theory and the center manifold theorem. In addition, the combination of numerical simulation results and theoretical calculation results indicates that time delays in p53-Mdm2 system are critical for p53 oscillations. The results may help us to better understand the biological functions of p53 pathway and provide clues for treatment of cancer.

A systematic approach to decipher crosstalk in the p53 signaling pathway using single cell dynamics

The transcription factors NF-κB and p53 are key regulators in the genotoxic stress response and are critical for tumor development. Although there is ample evidence for interactions between both networks, a comprehensive understanding of the crosstalk is lacking. Here, we developed a systematic approach to identify potential interactions between the pathways. We perturbed NF-κB signaling by inhibiting IKK2, a critical regulator of NF-κB activity, and monitored the altered response of p53 to genotoxic stress using single cell time lapse microscopy. Fitting subpopulation-specific computational p53 models to this time-resolved single cell data allowed to reproduce in a quantitative manner signaling dynamics and cellular heterogeneity for the unperturbed and perturbed conditions. The approach enabled us to untangle the integrated effects of IKK/ NF-κB perturbation on p53 dynamics and thereby derive potential interactions between both networks. Intriguingly, we find that a simultaneous perturbation of multiple processes is necessary to explain the observed changes in the p53 response. Specifically, we show interference with the activation and degradation of p53 as well as the degradation of Mdm2. Our results highlight the importance of the crosstalk and its potential implications in p53-dependent cellular functions.

Cells can respond to external and internal inputs by transducing information to the nucleus where transcription factors initiate corresponding cellular responses. Cellular signaling is mediated by several pathways; molecular networks that can interact with each other, which alters signal processing and modulates cellular responses. As deregulated signaling can lead to the development of tumors it is important to understand not only how signaling pathways function but also the contribution of their interaction on the signaling dynamics. Here, we analyzed the interplay of the IKK/ NF-κB and p53 pathway, which are both critical for the cellular response to DNA damage and have been implicated in tumor development. To systematically identify interaction points between both pathways, we perturbed IKK/ NF-κB signaling and tracked the changes in the response of p53 to DNA damage. Using computational methods, we show that several reactions in the p53 pathway are simultaneously affected by NF-κB signaling and that this combined action is necessary to explain altered behaviour of the p53 pathway. Hence, our results provide new insights into the interplay between the NF-κB and p53 pathway and help to gain a more comprehensive understanding of the crosstalk.

Frequency switching between oscillatory homeostats and the regulation of p53

Homeostasis is an essential concept to understand the stability of organisms and their adaptive behaviors when coping with external and internal assaults. Many hormones that take part in homeostatic control come in antagonistic pairs, such as glucagon and insulin reflecting the inflow and outflow compensatory mechanisms to control a certain internal variable, such as blood sugar levels. By including negative feedback loops homeostatic controllers can exhibit oscillations with characteristic frequencies. In this paper we demonstrate the associated frequency changes in homeostatic systems when individual controllers -in a set of interlocked feedback loops- gain control in response to environmental changes. Taking p53 as an example, we show how Per2, ATM and Mdm2 feedback loops -interlocked with p53- gain individual control in dependence to the level of DNA damage, and how each of these controllers provide certain functionalities in their regulation of p53. In unstressed cells, the circadian regulator Per2 ensures a basic p53 level to allow its rapid up-regulation in case of DNA damage. When DNA damage occurs the ATM controller increases the level of p53 and defends it towards uncontrolled degradation, which despite DNA damage, would drive p53 to lower values and p53 dysfunction. Mdm2 on its side keeps p53 at a high but sub-apoptotic level to avoid premature apoptosis. However, with on-going DNA damage the Mdm2 set-point is increased by HSP90 and other p53 stabilizers leading finally to apoptosis. An emergent aspect of p53 upregulation during cell stress is the coordinated inhibition of ubiquitin-independent and ubiquitin-dependent degradation reactions. Whether oscillations serve a function or are merely a by-product of the controllers are discussed in view of the finding that homeostatic control of p53, as indicated above, does in principle not require oscillatory homeostats.

Multi-scale modeling identifies the role of p53-Gys2 negative feedback loop in cellular homeostasis

The transcription factor p53 is a tumor suppressor and strictly controlled p53 protein abundance coordinates cellular outcomes in response to various stresses. The glycogen synthase 2 (Gys2) and p53 generates a novel negative feedback circuit in which p53 represses Gys2 expression whereas Gys2 can stabilize p53 by competitive binding with MDM2. However, the dynamic role of p53-Gys2 negative feedback is still elusive. In current work, we recapitulated the main experimental findings using multi-scale modeling and emphasized the pivotal role of p53-Gys2 negative feedback loop to main cellular homeostasis. The multi-scale modeling strategy was used to simulate both in vitro and in vivo experimental findings. We found that expression of a key oncoprotein HBx may facilitate cancer progression. Gys2 overexpression can inhibit hepatocellular carcinoma progression whereas Gys2 knockdown advanced cancer development. We also applied oscillatory and impulse disturbance to p53 signaling pathway and the results showed that optimal p53-Gys2 negative feedback loop was highly resistant to oscillatory or impulse disturbances. Instead, the canonical p53-MDM2 negative feedback circuit can significantly affect the dynamics of p53 and therefore effectively shaped pulsatile patterns. Therefore, the dual negative feedback loops in p53 signaling can provide features of both robustness and tunability. These dynamic features are critical for cellular homeostasis against tumor progression in p53 signaling pathway.

Sporadic activation of an oxidative stress-dependent NRF2-p53 signaling network in breast epithelial spheroids and premalignancies

Breast and mammary epithelial cells experience different local environments during tissue development and tumorigenesis. Microenvironmental heterogeneity gives rise to distinct cell regulatory states whose identity and importance are just beginning to be appreciated. Cellular states diversify when clonal three-dimensional (3D) spheroids are cultured in basement membrane, and one such state is associated with stress tolerance and poor response to anticancer therapeutics. Here, we found that this state was jointly coordinated by the NRF2 and p53 pathways, which were costabilized by spontaneous oxidative stress within 3D cultures. Inhibition of NRF2 or p53 individually disrupted some of the transcripts defining the regulatory state but did not yield a notable phenotype in nontransformed breast epithelial cells. In contrast, combined perturbation prevented 3D growth in an oxidative stress-dependent manner. By integrating systems models of NRF2 and p53 signaling in a single oxidative stress network, we recapitulated these observations and made predictions about oxidative stress profiles during 3D growth. NRF2 and p53 signaling were similarly coordinated in normal breast epithelial tissue and hormone-negative ductal carcinoma in situ lesions but were uncoupled in triple-negative breast cancer (TNBC), a subtype in which p53 is usually mutated. Using the integrated model, we correlated the extent of this uncoupling in TNBC cell lines with the importance of NRF2 in the 3D growth of these cell lines and their predicted handling of oxidative stress. Our results point to an oxidative stress tolerance network that is important for single cells during glandular development and the early stages of breast cancer.

Coupling of cell fate selection model enhances DNA damage response and may underlie BE phenomenon

Double-strand break-induced (DSB) cells send signal that induces DSBs in neighbour cells, resulting in the interaction among cells sharing the same medium. Since p53 network gives oscillatory response to DSBs, such interaction among cells could be modelled as an excitatory coupling of p53 network oscillators. This study proposes a plausible coupling model of three-mode two-dimensional oscillators, which models the p53-mediated cell fate selection in globally coupled DSB-induced cells. The coupled model consists of ATM and Wip1 proteins as variables. The coupling mechanism is realised through ATM variable via a mean-field modelling the bystander signal in the intercellular medium. Investigation of the model reveals that the coupling generates more sensitive DNA damage response by affecting cell fate selection. Additionally, the authors search for the cause-effect relationship between coupled p53 network oscillators and bystander effect (BE) endpoints. For this, they search for the possible values of uncertain parameters that may replicate BE experiments' results. At certain parametric regions, there is a correlation between the outcomes of cell fate and endpoints of BE, suggesting that the intercellular coupling of p53 network may manifest itself as the form of observed BEs.

Cell Cycle and DNA Repair Regulation in the Damage Response: Protein Phosphatases Take Over the Reins

Cells are constantly suffering genotoxic stresses that affect the integrity of our genetic material. Genotoxic insults must be repaired to avoid the loss or inappropriate transmission of the genetic information, a situation that could lead to the appearance of developmental abnormalities and tumorigenesis. To combat this threat, eukaryotic cells have evolved a set of sophisticated molecular mechanisms that are collectively known as the DNA damage response (DDR). This surveillance system controls several aspects of the cellular response, including the detection of lesions, a temporary cell cycle arrest, and the repair of the broken DNA. While the regulation of the DDR by numerous kinases has been well documented over the last decade, the complex roles of protein dephosphorylation have only recently begun to be investigated. Here, we review recent progress in the characterization of DDR-related protein phosphatases during the response to a DNA lesion, focusing mainly on their ability to modulate the DNA damage checkpoint and the repair of the damaged DNA. We also discuss their protein composition and structure, target specificity, and biochemical regulation along the different stages encompassed in the DDR. The compilation of this information will allow us to better comprehend the physiological significance of protein dephosphorylation in the maintenance of genome integrity and cell viability in response to genotoxic stress.

Inferring Leading Interactions in the p53/Mdm2/Mdmx Circuit through Live-Cell Imaging and Modeling

The tumor-suppressive transcription factor p53 is a master regulator of stress responses. In non-stressed conditions, p53 is maintained at low levels by the ubiquitin ligase Mdm2 and its binding partner Mdmx. Mdmx depletion leads to a biphasic p53 response, with an initial post-mitotic pulse followed by oscillations. The mechanism underlying this dynamical behavior is unknown. Two different roles for Mdmx have been proposed: enhancing p53 ubiquitination by Mdm2 and inhibiting p53 activity. Here, we developed a mathematical model of the p53/Mdm2/Mdmx network to investigate which Mdmx functions quantitatively affect specific features of p53 dynamics under various conditions. We found that enhancement of Mdm2 activity was the most critical role of Mdmx under unstressed conditions. The model also accurately predicted p53 dynamics in Mdmx-depleted cells following DNA damage. This work outlines a strategy for rapidly testing possible network interactions to reveal those most impactful in regulating the dynamics of key proteins.

Stochastic transcription in the p53-mediated response to DNA damage is modulated by burst frequency

Discontinuous transcription has been described for different mammalian cell lines and numerous promoters. However, our knowledge of how the activity of individual promoters is adjusted by dynamic signaling inputs from transcription factors is limited. To address this question, we characterized the activity of selected target genes that are regulated by pulsatile accumulation of the tumor suppressor p53 in response to ionizing radiation. We performed time-resolved measurements of gene expression at the single-cell level by smFISH and used the resulting data to inform a mathematical model of promoter activity. We found that p53 target promoters are regulated by frequency modulation of stochastic bursting and can be grouped along three archetypes of gene expression. The occurrence of these archetypes cannot solely be explained by nuclear p53 abundance or promoter binding of total p53. Instead, we provide evidence that the time-varying acetylation state of p53's C-terminal lysine residues is critical for gene-specific regulation of stochastic bursting.

Bifurcation analysis of bistable and oscillatory dynamics in biological networks using the root-locus method

Most of the biological systems including gene regulatory networks can be described well by ordinary differential equation models with rational non-linearities. These models are derived either based on the reaction kinetics or by curve fitting to experimental data. This study demonstrates the applicability of the root-locus-based bifurcation analysis method for studying the complex dynamics of such models. The effectiveness of the bifurcation analysis in determining the exact parameter regions in each of which the system shows a certain dynamical behaviour, such as bistability, oscillation, and asymptotically equilibrium dynamics is shown by considering two mostly studied gene regulatory networks, namely Gardner's genetic toggle switch and p53 gene network possessing two-phase (mono-stable/oscillation) dynamics.

Oscillation induced by Hopf bifurcation in the p53-Mdm2 feedback module

This study develops an integrated model of the p53-Mdm2 interaction composed of five basic components and time delay in the DNA damage response based on the existing research work. Some critical factors, including time delay, system parameters, and their interactions in the p53-Mdm2 system are investigated to examine their effects on the oscillatory behaviour induced by Hopf bifurcation. It is shown that the positive feedback formed between p53 and the activity of Mdm2 in the cytoplasm can cause a slight decrease in the amplitude of the p53 oscillation. The length of the time delay plays an important role in determining the amplitude and period of the oscillation and can significantly extend the parameter range for the system to demonstrate oscillatory behaviour. The numerical simulation results are found to be in good agreement with the published experimental observation. It is expected that the results of this research would be helpful to better understand the biological functions of p53 pathway and provide some clues in the treatment of cancer.

Life or Death after a Break: What Determines the Choice?

DNA double-strand breaks (DSBs) pose a constant threat to genomic integrity. Such DSBs need to be repaired to preserve homeostasis at both the cellular and organismal levels. Hence, the DNA damage response (DDR) has evolved to repair these lesions and limit their toxicity. The initiation of DNA repair depends on the activation of the DDR, and we know that the strength of DDR signaling may differentially affect cellular viability. However, we do not fully understand what determines the cytotoxicity of a DSB. Recent work has identified genomic location, (in)correct DNA repair pathway usage, and cell-cycle position as contributors to DSB-induced cytotoxicity. In this review, we discuss how these determinants affect cytotoxicity, highlight recent discoveries, and identify open questions that could help to improve our understanding about cell fate decisions after a DNA DSB.

Cell Membrane Composition Drives Selectivity and Toxicity of Designed Cyclic Helix-Loop-Helix Peptides with Cell Penetrating and Tumor Suppressor Properties

The tumor suppressor protein p53 is inactive in a large number of cancers, including some forms of sarcoma, breast cancer, and leukemia, due to overexpression of its intrinsic inhibitors MDM2 and MDMX. Reactivation of p53 tumor suppressor activity, via disruption of interactions between MDM2/X and p53 in the cytosol, is a promising strategy to treat cancer. Peptides able to bind MDM2 and/or MDMX were shown to prevent MDM2/X:p53 interactions, but most possess low cell penetrability, low stability, and/or high toxicity to healthy cells. Recently, the designed peptide cHLH-p53-R was reported to possess high affinity for MDM2, resistance toward proteases, cell-penetrating properties, and toxicity toward cancer cells. This peptide uses a stable cyclic helix-loop-helix (cHLH) scaffold, which includes two helices connected with a Gly loop and cyclized to improve stability. In the current study, we were interested in examining the cell selectivity of cHLH-p53-R, its cellular internalization, and ability to reactivate the p53 pathway. We designed analogues of cHLH-p53-R and employed biochemical and biophysical methodologies using in vitro model membranes and cell-based assays to compare their structure, activity, and mode-of-action. Our studies show that cHLH is an excellent scaffold to stabilize and constrain p53-mimetic peptides with helical conformation, and reveal that anticancer properties of cHLH-p53-R are mediated by its ability to selectively target, cross, and disrupt cancer cell membranes, and not by activation of the p53 pathway. These findings highlight the importance of examining the mode-of-action of designed peptides to fully exploit their potential to develop targeted therapies.

TRIM66 reads unmodified H3R2K4 and H3K56ac to respond to DNA damage in embryonic stem cells

Recognition of specific chromatin modifications by distinct structural domains within “reader” proteins plays a critical role in the maintenance of genomic stability. However, the specific mechanisms involved in this process remain unclear. Here we report that the PHD-Bromo tandem domain of tripartite motif-containing 66 (TRIM66) recognizes the unmodified H3R2-H3K4 and acetylated H3K56. The aberrant deletion of Trim66 results in severe DNA damage and genomic instability in embryonic stem cells (ESCs). Moreover, we find that the recognition of histone modification by TRIM66 is critical for DNA damage repair (DDR) in ESCs. TRIM66 recruits Sirt6 to deacetylate H3K56ac, negatively regulating the level of H3K56ac and facilitating the initiation of DDR. Importantly, Trim66-deficient blastocysts also exhibit higher levels of H3K56ac and DNA damage. Collectively, the present findings indicate the vital role of TRIM66 in DDR in ESCs, establishing the relationship between histone readers and maintenance of genomic stability.

TRIM66 protein has an N-terminal tripartite motif and a C-terminal PHD Bromodomain. Here the authors show the specific histone modification recognition of TRIM66-PHD-Bromodomain through crystallography and biochemistry assay, and further reveal that TRIM66 recognition of certain histone modification is important for DNA damage repair in ESCs.

p53 pulse modulation differentially regulates target gene promoters to regulate cell fate decisions

The p53 tumor suppressor regulates distinct responses to cellular stresses. Although different stresses generate different p53 dynamics, the mechanisms by which cells decode p53 dynamics to differentially regulate target genes are not well understood. Here, we determined in individual cells how canonical p53 target gene promoters vary in responsiveness to features of p53 dynamics. Employing a chemical perturbation approach, we independently modulated p53 pulse amplitude, duration, or frequency, and we then monitored p53 levels and target promoter activation in individual cells. We identified distinct signal processing features-thresholding in response to amplitude modulation, a refractory period in response to duration modulation, and dynamic filtering in response to frequency modulation. We then showed that the signal processing features not only affect p53 target promoter activation, they also affect p53 regulation and downstream cellular functions. Our study shows how different promoters can differentially decode features of p53 dynamics to generate distinct responses, providing insight into how perturbing p53 dynamics can be used to generate distinct cell fates.

The DNA-damage response and nuclear events as regulators of nonapoptotic forms of cell death

The maintenance of genome stability is essential for the cell as the integrity of genomic information guaranties reproduction of a whole organism. DNA damage occurring in response to different natural and nonnatural stimuli (errors in DNA replication, UV radiation, chemical agents, etc.) is normally detected by special cellular machinery that induces DNA repair. However, further accumulation of genetic lesions drives the activation of cell death to eliminate cells with defective genome. This particular feature is used for targeting fast-proliferating tumor cells during chemo-, radio-, and immunotherapy. Among different cell death modalities induced by DNA damage, apoptosis is the best studied. Nevertheless, nonapoptotic cell death and adaptive stress responses are also activated following genotoxic stress and play a crucial role in the outcome of anticancer therapy. Here, we provide an overview of nonapoptotic cell death pathways induced by DNA damage and discuss their interplay with cellular senescence, mitotic catastrophe, and autophagy.

Abrogation of FBW7α-dependent p53 degradation enhances p53's function as a tumor suppressor

The gene encoding the tumor suppressor p53 is mutated in most cancers. p53 expression is known to be tightly controlled by several E3 ligases. Here, we show that F-box and WD repeat domain-containing 7α (FBW7α), the substrate-recognition component of the SCF^FBW7 multiprotein E3 ligase complex, targets both WT and tumor-derived mutants of p53 for proteasomal degradation in multiple human cancer cell lines (HCT116 and U2OS). We found that lack of FBW7α stabilizes p53 levels, thereby increasing its half-life. p53 ubiquitylation and subsequent degradation require the F-box and the C-terminal WD40 repeats in FBW7α. The polyubiquitylation of p53 occurred via Lys-48 linkage and involved phosphorylation on p53 at Ser-33 and Ser-37 by glycogen synthase kinase 3β (GSK3β) and DNA-dependent protein kinase (DNA-PK), respectively. These phosphorylation events created a phosphodegron that enhanced p53 binding to FBW7α, allowing for the attachment of polyubiquitin moieties at Lys-132 in p53. FBW7α-dependent p53 polyubiquitylation apparently occurred during and immediately after DNA double-strand breaks induced by either doxorubicin or ionizing radiation. Accordingly, in cells lacking FBW7α, p53 induction was enhanced after DNA damage. Phosphodegron-mediated polyubiquitylation of p53 on Lys-132 had functional consequences, with cells in which FBW7α-mediated p53 degradation was abrogated exhibiting enhancement of their tumorigenic potential. We conclude that p53, which previously has been reported to transactivate FBW7, is also targeted by the same E3 ligase for degradation, suggesting the presence of a regulatory feedback loop that controls p53 levels and functions during DNA damage.

Helping the Released Guardian: Drug Combinations for Supporting the Anticancer Activity of HDM2 (MDM2) Antagonists

The protein p53, known as the "Guardian of the Genome", plays an important role in maintaining DNA integrity, providing protection against cancer-promoting mutations. Dysfunction of p53 is observed in almost every cancer, with 50% of cases bearing loss-of-function mutations/deletions in the TP53 gene. In the remaining 50% of cases the overexpression of HDM2 (mouse double minute 2, human homolog) protein, which is a natural inhibitor of p53, is the most common way of keeping p53 inactive. Disruption of HDM2-p53 interaction with the use of HDM2 antagonists leads to the release of p53 and expression of its target genes, engaged in the induction of cell cycle arrest, DNA repair, senescence, and apoptosis. The induction of apoptosis, however, is restricted to only a handful of p53wt cells, and, generally, cancer cells treated with HDM2 antagonists are not efficiently eliminated. For this reason, HDM2 antagonists were tested in combinations with multiple other therapeutics in a search for synergy that would enhance the cancer eradication. This manuscript aims at reviewing the recent progress in developing strategies of combined cancer treatment with the use of HDM2 antagonists.

Contribution of time delays to p53 oscillation in DNA damage response

Although the oscillatory dynamics of the p53 network have been extensively studied, the understanding of the mechanism of delay-induced oscillations is still limited. In this paper, a comprehensive mathematical model of p53 network is studied, which contains two delayed negative feedback loops. By studying the model with and without explicit delays, the results indicate that the time delay of Mdm2 protein synthesis can well control the pulse shape but cannot induce p53 oscillation alone, while the time delay required for Wip1 protein synthesis induces a Hopf bifurcation to drive p53 oscillation. In addition, the synergy of the two delays will cause the p53 network to oscillate in advance, indicating that p53 begins the repair process earlier in the damaged cell. Furthermore, the stability and bifurcation of the model are addressed, which may highlight the role of time delay in p53 oscillations.

Exploring the Influence of Parameters on the p53 Response When Single-Stranded Breaks and Double-Stranded Breaks Coexist

The p53 response to DNA damage is closely related to cell fate decisions. P53 preferentially responds to single-stranded breaks (SSBs) exhibiting a graded response when single-stranded breaks (SSBs) and double-stranded breaks (DSBs) coexist. However, how p53 natural preferential response is affected by kinetic parameters remains to be elucidated. Here, based on the hybrid model I, we computationally searched all the parameters and parameter combinations in the parameter space to identify those that could alter the natural preferential response of p53 when SSBs and DSBs coexist. Firstly, when a single parameter is changed, the parameters that can alter graded response to produce p53 pulse response are production rate of ATM- and Rad3-related kinase(ATR) (beta2), ATR degradation rate (alf2) and ATR-dependent p53 production rate (beta31). Secondly, when double parameters are changed, the combinations of beta2/alf2/beta31 and any other parameters are capable of altering the p53 natural preferential response, and the combination of ataxia-telangiectasia mutated kinase (ATM)-dependent p53 production rate (beta3) and Wip1-dependent p53 degradation rate (alf35) is also capable of altering the p53 natural preferential response. Thirdly, we analyzed the sensitivity of both pulse amplitude and apoptosis to kinetic parameters. We find that pulse amplitude is most sensitive to ATM-dependent p53 production rate (beta3), and apoptosis is more sensitive to damage-dependent ATM production rate (beta1), wip1-dependent ATM degradation rate (alf15), wip1 production rate (beta5) and wip1 degradation rate (alf5). What is more, the smaller the value of alf15/beta5 or the larger the value of beta1/alf5, the more susceptible the cells are to apoptosis. These results provide clues to design more effective and less toxic targeted treatments for cancer.

Oscillatory Dynamics of p53 Genetic Network Induced by Feedback Loops and Time Delays

DNA damage caused by γ -irradiation initiates oscillatory expression of the p53 genetic network. Although many studies revealed the effects of the p53-Mdm2 circuit on p53 dynamics, a few studies explored the contribution of upstream kinases to p53 oscillation. In this paper, an integrated mathematical model of the p53 network in response to γ -irradiation is studied, which consists of five basic components, two ubiquitous time delays, and two negative feedback loops. It is found that recurrent p53 pulses are externally initiated by ataxia telangiectasia mutated (ATM), and the negative feedback loop formed between ATM and p53, via Wip1, plays a dominant role in generating p53 oscillation. In addition, p53 oscillation requires not only an appropriate Mdm2 negative strength but also a threshold level of Wip1 negative strength. Furthermore, the time delays required for transcription and translation of Mdm2 and Wip1 proteins are essential for p53 oscillation. In particular, the critical value of time delay for inducing oscillation and the properties of delay-driven Hopf bifurcation are theoretically analyzed. As expected, the results are clearly in consistence with biological experiments and observations.

Integrated Regulation of HuR by Translation Repression and Protein Degradation Determines Pulsatile Expression of p53 Under DNA Damage

Expression of tumor suppressor p53 is regulated at multiple levels, disruption of which often leads to cancer. We have adopted an approach combining computational systems modeling with experimental validation to elucidate the translation regulatory network that controls p53 expression post DNA damage. The RNA-binding protein HuR activates p53 mRNA translation in response to UVC-induced DNA damage in breast carcinoma cells. p53 and HuR levels show pulsatile change post UV irradiation. The computed model fitted with the observed pulse of p53 and HuR only when hypothetical regulators of synthesis and degradation of HuR were incorporated. miR-125b, a UV-responsive microRNA, was found to represses the translation of HuR mRNA. Furthermore, UV irradiation triggered proteasomal degradation of HuR mediated by an E3-ubiquitin ligase tripartite motif-containing 21 (TRIM21). The integrated action of miR-125b and TRIM21 constitutes an intricate control system that regulates pulsatile expression of HuR and p53 and determines cell viability in response to DNA damage.

Apoptosis: A mammalian cell bioprocessing perspective

Apoptosis is a form of programmed and controlled cell death that accounts for the majority of cellular death in bioprocesses. Cell death affects culture longevity and product quality; it is instigated by several stresses experienced by the cells within a bioreactor. Understanding the factors that cause apoptosis as well as developing strategies that can protect cells is crucial for robust bioprocess development. This review aims to a) address apoptosis from a bioprocess perspective; b) describe the significant apoptotic mechanisms linking them to the most relevant stresses encountered in bioreactors; c) discuss the design of operating conditions in order to avoid cell death; d) focus on industrially relevant cell lines; and e) present anti-apoptosis strategies including cell engineering and model-based optimization of bioprocesses. In addition, the importance of apoptosis in quality-by-design bioprocess development from clone screening to production scale are highlighted.

Tuning of mRNA stability through altering 3'-UTR sequences generates distinct output expression in a synthetic circuit driven by p53 oscillations

Synthetic biological circuits that can generate outputs with distinct expression dynamics are useful for a variety of biomedical and industrial applications. We present a method to control output dynamics by altering output mRNA decay rates. Using oscillatory expression of the transcription factor p53 as the circuit regulator, we use two approaches for controlling target gene transcript degradation rates based on the output gene's 3'-untranslated region (3'-UTR): introduction of copies of destabilizing AU-rich elements into the 3'-UTR or swapping in naturally occurring 3'-UTRs conferring different transcript stabilities. As a proof of principle, we apply both methods to control the expression dynamics of a fluorescent protein and visualize the circuit output dynamics in single living cells. We then use the naturally occurring 3'-UTR approach to restore apoptosis in a tunable manner in a cancer cell line deficient for caspase-3 expression. Our method can be readily adapted to regulate multiple outputs each with different expression dynamics under the control of a single naturally occurring or synthetically constructed biological oscillator.

Roles of cellular heterogeneity, intrinsic and extrinsic noise in variability of p53 oscillation

The p53 protein is a key mediator of the cellular response to various stress signals. In response to DNA damage, the concentration of p53 can temporally oscillate with fluctuations in both the amplitude and period. The underlying mechanism for p53 variability is not fully understood. Here, we construct a core regulatory network of p53 dynamics comprising the ATM-p53-Wip1 and p53-Mdm2 negative feedback loops. We dissect the contributions of cellular heterogeneity, intrinsic noise, and multiple forms of extrinsic noise to p53 variability in terms of the coefficients of variation of four quantities. Cellular heterogeneity greatly determines the fraction of oscillating cells among a population of isogenic cells. Intrinsic noise-fluctuation in biochemical reactions-has little impact on p53 variability given large amounts of molecules, whereas extrinsic colored noise with proper strength and correlation time contributes much to oscillatory variability in individual cells. With the three sources of noise combined, our results reproduce the experimental observations, suggesting that the long correlation time of colored noise is essential to p53 variability. Compared with previous studies, the current work reveals both the individual and integrated effects of distinct noise sources on p53 variability. This study provides a framework for exploring the variability in oscillations in cellular signaling pathways.

Live-Cell Imaging and Analysis with Multiple Genetically Encoded Reporters

Genetically encoded live-cell reporters measure signaling pathway activity at the cellular level with high temporal resolution, often revealing a high degree of cell-to-cell heterogeneity. By using multiple spectrally distinct reporters within the same cell, signal transmission from one node to another within a signaling pathway can be analyzed to quantify factors such as signaling efficiency and delay. With other reporter configurations, correlation between different signaling pathways can be quantified. Such analyses can be used to establish the mechanisms and consequences of cell-to-cell heterogeneity and can inform new models of the functional properties of signaling pathways. In this unit, we describe an approach for designing and executing live-cell multiplexed reporter experiments. We also describe approaches for analyzing the resulting time-course data to quantify correlations and trends between the measured parameters at the single-cell level. © 2018 by John Wiley & Sons, Inc.

Sulforaphane Cannot Protect Human Fibroblasts From Repeated, Short and Sublethal Treatments with Hydrogen Peroxide

A delicate balance of reactive oxygen species (ROS) exists inside the cell: when the mechanisms that control the level of ROS fail, the cell is in an oxidative stress state, a condition that can accelerate aging processes. To contrast the pro-aging effect of ROS, the supplementation of antioxidants has been recently proposed. Sulforaphane (SFN) is an isothiocyanate isolated from Brassica plants that has been shown to modulate many critical factors inside the cells helping to counteract aging processes. In the present work, we exposed human dermal fibroblast to short, sublethal and repeated treatments with hydrogen peroxide for eight days, without or in combination with low concentration of SFN. Hydrogen peroxide treatments did not affect the oxidative status of the cells, without any significant change of the intracellular ROS levels or the number of mitochondria or thiols in total proteins. However, our regime promoted cell cycle progression and cell viability, increased the anti-apoptotic factor survivin and increased DNA damage, measured as number of foci positive for γ -H2AX. On the other hand, the treatment with SFN alone seemed to exert a protective effect, increasing the level of p53, which can block the expansion of possible DNA damaged cells. However, continued exposure to SFN at this concentration could not protect the cells from stress induced by hydrogen peroxide.

Protein stability of p53 targets determines their temporal expression dynamics in response to p53 pulsing

Oscillations in p53 expression are critical for regulating the cellular response to DNA damage. Hanson et al. show that the relationship between p53 pulse frequency and target mRNA and protein decay rates regulates stress response pathway dynamics and function.

In response to DNA damage, the transcription factor p53 accumulates in a series of pulses. While p53 dynamics play a critical role in regulating stress responses, how p53 pulsing affects target protein expression is not well understood. Recently, we showed that p53 pulses generate diversity in target mRNA expression dynamics; however, given that mRNA and protein expression are not necessarily well correlated, it remains to be determined how p53 pulses impact target protein expression. Using computational and experimental approaches, we show that target protein decay rates filter p53 pulses: Distinct target protein expression dynamics are generated depending on the relationship between p53 pulse frequency and target mRNA and protein stability. Furthermore, by mutating the targets MDM2 and PUMA to alter their stabilities, we show that downstream pathways are sensitive to target protein decay rates. This study delineates the mechanisms by which p53 dynamics play a crucial role in orchestrating the timing of events in the DNA damage response network.

Modulation of Phase Shift between Wnt and Notch Signaling Oscillations Controls Mesoderm Segmentation

How signaling dynamics encode information is a central question in biology. During vertebrate development, dynamic Notch signaling oscillations control segmentation of the presomitic mesoderm (PSM). In mouse embryos, this molecular clock comprises signaling oscillations of several pathways, i.e., Notch, Wnt, and FGF signaling. Here, we directly address the role of the relative timing between Wnt and Notch signaling oscillations during PSM patterning. To this end, we developed a new experimental strategy using microfluidics-based entrainment that enables specific control of the rhythm of segmentation clock oscillations. Using this approach, we find that Wnt and Notch signaling are coupled at the level of their oscillation dynamics. Furthermore, we provide functional evidence that the oscillation phase shift between Wnt and Notch signaling is critical for PSM segmentation. Our work hence reveals that dynamic signaling, i.e., the relative timing between oscillatory signals, encodes essential information during multicellular development.

•

Wnt and Notch signaling wave dynamics differ within segmenting mouse mesoderm

•

Entraining oscillations by microfluidics allows external control of the dynamics

•

Oscillatory Wnt and Notch signaling networks are coupled at the level of dynamics

•

Relative timing of Wnt and Notch signaling oscillations is critical for segmentation

Theoretical study on the oscillation mechanism of p53-Mdm2 network

In this paper, a delayed mathematical model was developed based on experimental data to understand how the time delays required for transcription and translation in Mdm2 gene expression affect the kinetic behavior of the p53-Mdm2 network. Taking the time delays as the main research parameters, the stability of the system at the positive equilibrium was studied by using the theoretical method of delay differential equation. We found that such delays can induce oscillations by undergoing a supercritical Hopf bifurcation. Then, we used the normal form theory and the center manifold reduction to study the direction and stability of the bifurcation in detail. Furthermore, we also studied the effects of the length of time delays and the model parameters by numerical simulations. We found that time delays in Mdm2 synthesis are required for p53 oscillations and the length of such delays can determine the amplitude and period of the oscillations. In addition, the model parameters can also change the stability of the system. These results illustrate that the repair process after DNA damage can be regulated by varying time delays and the model parameters.

Tissue-Specific Chk1 Activation Determines Apoptosis by Regulating the Balance of p53 and p21

The DNA damage response (DDR) protects cells against genomic instability. Surprisingly, little is known about the differences in DDR across tissues, which may affect cancer evolutionary trajectories and chemotherapy response. Using mathematical modeling and quantitative experiments, we found that the DDR is regulated differently in human breast and lung primary cells. Equal levels of cisplatin-DNA lesions caused stronger Chk1 activation in lung cells, leading to resistance. In contrast, breast cells were more resistant and showed more Chk2 activation in response to doxorubicin. Further analyses indicate that Chk1 activity played a regulatory role in p53 phosphorylation, whereas Chk2 activity was essential for p53 activation and p21 expression. We propose a novel “friction model,” in which the balance of p53 and p21 levels contributes to the apoptotic response in different tissues. Our results suggest that modulating the balance of p53 and p21 dynamics could optimize the response to chemotherapy.

•

Breast and lung cells show different sensitivities to chemotherapeutic drugs

•

Lung cells activate Chk1 more strongly than breast cells with chemotherapeutic drugs

•

Active Chk1 plays a regulatory role in p53 activation and apoptosis responses

•

The balance of p53 and p21 dynamics drives the apoptosis response to DNA damage

Signal Distortion: How Intracellular Pathogens Alter Host Cell Fate by Modulating NF-κB Dynamics

By uncovering complex dynamics in the expression or localization of transcriptional regulators in single cells that were otherwise hidden at the population level, live cell imaging has transformed our understanding of how cells sense and orchestrate appropriate responses to changes in their internal state or extracellular environment. This has proved particularly true for the nuclear factor-kappaB (NF-κB) family of transcription factors, key regulators of the inflammatory response and innate immune function, which are capable of encoding information about the mode and intensity of stimuli in the dynamics of NF-κB nuclear accumulation and loss. While live cell imaging continues to serve as a useful tool in ongoing efforts to characterize the feedbacks that shape these dynamics and to connect dynamics to downstream gene expression, it is also proving invaluable for recent studies that seek to determine how intracellular pathogens subvert NF-κB signaling to survive and replicate within host cells by providing quantitative information about the pathogen and changes in NF-κB activity during different stages of an infection. Here, we provide a brief overview of NF-κB signaling in innate immune cells and review recent literature that uses live imaging to investigate the mechanisms by which bacterial and yeast pathogens modulate NF-κB in a variety of different host cell types to evade destruction or maintain the viability of an intracellular growth niche.

Cell type-dependent bimodal p53 activation engenders a dynamic mechanism of chemoresistance

Studies of drug resistance mostly characterize genetic mutation, and we know much less about phenotypic mechanisms of drug resistance, especially at a quantitative level. p53 is an important mediator of cellular response to chemotherapy, but even p53 wild-type cells vary in drug sensitivity for unclear reasons. Here, we elucidated a new resistance mechanism to a DNA-damaging chemotherapeutic through bimodal modulation of p53 activation dynamics. By combining single-cell imaging with computational modeling, we characterized a four-component regulatory module, which generates bimodal p53 dynamics through coupled feed-forward and feedback, and found that the inhibitory strength between ATM and Mdm2 determined the differential modular output between drug-sensitive and drug-resistant cancer cell lines. We further showed that the combinatorial inhibition of Mdm2 and Wip1 was an effective strategy to alter p53 dynamics in resistant cancer cells and sensitize their apoptotic response. Our results point to p53 pulsing as a potentially druggable mechanism that mediates chemoresistance.

Mathematical model identifies effective P53 accumulation with target gene binding affinity in DNA damage response for cell fate decision

Functional p53 signaling is essential for appropriate responses to diverse stimuli. P53 dynamics employs the information from the stimulus leading to selective gene expression and cell fate decision. However, the decoding mechanism of p53 dynamics under DNA damage challenge remains poorly understood. Here we mathematically modeled the recently dual-phase p53 dynamics under doxorubicin treatment. We found that p53 could perform sequential pulses followed by a high-amplitude terminal pulse at relatively low doxorubicin treatment, whereas p53 became steadily accumulated when damage level was high. The effective p53 integral above a threshold but not the absolute accumulation of p53 precisely discriminated survival and death. Silencing negative regulators in p53 network might promote the occurrence of terminal pulse. Furthermore, lower binding affinity and degradation rate of p53 target genes could favorably discriminate high and low dose doxorubicin treatment. Grouping by temporal profiles suggested that the p53 dynamics rather than the doxorubicin doses could better discriminate cellular outcomes and confer less variation for effective p53 integral reemphasizing the importance of p53 level regulation. Our model has established a theoretical framework that p53 dynamics can work cooperatively with its binding affinity to target genes leading to cell fate choice, providing new clues of optimized clinical design by manipulating p53 dynamics.

Regulation of Tip60-dependent p53 acetylation in cell fate decision

The acetylation of p53 plays an essential role in regulating its transcriptional activity. Nevertheless, how p53 acetylation is modulated in cell fate decision is less understood. We developed a network model to reveal how Tip60-dependent p53 acetylation is regulated to modulate cellular outcome. We proposed that p53 is progressively activated and exhibits distinct dynamics depending on the severity of DNA damage. For mild damage, p53 is primarily activated to trigger cell cycle arrest by transactivating p21, with its concentration showing pulses. For severe damage, p53 is acetylated at Lysine 120 (K120) by Tip60 to trigger apoptosis by inducing PUMA, and its concentration increases to high levels. Several p53-centered feedback loops coordinate to regulate its acetylation status, ensuring a robust decision on cell fate.

Bi-functional biochemical networks

Understanding the relationship between the topology of a network and its function remains an important question in biological physics. However, this is not a one-to-one mapping. Often the behavior of a signaling system varies with the input signal it receives. For example, some biological systems show adaptation when they receive a low input signal while they show oscillation with a high input signal. We therefore set out to find all possible two-node and three-node networks that can perform both adaptation and oscillation with transcriptional regulation and enzymatic reactions. For two-node networks, we identified all bi-functional topologies by analyzing the Jacobean matrix. For three-node networks, they were identified by enumeration. We further investigated how the system can be transformed between these two functions. We found that the switching of functions can be achieved through changing anyone of the several key parameters, including the input signal level.

Defenses against Pro-oxidant Forces - Maintenance of Cellular and Genomic Integrity and Longevity

There has been enormous recent progress in understanding how human cells respond to oxidative stress, such as that caused by exposure to ionizing radiation. We have witnessed a significant deciphering of the events that underlie how antioxidant responses counter pro-oxidant damage to key biological targets in all cellular compartments, including the genome and mitochondria. These cytoprotective responses include: 1. The basal cellular repertoire of antioxidant capabilities and its supporting cast of facilitator enzymes; and 2. The inducible phase of the antioxidant response, notably that mediated by the Nrf2 transcription factor. There has also been frenetic progress in defining how reactive electrophilic species swamp existing protective mechanisms to augment DNA damage, events that are embodied in the cellular "DNA-damage response", including cell cycle checkpoint activation and DNA repair, which occur on a time scale of hours to days, as well as the implementation of cellular responses such as apoptosis, autophagy, senescence and reprograming that extend the time period of damage sensing and response into weeks, months and years. It has become apparent that, in addition to the initial oxidative insult, cells typically undergo further waves of secondary reactive oxygen/nitrogen species generation, DNA damage and signaling and that these may reemerge long after the initial events have subsided, probably being driven, at least in part, by persisting DNA damage. These reactive oxygen/nitrogen species are an integral part of the pathological consequences of radiation exposure and may persist across multiple cell divisions. Because of the pervasive nature of oxidative stress, a cell will manifest different responses in different subcellular compartments and to different levels of stress injury. Aspects of these compartmentalized responses can involve the same proteins (such as ATM, p53 and p21) but in different functional guises, e.g., in cytoplasmic versus nuclear responses or in early- versus late-phase events. Many of these responses involve gene activation and new protein synthesis as well as a plethora of post-translational modifications of both basal and induced response proteins. It is these responses that we focus on in this review.

Modeling Procedures

Being concerned by the understanding of the mechanism underlying chronic degenerative diseases , we presented in the previous chapter the medical systems biology conceptual framework that we present for that purpose in this volume. More specifically, we argued there the clear advantages offered by a state-space perspective when applied to the systems-level description of the biomolecular machinery that regulates complex degenerative diseases. We also discussed the importance of the dynamical interplay between the risk factors and the network of interdependencies that characterizes the biochemical, cellular, and tissue-level biomolecular reactions that underlie the physiological processes in health and disease. As we pointed out in the previous chapter, the understanding of this interplay (articulated around cellular phenotypic plasticity properties, regulated by specific kinds of gene regulatory networks) is necessary if prevention is chosen as the human-health improvement strategy (potentially involving the modulation of the patient's lifestyle). In this chapter we provide the medical systems biology mathematical and computational modeling tools required for this task.

Two-dimensional polynomial type canonical relaxation oscillator model for p53 dynamics

p53 network, which is responsible for DNA damage response of cells, exhibits three distinct qualitative behaviours; low state, oscillation and high state, which are associated with normal cell cycle progression, cell cycle arrest and apoptosis, respectively. The experimental studies demonstrate that these dynamics of p53 are due to the ATM and Wip1 interaction. This paper proposes a simple two-dimensional canonical relaxation oscillator model based on the identified topological structure of ATM and Wip1 interaction underlying these qualitative behaviours of p53 network. The model includes only polynomial terms that have the interpretability of known ATM and Wip1 interaction. The introduced model is useful for understanding relaxation oscillations in gene regulatory networks. Through mathematical analysis, we investigate the roles of ATM and Wip1 in forming of these three essential behaviours, and show that ATM and Wip1 constitute the core mechanism of p53 dynamics. In agreement with biological findings, we show that Wip1 degradation term is a highly sensitive parameter, possibly related to mutations. By perturbing the corresponding parameters, our model characterizes some mutations such as ATM deficiency and Wip1 overexpression. Finally, we provide intervention strategies considering our observation that Wip1 seems to be an important target to conduct therapies for these mutations.