Cell-to-Cell Variation in p53 Dynamics Leads to Fractional Killing.

Many chemotherapeutic drugs kill only a fraction of cancer cells, limiting their efficacy. We used live-cell imaging to investigate the role of p53 dynamics in fractional killing of colon cancer cells in response to chemotherapy. We found that both surviving and dying cells reach similar levels of p53, indicating that cell death is not determined by a fixed p53 threshold. Instead, a cell's probability of death depends on the time and levels of p53. Cells must reach a threshold level of p53 to execute apoptosis, and this threshold increases with time. The increase in p53 apoptotic threshold is due to drug-dependent induction of anti-apoptotic genes, predominantly in the inhibitors of apoptosis (IAP) family. Our study underlines the importance of measuring the dynamics of key players in response to chemotherapy to determine mechanisms of resistance and optimize the timing of combination therapy.

The p21CIP1-CDK4-DREAM axis is a master regulator of genotoxic stress-induced cellular senescence.

Cellular senescence, a major driver of aging, can be stimulated by DNA damage, and is counteracted by the DNA repair machinery. Here we show that in p16INK4a-deficient cells, senescence induction by the environmental genotoxin B[a]P or ionizing radiation (IR) completely depends on p21CIP1. Immunoprecipitation-based mass spectrometry interactomics data revealed that during senescence induction and maintenance, p21CIP1 specifically inhibits CDK4 and thereby activates the DREAM complex. Genome-wide transcriptomics revealed striking similarities in the response induced by B[a]P and IR. Among the top 100 repressed genes 78 were identical between B[a]P and IR and 76 were DREAM targets. The DREAM complex transcriptionally silences the main proliferation-associated transcription factors E2F1, FOXM1 and B-Myb as well as multiple DNA repair factors. Knockdown of p21CIP1, E2F4 or E2F5 diminished both, repression of these factors and senescence. The transcriptional profiles evoked by B[a]P and IR largely overlapped with the profile induced by pharmacological CDK4 inhibition, further illustrating the role of CDK4 inhibition in genotoxic stress-induced senescence. Moreover, data obtained by live-cell time-lapse microscopy suggest the inhibition of CDK4 by p21CIP1 is especially important for arresting cells which slip through mitosis. Overall, we identified the p21CIP1/CDK4/DREAM axis as a master regulator of genotoxic stress-induced senescence.

State- and stimulus-specific dynamics of SMAD signaling determine fate decisions in individual cells.

SMAD-mediated signaling regulates apoptosis, cell cycle arrest, and epithelial-to-mesenchymal transition to safeguard tissue homeostasis. However, it remains elusive how the relatively simple pathway can determine such a broad range of cell fate decisions and how it differentiates between varying ligands. Here, we systematically investigate how SMAD-mediated responses are modulated by various ligands of the transforming growth factor ß (TGFß) family and compare these ligand responses in quiescent and proliferating MCF10A cells. We find that the nature of the phenotypic response is mainly determined by the proliferation status, with migration and cell cycle arrest being dominant in proliferating cells for all tested TGFß family ligands, whereas cell death is the major outcome in quiescent cells. In both quiescent and proliferating cells, the identity of the ligand modulates the strength of the phenotypic response proportional to the dynamics of induced SMAD nuclear-to-cytoplasmic translocation and, as a consequence, the corresponding gene expression changes. Interestingly, the proliferation state of a cell has little impact on the set of genes induced by SMAD signaling; instead, it modulates the relative cellular sensitivity to TGFß superfamily members. Taken together, diversity of SMAD-mediated responses is mediated by differing cellular states, which determine ligand sensitivity and phenotypic effects, while the pathway itself merely serves as a quantitative relay from the cell membrane to the nucleus.

Basal dynamics of p53 reveal transcriptionally attenuated pulses in cycling cells.

The tumor suppressor p53 is activated by stress and leads to cellular outcomes such as apoptosis and cell-cycle arrest. Its activation must be highly sensitive to ensure that cells react appropriately to damage. However, proliferating cells often encounter transient damage during normal growth, where cell-cycle arrest or apoptosis may be unfavorable. How does the p53 pathway achieve the right balance between high sensitivity and tolerance to intrinsic damage? Using quantitative time-lapse microscopy of individual human cells, we found that proliferating cells show spontaneous pulses of p53, which are triggered by an excitable mechanism during cell-cycle phases associated with intrinsic DNA damage. However, in the absence of sustained damage, posttranslational modifications keep p53 inactive, preventing it from inducing p21 expression and cell-cycle arrest. Our approach of quantifying basal dynamics in individual cells can now be used to study how other pathways in human cells achieve sensitivity in noisy environments.

Data-based stochastic modeling reveals sources of activity bursts in single-cell TGF-ß signaling.

Cells sense their surrounding by employing intracellular signaling pathways that transmit hormonal signals from the cell membrane to the nucleus. TGF-ß/SMAD signaling encodes various cell fates, controls tissue homeostasis and is deregulated in diseases such as cancer. The pathway shows strong heterogeneity at the single-cell level, but quantitative insights into mechanisms underlying fluctuations at various time scales are still missing, partly due to inefficiency in the calibration of stochastic models that mechanistically describe signaling processes. In this work we analyze single-cell TGF-ß/SMAD signaling and show that it exhibits temporal stochastic bursts which are dose-dependent and whose number and magnitude correlate with cell migration. We propose a stochastic modeling approach to mechanistically describe these pathway fluctuations with high computational efficiency. Employing high-order numerical integration and fitting to burst statistics we enable efficient quantitative parameter estimation and discriminate models that assume noise in different reactions at the receptor level. This modeling approach suggests that stochasticity in the internalization of TGF-ß receptors into endosomes plays a key role in the observed temporal bursting. Further, the model predicts the single-cell dynamics of TGF-ß/SMAD signaling in untested conditions, e.g., successfully reflects memory effects of signaling noise and cellular sensitivity towards repeated stimulation. Taken together, our computational framework based on burst analysis, noise modeling and path computation scheme is a suitable tool for the data-based modeling of complex signaling pathways, capable of identifying the source of temporal noise.

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

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

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

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

Quantitative live cell imaging reveals a gradual shift between DNA repair mechanisms and a maximal use of HR in mid S phase.

DNA double-strand breaks are repaired by two main pathways: nonhomologous end joining (NHEJ) and homologous recombination (HR). The choice between these pathways depends on cell-cycle phase; however the continuous effect of cell cycle on the balance between them is still unclear. We used live cell imaging and fluorescent reporters for 53BP1, Rad52, and cell cycle to quantify the relative contribution of NHEJ and HR at different points of the cell cycle in single cells. We found that NHEJ is the dominant repair pathway in G1 and G2 even when both repair pathways are functional. The shift from NHEJ to HR is gradual, with the highest proportion of breaks repaired by HR in mid S, where the amount of DNA replication is highest. Higher proportions of HR also strongly correlate with slower rates of repair. Our study shows that the choice of repair mechanism is continuously adjusted throughout the cell cycle and suggests that the extent of active replication, rather than the presence of a sister chromatid influences the balance between the two repair pathways in human cells.

Comparing Agent-Based Delivery of DNA and PNA Forced Intercalation (FIT) Probes for Multicolor mRNA Imaging.

Fluorogenic oligonucleotide probes allow mRNA imaging in living cells. A key challenge is the cellular delivery of probes. Most delivery agents, such as cell-penetrating peptides (CPPs) and pore-forming proteins, require interactions with the membrane. Charges play an important role. To explore the influence of charge on fluorogenic properties and delivery efficiency, we compared peptide nucleic acid (PNA)- with DNA-based forced intercalation (FIT) probes. Perhaps counterintuitively, fluorescence signaling by charged DNA FIT probes proved tolerant to CPP conjugation, whereas CPP-FIT PNA conjugates were affected. Live-cell imaging was performed with a genetically engineered HEK293 cell line to allow the inducible expression of a specific mRNA target. Blob-like features and high background were recurring nuisances of the tested CPP and lipid conjugates. By contrast, delivery by streptolysin-O provided high enhancements of the fluorescence of the FIT probe upon target induction. Notably, DNA-based FIT probes were brighter and more responsive than PNA-based FIT probes. Optimized conditions enabled live-cell multicolor imaging of three different mRNA target sequences.

Cell-specific responses to the cytokine TGFß are determined by variability in protein levels.

The cytokine TGFß provides important information during embryonic development, adult tissue homeostasis, and regeneration. Alterations in the cellular response to TGFß are involved in severe human diseases. To understand how cells encode the extracellular input and transmit its information to elicit appropriate responses, we acquired quantitative time-resolved measurements of pathway activation at the single-cell level. We established dynamic time warping to quantitatively compare signaling dynamics of thousands of individual cells and described heterogeneous single-cell responses by mathematical modeling. Our combined experimental and theoretical study revealed that the response to a given dose of TGFß is determined cell specifically by the levels of defined signaling proteins. This heterogeneity in signaling protein expression leads to decomposition of cells into classes with qualitatively distinct signaling dynamics and phenotypic outcome. Negative feedback regulators promote heterogeneous signaling, as a SMAD7 knock-out specifically affected the signal duration in a subpopulation of cells. Taken together, we propose a quantitative framework that allows predicting and testing sources of cellular signaling heterogeneity.

Circular RNAs are a large class of animal RNAs with regulatory potency.

Circular RNAs (circRNAs) in animals are an enigmatic class of RNA with unknown function. To explore circRNAs systematically, we sequenced and computationally analysed human, mouse and nematode RNA. We detected thousands of well-expressed, stable circRNAs, often showing tissue/developmental-stage-specific expression. Sequence analysis indicated important regulatory functions for circRNAs. We found that a human circRNA, antisense to the cerebellar degeneration-related protein 1 transcript (CDR1as), is densely bound by microRNA (miRNA) effector complexes and harbours 63 conserved binding sites for the ancient miRNA miR-7. Further analyses indicated that CDR1as functions to bind miR-7 in neuronal tissues. Human CDR1as expression in zebrafish impaired midbrain development, similar to knocking down miR-7, suggesting that CDR1as is a miRNA antagonist with a miRNA-binding capacity ten times higher than any other known transcript. Together, our data provide evidence that circRNAs form a large class of post-transcriptional regulators. Numerous circRNAs form by head-to-tail splicing of exons, suggesting previously unrecognized regulatory potential of coding sequences.

Disorder and residual helicity alter p53-Mdm2 binding affinity and signaling in cells.

Levels of residual structure in disordered interaction domains determine in vitro binding affinities, but whether they exert similar roles in cells is not known. Here, we show that increasing residual p53 helicity results in stronger Mdm2 binding, altered p53 dynamics, impaired target gene expression and failure to induce cell cycle arrest upon DNA damage. These results establish that residual structure is an important determinant of signaling fidelity in cells.

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

DNA damage initiates a series of p53 pulses. Although much is known about the interactions surrounding p53, little is known about which interactions contribute to p53's dynamical behavior. The simplest explanation is that these pulses are oscillations intrinsic to the p53/Mdm2 negative feedback loop. Here we present evidence that this simple mechanism is insufficient to explain p53 pulses; we show that p53 pulses are externally driven by pulses in the upstream signaling kinases, ATM and Chk2, and that the negative feedback between p53 and ATM, via Wip1, is essential for maintaining the uniform shape of p53 pulses. We propose that p53 pulses result from repeated initiation by ATM, which is reactivated by persistent DNA damage. Our study emphasizes the importance of collecting quantitative dynamic information at high temporal resolution for understanding the regulation of signaling pathways and opens new ways to manipulate p53 pulses to ask questions about their function in response to DNA damage.

The p53 response in single cells is linearly correlated to the number of DNA breaks without a distinct threshold.

BACKGROUND: The tumor suppressor protein p53 is activated by cellular stress. DNA double strand breaks (DSBs) induce the activation of the kinase ATM, which stabilizes p53 and activates its transcriptional activity. Single cell analysis revealed that DSBs induced by gamma irradiation trigger p53 accumulation in a series of pulses that vary in number from cell to cell. Higher levels of irradiation increase the number of p53 pulses suggesting that they arise from periodic examination of the damage by ATM. If damage persists, additional pulses of p53 are triggered. The threshold of damage required for activating a p53 pulse is unclear. Previous studies that averaged the response across cell populations suggested that one or two DNA breaks are sufficient for activating ATM and p53. However, it is possible that by averaging over a population of cells important features of the dependency between DNA breaks and p53 dynamics are missed. RESULTS: Using fluorescent reporters we developed a system for following in individual cells the number of DSBs, the kinetics of repair and the p53 response. We found a large variation in the initial number of DSBs and the rate of repair between individual cells. Cells with higher number of DSBs had higher probability of showing a p53 pulse. However, there was no distinct threshold number of breaks for inducing a p53 pulse. We present evidence that the decision to activate p53 given a specific number of breaks is not entirely stochastic, but instead is influenced by both cell-intrinsic factors and previous exposure to DNA damage. We also show that the natural variations in the initial amount of p53, rate of DSB repair and cell cycle phase do not affect the probability of activating p53 in response to DNA damage. CONCLUSIONS: The use of fluorescent reporters to quantify DNA damage and p53 levels in live cells provided a quantitative analysis of the complex interrelationships between both processes. Our study shows that p53 activation differs even between cells that have a similar number of DNA breaks. Understanding the origin and consequences of such variability in normal and cancerous cells is crucial for developing efficient and selective therapeutic interventions.

Delivery to, and Reactivation of, the p53 Pathway in Cancer Cells Using a Grafted Cyclotide Conjugated with a Cell-Penetrating Peptide.

Peptides are promising drug modalities that can modulate protein-protein interactions, but their application is hampered by their limited ability to reach intracellular targets. Here, we improved the cytosolic delivery of a peptide blocking p53:MDM2/X interactions using a cyclotide as a stabilizing scaffold. We applied several design strategies to improve intracellular delivery and found that the conjugation of the lead cyclotide to the cyclic cell-penetrating peptide cR10 was the most effective. Conjugation allowed cell internalization at micromolar concentration and led to elevated intracellular p53 levels in A549, MCF7, and MCF10A cells, as well as inducing apoptosis in A549 cells without causing membrane disruption. The lead peptide had >35-fold improvement in inhibitory activity and increased cellular uptake compared to a previously reported cyclotide p53 activator. In summary, we demonstrated the delivery of a large polar cyclic peptide in the cytosol and confirmed its ability to modulate intracellular protein-protein interactions involved in cancer.

A synthetic-natural hybrid oscillator in human cells.

Recent studies have shown that many cell-signaling networks contain interactions and feedback loops that give rise to complex dynamics. Synthetic biology has allowed researchers to construct and analyze well-defined signaling circuits exhibiting behavior that can be predicted and quantitatively understood. Combining these approaches--wiring natural network components together with engineered interactions--has the potential to precisely modulate the dynamics of endogenous signaling processes and control the cell decisions they influence. Here, we focus on the p53 signaling pathway as a template for constructing a tunable oscillator comprised of both natural and synthetic components in mammalian cells. We find that a reduced p53 circuit implementing a single feedback loop preserves some features of the full network's dynamics, exhibiting pulses of p53 with tightly controlled timing. However, in contrast to the full natural p53 network, these pulses are damped in individual cells, with amplitude that depends on the input strength. Guided by a computational model of the reduced circuit, we constructed and analyzed circuit variants supplemented with synthetic positive and negative feedback loops and subjected to chemical perturbation. Our work demonstrates that three important features of oscillator dynamics--amplitude, period, and the rate of damping--can be controlled by manipulating stimulus level, interaction strength, and feedback topology. The approaches taken here may be useful for the rational design of synthetic networks with defined dynamics, and for identifying perturbations that control dynamics in natural biological circuits for research or therapeutic purposes.

Virtual stress plays tricks on cells.

Optogenetics enables the induction of virtual stress, which separates stress signaling from cellular damage. This provides new insights into the dynamics of the integrated stress response and reveals the mechanisms through which cells form memories of past stress events to guide their response to acute stress.

Resolving Crosstalk Between Signaling Pathways Using Mathematical Modeling and Time-Resolved Single Cell Data.

Crosstalk between signaling pathways can modulate the cellular response to stimuli and is therefore an important part of signal transduction. For a comprehensive understanding of cellular responses, identifying points of interaction between the underlying molecular networks is essential. Here, we present an approach that allows the systematic prediction of such interactions by perturbing one pathway and quantifying the concomitant alterations in the response of a second pathway. As the observed alterations contain information about the crosstalk, we use an ordinary differential equation-based model to extract this information by linking altered dynamics to individual processes. Consequently, we can predict the interaction points between two pathways. As an example, we employed our approach to investigate the crosstalk between the NF-κB and p53 signaling pathway. We monitored the response of p53 to genotoxic stress using time-resolved single cell data and perturbed NF-κB signaling by inhibiting the kinase IKK2. Employing a subpopulation-based modeling approach enabled us to identify multiple interaction points that are simultaneously affected by perturbation of NF-κB signaling. Hence, our approach can be used to analyze crosstalk between two signaling pathways in a systematic manner.

Stimulus-dependent dynamics of p53 in single cells.

Many biological networks respond to various inputs through a common signaling molecule that triggers distinct cellular outcomes. One potential mechanism for achieving specific input-output relationships is to trigger distinct dynamical patterns in response to different stimuli. Here we focused on the dynamics of p53, a tumor suppressor activated in response to cellular stress. We quantified the dynamics of p53 in individual cells in response to UV and observed a single pulse that increases in amplitude and duration in proportion to the UV dose. This graded response contrasts with the previously described series of fixed pulses in response to γ-radiation. We further found that while γ-triggered p53 pulses are excitable, the p53 response to UV is not excitable and depends on continuous signaling from the input-sensing kinases. Using mathematical modeling and experiments, we identified feedback loops that contribute to specific features of the stimulus-dependent dynamics of p53, including excitability and input-duration dependency. Our study shows that different stresses elicit different temporal profiles of p53, suggesting that modulation of p53 dynamics might be used to achieve specificity in this network.

Distinct mechanisms act in concert to mediate cell cycle arrest.

In response to DNA damage, cells arrest at specific stages in the cell cycle. This arrest must fulfill at least 3 requirements: it must be activated promptly; it must be sustained as long as damage is present to prevent loss of genomic information; and after the arrest, cells must re-enter into the appropriate cell cycle phase to ensure proper ploidy. Multiple molecular mechanisms capable of arresting the cell cycle have been identified in mammalian cells; however, it is unknown whether each mechanism meets all 3 requirements or whether they act together to confer specific functions to the arrest. To address this question, we integrated mathematical models describing the cell cycle and the DNA damage signaling networks and tested the contributions of each mechanism to cell cycle arrest and re-entry. Predictions from this model were then tested with quantitative experiments to identify the combined action of arrest mechanisms in irradiated cells. We find that different arrest mechanisms serve indispensable roles in the proper cellular response to DNA damage over time: p53-independent cyclin inactivation confers immediate arrest, whereas p53-dependent cyclin downregulation allows this arrest to be sustained. Additionally, p21-mediated inhibition of cyclin-dependent kinase activity is indispensable for preventing improper cell cycle re-entry and endoreduplication. This work shows that in a complex signaling network, seemingly redundant mechanisms, acting in a concerted fashion, can achieve a specific cellular outcome.