Mobilizing phospholipids on tumor plasma membrane implicates phosphatidylserine externalization blockade for cancer immunotherapy.

In "healthy" tumor cells, phosphatidylserine (PS) is predominately localized in the inner plasma membrane leaflet. During apoptosis, PS relocates to the outer leaflet. Herein, we established PSout tumor models with tumor cells lacking PS flippase component CDC50A, constantly exposing PS but alive. PSout tumors developed bigger than wild-type (WT) tumors, featuring M2 polarized tumor-associated macrophages (TAMs) and fewer tumor-antigen-specific T cells. The PS receptor TIM-3 is responsible for PS recognition. Employing an opposite tumor model, PSin, with tumor cells lacking the PS scramblase Xkr8 and unable to expose PS during otherwise normal apoptosis, we find that the accumulated apoptotic tumor cells produce and release cyclic GAMP (cGAMP) to immune cells to activate the STING pathway, leading to TAM M1 polarization, suppressed interleukin (IL)-10 secretion, and natural killer (NK) cell cytotoxicity. Silencing Xkr8 in vivo by either short hairpin RNA (shRNA) or small interfering RNA (siRNA) to achieve a PS externalization blockade provides robust therapeutic anti-tumor efficiency.

Telomeric 8-oxo-guanine drives rapid premature senescence in the absence of telomere shortening.

Oxidative stress is a primary cause of cellular senescence and contributes to the etiology of numerous human diseases. Oxidative damage to telomeric DNA has been proposed to cause premature senescence by accelerating telomere shortening. Here, we tested this model directly using a precision chemoptogenetic tool to produce the common lesion 8-oxo-guanine (8oxoG) exclusively at telomeres in human fibroblasts and epithelial cells. A single induction of telomeric 8oxoG is sufficient to trigger multiple hallmarks of p53-dependent senescence. Telomeric 8oxoG activates ATM and ATR signaling, and enriches for markers of telomere dysfunction in replicating, but not quiescent cells. Acute 8oxoG production fails to shorten telomeres, but rather generates fragile sites and mitotic DNA synthesis at telomeres, indicative of impaired replication. Based on our results, we propose that oxidative stress promotes rapid senescence by producing oxidative base lesions that drive replication-dependent telomere fragility and dysfunction in the absence of shortening and shelterin loss.

p53 Promotes Cytokine Expression in Melanoma to Regulate Drug Resistance and Migration.

The transcription factor p53 is frequently lost during tumor development in solid tumors; however, most melanomas retain a wild type p53 protein. The presence of wild type p53 in melanoma has fueled speculation that p53 may play a neutral or pro-tumorigenic role in this disease. Here we show that p53 is functional in human melanoma cell lines, and that loss of p53 results in a general reduction in basal NF-kB regulated cytokine production. The reduced cytokine expression triggered by p53 loss is broad and includes key inflammatory chemokines, such as CXCL1, CXCL8, and the IL6 class cytokine LIF, resulting in a reduced ability to induce chemotactic-dependent migration of tumor cells and immune cells and increased sensitivity to BRAF inhibition. Taken together, this result indicates that wild type p53 regulates cytokine expression and induces cytokine-dependent phenotype on melanoma.

Sensing plasma membrane pore formation induces chemokine production in survivors of regulated necrosis.

Although overwhelming plasma membrane integrity loss leads to cell lysis and necrosis, cells can tolerate a limited level of plasma membrane damage, undergo ESCRT-III-mediated repair, and survive. Here, we find that cells which undergo limited plasma membrane damage from the pore-forming actions of MLKL, GSDMD, perforin, or detergents experience local activation of PKCs through Ca2+ influx at the damage sites. S660-phosphorylated PKCs subsequently activate the TAK1/IKKs axis and RelA/Cux1 complex to trigger chemokine expressions. We observe that in late-stage cancers, cells with active MLKL show expression of CXCL8. Similar expression induction is also found in ischemia-injured kidneys. Chemokines generated in this manner are also indispensable for recruiting immune cells to the dead and dying cells. This plasma membrane integrity-sensing pathway is similar to the well-established yeast cell wall integrity signaling pathway at molecular level, and this suggests an evolutionary conserved mechanism to respond to the cellular barrier damage.

p53 dynamics vary between tissues and are linked with radiation sensitivity.

Radiation sensitivity varies greatly between tissues. The transcription factor p53 mediates the response to radiation; however, the abundance of p53 protein does not correlate well with the extent of radiosensitivity across tissues. Given recent studies showing that the temporal dynamics of p53 influence the fate of cultured cells in response to irradiation, we set out to determine the dynamic behavior of p53 and its impact on radiation sensitivity in vivo. We find that radiosensitive tissues show prolonged p53 signaling after radiation, while more resistant tissues show transient p53 activation. Sustaining p53 using a small molecule (NMI801) that inhibits Mdm2, a negative regulator of p53, reduced viability in cell culture and suppressed tumor growth. Our work proposes a mechanism for the control of radiation sensitivity and suggests tools to alter the dynamics of p53 to enhance tumor clearance. Similar approaches can be used to enhance killing of cancer cells or reduce toxicity in normal tissues following genotoxic therapies.

Acetylation-mediated remodeling of the nucleolus regulates cellular acetyl-CoA responses.

The metabolite acetyl-coenzyme A (acetyl-CoA) serves as an essential element for a wide range of cellular functions including adenosine triphosphate (ATP) production, lipid synthesis, and protein acetylation. Intracellular acetyl-CoA concentrations are associated with nutrient availability, but the mechanisms by which a cell responds to fluctuations in acetyl-CoA levels remain elusive. Here, we generate a cell system to selectively manipulate the nucleo-cytoplasmic levels of acetyl-CoA using clustered regularly interspaced short palindromic repeat (CRISPR)-mediated gene editing and acetate supplementation of the culture media. Using this system and quantitative omics analyses, we demonstrate that acetyl-CoA depletion alters the integrity of the nucleolus, impairing ribosomal RNA synthesis and evoking the ribosomal protein-dependent activation of p53. This nucleolar remodeling appears to be mediated through the class IIa histone deacetylases (HDACs). Our findings highlight acetylation-mediated control of the nucleolus as an important hub linking acetyl-CoA fluctuations to cellular stress responses.

Optogenetic Control Reveals Differential Promoter Interpretation of Transcription Factor Nuclear Translocation Dynamics.

Gene expression is thought to be affected not only by the concentration of transcription factors (TFs) but also the dynamics of their nuclear translocation. Testing this hypothesis requires direct control of TF dynamics. Here, we engineer CLASP, an optogenetic tool for rapid and tunable translocation of a TF of interest. Using CLASP fused to Crz1, we observe that, for the same integrated concentration of nuclear TF over time, changing input dynamics changes target gene expression: pulsatile inputs yield higher expression than continuous inputs, or vice versa, depending on the target gene. Computational modeling reveals that a dose-response saturating at low TF input can yield higher gene expression for pulsatile versus continuous input, and that multi-state promoter activation can yield the opposite behavior. Our integrated tool development and modeling approach characterize promoter responses to Crz1 nuclear translocation dynamics, extracting quantitative features that may help explain the differential expression of target genes.

Quantifying the Central Dogma in the p53 Pathway in Live Single Cells.

Transcription factors (TFs) integrate signals to regulate target gene expression, but we generally lack a quantitative understanding of how changes in TF levels regulate mRNA and protein production. Here, we established a system to simultaneously monitor the levels of p53, a TF that shows oscillations following DNA damage, and the transcription and protein levels of its target p21 in individual cells. p21 transcription tracked p53 dynamics, while p21 protein steadily accumulated. p21 transcriptional activation showed bursts of mRNA production, with p53 levels regulating the probability but not magnitude of activation. Variations in p53 levels between cells contributed to heterogeneous p21 transcription while independent p21 alleles exhibited highly correlated behaviors. Pharmacologically elevating p53 increased the probability of p21 transcription with minor effects on its magnitude, leading to a strong increase in p21 protein levels. Our results reveal quantitative mechanisms by which TF dynamics can regulate protein levels of its targets. A record of this paper's transparent peer review process is included in the Supplemental Information.

Identification of universal and cell-type specific p53 DNA binding.

BACKGROUND: The tumor suppressor p53 is a major regulator of the DNA damage response and has been suggested to selectively bind and activate cell-type specific gene expression programs. However recent studies and meta-analyses of genomic data propose largely uniform, and condition independent p53 binding and thus question the selective and cell-type dependent function of p53. RESULTS: To systematically assess the cell-type specificity of p53, we measured its association with DNA in 12 p53 wild-type cancer cell lines, from a range of epithelial linages, in response to ionizing radiation. We found that the majority of bound sites were occupied across all cell lines, however we also identified a subset of binding sites that were specific to one or a few cell lines. Unlike the shared p53-bound genome, which was not dependent on chromatin accessibility, the association of p53 with these atypical binding sites was well explained by chromatin accessibility and could be modulated by forcing cell state changes such as the epithelial-to-mesenchymal transition. CONCLUSIONS: Our study reconciles previous conflicting views in the p53 field, by demonstrating that although the majority of p53 DNA binding is conserved across cell types, there is a small set of cell line specific binding sites that depend on cell state.

The effects of proliferation status and cell cycle phase on the responses of single cells to chemotherapy.

DNA-damaging chemotherapeutics are widely used in cancer treatments, but for solid tumors they often leave a residual tumor-cell population. Here we investigated how cellular states might affect the response of individual cells in a clonal population to cisplatin, a DNA-damaging chemotherapeutic agent. Using a live-cell reporter of cell cycle phase and long-term imaging, we monitored single-cell proliferation before, at the time of, and after treatment. We found that in response to cisplatin, cells either arrested or died, and the ratio of these outcomes depended on the dose. While we found that the cell cycle phase at the time of cisplatin addition was not predictive of outcome, the proliferative history of the cell was: highly proliferative cells were more likely to arrest than to die, whereas slowly proliferating cells showed a higher probability of death. Information theory analysis revealed that the dose of cisplatin had the greatest influence on the cells' decisions to arrest or die, and that the proliferation status interacted with the cisplatin dose to further guide this decision. These results show an unexpected effect of proliferation status in regulating responses to cisplatin and suggest that slowly proliferating cells within tumors may be acutely vulnerable to chemotherapy.

RNA-seq reveals novel mechanistic targets of withaferin A in prostate cancer cells.

Withaferin A (WA) is a promising phytochemical exhibiting in vitro and in vivo anticancer activities against prostate and other cancers, but the mechanism of its action is not fully understood. In this study, we performed RNA-seq analysis using 22Rv1 human prostate cancer cell line to identify mechanistic targets of WA. Kyoto Encyclopedia of Genes and Genomes pathway analysis of the differentially expressed genes showed most significant enrichment of genes associated with metabolism. These results were validated using LNCaP and 22Rv1 human prostate cancer cells and Hi-Myc transgenic mice as models. The intracellular levels of acetyl-CoA, total free fatty acids and neutral lipids were decreased significantly following WA treatment in both cells, which was accompanied by downregulation of mRNA (confirmed by quantitative reverse transcription-polymerase chain reaction) and protein levels of key fatty acid synthesis enzymes, including ATP citrate lyase, acetyl-CoA carboxylase 1, fatty acid synthase and carnitine palmitoyltransferase 1A. Ectopic expression of c-Myc, but not constitutively active Akt, conferred a marked protection against WA-mediated suppression of acetyl-CoA carboxylase 1 and fatty acid synthase protein expression, and clonogenic cell survival. WA was a superior inhibitor of cell proliferation and fatty acid synthesis in comparison with known modulators of fatty acid metabolism including cerulenin and etomoxir. Intraperitoneal WA administration to Hi-Myc transgenic mice (0.1 mg/mouse, three times/week for 5 weeks) also resulted in a significant decrease in circulating levels of total free fatty acids and phospholipids, and expression of ATP citrate lyase, acetyl-CoA carboxylase 1, fatty acid synthase and carnitine palmitoyltransferase 1A proteins in the prostate in vivo.

Fluctuations in p53 Signaling Allow Escape from Cell-Cycle Arrest.

Biological signals need to be robust and filter small fluctuations yet maintain sensitivity to signals across a wide range of magnitudes. Here, we studied how fluctuations in DNA damage signaling relate to maintenance of long-term cell-cycle arrest. Using live-cell imaging, we quantified division profiles of individual human cells in the course of 1 week after irradiation. We found a subset of cells that initially establish cell-cycle arrest and then sporadically escape and divide. Using fluorescent reporters and mathematical modeling, we determined that fluctuations in the oscillatory pattern of the tumor suppressor p53 trigger a sharp switch between p21 and CDK2, leading to escape from arrest. Transient perturbation of p53 stability mimicked the noise in individual cells and was sufficient to trigger escape from arrest. Our results show that the self-reinforcing circuitry that mediates cell-cycle transitions can translate small fluctuations in p53 signaling into large phenotypic changes.

Conservation and Divergence of p53 Oscillation Dynamics across Species.

The tumor-suppressing transcription factor p53 is highly conserved at the protein level and plays a key role in the DNA damage response. One important aspect of p53 regulation is its dynamics in response to DNA damage, which include oscillations. Here, we observe that, while the qualitative oscillatory nature of p53 dynamics is conserved across cell lines derived from human, monkey, dog, mouse, and rat, the oscillation period is variable. Specifically, rodent cells exhibit rapid p53 oscillations, whereas dog, monkey, and human cells show slower oscillations. Computational modeling and experiments identify stronger negative feedback between p53 and MDM2 as the driver of faster oscillations in rodents, suggesting that the period of oscillation is a network-level property. In total, our study shows that despite highly conserved signaling, the quantitative features of p53 oscillations can diverge across evolution. We caution that strong amino acid conservation of proteins and transcriptional network similarity do not necessarily imply conservation of time dynamics.

Integrating genomic information and signaling dynamics for efficient cancer therapy.

The field of cancer systems biology has made great strides in understanding oncogenic pathway signaling and enumerating mutations involved in oncogenesis. However, application of these datasets to patient stratification, and to the design of personalized therapy, is in its infancy. We review BRAF and BRCA mutant targeted therapy, where patient stratification has had critical, albeit mixed success. We contrast the work on genomic targeted therapy with orthogonal studies on the dynamics of signaling pathways for designing optimal treatment schedules. We suggest that an integrated approach, combining genomic data and the dynamics of signaling pathways, is required for developing pathway specific computational models and for systematic deployment of targeted combination regimes. The field of cancer systems biology has made great strides in understanding oncogenic pathway signaling and enumerating mutations involved in oncogenesis. However, application of the existing approaches and datasets to patient stratification, and to the design of personalized therapy, is in its infancy. Here we discuss an integrated approach combining genomic data and the dynamics of signaling pathway to develop pathway specific computational models and the systematic deployment of targeted combination regimes. We review recent studies and existing datasets in the field of cancer systems biology and highlight potentially fruitful synergies between the different strands of this discipline.

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.

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.

Model-guided optogenetic study of PKA signaling in budding yeast.

In eukaryotes, protein kinase A (PKA) is a master regulator of cell proliferation and survival. The activity of PKA is subject to elaborate control and exhibits complex time dynamics. To probe the quantitative attributes of PKA dynamics in the yeast Saccharomyces cerevisiae, we developed an optogenetic strategy that uses a photoactivatable adenylate cyclase to achieve real-time regulation of cAMP and the PKA pathway. We capitalize on the precise and rapid control afforded by this optogenetic tool, together with quantitative computational modeling, to study the properties of feedback in the PKA signaling network and dissect the nonintuitive dynamic effects that ensue from perturbing its components. Our analyses reveal that negative feedback channeled through the Ras1/2 GTPase is delayed, pinpointing its time scale and its contribution to the dynamic features of the cAMP/PKA signaling network.

Dynamics of CDKN1A in Single Cells Defined by an Endogenous Fluorescent Tagging Toolkit.

Observing the endogenous abundance, localization, and dynamics of proteins in mammalian cells is crucial to understanding their function and behavior. Currently, there is no systematic approach for the fluorescent tagging of endogenous loci. Here, we used Cas9-catalyzed DNA breaks, short homology arms, and a family of donor plasmids to establish endogenous Fluorescent tagging (eFlut): a low-cost and efficient approach to generating endogenous proteins with fluorescent labels. We validated this protocol on multiple proteins in several cell lines and species and applied our tools to study the cell-cycle inhibitor CDKN1A in single cells. We uncover heterogeneity in the timing and rate of CDKN1A induction post-DNA damage and show that this variability is post-transcriptionally regulated, depends on cell-cycle position, and has long-term consequences for cellular proliferation. The tools developed in this study should support widespread study of the dynamics and localization of diverse proteins in mammalian cells.

Single-cell analysis of circadian dynamics in tissue explants.

Tracking molecular dynamics in single cells in vivo is instrumental to understanding how cells act and interact in tissues. Current tissue imaging approaches focus on short-term observation and typically nonendogenous or implanted samples. Here we develop an experimental and computational setup that allows for single-cell tracking of a transcriptional reporter over a period of >1 wk in the context of an intact tissue. We focus on the peripheral circadian clock as a model system and measure the circadian signaling of hundreds of cells from two tissues. The circadian clock is an autonomous oscillator whose behavior is well described in isolated cells, but in situ analysis of circadian signaling in single cells of peripheral tissues is as-yet uncharacterized. Our approach allowed us to investigate the oscillatory properties of individual clocks, determine how these properties are maintained among different cells, and assess how they compare to the population rhythm. These experiments, using a wide-field microscope, a previously generated reporter mouse, and custom software to track cells over days, suggest how many signaling pathways might be quantitatively characterized in explant models.