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

Encoding Growth Factor Identity in the Temporal Dynamics of FOXO3 under the Combinatorial Control of ERK and AKT Kinases

Extracellular growth factors signal to transcription factors via a limited number of cytoplasmic kinase cascades. It remains unclear how such cascades encode ligand identities and concentrations. In this paper, we use live-cell imaging and statistical modeling to study FOXO3, a transcription factor regulating diverse aspects of cellular physiology that is under combinatorial control. We show that FOXO3 nuclear-to-cytosolic translocation has two temporally distinct phases varying in magnitude with growth factor identity and cell type. These phases comprise synchronous translocation soon after ligand addition followed by an extended back-and-forth shuttling; this shuttling is pulsatile and does not have a characteristic frequency, unlike a simple oscillator. Early and late dynamics are differentially regulated by Akt and ERK and have low mutual information, potentially allowing the two phases to encode different information. In cancer cells in which ERK and Akt are dysregulated by oncogenic mutation, the diversity of states is lower.

Cidofovir is active against human papillomavirus positive and negative head and neck and cervical tumor cells by causing DNA damage as one of its working mechanisms

Human papillomavirus (HPV) causes cervical cancer and a large fraction of head and neck squamous cell carcinomas (HNSCC). Cidofovir (CDV) proved efficacious in the treatment of several HPV-induced benign and malignant hyper proliferations. To provide a better insight into how CDV selectively eradicates transformed cells, HPV+ and HPV- cervical carcinoma and HNSCC cell lines were compared to normal cells for antiproliferative effects, CDV metabolism, drug incorporation into cellular DNA, and DNA damage. Incorporation of CDV into cellular DNA was higher in tumor cells than in normal cells and correlated with CDV antiproliferative effects, which were independent of HPV status. Increase in phospho-ATM levels was detected following CDV exposure and higher levels of γ-H2AX (a quantitative marker of double-strand breaks) were measured in tumor cells compared to normal cells. A correlation between DNA damage and CDV incorporation into DNA was found but not between DNA damage and CDV antiproliferative effects. These data indicate that CDV antiproliferative effects result from incorporation of the drug into DNA causing DNA damage. However, the anti-tumor effects of CDV cannot be exclusively ascribed to DNA damage. Furthermore, CDV can be considered a promising broad spectrum anti-cancer agent, not restricted to HPV+ lesions.

Optimal Fault Detection and Diagnosis in Transcriptional Circuits Using Next-Generation Sequencing

We propose a methodology for model-based fault detection and diagnosis for stochastic Boolean dynamical systems indirectly observed through a single time series of transcriptomic measurements using Next Generation Sequencing (NGS) data. The fault detection consists of an innovations filter followed by a fault certification step, and requires no knowledge about the possible system faults. The innovations filter uses the optimal Boolean state estimator, called the Boolean Kalman Filter (BKF). In the presence of knowledge about the possible system faults, we propose an additional step of fault diagnosis based on a multiple model adaptive estimation (MMAE) method consisting of a bank of BKFs running in parallel. Performance is assessed by means of false detection and misdiagnosis rates, as well as average times until correct detection and diagnosis. The efficacy of the proposed methodology is demonstrated via numerical experiments using a p53-MDM2 negative feedback loop Boolean network with stuck-at faults that model molecular events commonly found in cancer.

Understanding non-linear effects from Hill-type dynamics with application to decoding of p53 signaling

Analytical equations are derived depicting four possible scenarios resulting from pulsed signaling of a system subject to Hill-type dynamics. Pulsed Hill-type dynamics involves the binding of multiple signal molecules to a receptor and occurs e.g., when transcription factor p53 orchestrates cancer prevention, during calcium signaling, and during circadian rhythms. The scenarios involve: (i) enhancement of high-affinity binders compared to low-affinity ones, (ii) slowing reactions involving high-affinity binders, (iii) transfer of the clocking of low-affinity binders from the signal molecule to the products, and (iv) a unique clocking process that produces incremental increases in the activity of high-affinity binders with each signal pulse. In principle, these mostly non-linear effects could control cellular outcomes. An applications to p53 signaling is developed, with binding to most gene promoters identified as category (iii) responses. However, currently unexplained enhancement of high-affinity promoters such as CDKN1a (p21) by pulsed signaling could be an example of (i). In general, provision for all possible scenarios is required in the design of mathematical models incorporating pulsed Hill-type signaling as some aspect.

Revealing determinants of two-phase dynamics of P53 network under gamma irradiation based on a reduced 2D relaxation oscillator model

This study proposes a two-dimensional (2D) oscillator model of p53 network, which is derived via reducing the multidimensional two-phase dynamics model into a model of ataxia telangiectasia mutated (ATM) and Wip1 variables, and studies the impact of p53-regulators on cell fate decision. First, the authors identify a 6D core oscillator module, then reduce this module into a 2D oscillator model while preserving the qualitative behaviours. The introduced 2D model is shown to be an excitable relaxation oscillator. This oscillator provides a mechanism that leads diverse modes underpinning cell fate, each corresponding to a cell state. To investigate the effects of p53 inhibitors and the intrinsic time delay of Wip1 on the characteristics of oscillations, they introduce also a delay differential equation version of the 2D oscillator. They observe that the suppression of p53 inhibitors decreases the amplitudes of p53 oscillation, though the suppression increases the sustained level of p53. They identify Wip1 and P53DINP1 as possible targets for cancer therapies considering their impact on the oscillator, supported by biological findings. They model some mutations as critical changes of the phase space characteristics. Possible cancer therapeutic strategies are then proposed for preventing these mutations' effects using the phase space approach.

p53 elevation in human cells halt SV40 infection by inhibiting T-ag expression

SV40 large T-antigen (T-ag) has been known for decades to inactivate the tumor suppressor p53 by sequestration and additional mechanisms. Our present study revealed that the struggle between p53 and T-ag begins very early in the infection cycle. We found that p53 is activated early after SV40 infection and defends the host against the infection. Using live cell imaging and single cell analyses we found that p53 dynamics are variable among individual cells, with only a subset of cells activating p53 immediately after SV40 infection. This cell-to-cell variabilty had clear consequences on the outcome of the infection. None of the cells with elevated p53 at the beginning of the infection proceeded to express T-ag, suggesting a p53-dependent decision between abortive and productive infection. In addition, we show that artificial elevation of p53 levels prior to the infection reduces infection efficiency, supporting a role for p53 in defending against SV40. We further found that the p53-mediated host defense mechanism against SV40 is not facilitated by apoptosis nor via interferon-stimulated genes. Instead p53 binds to the viral DNA at the T-ag promoter region, prevents its transcriptional activation by Sp1, and halts the progress of the infection. These findings shed new light on the long studied struggle between SV40 T-ag and p53, as developed during virus-host coevolution. Our studies indicate that the fate of SV40 infection is determined as soon as the viral DNA enters the nucleus, before the onset of viral gene expression.

The relationship between DNA single-stranded damage response and double-stranded damage response

The damage response of DNA single-stranded breaks(SSBs) and double-stranded breaks(DSBs) are two relatively independent processes involving different signaling pathways and protein factors, but there are still many overlapping parts. All of them can activate p53 protein, then the activated p53 regulates the damage response of single-stranded breaks or double-stranded breaks in transcriptional regulation and non-transcriptional regulation. Especially, the two types of damage would compete for RPA and ATR resources in damage repair process. The research has been focused on damage response of DNA single-stranded breaks or DNA double-stranded breaks. However, when single-stranded breaks and double-stranded breaks exist simultaneously, the DNA damage response remains to be elucidated. Here, we present a hybrid numerical model of p53 response and a hybrid numerical model of DNA damage repair exploring DNA damage repair and apoptosis mechanisms when DNA single-stranded breaks and DNA double-stranded breaks exist simultaneously. Firstly, when two kinds of damage are present at the same time, the response of p53 is graded, it means that p53 responds to single-stranded breaks preferentially; Secondly, DNA single-stranded breaks are repaired preferentially, and single-stranded breaks and double-stranded breaks can be repaired simultaneously after most of single-stranded breaks having been repaired; Moreover, single-stranded breaks are more likely to cause apoptosis, because the accumulation of p53 in DNA single-stranded breaks is faster than it in DNA double-stranded breaks and single-stranded breaks has lower threshold of apoptosis.

Information processing in the NF-κB pathway

The NF-κB pathway is known to transmit merely 1 bit of information about stimulus level. We combined experimentation with mathematical modeling to elucidate how information about TNF concentration is turned into a binary decision. Using Kolmogorov-Smirnov distance, we quantified the cell’s ability to discern 8 TNF concentrations at each step of the NF-κB pathway, to find that input discernibility decreases as signal propagates along the pathway. Discernibility of low TNF concentrations is restricted by noise at the TNF receptor level, whereas discernibility of high TNF concentrations it is restricted by saturation/depletion of downstream signaling components. Consequently, signal discernibility is highest between 0.03 and 1 ng/ml TNF. Simultaneous exposure to TNF or LPS and a translation inhibitor, cycloheximide, leads to prolonged NF-κB activation and a marked increase of transcript levels of NF-κB inhibitors, IκBα and A20. The impact of cycloheximide becomes apparent after the first peak of nuclear NF-κB translocation, meaning that the NF-κB network not only relays 1 bit of information to coordinate the all-or-nothing expression of early genes, but also over a longer time course integrates information about other stimuli. The NF-κB system should be thus perceived as a feedback-controlled decision-making module rather than a simple information transmission channel.

p53 dynamics orchestrates with binding affinity to target genes for cell fate decision

Emerging evidence support that temporal dynamics is pivotal for signaling molecules in orchestrating smart responses to diverse stimuli. p53 is such a signaling molecule that employs temporal dynamics for the selective activation of downstream target genes and ultimately for cell fate decision. Yet how this fine-tuned p53 machinery is quantitatively decoded remains largely unclear. Here we report a quantitative mechanism defining how p53 dynamics orchestrates with binding affinity to target genes for cell fate decision. Treating cells with a genotoxic drug doxorubicin at various doses and durations, we found that a mild and prolonged challenge triggered sequential p53 pulses and ultimately resulted in a terminal pulse enacting apoptosis in a comparable rate with that induced by an acute and high-dose treatment. To transactivate proapoptotic genes and thereafter executing apoptosis, p53 must exceed a certain threshold and accumulate for sufficient time at levels above it. Effective cumulative levels above the threshold, defined as E∫p53, but not the total accumulation levels of p53, precisely discriminate survival and apoptotic cells. p53 accumulation below this threshold, even with prolonging time to reach a total level comparable to that from the accumulation over the threshold, could not transactivate proapoptotic genes to which the binding affinity of p53 is lower than that of proarrest genes, and this property is independent of dynamic features. Our findings indicate that the dynamic feature per se does not directly control cell fate, but rather it orchestrates with the binding affinity to target genes to confer an appropriate time window for cell fate choice. Our study provides a quantitative mechanism unifying p53 dynamics and binding affinity to target genes, providing novel insights to understand how p53 can respond quantitatively to chemotherapeutic drugs, and guiding the design of metronomic regimens for chemotherapeutic drugs.

MicroRNA-16 feedback loop with p53 and Wip1 can regulate cell fate determination between apoptosis and senescence in DNA damage response

Cell fate regulation is an open problem whose comprehension impacts several areas of the biosciences. DNA damage induces cell cycle checkpoints that activate the p53 pathway to regulate cell fate mechanisms such as apoptosis or senescence. Experiments with different cell types show that the p53 pathway regulates cell fate through a switch behavior in its dynamics. For low DNA damage the pathway presents an oscillatory pattern associated with intense DNA damage repair while for high damage there are no oscillations and either p53 concentration increases inducing apoptosis or the cell enters a senescence state. Apoptosis and senescence phenotypes seem to have compensatory functions in tissues and the microRNA 16-1 (miR-16) is involved in the regulation of the fate between both phenotypes in cancer cells. To investigate the regulation of cell fate we developed a logical model of the G1/S checkpoint in DNA damage response that takes into account different levels of damage and contemplates the influence of miR-16 through its positive feedback loop formed with p53 and Wip1. The model reproduces the observed cellular phenotypes in experiments: oscillatory (for low DNA damage) regulated by negative feedback loops involving mainly p53 and Mdm2 and apoptotic or senescent (for high DNA damage) regulated by the positive p53/Wip1/miR-16 feedback loop. We find good agreement between the level of DNA damage and the probability of the phenotype produced according to experiments. We also find that this positive feedback makes senescent and apoptotic phenotypes to be determined stochastically (bistable), however controlling the expression level of miR-16 allows the control of fate determination as observed experimentally.

p53 pulses lead to distinct patterns of gene expression albeit similar DNA-binding dynamics

The dynamics of transcription factors play important roles in a variety of biological systems. However, the mechanisms by which these dynamics are decoded into different transcriptional responses are not well understood. Here we focus on the dynamics of the tumor-suppressor protein p53, which exhibits a series of pulses in response to DNA damage. We performed time course RNA sequencing (RNA-seq) and chromatin immunoprecipitation sequencing (ChIP-seq) measurements to determine how p53 oscillations are linked with gene expression genome wide. We discovered multiple distinct patterns of gene expression in response to p53 pulses. Surprisingly, p53-binding dynamics were uniform across all genomic loci, even for genes that exhibited distinct mRNA dynamics. Using a mathematical model, supported by additional experimental measurements in response to sustained p53 input, we determined that p53 binds to and activates transcription of its target genes uniformly, whereas post-transcriptional mechanisms are responsible for the differences in gene expression dynamics.

Global Inhibition with Specific Activation: How p53 and MYC Redistribute the Transcriptome in the DNA Double-Strand Break Response

In response to stresses, cells often halt normal cellular processes, yet stress-specific pathways must bypass such inhibition to generate effective responses. We investigated how cells redistribute global transcriptional activity in response to DNA damage. We show that an oscillatory increase of p53 levels in response to double-strand breaks drives a counter-oscillatory decrease of MYC levels. Using RNA sequencing (RNA-seq) of newly synthesized transcripts, we found that p53-mediated reduction of MYC suppressed general transcription, with the most highly expressed transcripts reduced to a greater extent. In contrast, upregulation of p53 targets was relatively unaffected by MYC suppression. Reducing MYC during the DNA damage response was important for cell-fate regulation, as counteracting MYC repression reduced cell-cycle arrest and elevated apoptosis. Our study shows that global inhibition with specific activation of transcriptional pathways is important for the proper response to DNA damage; this mechanism may be a general principle used in many stress responses.

Effect of HPV E6/E7 siRNA with Chemotherapeutic Agents on the Regulation of TP53/E2F Dynamic Behavior for Cell Fate Decisions

Toxicity and resistance remain major challenges for advanced or recurrent cervical cancer therapies, as treatment requires high doses of chemotherapeutic agents. Restoration of TP53 and hypophosphorylated-retinoblastoma (pRB) proteins by human papillomavirus (HPV) E6/E7 siRNA sensitizes HPV-positive cervical cancer cells toward chemotherapeutic agents. Here, we investigated the therapeutic effects of E6/E7 siRNA on the dynamic behavior of TP53 and RB/E2F signaling networks in deciding the cell fate. The synergistic effect of HPV E6/E7 siRNA pool (SP) with chemotherapeutic agents on TP53 and RB/E2F signaling, proliferation, and apoptosis was analyzed in vitro and in vivo. Compared to the E6/E7 SP alone, E6/E7 SP with cisplatin treatment effectively restored TP53 and RB/E2F signaling and contributes to differences in cell fate, such as apoptosis or cell cycle arrest. We also developed a cellular dynamics model that includes TP53-RB/E2F dynamics and cell proliferation profiles, and confirmed its utility for investigating E6/E7 siRNA-based combination regimens. Using a dual reporter system, we further confirmed the cross talk between TP53 and RB/E2F signaling mechanisms. Treatment of E6/E7 SP cationic liposome (i.v.) with cisplatin and paclitaxel (i.p.) potentially inhibited tumor growth in BALB/c-nude mice. Altogether, our findings suggest that stabilization of TP53 and the RB/E2F repressor complex by E6/E7 SP combined with low-dose chemotherapy can effectively suppress tumor growth.

Census and evaluation of p53 target genes

The tumor suppressor p53 functions primarily as a transcription factor. Mutation of the TP53 gene alters its response pathway, and is central to the development of many cancers. The discovery of a large number of p53 target genes, which confer p53's tumor suppressor function, has led to increasingly complex models of p53 function. Recent meta-analysis approaches, however, are simplifying our understanding of how p53 functions as a transcription factor. In the survey presented here, a total set of 3661 direct p53 target genes is identified that comprise 3509 potential targets from 13 high-throughput studies, and 346 target genes from individual gene analyses. Comparison of the p53 target genes reported in individual studies with those identified in 13 high-throughput studies reveals limited consistency. Here, p53 target genes have been evaluated based on the meta-analysis data, and the results show that high-confidence p53 target genes are involved in multiple cellular responses, including cell cycle arrest, DNA repair, apoptosis, metabolism, autophagy, mRNA translation and feedback mechanisms. However, many p53 target genes are identified only in a small number of studies and have a higher likelihood of being false positives. While numerous mechanisms have been proposed for mediating gene regulation in response to p53, recent advances in our understanding of p53 function show that p53 itself is solely an activator of transcription, and gene downregulation by p53 is indirect and requires p21. Taking into account the function of p53 as an activator of transcription, recent results point to an unsophisticated means of regulation.

Significance of Wild-Type p53 Signaling in Suppressing Apoptosis in Response to Chemical Genotoxic Agents: Impact on Chemotherapy Outcome

Our genomes are subject to potentially deleterious alterations resulting from endogenous sources (e.g., cellular metabolism, routine errors in DNA replication and recombination), exogenous sources (e.g., radiation, chemical agents), and medical diagnostic and treatment applications. Genome integrity and cellular homeostasis are maintained through an intricate network of pathways that serve to recognize the DNA damage, activate cell cycle checkpoints and facilitate DNA repair, or eliminate highly injured cells from the proliferating population. The wild-type p53 tumor suppressor and its downstream effector p21^WAF1 (p21) are key regulators of these responses. Although extensively studied for its ability to control cell cycle progression, p21 has emerged as a multifunctional protein capable of downregulating p53, suppressing apoptosis, and orchestrating prolonged growth arrest through stress-induced premature senescence. Studies with solid tumors and solid tumor-derived cell lines have revealed that such growth-arrested cancer cells remain viable, secrete growth-promoting factors, and can give rise to progeny with stem-cell-like properties. This article provides an overview of the mechanisms by which p53 signaling suppresses apoptosis following genotoxic stress, facilitating repair of genomic injury under physiological conditions but having the potential to promote tumor regrowth in response to cancer chemotherapy.

p53 dynamics in response to DNA damage vary across cell lines and are shaped by efficiency of DNA repair and activity of the kinase ATM

Cellular systems show a wide range of signaling dynamics. Many of these dynamics are highly stereotyped, such as oscillations at a fixed frequency. However, most studies looking at the role of signaling dynamics focus on one or a few cell lines, leaving the diversity of dynamics across tissues or cell lines a largely unexplored question. We focused on the dynamics of the tumor suppressor protein p53, which regulates cell cycle arrest and apoptosis in response to DNA damage. We established live-cell reporters for 12 cancer cell lines expressing wild-type p53 and quantified p53 dynamics in response to double-strand break-inducing DNA damage. In many of the tested cell lines, we found that p53 abundance oscillated in response to ionizing radiation or the DNA-damaging chemotherapeutic neocarzinostatin and that the periodicity of the oscillations was fixed. In other cell lines, p53 abundance dynamically changed in different ways, such as a single broad pulse or a continuous induction. By combining single-cell assays of p53 signaling dynamics, small-molecule screening in live cells, and mathematical modeling, we identified molecules that perturbed p53 dynamics and determined that cell-specific variation in the efficiency of DNA repair and the activity of the kinase ATM (ataxia-telangiectasia mutated) controlled the signaling landscape of p53 dynamics. Because the dynamics of wild-type p53 varied substantially between cell lines, our study highlights the limitation of using one line as a model system and emphasizes the importance of studying the dynamics of other signaling pathways across different cell lines and genetic backgrounds.

Excitability in the p53 network mediates robust signaling with tunable activation thresholds in single cells

Cellular signaling systems precisely transmit information in the presence of molecular noise while retaining flexibility to accommodate the needs of individual cells. To understand design principles underlying such versatile signaling, we analyzed the response of the tumor suppressor p53 to varying levels of DNA damage in hundreds of individual cells and observed a switch between distinct signaling modes characterized by isolated pulses and sustained oscillations of p53 accumulation. Guided by dynamic systems theory we show that this requires an excitable network structure comprising positive feedback and provide experimental evidence for its molecular identity. The resulting data-driven model reproduced all features of measured signaling responses and is sufficient to explain their heterogeneity in individual cells. We present evidence that heterogeneity in the levels of the feedback regulator Wip1 sets cell-specific thresholds for p53 activation, providing means to modulate its response through interacting signaling pathways. Our results demonstrate how excitable signaling networks can provide high specificity, sensitivity and robustness while retaining unique possibilities to adjust their function to the physiology of individual cells.

Recent progress and open challenges in modeling p53 dynamics in single cells

In mammalian cells, the tumor suppressor p53 is activated upon a variety of cellular stresses and ensures an appropriate response ranging from arrest and repair to the induction of senescence and apoptosis. Quantitative measurements in individual living cells showed stimulus-dependent dynamics of p53 accumulation upon stress induction. Due to the complexity of the underlying biochemical interactions, mathematical models were indispensable for understanding the topology of the network regulating p53 dynamics. Recent work provides furhter insights into the causes of heterogeneous responses in individual cells, the rewiring of the network in response to different inputs and the role of the downstream processes in determining the cellular fate upon stress.

Modeling amplified p53 responses under DNA-PK inhibition in DNA damage response

During DNA double strand breaks (DSBs) repair, coordinated activation of phosphatidylinositol 3-kinase (PI3K)-like kinases can activate p53 signaling pathway. Recent findings have identified novel interplays among these kinases demonstrating amplified first p53 pulses under DNA-PK inhibition. However, no theoretical model has been developed to characterize such dynamics. In current work, we modeled the prolonged p53 pulses with DNA-PK inhibitor. We could identify a dose-dependent increase in the first pulse amplitude and width. Meanwhile, weakened DNA-PK mediated ATM inhibition was insufficient to reproduce such dynamic behavior. Moreover, the information flow was shifted predominantly to the first pulse under DNA-PK inhibition. Furthermore, the amplified p53 responses were relatively robust. Taken together, our model can faithfully replicate amplified p53 responses under DNA-PK inhibition and provide insights into cell fate decision by manipulating p53 dynamics.

Single-cell Gene Expression Profiling Using FACS and qPCR with Internal Standards

Gene expression measurements from bulk populations of cells can obscure the considerable transcriptomic variation of individual cells within those populations. Single-cell gene expression measurements can help assess the role of noise in gene expression, identify correlations in the expression of pairs of genes, and reveal subpopulations of cells that respond differently to a stimulus. Here, we describe a procedure to measure the expression of up to 96 genes in single mammalian cells isolated from a population growing in tissue culture. Cells are sorted into lysis buffer by fluorescence-activated cell sorting (FACS), and the mRNA species of interest are reverse-transcribed and amplified. Gene expression is then measured using a microfluidic real-time PCR machine, which performs up to 96 qPCR assays on up to 96 samples at a time. We also describe the generation and use of PCR amplicon standards to enable the estimation of the absolute number of each transcript. Compared with other methods of measuring gene expression in single cells, this approach allows for the quantification of more distinct transcripts than RNA FISH at a lower cost than RNA-Seq.

Acetylation-dependent regulation of MDM2 E3 ligase activity dictates its oncogenic function

Abnormal activation of the oncogenic E3 ubiquitin ligase murine double minute 2 (MDM2) is frequently observed in human cancers. By ubiquitinating the tumor suppressor p53 protein, which leads to its proteasome-mediated destruction, MDM2 limits the tumor-suppressing activity of p53. On the other hand, by ubiquitinating itself, MDM2 targets itself for destruction and promotes the p53 tumor suppressor pathway, a process that can be antagonized by the deubiquitinase herpesvirus-associated ubiquitin-specific protease (HAUSP). We investigated the regulation of MDM2 substrate specificity and found that acetyltransferase p300-mediated acetylation and stabilization of MDM2 are molecular switches that block self-ubiquitination, thereby shifting its E3 ligase activity toward p53. In vitro and in cancer cell lines, p300-mediated acetylation of MDM2 on Lys¹⁸² and Lys¹⁸⁵ enabled HAUSP to bind, presumably deubiquitinate, and stabilize MDM2. This acetylation within the nuclear localization signal domain decreased its interaction with the acidic domain, subsequently increased the interaction between the acidic domain and RING domain in MDM2, enabled the binding of HAUSP to the acidic domain in MDM2, and shifted MDM2 activity from autoubiquitination to p53 ubiquitination. However, upon genotoxic stress through exposure to etoposide, the deacetylase sirtuin 1 (SIRT1) deacetylated MDM2 at Lys¹⁸² and Lys¹⁸⁵, thereby promoting self-ubiquitination and less ubiquitination and subsequent degradation of p53, thus increasing p53-dependent apoptosis. Therefore, this study indicates that dynamic acetylation is a molecular switch in the regulation of MDM2 substrate specificity, revealing further insight into the posttranslational regulation of the MDM2/p53 cell survival axis.

Dynamics of p53: A Master Decider of Cell Fate

Cellular stress-induced temporal alterations—i.e., dynamics—are typically exemplified by the dynamics of p53 that serve as a master to determine cell fate. p53 dynamics were initially identified as the variations of p53 protein levels. However, a growing number of studies have shown that p53 dynamics are also manifested in variations in the activity, spatial location, and posttranslational modifications of p53 proteins, as well as the interplay among all p53 dynamical features. These are essential in determining a specific outcome of cell fate. In this review, we discuss the importance of the multifaceted features of p53 dynamics and their roles in the cell fate decision process, as well as their potential applications in p53-based cancer therapy. The review provides new insights into p53 signaling pathways and their potentials in the development of new strategies in p53-based cancer therapy.

p53 and RAD9, the DNA Damage Response, and Regulation of Transcription Networks

The way cells respond to DNA damage is important since inefficient repair or misrepair of lesions can have deleterious consequences, including mutation, genomic instability, neurodegenerative disorders, premature aging, cancer or death. Whether damage occurs spontaneously as a byproduct of normal metabolic processes, or after exposure to exogenous agents, cells muster a coordinated, complex DNA damage response (DDR) to mitigate potential harmful effects. A variety of activities are involved to promote cell survival, and include DNA repair, DNA damage tolerance, as well as transient cell cycle arrest to provide time for repair before entry into critical cell cycle phases, an event that could be lethal if traversal occurs while damage is present. When such damage is prolonged or not repairable, senescence, apoptosis or autophagy is induced. One major level of DDR regulation occurs via the orchestrated transcriptional control of select sets of genes encoding proteins that mediate the response. p53 is a transcription factor that transactivates specific DDR downstream genes through binding DNA consensus sequences usually in or near target gene promoter regions. The profile of p53-regulated genes activated at any given time varies, and is dependent upon type of DNA damage or stress experienced, exact composition of the consensus DNA binding sequence, presence of other DNA binding proteins, as well as cell context. RAD9 is another protein critical for the response of cells to DNA damage, and can also selectively regulate gene transcription. The limited studies addressing the role of RAD9 in transcription regulation indicate that the protein transactivates at least one of its target genes, p21/waf1/cip1, by binding to DNA sequences demonstrated to be a p53 response element. NEIL1 is also regulated by RAD9 through a similar DNA sequence, though not yet directly verified as a bonafide p53 response element. These findings suggest a novel pathway whereby p53 and RAD9 control the DDR through a shared mechanism involving an overlapping network of downstream target genes. Details and unresolved questions about how these proteins coordinate or compete to execute the DDR through transcriptional reprogramming, as well as biological implications, are discussed.

Phosphatase Wip1 in Immunity: An Overview and Update

Wild-type p53-induced phosphatase 1 (Wip1) is a newly identified serine/threonine phosphatase, which belongs to the PP2C family. Due to its involvement in stress-induced networks and overexpression in human tumors, primary studies have mainly focused on the role of Wip1 in tumorigenesis. It now has also been implicated in regulating several other physiological processes such as organism aging and neurogenesis. Recent evidence highlights a new role of Wip1 in controlling immune response through regulating immune cell development and function, as well as through the interplay with inflammatory signaling pathways such NF-κB and p38 mitogen-activated protein kinase. In this short review, we will give an overview of Wip1 in immunity to better understand this important phosphatase.

On p53 revival using system oriented drug dosage design

We propose a new paradigm in the drug design for the revival of the p53 pathway in cancer cells. It is shown that the current strategy of using small molecule based Mdm2 inhibitors is not enough to adequately revive p53 in cancerous cells, especially when it comes to the extracting pulsating behavior of p53. This fact has come to notice when a novel method for the drug dosage design is introduced using system oriented concepts. As a test case, small molecule drug Mdm2 repressor Nutlin 3a is considered. The proposed method determines the dose of Nutlin to revive p53 pathway functionality. For this purpose, PBK dynamics of Nutlin have also been integrated with p53 pathway model. The p53 pathway is the focus of researchers for the last thirty years for its pivotal role as a frontline cancer suppressant protein due to its effect on cell cycle checkpoints and cell apoptosis in response to a DNA strand break. That is the reason for finding p53 being absent in more than 50% of tumor cancers. Various drugs have been proposed to revive p53 in cancer cells. Small molecule based drugs are at the foremost and are the subject of advanced clinical trials. The dosage design of these drugs is an important issue. We use control systems concepts to develop the drug dosage so that the cancer cells can be treated in appropriate time. We investigate by using a computational model how p53 protein responds to drug Nutlin 3a, an agent that interferes with the MDM2-mediated p53 regulation. The proposed integrated model describes in some detail the regulation network of p53 including the negative feedback loop mediated by MDM2 and the positive feedback loop mediated by Mdm2 mRNA as well as the reversible represses of MDM2 caused by Nutlin. The reported PBK dynamics of Nutlin 3a are also incorporated to see the full effect. It has been reported that p53 response to stresses in two ways. Either it has a sustained (constant) p53 response, or there are oscillations in p53 concentration. The claimed dosage strategy achieves the p53 response in the first case. However, for the induction of oscillations, it is shown through bifurcation analysis that to achieve oscillating behavior of p53 inhibition of Mdm2 is not enough, rather antirepression of the p53-Mdm2 complex is also needed which leads to the need of a new drug design paradigm.

Model-driven experimental approach reveals the complex regulatory distribution of p53 by the circadian factor Period 2

The circadian clock and cell cycle networks are interlocked on the molecular level, with the core clock loop exerting a multilevel regulatory role over cell cycle components. This is particularly relevant to the circadian factor Period 2 (Per2), which modulates the stability of the tumor suppressor p53 in unstressed cells and transcriptional activity in response to genotoxic stress. Per2 binding prevents Mdm2-mediated ubiquitination of p53 and, therefore, its degradation, and oscillations in the peaks of Per2 and p53 were expected to correspond. However, our findings showed that Per2 and p53 rhythms were significantly out-of-phase relative to each other in cell lysates and in purified cytoplasmic fractions. These seemingly conflicting experimental data motivated the use of a combined theoretical and experimental approach focusing on the role played by Per2 in dictating the phase of p53 oscillations. Systematic modeling of all possible regulatory scenarios predicted that the observed phase relationship between Per2 and p53 could be simulated if (i) p53 was more stable in the nucleus than in the cytoplasm, (ii) Per2 associates to various ubiquitinated forms of p53, and (iii) Per2 mediated p53 nuclear import. These predictions were supported by a sevenfold increase in p53's half-life in the nucleus and by in vitro binding of Per2 to the various ubiquitinated forms of p53. Last, p53's nuclear shuttling was significantly favored by ectopic expression of Per2 and reduced because of Per2 down-regulation. Our combined theoretical/mathematical approach reveals how clock regulatory nodes can be inferred from oscillating time course data.

Estimation of delays in generalized asynchronous Boolean networks

A new generalized asynchronous Boolean network (GABN) model has been proposed in this paper. This continuous-time discrete-state model captures the biological reality of cellular dynamics without compromising the computational efficiency of the Boolean framework. The GABN synthesis procedure is based on the prior knowledge of the logical structure of the regulatory network, and the experimental transcriptional parameters. The novelty of the proposed methodology lies in considering different delays associated with the activation and deactivation of a particular protein (especially the transcription factors). A few illustrative examples of some well-studied network motifs have been provided to explore the scope of using the GABN model for larger networks. The GABN model of the p53-signaling pathway in response to γ-irradiation has also been simulated in the current paper to provide an indirect validation of the proposed schema.

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.

State feedback control design for Boolean networks

Background

Driving Boolean networks to desired states is of paramount significance toward our ultimate goal of controlling the progression of biological pathways and regulatory networks. Despite recent computational development of controllability of general complex networks and structural controllability of Boolean networks, there is still a lack of bridging the mathematical condition on controllability to real boolean operations in a network. Further, no realtime control strategy has been proposed to drive a Boolean network.

Results

In this study, we applied semi-tensor product to represent boolean functions in a network and explored controllability of a boolean network based on the transition matrix and time transition diagram. We determined the necessary and sufficient condition for a controllable Boolean network and mapped this requirement in transition matrix to real boolean functions and structure property of a network. An efficient tool is offered to assess controllability of an arbitrary Boolean network and to determine all reachable and non-reachable states. We found six simplest forms of controllable 2-node Boolean networks and explored the consistency of transition matrices while extending these six forms to controllable networks with more nodes. Importantly, we proposed the first state feedback control strategy to drive the network based on the status of all nodes in the network. Finally, we applied our reachability condition to the major switch of P53 pathway to predict the progression of the pathway and validate the prediction with published experimental results.

Conclusions

This control strategy allowed us to apply realtime control to drive Boolean networks, which could not be achieved by the current control strategy for Boolean networks. Our results enabled a more comprehensive understanding of the evolution of Boolean networks and might be extended to output feedback control design.

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.

Two is better than one; toward a rational design of combinatorial therapy

Drug combination is an appealing strategy for combating the heterogeneity of tumors and evolution of drug resistance. However, the rationale underlying combinatorial therapy is often not well established due to lack of understandings of the specific pathways responding to the drugs, and their temporal dynamics following each treatment. Here we present several emerging trends in harnessing properties of biological systems for the optimal design of drug combinations, including the type of drugs, specific concentration, sequence of addition and the temporal schedule of treatments. We highlight recent studies showing different approaches for efficient design of drug combinations including single-cell signaling dynamics, adaption and pathway crosstalk. Finally, we discuss novel and feasible approaches that can facilitate the optimal design of combinatorial therapy.

Hyperactivation of ATM upon DNA-PKcs inhibition modulates p53 dynamics and cell fate in response to DNA damage

Time-resolved single-cell analysis and pharmacological perturbations reveal a new regulatory interplay between PI3K-like kinases in response to DNA double-strand breaks: upon loss of DNA-PKcs activity, ATM is hyperactivated and induces an amplified p53 response, which sensitizes cells for damage-induced senescence.

A functional DNA damage response is essential for maintaining genome integrity in the presence of DNA double-strand breaks. It is mainly coordinated by the kinases ATM, ATR, and DNA-PKcs, which control the repair of broken DNA strands and relay the damage signal to the tumor suppressor p53 to induce cell cycle arrest, apoptosis, or senescence. Although many functions of the individual kinases have been identified, it remains unclear how they act in concert to ensure faithful processing of the damage signal. Using specific inhibitors and quantitative analysis at the single-cell level, we systematically characterize the contribution of each kinase for regulating p53 activity. Our results reveal a new regulatory interplay in which loss of DNA-PKcs function leads to hyperactivation of ATM and amplification of the p53 response, sensitizing cells for damage-induced senescence. This interplay determines the outcome of treatment regimens combining irradiation with DNA-PKcs inhibitors in a p53-dependent manner.

p53 Pulses Diversify Target Gene Expression Dynamics in an mRNA Half-Life-Dependent Manner and Delineate Co-regulated Target Gene Subnetworks

The transcription factor p53 responds to DNA double-strand breaks by increasing in concentration in a series of pulses of fixed amplitude, duration, and period. How p53 pulses influence the dynamics of p53 target gene expression is not understood. Here, we show that, in bulk cell populations, patterns of p53 target gene expression cluster into groups with stereotyped temporal behaviors, including pulsing and rising dynamics. These behaviors correlate statistically with the mRNA decay rates of target genes: short mRNA half-lives produce pulses of gene expression. This relationship can be recapitulated by mathematical models of p53-dependent gene expression in single cells and cell populations. Single-cell transcriptional profiling demonstrates that expression of a subset of p53 target genes is coordinated across time within single cells; p53 pulsing attenuates this coordination. These results help delineate how p53 orchestrates the complex DNA damage response and give insight into the function of pulsatile signaling pathways.

A suppressive role of ionizing radiation-responsive miR-29c in the development of liver carcinoma via targeting WIP1

Hepatocellular carcinoma (HCC) is the third most common cause of cancer-related deaths worldwide, and it has been linked to radiation exposure. As a well-defined oncogene, wild-type p53-induced phosphatase 1 (WIP1) plays an inhibitory role in several tumor suppressor pathways, including p53. WIP1 is amplified and overexpressed in many malignancies, including HCC. However, the underlying mechanisms remain largely unknown. Here, we show that low-dose ionizing radiation (IR) induces miR-29c expression in female mouse liver, while inhibiting its expression in HepG2, a human hepatocellular carcinoma cell line which is used as a model system in this study. miR-29c expression is downregulated in human hepatocellular carcinoma cells, which is inversely correlated with WIP1 expression. miR-29c attenuates luciferase activity of a reporter harboring the 3'UTR binding motif of WIP1 mRNA. Ectopic expression of miR-29c significantly represses cell proliferation and induces apoptosis and G1 arrest in HepG2. In contrast, the knockdown of miR-29c greatly enhances HepG2 cell proliferation and suppresses apoptosis. The biological effects of miR-29c may be mediated by its target WIP1 which regulates p53 activity via dephosphorylation at Ser-15. Finally, fluorescence in situ hybridization (FISH) and immunohistochemical analyses indicate that miR-29c is downregulated in 50.6% of liver carcinoma tissues examined, whereas WIP1 is upregulated in 45.4% of these tissues. The expression of miR-29c inversely correlates with that of WIP1 in HCC. Our results suggest that the IR-responsive miR-29c may function as a tumor suppressor that plays a crucial role in the development of liver carcinoma via targeting WIP1, therefore possibly representing a target molecule for therapeutic intervention for this disease.

Schedule-dependent interaction between anticancer treatments

The oncogene MDMX is overexpressed in many cancers, leading to suppression of the tumor suppressor p53. Inhibitors of the oncogene product MDMX therefore might help reactivate p53 and enhance the efficacy of DNA-damaging drugs. However, we currently lack a quantitative understanding of how MDMX inhibition affects the p53 signaling pathway and cell sensitivity to DNA damage. Live cell imaging showed that MDMX depletion triggered two distinct phases of p53 accumulation in single cells: an initial postmitotic pulse, followed by low-amplitude oscillations. The response to DNA damage was sharply different in these two phases; in the first phase, MDMX depletion was synergistic with DNA damage in causing cell death, whereas in the second phase, depletion of MDMX inhibited cell death. Thus a quantitative understanding of signal dynamics and cellular states is important for designing an optimal schedule of dual-drug administration.

Redefining Signaling Pathways with an Expanding Single-Cell Toolbox

Genetically identical cells respond heterogeneously to uniform environmental stimuli. Consequently, investigating the signaling networks that control these cell responses using 'average' bulk cell measurements can obscure underlying mechanisms and misses information emerging from cell-to-cell variability. Here we review recent technological advances including live-cell fluorescence imaging-based approaches and microfluidic devices that enable measurements of signaling networks, dynamics, and responses in single cells. We discuss how these single-cell tools have uncovered novel mechanistic insights for canonical signaling pathways that control cell proliferation (ERK), DNA-damage responses (p53), and innate immune and stress responses (NF-κB). Future improvements in throughput and multiplexing, analytical pipelines, and in vivo applicability will all significantly expand the biological information gained from single-cell measurements of signaling pathways.

Feedbacks, Bifurcations, and Cell Fate Decision-Making in the p53 System

The p53 transcription factor is a regulator of key cellular processes including DNA repair, cell cycle arrest, and apoptosis. In this theoretical study, we investigate how the complex circuitry of the p53 network allows for stochastic yet unambiguous cell fate decision-making. The proposed Markov chain model consists of the regulatory core and two subordinated bistable modules responsible for cell cycle arrest and apoptosis. The regulatory core is controlled by two negative feedback loops (regulated by Mdm2 and Wip1) responsible for oscillations, and two antagonistic positive feedback loops (regulated by phosphatases Wip1 and PTEN) responsible for bistability. By means of bifurcation analysis of the deterministic approximation we capture the recurrent solutions (i.e., steady states and limit cycles) that delineate temporal responses of the stochastic system. Direct switching from the limit-cycle oscillations to the “apoptotic” steady state is enabled by the existence of a subcritical Neimark—Sacker bifurcation in which the limit cycle loses its stability by merging with an unstable invariant torus. Our analysis provides an explanation why cancer cell lines known to have vastly diverse expression levels of Wip1 and PTEN exhibit a broad spectrum of responses to DNA damage: from a fast transition to a high level of p53 killer (a p53 phosphoform which promotes commitment to apoptosis) in cells characterized by high PTEN and low Wip1 levels to long-lasting p53 level oscillations in cells having PTEN promoter methylated (as in, e.g., MCF-7 cell line).

Cancers are diseases of signaling networks. Transcription factor p53 is a pivotal node of a network that integrates a variety of stress signals and governs critical processes of DNA repair, cell cycle arrest, and apoptosis. Somewhat paradoxically, despite the fact that carcinogenesis is prevalently caused by p53 network malfunction, most of our knowledge about p53 signaling is based on cancer or immortalized cell lines. In this paper, we construct a mathematical model of intact p53 network to understand dynamics of non-cancerous cells and then dynamics of cancerous cells by introducing perturbations to the regulatory system. Cell fate decisions are enabled by the presence of interlinked feedback loops which give rise to a rich repertoire of behaviors. We explain and analyze by means of numerical simulations how the dynamical structure of the regulatory system allows for generating unambiguous single-cell fate decisions, also in the case when the cell population splits into an apoptotic and a surviving subpopulation. Perturbation analysis provides an explanation why cancer cell lines known to have vastly diverse expression levels of p53 regulators can exhibit a broad spectrum of responses to DNA damage.

Modeling the heterogeneity of p53 dynamics in DNA damage response

The tumor suppressor p53 can be activated by DNA damage and exhibits undamped pulses. Recent reports have demonstrated a non-threshold mechanism for p53 dynamics. However, no related theoretical studies have been proposed. Here, we constructed a refined DNA damage repair model that incorporated both intrinsic and extrinsic DNA lesions. We proposed that the basal DNA damage may trigger significant fractions of p53 pulses. We also reproduced the heterogeneity of p53 dynamics in experiments. The number of p53 pulses showed no correlations with DNA damage. We also replicated the linear correlation between DNA damage and the probability of igniting a pulse. Our model has unraveled the heterogeneous p53 responses.

Inference of asynchronous Boolean network from biological pathways

Gene regulation is a complex process with multiple levels of interactions. In order to describe this complex dynamical system with tractable parameterization, the choice of the dynamical system model is of paramount importance. The right abstraction of the modeling scheme can reduce the complexity in the inference and intervention design, both computationally and experimentally. This article proposes an asynchronous Boolean network framework to capture the transcriptional regulation as well as the protein-protein interactions in a genetic regulatory system. The inference of asynchronous Boolean network from biological pathways information and experimental evidence are explained using an algorithm. The suitability of this paradigm for the variability of several reaction rates is also discussed. This methodology and model selection open up new research challenges in understanding gene-protein interactive system in a coherent way and can be beneficial for designing effective therapeutic intervention strategy.

MiR-192-Mediated Positive Feedback Loop Controls the Robustness of Stress-Induced p53 Oscillations in Breast Cancer Cells

The p53 tumor suppressor protein plays a critical role in cellular stress and cancer prevention. A number of post-transcriptional regulators, termed microRNAs, are closely connected with the p53-mediated cellular networks. While the molecular interactions among p53 and microRNAs have emerged, a systems-level understanding of the regulatory mechanism and the role of microRNAs-forming feedback loops with the p53 core remains elusive. Here we have identified from literature that there exist three classes of microRNA-mediated feedback loops revolving around p53, all with the nature of positive feedback coincidentally. To explore the relationship between the cellular performance of p53 with the microRNA feedback pathways, we developed a mathematical model of the core p53-MDM2 module coupled with three microRNA-mediated positive feedback loops involving miR-192, miR-34a, and miR-29a. Simulations and bifurcation analysis in relationship to extrinsic noise reproduce the oscillatory behavior of p53 under DNA damage in single cells, and notably show that specific microRNA abrogation can disrupt the wild-type cellular phenotype when the ubiquitous cell-to-cell variability is taken into account. To assess these in silico results we conducted microRNA-perturbation experiments in MCF7 breast cancer cells. Time-lapse microscopy of cell-population behavior in response to DNA double-strand breaks, together with image classification of single-cell phenotypes across a population, confirmed that the cellular p53 oscillations are compromised after miR-192 perturbations, matching well with the model predictions. Our study via modeling in combination with quantitative experiments provides new evidence on the role of microRNA-mediated positive feedback loops in conferring robustness to the system performance of stress-induced response of p53.

DNA damage triggered activities of the tumor suppressor protein p53 could be significantly dynamical. The functional role of p53 oscillations in cellular decision making during cancer development has been appreciated. A set of recent studies have revealed extensive crosstalk between the p53 network and microRNAs, but the specifics of the participation of microRNAs in the regulation of the p53 signaling pathway remains largely elusive. Here we investigated microRNAs that form feedback regulation with p53. We enumerated the molecular interactions among these microRNAs and the p53 core and developed a mathematical model to reproduce the DNA damage induced p53 oscillations in single cells. We performed computer simulations and system analysis in combination with experimental assessment to probe the behavior of p53 under microRNA-inhibited conditions. We show that the robust cellular performance of the stress response of p53 in a breast cancer cell line is controlled by miR-192, which forms positive feedback loops with p53.

Role of the nucleus in apoptosis: signaling and execution

Since their establishment in the early 1970s, the nuclear changes upon apoptosis induction, such as the condensation of chromatin, disassembly of nuclear scaffold proteins and degradation of DNA, were, and still are, considered as the essential steps and hallmarks of apoptosis. These are the characteristics of the execution phase of apoptotic cell death. In addition, accumulating data clearly show that some nuclear events can lead to the induction of apoptosis. In particular, if DNA lesions resulting from deregulation during the cell cycle or DNA damage induced by chemotherapeutic drugs or viral infection cannot be efficiently eliminated, apoptotic mechanisms, which enable cellular transformation to be avoided, are activated in the nucleus. The functional heterogeneity of the nuclear organization allows the tight regulation of these signaling events that involve the movement of various nuclear proteins to other intracellular compartments (and vice versa) to initiate and govern apoptosis. Here, we discuss how these events are coordinated to execute apoptotic cell death.

Bifurcation analysis and potential landscapes of the p53-Mdm2 module regulated by the co-activator programmed cell death 5

The dynamics of p53 play important roles in the regulation of cell fate decisions in response to various stresses, and programmed cell death 5 (PDCD5) functions as a co-activator of p53 that modulates p53 dynamics. In the present paper, we investigated how p53 dynamics are modulated by PDCD5 during the deoxyribose nucleic acid damage response using methods of bifurcation analysis and potential landscape. Our results revealed that p53 activities display rich dynamics under different PDCD5 levels, including monostability, bistability with two stable steady states, oscillations, and the coexistence of a stable steady state (or two states) and an oscillatory state. The physical properties of the p53 oscillations were further demonstrated by the potential landscape in which the potential force attracts the system state to the limit cycle attractor, and the curl flux force drives coherent oscillation along the cyclic trajectory. We also investigated the efficiency with which PDCD5 induced p53 oscillations. We show that Hopf bifurcation can be induced by increasing the PDCD5 efficiency and that the system dynamics exhibited clear transition features in both barrier height and energy dissipation when the efficiency was close to the bifurcation point.

An acetyl-methyl switch drives a conformational change in p53

Individual posttranslational modifications (PTMs) of p53 mediate diverse p53-dependent responses; however, much less is known about the combinatorial action of adjacent modifications. Here, we describe crosstalk between the early DNA damage response mark p53K382me2 and the surrounding PTMs that modulate binding of p53 cofactors, including 53BP1 and p300. The 1.8 Å resolution crystal structure of the tandem Tudor domain (TTD) of 53BP1 in complex with p53 peptide acetylated at K381 and dimethylated at K382 (p53K381acK382me2) reveals that the dual PTM induces a conformational change in p53. The α-helical fold of p53K381acK382me2 positions the side chains of R379, K381ac, and K382me2 to interact with TTD concurrently, reinforcing a modular design of double PTM mimetics. Biochemical and nuclear magnetic resonance analyses show that other surrounding PTMs, including phosphorylation of serine/threonine residues of p53, affect association with TTD. Our findings suggest a novel PTM-driven conformation switch-like mechanism that may regulate p53 interactions with binding partners.

Studying Cellular Signal Transduction with OMIC Technologies

In the gulf between genotype and phenotype exists proteins and, in particular, protein signal transduction systems. These systems use a relatively limited parts list to respond to a much longer list of extracellular, environmental, and/or mechanical cues with rapidity and specificity. Most signaling networks function in a highly non-linear and often contextual manner. Furthermore, these processes occur dynamically across space and time. Because of these complexities, systems and "OMIC" approaches are essential for the study of signal transduction. One challenge in using OMIC-scale approaches to study signaling is that the "signal" can take different forms in different situations. Signals are encoded in diverse ways such as protein-protein interactions, enzyme activities, localizations, or post-translational modifications to proteins. Furthermore, in some cases, signals may be encoded only in the dynamics, duration, or rates of change of these features. Accordingly, systems-level analyses of signaling may need to integrate multiple experimental and/or computational approaches. As the field has progressed, the non-triviality of integrating experimental and computational analyses has become apparent. Successful use of OMIC methods to study signaling will require the "right" experiments and the "right" modeling approaches, and it is critical to consider both in the design phase of the project. In this review, we discuss common OMIC and modeling approaches for studying signaling, emphasizing the philosophical and practical considerations for effectively merging these two types of approaches to maximize the probability of obtaining reliable and novel insights into signaling biology.