The p53 protein and its molecular network: modelling a missing link between DNA damage and cell fate

Circumventing glioblastoma resistance to temozolomide through optimal drug combinations designed by systems pharmacology and machine learning

BACKGROUND AND PURPOSE

Glioblastoma (GBM), the most frequent and aggressive brain tumour in adults, is associated with a dismal prognostic despite intensive treatment involving surgery, radiotherapy and temozolomide (TMZ)-based chemotherapy. The initial or acquired resistance of GBM to TMZ appeals for precision medicine approaches to the design of novel efficient combination pharmacotherapies. Such investigation needs to account for the overexpression of the O6-methylguanine-DNA methyl-transferase (MGMT) repair enzyme which is responsible for TMZ resistance in patients.

EXPERIMENTAL APPROACH

A comprehensive approach combining quantitative systems pharmacology (QSP) models and machine learning (ML) was undertaken to design TMZ-based drug combinations circumventing the initial resistance to the alkylating agent.

KEY RESULTS

A QSP model representing TMZ cellular pharmacokinetics-pharmacodynamics and dysregulated pathways in GBM was developed and validated using multi-type time- and dose-resolved datasets, available in control or MGMT-overexpressing cells. In silico drug screening and subsequent experimental validation identified a strategy to re-sensitise TMZ-resistant cells consisting in combining TMZ with inhibitors of the base excision repair and of homologous recombination. Using ML, functional signatures of response to such optimal multi-agent therapy were derived to assist decision-making in patients.

CONCLUSION AND IMPLICATIONS

We successfully demonstrated the relevance of combined QSP and ML to design efficient drug combinations re-sensitising glioblastoma cells initially resistant to TMZ. The developed framework may further serve to identify personalised therapies and administration schedules by extending it to account for additional patient-specific altered pathways and whole-body features.

Cell cycle deregulation in neurodegenerative diseases

Background: Cell cycle is critical for a wide range of cellular processes such as proliferation, differentiation and apoptosis in dividing cells. Neurons are postmitotic cells which have withdrawn from the cell division cycle. Recent data show us that inappropriate activation of cell cycle regulators including cyclins, cyclin dependent kinases (CDKs) and endogenous cyclin dependent kinase inhibitors (CDKIs) may take part in the aetiology of neurodegenerative diseases. However, the mechanisms for cell cycle reentry in neurodegenerative disease remain unclear.Methods: Electronic databases such as Pubmed, Science Direct, Directory of Open Access Journals, PLOS were searched for relevant articles.Conclusion: The present work reviews basic aspects of cell cycle mechanism, as well as the evidence showing the expression of cell cycle proteins in neurodegenerative disease. We provide a brief summary of these findings and hope to highlight the interaction between the cell cycle reentry and neurodegenerative diseases. Moreover, we outline the possible signaling pathways. However more understanding of the mechanism of cell cycle is of great importance. Because these represents an alternative target for therapeutic interventions, leading to novel treatments of neurodegenerative diseases.

Mechanistic approach of the therapeutic potential of mesenchymal stem cells on brain damage in irradiated mice: emphasis on anti-inflammatory and anti-apoptotic effects

BACKGROUND AND OBJECTIVES

Brain damage which has been induced by radiation generally occurs in radiotherapeutics patients. Stem cell transplantation represents a vital applicant for alleviating neurodegenerative disorders. This work aims at exploring the potential of bone marrow-derived mesenchymal stem cells (BM-MSCs) on brain injury induced by γ radiation in mice and the possible underlying mechanisms were elucidated.

MATERIALS AND METHODS

Mice were allocated into three groups; Group I (Control), Group II (Irradiated control) where mice submitted to 5 Gy of whole-body γ radiation, Group III (Irradiated + BM-MSCs) where mice were intravenously injected of BM-MSCs at a dose of 106 cells/mice 24 h following irradiation. Animals were sacrificed 28 d following exposure to γ radiation.

RESULTS

It was observed that BM-MSCs therapy provided a valuable tissue repair as evidenced by a reduction in inflammatory mediators including tumor necrosis factor alpha (TNF-α), interleukin-1β (IL-1β), nuclear factor kappa (NF-κβ), phosphorylated NF-κβ-p65 (P-NF-κβ-p65), interferon-gamma (IFNγ) and monocyte chemoattractant protein-1 (MCP-1) associated with decreased levels of transforming growth factor-β (TGF-β) and vascular endothelial growth factor (VEGF) in brain tissues of irradiated mice. Furthermore, neuronal apoptosis was declined in brain tissues of the BM-MSCs group as remarkable inhibition of caspase-3 and Bax accompanied by elevation of Bcl-2 proteins expression. These results were supported by histopathological investigation.

CONCLUSIONS

In conclusion, BM-MSCs could display a vital rule in alleviating brain injury in radio-therapeutic patients.

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.

Mesenchymal Stem Cell-Conditioned Medium Protects Hippocampal Neurons From Radiation Damage by Suppressing Oxidative Stress and Apoptosis

Objective

To investigate the effects of mesenchymal stem cell-conditioned medium (MSC-CM) on radiation-induced oxidative stress, survival and apoptosis in hippocampal neurons.

Methods

The following groups were defined: Control, radiation treatment (RT), RT+MSC-CM, MSC-CM, RT + N-Acetylcysteine (RT+NAC), and RT + MSC-CM + PI3 K inhibitor (LY294002). A cell Counting Kit-8 (CCK-8) was used to measure cell proliferation. Apoptosis was examined by AnnexinV/PI flow cytometric analyses. Intracellular reactive oxygen species (ROS) were detected by DCFH-DA. Intracellular glutathione (GSH), malondialdehyde (MDA) content, and superoxide dismutase (SOD) activity were detected by colorimetric assays. Protein levels of γ-H2AX, PI3K-AKT, P53, cleaved caspase-3, Bax, and BCl-2 were analyzed by Western blotting.

Results

The proliferation of HT22 cells was significantly inhibited in the RT group, but was significantly preserved in the RT + MSC-CM group (P < 0.01). Apoptosis was significantly higher in the RT group than in the RT+ MSC-CM group (P < 0.01). MSC-CM decreased intracellular ROS and MDA content after irradiation (P < 0.01). GSH level and SOD activity were higher in the RT + MSC-CM group than in the RT group, as was MMP (P < 0.01). MSC-CM decreased expression of γ-H2AX, P53, Bax, and cleaved-caspase-3, but increased Bcl-2 expression (P < 0.01).

Conclusion

MSC-CM attenuated radiation-induced hippocampal neuron cell line damage by alleviating oxidative stress and suppressing apoptosis.

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.

PD-L1 expression correlated with p53 expression in oral squamous cell carcinoma

Background

Programmed cell death ligand 1 (PD-L1) is an immune checkpoint molecule that attenuates the immune response. PD-L1 contributes to failed antitumor immunity; thereby, blockade of PD-L1 with monoclonal antibody enhances the immune response. Recently, it was reported that PD-L1 was regulated by protein 53 (p53). Besides, cytokeratin 17 (CK17) is thought to be a diagnostic marker of oral squamous cell carcinoma (OSCC). Our aim was to evaluate the correlation between the immunohistochemical expression of PD-L1, p53 and CK17 with clinicopathological characteristics and disease-specific survival in patients with OSCC.

Methods

A total of 48 patients with OSCC were included in this study. Immunohistochemical staining was performed to evaluate the correlation among the expressions of PD-L1, p53 and CK17, and furthermore the correlation among various clinicopathological factors, PD-L1, p53 and CK17.

Results

The positive rate of p53, CK17, PD-L1 (tumor cells) and PD-L1 (tumor-infiltrating lymphocytes) was 63.2%, 91.7%, 48.9% and 57.1%. A statistically significant correlation between p53 expression and T stage and TNM stage (p = 0.049, p = 0.03, respectively) was observed. Also, a statistically significant correlation between p53 and PD-L1 (TCs) expression (p = 0.0009) was observed. Five-year disease-specific survival rate was not significantly correlated with gender, TNM stage, p53 expression, PD-L1 expression and CK17 expression.

Conclusion

The expression of p53 and PD-L1 shows significantly positive correlation in oral squamous cell carcinoma in tumor cells. Also, a significant correlation between p53 expression and T stage and TNM stage was observed. No other significant correlation between PD-L1 staining or CK17 and clinical or pathologic characteristics was identified.

Mathematical Modeling of p53 Pathways

Cells have evolved balanced systems that ensure an appropriate response to stress. The systems elicit repair responses in temporary or moderate stress but eliminate irreparable cells via apoptosis in detrimental conditions of prolonged or severe stress. The tumor suppressor p53 is a central player in these stress response systems. When activated under DNA damage stress, p53 regulates hundreds of genes that are involved in DNA repair, cell cycle, and apoptosis. Recently, increasing studies have demonstrated additional regulatory roles of p53 in metabolism and mitochondrial physiology. Due to the inherent complexity of feedback loops between p53 and its target genes, the application of mathematical modeling has emerged as a novel approach to better understand the multifaceted functions and dynamics of p53. In this review, we discuss several mathematical modeling approaches in exploring the p53 pathways.

Simulation-assisted machine learning

MOTIVATION

In a predictive modeling setting, if sufficient details of the system behavior are known, one can build and use a simulation for making predictions. When sufficient system details are not known, one typically turns to machine learning, which builds a black-box model of the system using a large dataset of input sample features and outputs. We consider a setting which is between these two extremes: some details of the system mechanics are known but not enough for creating simulations that can be used to make high quality predictions. In this context we propose using approximate simulations to build a kernel for use in kernelized machine learning methods, such as support vector machines. The results of multiple simulations (under various uncertainty scenarios) are used to compute similarity measures between every pair of samples: sample pairs are given a high similarity score if they behave similarly under a wide range of simulation parameters. These similarity values, rather than the original high dimensional feature data, are used to build the kernel.

RESULTS

We demonstrate and explore the simulation-based kernel (SimKern) concept using four synthetic complex systems-three biologically inspired models and one network flow optimization model. We show that, when the number of training samples is small compared to the number of features, the SimKern approach dominates over no-prior-knowledge methods. This approach should be applicable in all disciplines where predictive models are sought and informative yet approximate simulations are available.

AVAILABILITY AND IMPLEMENTATION

The Python SimKern software, the demonstration models (in MATLAB, R), and the datasets are available at https://github.com/davidcraft/SimKern.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Alteration of gene expression profile in CD3+ T-cells after downregulating MALT1

Background

T cell immunodeficiency is a common feature in patients with different kinds of hematological disease such as T cell non-Hodgkin lymphoma (T-NHL), B cells NHL (B-NHL), NK/T cell NHL (NK/T-CL) and acute myeloid leukemia (AML). In our recent research, we found that significantly lower expression levels in MALT1 and NF-κB were related to suppression of T cell activation. Therefore, this study was conducted to further investigate the role of downregulating MALT1 in the development of immunodeficiency in T cells.

Methods

We induced activation inhibition in CD3⁺ T cells by MALT1 knockdown. Then we characterized the gene expression profile after MALT1 suppression by microarray analysis.

Result

The differentially expressed genes were ZAP-70, p65, MDM2, ATM, NFATC2 which participate in the NF-κB, p53, and NFAT pathways in CD3⁺ T cells after MALT1 downregulation.

Conclusion

MALT1 suppression may contribute to immunodeficiency in T cells via suppression of T cell activation and proliferation pathways. These data may help to explain some of the characteristics of immunodeficiency of T cells.

Vitexin compound 1, a novel extraction from a Chinese herb, suppresses melanoma cell growth through DNA damage by increasing ROS levels

BACKGROUND

Vitex negundo L (Verbenaceae) is an aromatic shrub that is abundant in Asian countries. A series of compounds from Vitex negundo have been used in traditional Chinese medicine for the treatment of various diseases. Cutaneous melanoma is one of the most aggressive malignancies. A significant feature of melanoma is its resistance to traditional chemotherapy and radiotherapy; therefore, there is an urgent need to develop novel treatments for melanoma.

METHODS

We first examined the effects of VB1 (vitexin compound 1) on cell viability by CCK-8 (cell counting kit) and Colony Formation Assay; And then, we analyzed the apoptosis and cell cycle by flow cytometry, verified apoptosis by Immunoblotting. The in vivo effect of VB1 was evaluated in xenograft mouse model. Potential mechanisms of VB1's antitumor effects were explored by RNA sequencing and the key differential expression genes were validated by real-time quantitative PCR. Finally, the intracellular reactive oxygen species (ROS) level was detected by flow cytometry, and the DNA damage was revealed by Immunofluorescence and Immunoblotting.

RESULTS

In this study, we show that VB1, which is a compound purified from the seed of the Chinese herb Vitex negundo, blocks melanoma cells growth in vitro and in vivo, arrests the cell cycle in G2/M phase and induces apoptosis in melanoma cell lines, whereas the effects are not significantly observed in normal cells. To study the details of VB1, we analyzed the alteration of gene expression profiles after treatment with VB1 in melanoma cells. The findings showed that VB1 can affect various pathways, including p53, apoptosis and the cell cycle pathway, in a variety of melanoma cell lines. Furthermore, we confirmed that VB1 restored the P53 pathway protein level, and then we demonstrated that VB1 significantly induced the accumulation of ROS, which resulted in DNA damage in melanoma cell lines. Interestingly, our results showed that VB1 also increased the ROS levels in BRAFi (BRAF inhibitor)-resistant melanoma cells, leading to DNA cytotoxicity, which caused G2/M phase arrest and apoptosis.

CONCLUSIONS

Taken together, our findings indicate that vitexin compound 1 might be a promising therapeutic Chinese medicine for melanoma treatment regardless of BRAFi resistance.

Involvement of Prokineticin 2-expressing Neutrophil Infiltration in 5-Fluorouracil-induced Aggravation of Breast Cancer Metastasis to Lung

Adjuvant chemotherapy is used for human breast cancer patients, even after curative surgery of primary tumor, to prevent tumor recurrence primarily as a form of metastasis. However, anticancer drugs can accelerate metastasis in several mouse metastasis models. Hence, we examined the effects of postsurgical administration with 5-fluorouracil (5-FU), doxorubicin, and cyclophosphamide, on lung metastasis process, which developed after the resection of the primary tumor arising from the orthotopic injection of a mouse triple-negative breast cancer cell line, 4T1. Only 5-FU markedly increased the numbers and sizes of lung metastasis foci, with enhanced tumor cell proliferation and angiogenesis as evidenced by increases in Ki67-positive cell numbers and CD31-positive areas, respectively. 5-FU-mediated augmented lung metastasis was associated with increases in intrapulmonary neutrophil numbers and expression of neutrophilic chemokines, Cxcl1 and Cxcl2 in tumor cells, with few effects on intrapulmonary T-cell or macrophage numbers. 5-FU enhanced Cxcl1 and Cxcl2 expression in 4T1 cells in a NFκB-dependent manner. Moreover, the administration of a neutrophil-depleting antibody or a Cxcr2 antagonist, SB225002, significantly attenuated 5-FU-mediated enhanced lung metastasis with depressed neutrophil infiltration. Furthermore, infiltrating neutrophils and 4T1 cells abundantly expressed prokineticin-2 (Prok2) and its receptor, Prokr1, respectively. Finally, the administration of 5-FU after the resection of the primary tumor failed to augment lung metastasis in the mice receiving Prokr1-deleted 4T1 cells. Collectively, 5-FU can enhance lung metastasis by inducing tumor cells to produce Cxcl1 and Cxcl2, which induced the migration of neutrophils expressing Prok2 with a capacity to enhance 4T1 cell proliferation. Mol Cancer Ther; 17(7); 1515-25. ©2018 AACR.

Bacterial Cytological Profiling Reveals the Mechanism of Action of Anticancer Metal Complexes

Target identification and mechanistic studies of cytotoxic agents are challenging processes that are both time-consuming and costly. Here we describe an approach to mechanism of action studies for potential anticancer compounds by utilizing the simple prokaryotic system, E. coli, and we demonstrate its utility with the characterization of a ruthenium polypyridyl complex [Ru(bpy)₂dmbpy²⁺]. Expression of the photoconvertible fluorescent protein Dendra2 facilitated both high throughput studies and single-cell imaging. This allowed for simultaneous ratiometric analysis of inhibition of protein production and phenotypic investigations. The profile of protein production, filament size and population, and nucleoid morphology revealed important differences between inorganic agents that damage DNA vs more selective inhibitors of transcription and translation. Trace metal analysis demonstrated that DNA is the preferred nucleic acid target of the ruthenium complex, but further studies in human cancer cells revealed altered cell signaling pathways compared to the commonly administrated anticancer agent cisplatin. This study demonstrates E. coli can be used to rapidly distinguish between compounds with disparate mechanisms of action and also for more subtle distinctions within in studies in mammalian cells.

System-based strategies for p53 recovery

The authors have proposed a systems theory-based novel drug design approach for the p53 pathway. The pathway is taken as a dynamic system represented by ordinary differential equations-based mathematical model. Using control engineering practices, the system analysis and subsequent controller design is performed for the re-activation of wild-type p53. p53 revival is discussed for both modes of operation, i.e. the sustained and oscillatory. To define the problem in control system paradigm, modification in the existing mathematical model is performed to incorporate the effect of Nutlin. Attractor point analysis is carried out to select the suitable domain of attraction. A two-loop negative feedback control strategy is devised to drag the system trajectories to the attractor point and to regulate cellular concentration of Nutlin, respectively. An integrated framework is constituted to incorporate the pharmacokinetic effects of Nutlin in the cancerous cells. Bifurcation analysis is also performed on the p53 model to see the conditions for p53 oscillation.

Circulating exosomal microRNAs reveal the mechanism of Fructus Meliae Toosendan-induced liver injury in mice

The toxicological mechanisms of liver injury caused by most traditional Chinese medicine (TCM) remain largely unknown. Due to the unique features, exosomal microRNAs (miRNAs) are currently attracting major interests to provide further insights into toxicological mechanisms. Thus, taking Fructus Meliae Toosendan as an example of hepatoxic TCM, this study aimed to elucidate its hepatotoxicity mechanisms through profiling miRNAs in circulating exosomes of Fructus Meliae Toosendan water extract (FMT)-exposed mice. Biological pathway analysis of the 64 differentially expressed exosomal miRNAs (DEMs) showed that hepatic dysfunction induced by FMT likely related to apoptosis, mitochondrial dysfunction, and cell cycle dysregulation. Integrated analysis of serum exosomal DEMs and hepatic differentially expressed mRNAs further enriched oxidative stress and apoptosis related pathways. In vitro validation studies for omics results suggested that FMT-induced DNA damage was mediated by generating intracellular reactive oxygen species, leading to cell apoptosis through p53-dependent mitochondrial damage and S-phase arrest. Nrf2-mediated antioxidant response was activated to protect liver cells. Moreover, serum exosomal miR-370-3p, the most down-regulated miRNA involving in these pathways, might be the momentous event in aggravating cytotoxic effect of FMT by elevating p21 and Cyclin E. In conclusion, circulating exosomal miRNAs profiling could contribute to deepen the understanding of TCM-induced hepatotoxicity.

The role of estradiol metabolism in urogenital schistosomiasis-induced bladder cancer

Urogenital schistosomiasis is a neglected tropical disease that can lead to bladder cancer. How urogenital schistosomiasis induces carcinogenesis remains unclear, although there is evidence that the human blood fluke Schistosoma haematobium, the infectious agent of urogenital schistosomiasis, releases estradiol-like metabolites. These kind of compounds have been implicated in other cancers. Aiming for enhanced understanding of the pathogenesis of the urogenital schistosomiasis-induced bladder cancer, here we review, interpret, and discuss findings of estradiol-like metabolites detected in both the parasite and in the human urine during urogenital schistosomiasis. Moreover, we predict pathways and enzymes that are involved in the production of these metabolites emphasizing their potential effects on the dysregulation of the tumor suppressor gene p53 expression during urogenital schistosomiasis. Enhanced understanding of these potential carcinogens may not only shed light on urogenital schistosomiasis-induced neoplasia of the bladder, but would also facilitate development of interventions and biomarkers for this and other infection-associated cancers at large.

Wip1: A candidate phosphatase for cancer diagnosis and treatment

The critical regulatory mechanisms in numerous cellular pathways including cell survival and DNA damage response mostly depend on phosphorylation and dephosphorylation of proteins. The serine/threonine phosphatase wild-type p53-induced phosphatase 1 (Wip1) is a growth-promoting phosphatase and its numerous downstream targets are important tumor suppressors. Here, we review the Wip1 activity and its relevance to cancer as an oncoprotein. Consecutive investigations about Wip1 and its relation to cancer is critical, as these studies ultimately contribute to the etiology of cancer. A number of innovative studies have recently investigated the importance of Wip1 as a new candidate for cancer diagnosis and prognosis. Accordingly, we discuss the present challenges of using Wip1 as a target for cancer treatment.

Spatio-temporal models of synthetic genetic oscillators

Signal transduction pathways play a major role in many important aspects of cellular function e.g. cell division, apoptosis. One important class of signal transduction pathways is gene regulatory networks (GRNs). In many GRNs, proteins bind to gene sites in the nucleus thereby altering the transcription rate. Such proteins are known as transcription factors. If the binding reduces the transcription rate there is a negative feedback leading to oscillatory behaviour in mRNA and protein levels, both spatially (e.g. by observing fluorescently labelled molecules in single cells) and temporally (e.g. by observing protein/mRNA levels over time). Recent computational modelling has demonstrated that spatial movement of the molecules is a vital component of GRNs and may cause the oscillations. These numerical findings have subsequently been proved rigorously i.e. the diffusion coefficient of the protein/mRNA acts as a bifurcation parameter and gives rise to a Hopf bifurcation. In this paper we first present a model of the canonical GRN (the Hes1 protein) and show the effect of varying the spatial location of gene and protein production sites on the oscillations. We then extend the approach to examine spatio-temporal models of synthetic gene regulatory networks e.g. n-gene repressilators and activator-repressor systems.

Unraveling cellular pathways contributing to drug-induced liver injury by dynamical modeling

INTRODUCTION

Drug-induced liver injury (DILI) is a significant threat to human health and a major problem in drug development. It is hard to predict due to its idiosyncratic nature and which does not show up in animal trials. Hepatic adaptive stress response pathway activation is generally observed in drug-induced liver injury. Dynamical pathway modeling has the potential to foresee adverse effects of drugs before they go in trial. Ordinary differential equation modeling can offer mechanistic insight, and allows us to study the dynamical behavior of stress pathways involved in DILI. Areas covered: This review provides an overview on the progress of the dynamical modeling of stress and death pathways pertinent to DILI, i.e. pathways relevant for oxidative stress, inflammatory stress, DNA damage, unfolded proteins, heat shock and apoptosis. We also discuss the required steps for applying such modeling to the liver. Expert opinion: Despite the strong progress made since the turn of the century, models of stress pathways have only rarely been specifically applied to describe pathway dynamics for DILI. We argue that with minor changes, in some cases only to parameter values, many of these models can be repurposed for application in DILI research. Combining both dynamical models with in vitro testing might offer novel screening methods for the harmful side-effects of drugs.

A novel mathematical model of ATM/p53/NF- κB pathways points to the importance of the DDR switch-off mechanisms

Background

Ataxia telangiectasia mutated (ATM) is a detector of double-strand breaks (DSBs) and a crucial component of the DNA damage response (DDR) along with p53 and NF- κB transcription factors and Wip1 phosphatase. Despite the recent advances in studying the DDR, the mechanisms of cell fate determination after DNA damage induction is still poorly understood.

Results

To investigate the importance of various DDR elements with particular emphasis on Wip1, we developed a novel mathematical model of ATM/p53/NF- κB pathways. Our results from in silico and in vitro experiments performed on U2-OS cells with Wip1 silenced to 25 % (Wip1-RNAi) revealed a strong dependence of cellular response to DNA damages on this phosphatase. Notably, Wip1-RNAi cells exhibited lower resistance to ionizing radiation (IR) resulting in smaller clonogenicity and higher apoptotic fraction.

Conclusions

In this article, we demonstrated that Wip1 plays a role as a gatekeeper of apoptosis and influences the pro-survival behaviour of cells – the level of Wip1 increases to block the apoptotic decision when DNA repair is successful. Moreover, we were able to verify the dynamics of proteins and transcripts, apoptotic fractions and cells viability obtained from stochastic simulations using in vitro approaches. Taken together, we demonstrated that the model can be successfully used in prediction of cellular behaviour after exposure to IR. Thus, our studies may provide further insights into key elements involved in the underlying mechanisms of the DDR.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0293-0) contains supplementary material, which is available to authorized users.

Differences in predictions of ODE models of tumor growth: a cautionary example

BACKGROUND

While mathematical models are often used to predict progression of cancer and treatment outcomes, there is still uncertainty over how to best model tumor growth. Seven ordinary differential equation (ODE) models of tumor growth (exponential, Mendelsohn, logistic, linear, surface, Gompertz, and Bertalanffy) have been proposed, but there is no clear guidance on how to choose the most appropriate model for a particular cancer.

METHODS

We examined all seven of the previously proposed ODE models in the presence and absence of chemotherapy. We derived equations for the maximum tumor size, doubling time, and the minimum amount of chemotherapy needed to suppress the tumor and used a sample data set to compare how these quantities differ based on choice of growth model.

RESULTS

We find that there is a 12-fold difference in predicting doubling times and a 6-fold difference in the predicted amount of chemotherapy needed for suppression depending on which growth model was used.

CONCLUSION

Our results highlight the need for careful consideration of model assumptions when developing mathematical models for use in cancer treatment planning.

Oleiferasaponin C6 from the seeds of Camellia oleifera Abel.: a novel compound inhibits proliferation through inducing cell-cycle arrest and apoptosis on human cancer cell lines in vitro

Oleiferasaponin C₆ was isolated from Camellia oleifera Abel. and inhibits proliferation through inducing cell-cycle arrest and apoptosis on cancer cell lines in vitro.

A Novel Naphthalimide Compound Restores p53 Function in Non-small Cell Lung Cancer by Reorganizing the Bak·Bcl-xl Complex and Triggering Transcriptional Regulation

p53 inactivation is a hallmark in non-small-cell lung cancer (NSCLC). It is therefore highly desirable to develop tumor-specific treatment for NSCLC therapy by restoring p53 function. Herein, a novel naphthalimide compound, NA-17, was identified as a promising drug candidate in view of both its anticancer activity and mechanism of action. NA-17 exhibited strong anticancer activity on a broad range of cancer cell lines but showed low toxicity to normal cell lines, such as HL-7702 and WI-38. Moreover, NA-17 showed p53-dependent inhibition selectivity in different NSCLC cell lines due to the activation state of endogenous p53 in the background level. Further studies revealed that NA-17 caused cell cycle arrest at the G1 phase, changed cell size, and induced apoptosis and cell death by increasing the proportion of sub-G1 cells. Molecular mechanism studies suggested that targeted accumulation of phospho-p53 in mitochondria and nuclei induced by NA-17 resulted in activation of Bak and direct binding of phospho-p53 to the target DNA sequences, thereby evoking cell apoptosis and cell cycle arrest and eventually leading to irreversible cancer cell inhibition. This work provided new insights into the molecular interactions and anticancer mechanisms of phospho-p53-dependent naphthalimide compounds.

PDL1 Regulation by p53 via miR-34

BACKGROUND

Although clinical studies have shown promise for targeting PD1/PDL1 signaling in non-small cell lung cancer (NSCLC), the regulation of PDL1 expression is poorly understood. Here, we show that PDL1 is regulated by p53 via miR-34.

METHODS

p53 wild-type and p53-deficient cell lines (p53(-/-) and p53(+/+) HCT116, p53-inducible H1299, and p53-knockdown H460) were used to determine if p53 regulates PDL1 via miR-34. PDL1 and miR-34a expression were analyzed in samples from patients with NSCLC and mutated p53 vs wild-type p53 tumors from The Cancer Genome Atlas for Lung Adenocarcinoma (TCGA LUAD). We confirmed that PDL1 is a direct target of miR-34 with western blotting and luciferase assays and used a p53(R172HΔ)g/+K-ras(LA1/+) syngeneic mouse model (n = 12) to deliver miR-34a-loaded liposomes (MRX34) plus radiotherapy (XRT) and assessed PDL1 expression and tumor-infiltrating lymphocytes (TILs). A two-sided t test was applied to compare the mean between different treatments.

RESULTS

We found that p53 regulates PDL1 via miR-34, which directly binds to the PDL1 3' untranslated region in models of NSCLC (fold-change luciferase activity to control group, mean for miR-34a = 0.50, SD = 0.2, P < .001; mean for miR-34b = 0.52, SD = 0.2, P = .006; and mean for miR-34c = 0.59, SD = 0.14, and P = .006). Therapeutic delivery of MRX34, currently the subject of a phase I clinical trial, promoted TILs (mean of CD8 expression percentage of control group = 22.5%, SD = 1.9%; mean of CD8 expression percentage of MRX34 = 30.1%, SD = 3.7%, P = .016, n = 4) and reduced CD8(+)PD1(+) cells in vivo (mean of CD8/PD1 expression percentage of control group = 40.2%, SD = 6.2%; mean of CD8/PD1 expression percentage of MRX34 = 20.3%, SD = 5.1%, P = .001, n = 4). Further, MRX34 plus XRT increased CD8(+) cell numbers more than either therapy alone (mean of CD8 expression percentage of MRX34 plus XRT to control group = 44.2%, SD = 8.7%, P = .004, n = 4). Finally, miR-34a delivery reduced the numbers of radiation-induced macrophages (mean of F4-80 expression percentage of control group = 52.4%, SD = 1.7%; mean of F4-80 expression percentage of MRX34 = 40.1%, SD = 3.5%, P = .008, n = 4) and T-regulatory cells.

CONCLUSIONS

We identified a novel mechanism by which tumor immune evasion is regulated by p53/miR-34/PDL1 axis. Our results suggest that delivery of miRNAs with standard therapies, such as XRT, may represent a novel therapeutic approach for lung cancer.

The reversal of antineoplastic drug resistance in cancer cells by β-elemene

Multidrug resistance (MDR), defined as the resistance of cancer cells to compounds with diverse structures and mechanisms of actions, significantly limits the efficacy of antitumor drugs. A major mechanism that mediates MDR in cancer is the overexpression of adenosine triphosphate (ATP)-binding cassette transporters. These transporters bind to their respective substrates and catalyze their efflux from cancer cells, thereby lowering the intracellular concentrations of the substrates and thus attenuating or even abolishing their efficacy. In addition, cancer cells can become resistant to drugs via mechanisms that attenuate apoptosis and cell cycle arrest such as alterations in the p53, check point kinase, nuclear factor kappa B, and the p38 mitogen-activated protein kinase pathway. In this review, we discuss the mechanisms by which β-elemene, a compound extracted from Rhizoma zedoariae that has clinical antitumor efficacy, overcomes drug resistance in cancer.

The Antitumor Activity of the Novel Compound Jesridonin on Human Esophageal Carcinoma Cells

Jesridonin, a small molecule obtained through the structural modification of Oridonin, has extensive antitumor activity. In this study, we evaluated both its in vitro activity in the cancer cell line EC109 and its in vivo effect on tumor xenografts in nude mice. Apoptosis induced by Jesridonin was determined using an MTT assay, Annexin-V FITC assay and Hoechest 33258 staining. Apoptosis via mitochondrial and death receptor pathways were confirmed by detecting the regulation of MDM2, p53, and Bcl-2 family members and by activation of caspase-3/-8/-9. In addition, vena caudalis injection of Jesridonin showed significant inhibition of tumor growth in the xenograft model, and Jesridonin-induced cell apoptosis in tumor tissues was determined using TUNEL. Biochemical serum analysis of alkaline phosphatase (ALP), alanine transaminase (ALT), aspartate transaminase (AST), gamma-glutamyl transferase (GGT), total protein (TP) and albumin (ALB) indicated no obvious effects on liver function. Histopathological examination of the liver, kidney, lung, heart and spleen revealed no signs of JD-induced toxicity. Taken together, these results demonstrated that Jesridonin exhibits antitumor activity in human esophageal carcinomas EC109 cells both in vitro and in vivo and demonstrated no adverse effects on major organs in nude mice. These studies provide support for new drug development.

Regulation of p53 in NIH3T3 mouse fibroblasts following hyperosmotic stress

The aim of this project was to analyze the regulation of p53 expression in NIH3T3 fibroblasts under the influence of increasing hyperosmotic stress. Expression of p53 showed a biphasic response pattern in NIH3T3 cells under increasing osmotic stress (337 mOsm to 737 mOsm) with a maximum at 587 mOsm. Under isotonic conditions p53 expression increased after addition of the proteasome inhibitor MG132 indicating that cellular p53 levels in unperturbed cells is kept low by proteasomal degradation. However, under hypertonic conditions p53 synthesis as well as p53 degradation were significantly reduced and it is demonstrated that the increase in p53 expression observed when tonicity is increased from 337 to 587 mOsm reflects that degradation is more inhibited than synthesis, whereas the decrease in p53 expression at higher tonicities reflects that synthesis is more inhibited than degradation. The activity of the p53 regulating proteins p38 MAP kinase and the ubiquitin ligase MDM2 were studied as a function of increasing osmolarity. MDM2 protein expression was unchanged at all osmolarities, whereas MDM2 phosphorylation (Ser¹⁶⁶) increased at osmolarities up to 537 mOsm and remained constant at higher osmolarities. Phosphorylation of p38 increased at osmolarities up to 687 mOsm which correlated with an increased phosphorylation of p53 (Ser¹⁵) and the decreased p53 degradation. Caspase-3 activity increased gradually with hypertonicity and at 737 mOsm both Caspase-3 activity and annexin V binding are high even though p53 expression and activity are low, indicating that initiation of apoptosis under severe hypertonic conditions is not strictly controlled by p53.

Involvement of dysregulated Wip1 in manganese-induced p53 signaling and neuronal apoptosis

Overexposure to manganese (Mn) has been known to induce neuronal death and neurodegenerative symptoms. However, the precise mechanisms underlying Mn neurotoxicity remain incompletely understood. In the present study, we established a Mn-exposed rat model and found that downregulation of wild type p53-induced phosphatase 1 (Wip1) might contribute to p53 activation and resultant neuronal apoptosis following Mn exposure. Western blot and immunohistochemical analyses revealed that the expression of Wip1 was markedly decreased following Mn exposure. In addition, immunofluorescence assay demonstrated that Mn exposure led to significant reduction in the number of Wip1-positive neurons. Accordingly, the expression of Mdm2 was progressively decreased, which was accompanied with markedly increased expression of p53, as well as the ratio of Bax/Bcl-xl. Furthermore, we showed that Mn exposure decreased the viability and induced apparent apoptosis in NFG-differentiated neuron-like PC12 cells. Importantly, the expression of Wip1 decreased progressively, whereas the level of cellular p53 and the ratio of Bax/Bcl-xl were elevated, which resembled the expression of the proteins in animal model studies. Depletion of p53 significantly ameliorated Mn-mediated cytotoxic effect in PC12 cells. In addition, ectopic expression of Wip1 attenuated Mn-induced p53 signaling as well as apoptosis in PC12 cells. Finally, we observed that depletion of Wip1 augmented Mn-induced apoptosis in PC12 cells. Collectively, these findings suggest that downregulated Wip1 expression plays an important role in Mn-induced neuronal death in the brain striatum via the modulation of p53 signaling.

Physiological role of taurine--from organism to organelle

Taurine is often referred to as a semi-essential amino acid as newborn mammals have a limited ability to synthesize taurine and have to rely on dietary supply. Taurine is not thought to be incorporated into proteins as no aminoacyl tRNA synthetase has yet been identified and is not oxidized in mammalian cells. However, taurine contributes significantly to the cellular pool of organic osmolytes and has accordingly been acknowledged for its role in cell volume restoration following osmotic perturbation. This review describes taurine homeostasis in cells and organelles with emphasis on taurine biophysics/membrane dynamics, regulation of transport proteins involved in active taurine uptake and passive taurine release as well as physiological processes, for example, development, lung function, mitochondrial function, antioxidative defence and apoptosis which seem to be affected by a shift in the expression of the taurine transporters and/or the cellular taurine content.

SC-III3, a novel scopoletin derivative, induces cytotoxicity in hepatocellular cancer cells through oxidative DNA damage and ataxia telangiectasia-mutated nuclear protein kinase activation

Background

Natural products from plants have been proven to be important resources of antitumor agents. In this study, we exploited the antitumor activity of (E)-3-(4-chlorophenyl)-N-(7-hydroxy-6-methoxy-2-oxo-2H-chromen-3-yl) acrylamide (SC-III3), a newly synthesized derivative of scopoletin, by in vitro and in vivo experiments.

Methods

Human hepatocellular carcinoma cell line HepG2 cells and xenograft of HepG2 cells in BALB/c nude mice were used to investigate the effects of SC-III3 on hepatocellular cancers. Cell cycle arrest and apoptosis were analyzed by flow cytometry. Cell cycle arrest, apoptosis and ATM-Chk pathway-related proteins were characterized by western blot.

Results

SC-III3 selectively inhibited the viability of HepG2 cells without significant cytotoxicity against human normal liver cells LO2. In mouse xenograft model of HepG2 cells, SC-III3 showed a marked inhibition of tumor growth in a dose-dependent manner. Cell cycle analysis revealed that SC-III3 induced cells to accumulate in S phase, which was accompanied by a marked decrease of the expressions of cyclin A, cyclin B, cyclin E and Cdk2 proteins, the crucial regulators of S phase cell cycle. SC-III3 treatment resulted in DNA breaks in HepG2 cells, which might contribute to its S phase arrest. The S arrest and the activation of ATM-Chk1/Chk2-Cdc25A-Cdk2 pathways induced by SC-III3 in HepG2 cells could be efficiently abrogated by pretreatments of either Ku55933 (an inhibitor of ATM) or UCN-01 (an inhibitor of Chk1/Chk2). The activation of p53-p21 pathway by SC-III3 was also reversed by Ku55933 treatment. SC-III3 led to significant accumulation of intracellular reactive oxygen species (ROS), a breaker of DNA strand, in HepG2 cells but not LO2 cells. Pretreatment with N-acetyl-l-cysteine (NAC), a ROS scavenger, could reverse SC-III3-caused ROS accumulation, DNA damage, activation of signal pathways relevant to DNA damage, S phase arrest and cell viability decrease in HepG2 cells.

Conclusion

SC-III3 is able to efficiently inhibit the growth of hepatocellular carcinoma through inducing the generation of intracellular ROS, DNA damage and consequent S phase arrest, but lack of significant cytotoxicity against normal liver cells. This compound deserves further studies as a candidate of anticancer drugs.

Knockdown of PSF1 expression inhibits cell proliferation in lung cancer cells in vitro

Partner of sld five 1 (PSF1) is a member of the heterotetrameric complex termed GINS. Previous studies have shown that PSF1 is unregulated in several cancer and associated with tumor malignant characters. However, the effects of PSF1 in lung cancer are still unclear. The goal of this study was to investigate the effects of PSF1 on the proliferation capacities of lung cancer. To start with, expression of PSF1 in 22 human lung cancer samples and adjacent non-tumor samples were detected by real-time RT-PCR and Western blotting. Our results showed that PSF1 was overexpressed in lung cancer samples compared to adjacent non-tumor samples. To achieve better insights of PSF1 functions in lung cancer cells, we used PSF1-specific small interfering RNA (siRNA) successfully inhibit the expression of PSF1 in messenger RNA (mRNA) and protein levels. In addition, we used lung cancer cell lines with different p53 gene background (p53 null and p53 wild-type). The results showed that knockdown of PSF1 inhibited cell proliferation and caused cell cycle arrest of lung cancer cells in a p53-independent manner. Our data indicated that PSF1 is functionally involved in lung cancer cell proliferation and is a potential target for lung cancer therapy.

Dynamic multiprotein assemblies shape the spatial structure of cell signaling

Cell signaling underlies critical cellular decisions. Coordination, efficiency as well as fail-safe mechanisms are key elements. How the cell ensures that these hallmarks are at play are important questions. Cell signaling is often viewed as taking place through discrete and cross-talking pathways; oftentimes these are modularized to emphasize distinct functions. While simple, convenient and clear, such models largely neglect the spatial structure of cell signaling; they also convey inter-modular (or inter-protein) spatial separation that may not exist. Here our thesis is that cell signaling is shaped by a network of multiprotein assemblies. While pre-organized, the assemblies and network are loose and dynamic. They contain transiently-associated multiprotein complexes which are often mediated by scaffolding proteins. They are also typically anchored in the membrane, and their continuum may span the cell. IQGAP1 scaffolding protein which binds proteins including Raf, calmodulin, Mek, Erk, actin, and tens more, with actin shaping B-cell (and likely other) membrane-anchored nanoclusters and allosterically polymerizing in dynamic cytoskeleton formation, and Raf anchoring in the membrane along with Ras, provides a striking example. The multivalent network of dynamic proteins and lipids, with specific interactions forming and breaking, can be viewed as endowing gel-like properties. Collectively, this reasons that efficient, productive and reliable cell signaling takes place primarily through transient, preorganized and cooperative protein-protein interactions spanning the cell rather than stochastic, diffusion-controlled processes.

Transient expression of RAD51 in the late G2-phase is required for cell cycle progression in synchronous Physarum cells

The homologous recombination factor RAD51 is highly conserved. This criterion enabled us to identify a RAD51 ortholog in Physarum polycephalum. We found that the Physarum protein presents a high homology to the human protein and cross-reacted with antibodies directed against the human RAD51. Taking advantage of the natural synchrony of millions of nuclei within a single cell of Physarum, we investigated the fluctuation of the amount of the PpRAD51 throughout the cell cycle. Our results showed that in the late G2-phase, RAD51 was transiently expressed in a large quantity. Furthermore, knocking-down RAD51 in the G2-phase abolished this transient expression before mitosis and affected cell cycle progression. These results support the idea that RAD51 plays a role in the progression of the cell cycle in the late G2-phase.

The dynamics of p53 in single cells: physiologically based ODE and reaction-diffusion PDE models

The intracellular signalling network of the p53 protein plays important roles in genome protection and the control of cell cycle phase transitions. Recently observed oscillatory behaviour in single cells under stress conditions has inspired several research groups to simulate and study the dynamics of the protein with the aim of gaining a proper understanding of the physiological meanings of the oscillations. We propose compartmental ODE and PDE models of p53 activation and regulation in single cells following DNA damage and we show that the p53 oscillations can be retrieved by plainly involving p53-Mdm2 and ATM-p53-Wip1 negative feedbacks, which are sufficient for oscillations experimentally, with no further need to introduce any delays into the protein responses and without considering additional positive feedback.

Inhibition of G9a induces DUSP4-dependent autophagic cell death in head and neck squamous cell carcinoma

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) is a common cancer worldwide. Emerging evidence indicates that alteration of epigenetics might be a key event in HNSCC progression. Abnormal expression of histone methyltransferase G9a, which contributes to transcriptional repression of tumor suppressors, has been implicated in promoting cancerous malignancies. However, its role in HNSCC has not been previously characterized. In this study, we elucidate the function of G9a and its downstream mechanism in HNSCC.

METHODS

We investigated the clinical relevance of G9a in HNSCC using immunohistochemistry (IHC) staining. In vitro cell proliferation and tumorigenesis ability of G9a-manipulated HNSCC cells were examined with MTT assays, clonogenic assays, and soft agar assays. We examined different routes of cell death in HNSCC cells induced by G9a-depletion or enzymatic inhibition by immunoblot, flow cytometry, fluorescent and transmission electron microscopy analysis. Specific targets of G9a were identified by affymetrix microarray and quantitative reverse transcription-polymerase chain reaction (qRT-PCR). Lastly, functions of G9a in vivo were confirmed with a xenograft tumor model.

RESULTS

G9a expression is positively correlated to proliferation marker Ki-67 and to poor prognosis in HNSCC patients. Genetic or pharmacological inhibition of G9a reduced cell proliferation without inducing necrosis or apoptosis. Instead, autophagic cell death was the major consequence, and our investigation of mechanisms suggested it is mediated via the dual specificity phosphatase-4 (DUSP4) dependent ERK inactivation pathway. An orthotopic tumor model further confirmed the growth inhibiting effect and induction of autophagy that followed suppression of G9a.

CONCLUSIONS

In this study, we provide evidence that G9a confers the survival advantage of HNSCC. Genetic or pharmacological inhibition of G9a induces autophagic cell death; this finding provides a basis for new therapeutic targets for treating HNSCC.

Reaction-diffusion systems for spatio-temporal intracellular protein networks: A beginner's guide with two examples

Spatio-temporal dynamics of a variety of proteins is, among other things, regulated by post-translational modifications of these proteins. Such modifications can thus influence stability and biochemical activities of the proteins, activity and stability of their upstream targets within specific signalling pathways. Commonly used mathematical tools for such protein–protein (and/or protein-mRNA) interactions in single cells, namely, Michaelis–Menten and Hill kinetics, yielding a system of ordinary differential equations, are extended here into (non-linear) partial differential equations by taking into account a more realistic spatial representation of the environment where these reactions occur. In the modelling framework under consideration, all interactions occur in a cell divided into two compartments, the nucleus and the cytoplasm, connected by the semipermeable nuclear membrane and bounded by the impermeable cell membrane. Passive transport mechanism, modelled by the so-called Kedem–Katchalsky boundary conditions, is used here to represent migration of species throughout the nuclear membrane. Nonlinear systems of partial differential equations are solved by the semi-implicit Rothe method. Examples of two spatial oscillators are shown. Namely, these are the circadian rhythm for concentration of the FRQ protein in Neurospora crassa and oscillatory dynamics observed in the activation and regulation of the p53 protein following DNA damage in mammalian cells.

Benzo(a)pyrene-7,8-diol-9,10-epoxide induced p53-independent necrosis via the mitochondria-associated pathway involving Bax and Bak activation

Benzo(a)pyrene-7,8-diol-9,10-epoxide (BPDE) is a highly reactive DNA damage agent and can induce cell death through both p53-independent and -dependent pathways. However, little is known about the molecular mechanisms of p53-independent pathways in BPDE-induced cell death. To understand the p53-independent mechanisms, we have now examined BPDE-induced cytotoxicity in p53-deficient baby mouse kidney (BMK) cells. The results showed that BPDE could induce Bax and Bak activation, cytochrome c release, caspases activation, and necrotic cell death in the BMK cells. Bax and Bak, two key molecules of mitochondrial permeability transition pore, were interdependently activated by BPDE, with Bax and Bak translocation to and Bax/Bak homo-oligomerization in mitochondria, release of cytochrome c was induced. Importantly, cytochrome c release and necrotic cell death were diminished in BMK cells (Bax−/−), BMK cells (Bak−/−), and BMK cells (Bax−/−/Bak−/−). Furthermore, overexpression of Bcl-2 could ameliorate BPDE-induced cytochrome c release and necrosis. Together the findings suggested that BPDE-induced necrosis was modulated by the p53-independent pathway, which was related to the translocation of Bax and Bak to mitochondria, release of cytochrome c, and activation of caspases.

Dynamics of posttranslational modifications of p53

The latest experimental evidence indicates that acetylation of p53 at K164 (lysine 164) and K120 may induce directly cell apoptosis under severe DNA damage. However, previous cell apoptosis models only studied the effects of active and/or inactive p53, that is, phosphorylation/dephosphorylation of p53. In the present paper, based partly on Geva-Zatorsky et al. (2006) and Batchelor et al. (2008), we propose a new cell apoptosis network, in which p53 has three statuses, that is, unphosphorylated p53, phosphorylated p53, and acetylated p53. The time delay differential equations (DDEs) are formulated based on our network to investigate the dynamical insights of p53-induced cell apoptosis. In agreement with experiments (Loewer et al. (2010)), our simulations indicate that acetylated p53 accumulates gradually and then induces the proapoptotic protein Bax under enough DNA damage. Moreover, phosphorylated p53 oscillates and initiates cell repair during DNA damage.