Probing the functional impact of sequence variation on p53-DNA interactions using a novel microsphere assay for protein-DNA binding with human cell extracts

ΔNp63 is a pioneer factor that binds inaccessible chromatin and elicits chromatin remodeling

BACKGROUND

ΔNp63 is a master transcriptional regulator playing critical roles in epidermal development and other cellular processes. Recent studies suggest that ΔNp63 functions as a pioneer factor that can target its binding sites within inaccessible chromatin and induce chromatin remodeling.

METHODS

In order to examine if ΔNp63 can bind to inaccessible chromatin and to determine if specific histone modifications are required for binding, we induced ΔNp63 expression in two p63-naïve cell lines. ΔNp63 binding was then examined by ChIP-seq and the chromatin at ΔNp63 targets sites was examined before and after binding. Further analysis with competitive nucleosome binding assays was used to determine how ΔNp63 directly interacts with nucleosomes.

RESULTS

Our results show that before ΔNp63 binding, targeted sites lack histone modifications, indicating ΔNp63's capability to bind at unmodified chromatin. Moreover, the majority of the sites that are bound by ectopic ΔNp63 expression exist in an inaccessible state. Once bound, ΔNp63 induces acetylation of the histone and the repositioning of nucleosomes at its binding sites. Further analysis with competitive nucleosome binding assays reveal that ΔNp63 can bind directly to nucleosome edges with significant binding inhibition occurring within 50 bp of the nucleosome dyad.

CONCLUSION

Overall, our results demonstrate that ΔNp63 is a pioneer factor that binds nucleosome edges at inaccessible and unmodified chromatin sites and induces histone acetylation and nucleosome repositioning.

Revealing a human p53 universe

p53 transcriptional networks are well-characterized in many organisms. However, a global understanding of requirements for in vivo p53 interactions with DNA and relationships with transcription across human biological systems in response to various p53 activating situations remains limited. Using a common analysis pipeline, we analyzed 41 data sets from genome-wide ChIP-seq studies of which 16 have associated gene expression data, including our recent primary data with normal human lymphocytes. The resulting extensive analysis, accessible at p53 BAER hub via the UCSC browser, provides a robust platform to characterize p53 binding throughout the human genome including direct influence on gene expression and underlying mechanisms. We establish the impact of spacers and mismatches from consensus on p53 binding in vivo and propose that once bound, neither significantly influences the likelihood of expression. Our rigorous approach revealed a large p53 genome-wide cistrome composed of >900 genes directly targeted by p53. Importantly, we identify a core cistrome signature composed of genes appearing in over half the data sets, and we identify signatures that are treatment- or cell-specific, demonstrating new functions for p53 in cell biology. Our analysis reveals a broad homeostatic role for human p53 that is relevant to both basic and translational studies.

Defining TP53 pioneering capabilities with competitive nucleosome binding assays

Accurate gene expression requires the targeting of transcription factors (TFs) to regulatory sequences often occluded within nucleosomes. The ability to target a TF binding site (TFBS) within a nucleosome has been the defining characteristic for a special class of TFs known as pioneer factors. Recent studies suggest TP53 functions as a pioneer factor that can target its TFBS within nucleosomes, but it remains unclear how TP53 binds to nucleosomal DNA. To comprehensively examine TP53 nucleosome binding, we competitively bound TP53 to multiple in vitro–formed nucleosomes containing a high- or low-affinity TP53 TFBS located at differing translational and rotational positions within the nucleosome. Stable TP53–nucleosome complexes were isolated and quantified using next-generation sequencing. Our results demonstrate TP53 binding is limited to nucleosome edges with significant binding inhibition occurring within 50 bp of the nucleosome dyad. Binding site affinity only affects TP53 binding for TFBSs located at the same nucleosomal positions; otherwise, nucleosome position takes precedence. Furthermore, TP53 has strong nonspecific nucleosome binding facilitating its interaction with chromatin. Our in vitro findings were confirmed by examining TP53-induced binding in a cell line model, showing induced binding at nucleosome edges flanked by a nucleosome-free region. Overall, our results suggest that the pioneering capabilities of TP53 are driven by nonspecific nucleosome binding with specific binding at nucleosome edges.

A novel molecular rotor facilitates detection of p53-DNA interactions using the Fluorescent Intercalator Displacement Assay

We have investigated the use of fluorescent molecular rotors as probes for detection of p53 binding to DNA. These are a class of fluorophores that undergo twisted intramolecular charge transfer (TICT). They are non-fluorescent in a freely rotating conformation and experience a fluorescence increase when restricted in the planar conformation. We hypothesized that intercalation of a molecular rotor between DNA base pairs would result in a fluorescence turn-on signal. Upon displacement by a DNA binding protein, measurable loss of signal would facilitate use of the molecular rotor in the fluorescent intercalator displacement (FID) assay. A panel of probes was interrogated using the well-established p53 model system across various DNA response elements. A novel, readily synthesizable molecular rotor incorporating an acridine orange DNA intercalating group (AO-R) outperformed other conventional dyes in the FID assay. It enabled relative measurement of p53 sequence-specific DNA interactions and study of the dominant-negative effects of cancer-associated p53 mutants. In a further application, AO-R also proved useful for staining apoptotic cells in live zebrafish embryos.

True equilibrium measurement of transcription factor-DNA binding affinities using automated polarization microscopy

The complex patterns of gene expression in metazoans are controlled by selective binding of transcription factors (TFs) to regulatory DNA. To improve the quantitative understanding of this process, we have developed a novel method that uses fluorescence anisotropy measurements in a controlled delivery system to determine TF-DNA binding energies in solution with high sensitivity and throughput. Owing to its large dynamic range, the method, named high performance fluorescence anisotropy (HiP-FA), allows for reliable quantification of both weak and strong binding; binding specificities are calculated on the basis of equilibrium constant measurements for mutational DNA variants. We determine the binding preference landscapes for 26 TFs and measure high absolute affinities, but mostly lower binding specificities than reported by other methods. The revised binding preferences give rise to improved predictions of in vivo TF occupancy and enhancer expression. Our approach provides a powerful new tool for the systems-biological analysis of gene regulation.

Methods to measure selective transcription factor-DNA binding often lack sensitivity and are not performed in solution. Here the authors develop a method to perform fluorescence anisotropy measurements of transcription factor-DNA binding energies with high sensitivity and throughput.

Functional characterization of p53 pathway components in the ancient metazoan Trichoplax adhaerens

The identification of genes encoding a p53 family member and an Mdm2 ortholog in the ancient placozoan Trichoplax adhaerens advocates for the evolutionary conservation of a pivotal stress-response pathway observed in all higher eukaryotes. Here, we recapitulate several key functionalities ascribed to this known interacting protein pair by analysis of the placozoan proteins (Tap53 and TaMdm2) using both in vitro and cellular assays. In addition to interacting with each other, the Tap53 and TaMdm2 proteins are also able to respectively bind human Mdm2 and p53, providing strong evidence for functional conservation. The key p53-degrading function of Mdm2 is also conserved in TaMdm2. Tap53 retained DNA binding associated with p53 transcription activation function. However, it lacked transactivation function in reporter genes assays using a heterologous cell line, suggesting a cofactor incompatibility. Overall, the data supports functional roles for TaMdm2 and Tap53, and further defines the p53 pathway as an evolutionary conserved fulcrum mediating cellular response to stress.

Polymorphisms in cancer-related pathway genes and lung cancer

We evaluated the associations between potentially functional variants in a comprehensive list of cancer-related genes and lung cancer in a Korean population.

A total of 1969 potentially functional single nucleotide polymorphisms (SNPs) of 1151 genes involved in carcinogenesis were evaluated using an Affymetrix custom-made GeneChip in 610 nonsmall cell lung cancer patients and 610 healthy controls. A replication study was conducted in an independent set of 490 cases and 486 controls. 68 SNPs were significantly associated with lung cancer in the discovery set and tested for replication.

Among the 68 SNPs, three SNPs (corepressor interacting with RBPJ 1 (CIR1) rs13009079T>C, ribonucleotide reductase M1 (RRM1) rs1465952T>C and solute carrier family 38, member 4 (SLC38A4) rs2429467C>T) consistantly showed significant associations with lung cancer in the replication study. In combined analysis, adjusted odds ratio for CIR1 rs13009079T>C, RRM1 rs1465952T>C and SLC38A4 rs2429467C>T were 0.69, 0.71 and 0.73, respectively (p=4×10−5, 0.01 and 0.001, respectively) under the dominant model. The relative mRNA expression level of CIR1 was significantly associated with rs13009079T>C genotypes in normal lung tissues (ptrend=0.03).

These results suggest that the three SNPs, particularly CIR1 rs13009079T>C, may play a role in the pathogenesis of lung cancer.

Basal and copper-induced expression of metallothionein isoform 1,2 and 3 genes in epithelial cancer cells: The role of tumor suppressor p53

Metallothioneins (MTs) are a ubiquitous low-molecular weight, cysteine rich proteins with a high affinity for metal ions. The expression and induction of MTs have been associated with protection against DNA damage, oxidative stress, and apoptosis. Our past research had shown that p53 is an important factor in metal regulation of MTs. The present study was undertaken to explore further the interrelationship between p53 and MTs. We investigated whether silencing of p53 could affect expression pattern of basal and copper induced metallothioneins. The silencing of wild-type p53 (wt-p53) in epithelial breast cancer MCF7 cells affected the basal level of MT-2A RNA, whereas the levels of MT-1A and MT-1X RNA remained largely unchanged. The expression of MT-3 was undetectable in MCF7 with either functional or silenced p53. MCF7 cells with silenced wt-p53 failed to upregulate MT-2A in response to copper and showed a reduced sensitivity toward copper induced cell apoptotic death. Similarly in MCF7-E6 and MDA-MB-231 cells, the presence of inactive/mutated p53 halted MT-1A and MT-2A gene expression in response to copper. Constitutive expression of MT-3 RNA was detectable in the presence of mutated p53 (mtp53). Transient transfection of MDA-MB-231 cells with wt-p53 enabled copper induced upregulation of both MT-1A and MT-2A but not basal level of MT-2A, MT-1E, MT-1X and MT-3. Inactivation of p53 in HepG2 cells amplified the basal expression of studied MT isoforms, including MT-3, as well as copper-induced mRNA expression of MTs except MT-1H and MT-3. Presented data demonstrate a direct relation between p53 and MT-1A and MT-2A and they also indicate that wt-p53 might be a negative regulator of MT-3 in epithelial cancer cells.

An African-specific polymorphism in the TP53 gene impairs p53 tumor suppressor function in a mouse model

In this study, Jennis et al. characterize the first mouse model of an African-specific naturally occurring coding region variant at codon 47 of the p53 tumor suppressor gene (S47). They show that homozygous S47 mice are markedly tumor-prone and that the S47 variant impairs not only p53-mediated cell death but also the ability of p53 to transactivate a subset of genes involved in metabolism and ferroptosis.

A nonsynonymous single-nucleotide polymorphism at codon 47 in TP53 exists in African-descent populations (P47S, rs1800371; referred to here as S47). Here we report that, in human cell lines and a mouse model, the S47 variant exhibits a modest decrease in apoptosis in response to most genotoxic stresses compared with wild-type p53 but exhibits a significant defect in cell death induced by cisplatin. We show that, compared with wild-type p53, S47 has nearly indistinguishable transcriptional function but shows impaired ability to transactivate a subset of p53 target genes, including two involved in metabolism: Gls2 (glutaminase 2) and Sco2. We also show that human and mouse cells expressing the S47 variant are markedly resistant to cell death by agents that induce ferroptosis (iron-mediated nonapoptotic cell death). We show that mice expressing S47 in homozygous or heterozygous form are susceptible to spontaneous cancers of diverse histological types. Our data suggest that the S47 variant may contribute to increased cancer risk in individuals of African descent, and our findings highlight the need to assess the contribution of this variant to cancer risk in these populations. These data also confirm the potential relevance of metabolism and ferroptosis to tumor suppression by p53.

Regulation of carboxylesterase-2 expression by p53 family proteins and enhanced anti-cancer activities among 5-fluorouracil, irinotecan and doxazolidine prodrug

BACKGROUND AND PURPOSE

For four decades, 5-fluorouracil (5-FU) has been a major anti-cancer medicine. This drug is increasingly used with other anti-cancer agents such as irinotecan. Irinotecan and many others such as PPD (pentyl carbamate of p-aminobenzyl carbamate of doxazolidine) require activation by carboxylesterase-2 (CES2). 5-FU, on the other hand, reportedly induces CES2 in colorectal tumour lines. The aims of this study were to determine the molecular basis for the induction and to ascertain interactive cell-killing activity between 5-FU and ester prodrugs.

EXPERIMENTAL APPROACH

Colorectal and non-colorectal lines and xenografts were treated with 5-FU and the expression of CES2 was determined. Cell-killing activity of irinotecan and PPD were determined in the presence or absence of CES2 inhibitor. Several molecular experiments were used to determine the molecular basis for the induction.

KEY RESULTS

Without exceptions, robust induction was detected in cell lines expressing functional p53. High-level induction was also detected in xenografts. 5-FU pretreatment significantly increased cell-killing activity of irinotecan and PPD. Molecular experiments established that the induction was achieved by both transactivation and increased mRNA stability through p53. Either p63 or p73, functionally related to p53, did not support the transactivation.

CONCLUSIONS AND IMPLICATIONS

The results in this study suggest that FOLFIRI, a common regimen combining irinotecan and 5-FU, should switch the dosing sequence, namely from 5-FU to irinotecan, to enhance hydrolytic activation of irinotecan. This modified order likely reduces the dose of anti-cancer agents, thus minimizing overall toxicity. The results also conclude that p53 family members act differently in regulating gene expression.

Whole-genome cartography of p53 response elements ranked on transactivation potential

Background

Many recent studies using ChIP-seq approaches cross-referenced to trascriptome data and also to potentially unbiased in vitro DNA binding selection experiments are detailing with increasing precision the p53-directed gene regulatory network that, nevertheless, is still expanding. However, most experiments have been conducted in established cell lines subjected to specific p53-inducing stimuli, both factors potentially biasing the results.

Results

We developed p53retriever, a pattern search algorithm that maps p53 response elements (REs) and ranks them according to predicted transactivation potentials in five classes. Besides canonical, full site REs, we developed specific pattern searches for non-canonical half sites and 3/4 sites and show that they can mediate p53-dependent responsiveness of associated coding sequences. Using ENCODE data, we also mapped p53 REs in about 44,000 distant enhancers and identified a 16-fold enrichment for high activity REs within those sites in the comparison with genomic regions near transcriptional start sites (TSS). Predictions from our pattern search were cross-referenced to ChIP-seq, ChIP-exo, expression, and various literature data sources. Based on the mapping of predicted functional REs near TSS, we examined expression changes of thirteen genes as a function of different p53-inducing conditions, providing further evidence for PDE2A, GAS6, E2F7, APOBEC3H, KCTD1, TRIM32, DICER, HRAS, KITLG and TGFA p53-dependent regulation, while MAP2K3, DNAJA1 and potentially YAP1 were identified as new direct p53 target genes.

Conclusions

We provide a comprehensive annotation of canonical and non-canonical p53 REs in the human genome, ranked on predicted transactivation potential. We also establish or corroborate direct p53 transcriptional control of thirteen genes. The entire list of identified and functionally classified p53 REs near all UCSC-annotated genes and within ENCODE mapped enhancer elements is provided. Our approach is distinct from, and complementary to, existing methods designed to identify p53 response elements. p53retriever is available as an R package at: http://tomateba.github.io/p53retriever.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1643-9) contains supplementary material, which is available to authorized users.

Interactions of chromatin context, binding site sequence content, and sequence evolution in stress-induced p53 occupancy and transactivation

Cellular stresses activate the tumor suppressor p53 protein leading to selective binding to DNA response elements (REs) and gene transactivation from a large pool of potential p53 REs (p53REs). To elucidate how p53RE sequences and local chromatin context interact to affect p53 binding and gene transactivation, we mapped genome-wide binding localizations of p53 and H3K4me3 in untreated and doxorubicin (DXR)-treated human lymphoblastoid cells. We examined the relationships among p53 occupancy, gene expression, H3K4me3, chromatin accessibility (DNase 1 hypersensitivity, DHS), ENCODE chromatin states, p53RE sequence, and evolutionary conservation. We observed that the inducible expression of p53-regulated genes was associated with the steady-state chromatin status of the cell. Most highly inducible p53-regulated genes were suppressed at baseline and marked by repressive histone modifications or displayed CTCF binding. Comparison of p53RE sequences residing in different chromatin contexts demonstrated that weaker p53REs resided in open promoters, while stronger p53REs were located within enhancers and repressed chromatin. p53 occupancy was strongly correlated with similarity of the target DNA sequences to the p53RE consensus, but surprisingly, inversely correlated with pre-existing nucleosome accessibility (DHS) and evolutionary conservation at the p53RE. Occupancy by p53 of REs that overlapped transposable element (TE) repeats was significantly higher (p<10⁻⁷) and correlated with stronger p53RE sequences (p<10⁻¹¹⁰) relative to nonTE-associated p53REs, particularly for MLT1H, LTR10B, and Mer61 TEs. However, binding at these elements was generally not associated with transactivation of adjacent genes. Occupied p53REs located in L2-like TEs were unique in displaying highly negative PhyloP scores (predicted fast-evolving) and being associated with altered H3K4me3 and DHS levels. These results underscore the systematic interaction between chromatin status and p53RE context in the induced transactivation response. This p53 regulated response appears to have been tuned via evolutionary processes that may have led to repression and/or utilization of p53REs originating from primate-specific transposon elements.

It is well established that p53 binds DNA elements near p53 target genes to regulate the response to cellular stress. To assess factors influencing binding to response elements and subsequent gene expression, we have analyzed 2932 p53-occupied response elements (p53REs) in the context of genome-wide chromatin state, DNA accessibility and dynamics, and considered roles for binding-sequence specificity and evolutionary conservation. While p53 occupancy level shows little apparent direct relationship to gene expression change, after grouping expressed genes by their chromatin status at baseline, a relationship between occupancy of p53REs and gene expression change emerged. Analysis of p53RE sequences demonstrated that p53 occupancy was strongly correlated with sequence similarity to p53RE consensus, but surprisingly, was inversely correlated with nucleosome accessibility (DHS) and evolutionary conservation. These data revealed a systematic interaction between p53RE content and chromatin context that affects both quantitative p53 occupancy and the induced transactivation response to exposure. Moreover, this interaction appears to have been tuned via evolutionary events involving transposable elements, which strongly bind p53, but in only a few instances affect gene expression levels. Models of p53-regulated gene expression response that consider both chromatin state and sequence context may prove useful in guiding strategies for cancer prevention or therapy.

Linking polymorphic p53 response elements with gene expression in airway epithelial cells of smokers and cancer risk

Chronic cigarette smoking exposes airway epithelial cells to thousands of carcinogens, oxidants and DNA-damaging agents, creating a field of molecular injury in the airway and altering gene expression. Studies of cytologically normal bronchial epithelial cells from smokers have identified transcription-based biomarkers that may prove useful in early diagnosis of lung cancer, including a number of p53-regulated genes. The ability of p53 to regulate transcription is critical for tumor suppression, and this suggests that single-nucleotide polymorphisms (SNPs) in functional p53 binding sites (p53 response elements, or p53REs) that affect gene expression could influence susceptibility to cancer. To connect p53RE SNP genotype with gene expression and cancer risk, we identified a set of 204 SNPs in putative p53REs, and performed cis expression quantitative trait loci (eQTL) analysis, assessing associations between SNP genotypes and mRNA levels of adjacent genes in bronchial epithelial cells obtained from 44 cigarette smokers. To further test and validate these genotype-expression associations, we searched published eQTL studies from independent populations and determined that 53% (39/74) of the bronchial epithelial eQTLs were observed in at least one of other studies. SNPs in p53REs were also evaluated for effects on p53-DNA binding using a quantitative in vitro protein-DNA binding assay. Last, based on linkage disequilibrium, we found 6 p53RE SNPs associated with gene expression were identified as cancer risk SNPs by either genome-wide association studies or candidate gene studies. We provide an approach for identifying and evaluating potentially functional SNPs that may modulate the airway gene expression response to smoking and may influence susceptibility to cancers.

Rapid and Sensitive Detection of p53 Based on DNA-Protein Binding Interactions Using Silver Nanoparticle Films and Microwave Heating

Tumor detection can be carried out via the detection of proteins, such as p53, which is known to play vital role in more than 50% of all cancers affecting humans. Early diagnosis of tumor detection can be achieved by decreasing the lower detection limit of p53 bioassays. Microwave-accelerated bioassay (MAB) technique, which is based on the use of circular bioassay platforms in combination with microwave heating, is employed for the rapid and sensitive detection of p53 protein. Direct sandwich ELISA was constructed on our circular bioassay platforms based on DNA-protein binding interactions. Colorimetric and fluorescence based detection methods were used for room temperature bioassay (control bioassay; total bioassay time is 27 hours) and bioassay using microwave heating (i.e., the MAB technique; total bioassay time is 10 minutes). In the colorimetric based detection, a very high background signal due to the non-specific binding of proteins for the bioassay carried out at room temperature and a LLOD of 0.01 ng/mL for p53 was observed using the MAB technique. The LLOD for the fluorescence-based detection using the MAB technique was found to be 0.01 ng/mL. The use of circular bioassay platforms in the MAB technique results in microwave-induced temperature gradient, where the specific protein binding interactions are significantly accelerated; thereby reducing the background signal and the lower limit of detection of p53 protein.

A functional p53 responsive polymorphism in KITLG, rs4590952, does not affect the risk of breast cancer

Recently, a functional polymorphism in KITLG, rs4590952, was identified to be associated with testicular cancer susceptibility through increasing the p53-dependent KITLG expression and disrupting the function of p53. We performed a hospital-based case-control study, including 1241 breast cancer (BC) patients and 1259 cancer-free controls, to investigate the role of this polymorphism in the risk of BC in Chinese Han population. However, no significant association between rs4590952 and BC risk was identified in allelic model with the odds ratio (OR) of 1.04 (95% confidence interval (CI) = 0.73–1.46, P = 0.839) or in any other genetic models. When performed stratified analysis according to the Estrogen Receptor (ER) and Progesterone Receptor (PR) status, rs4590952 was neither associated with ER+/PR+ nor ER−PR− subgroups. Our results suggested that rs4590952 was not associated with the risk of BC in Chinese population, implying that heterogeneous distinct mechanisms might exist in the etiology of different cancers.

A polymorphic p53 response element in KIT ligand influences cancer risk and has undergone natural selection

The ability of p53 to regulate transcription is crucial for tumor suppression and implies that inherited polymorphisms in functional p53-binding sites could influence cancer. Here, we identify a polymorphic p53 responsive element and demonstrate its influence on cancer risk using genome-wide data sets of cancer susceptibility loci, genetic variation, p53 occupancy, and p53-binding sites. We uncover a single-nucleotide polymorphism (SNP) in a functional p53-binding site and establish its influence on the ability of p53 to bind to and regulate transcription of the KITLG gene. The SNP resides in KITLG and associates with one of the largest risks identified among cancer genome-wide association studies. We establish that the SNP has undergone positive selection throughout evolution, signifying a selective benefit, but go on to show that similar SNPs are rare in the genome due to negative selection, indicating that polymorphisms in p53-binding sites are primarily detrimental to humans.

Diverse stresses dramatically alter genome-wide p53 binding and transactivation landscape in human cancer cells

The effects of diverse stresses on promoter selectivity and transcription regulation by the tumor suppressor p53 are poorly understood. We have taken a comprehensive approach to characterizing the human p53 network that includes p53 levels, binding, expression and chromatin changes under diverse stresses. Human osteosarcoma U2OS cells treated with anti-cancer drugs Doxorubicin (DXR) or Nutlin-3 (Nutlin) led to strikingly different p53 gene binding patterns based on chromatin immunoprecipitation with high-throughput sequencing experiments. Although two contiguous RRRCWWGYYY decamers is the consensus binding motif, p53 can bind a single decamer and function in vivo. Although the number of sites bound by p53 was six times greater for Nutlin than DXR, expression changes induced by Nutlin were much less dramatic compared with DXR. Unexpectedly, the solvent dimethylsulphoxide (DMSO) alone induced p53 binding to many sites common to DXR; however, this binding had no effect on target gene expression. Together, these data imply a two-stage mechanism for p53 transactivation where p53 binding only constitutes the first stage. Furthermore, both p53 binding and transactivation were associated with increased active histone modification histone H3 lysine 4 trimethylation. We discovered 149 putative new p53 target genes including several that are relevant to tumor suppression, revealing potential new targets for cancer therapy and expanding our understanding of the p53 regulatory network.

In vitro selection of mutant HDM2 resistant to Nutlin inhibition

HDM2 binds to the p53 tumour suppressor and targets it for proteosomal degradation. Presently in clinical trials, the small molecule Nutlin-3A competitively binds to HDM2 and abrogates its repressive function. Using a novel in vitro selection methodology, we simulated the emergence of resistance by evolving HDM2 mutants capable of binding p53 in the presence of Nutlin concentrations that inhibit the wild-type HDM2-p53 interaction. The in vitro phenotypes were recapitulated in ex vivo assays measuring both p53 transactivation function and the direct p53-HDM2 interaction in the presence of Nutlin. Mutations conferring drug resistance were not confined to the N-terminal p53/Nutlin–binding domain, and were additionally seen in the acidic, zinc finger and RING domains. Mechanistic insights gleaned from this broad spectrum of mutations will aid in future drug design and further our understanding of the complex p53-HDM2 interaction.

Low-level p53 expression changes transactivation rules and reveals superactivating sequences

Transcriptional activation by the tumor suppressor p53 is considered to depend on cellular level, although there are few systems where this dependence on cellular level of p53 has been directly addressed. Previously, we reported that transactivation from p53 targets was sensitive to both p53 amount and DNA sequence, with some sequences being responsive to much lower p53 levels than others when examined in yeast model systems or human cells. Because p53 is normally present at low levels and perturbations might lead to small increases, we examined transactivation under limiting p53. Unlike the positive relationship between transactivation and binding affinity from target sequences at high cellular levels of human p53 in yeast, no such relationship was found at low levels. However, transactivation in the yeast system and the torsional flexibility of target sequences were highly correlated, revealing a unique structural relationship between transcriptional function and sequence. Surprisingly, a few sequences supported high transactivation at low p53 levels in yeast or when transfected into human cells. On the basis of kinetic and flexibility analyses the "supertransactivation" property was due to low binding off rates of flexible target sites. Interestingly, a supertransactivation response element can differentiate transcriptional capacities of many breast cancer-associated p53 mutants. Overall, these studies, which are relevant to other transcription factors, address the extent to which transactivation properties of p53 target sequences are determined by their intrinsic physical properties and reveal unique rules of engagement of target sequences at low p53 levels.

Analysis of p53 binding to DNA by fluorescence imaging microscopy

Transcription factors play a central role in cell biology through binding to target DNA elements and regulating gene expression. In this study, we used the p53 tumour suppressor as a model transcription factor to develop an imaging based assay to measure DNA binding. The assay utilizes fluorescence imaging microscopy to detect labelled p53 bound to DNA coated on microbeads. We demonstrate the ability to multiplex the assay by interrogating simultaneous binding to variant DNA sequences present on tractable beads. Additionally, the assay measures activation of p53 for increased DNA binding by a known peptide in addition to reactivation of mutant p53 by a small molecule. It may therefore be adaptable to a high-content imaging screen for compounds capable of restoring the function of mutant p53 associated with cancer.

A genome-wide screen for genetic variants that modify the recruitment of REST to its target genes

Increasing numbers of human diseases are being linked to genetic variants, but our understanding of the mechanistic links leading from DNA sequence to disease phenotype is limited. The majority of disease-causing nucleotide variants fall within the non-protein-coding portion of the genome, making it likely that they act by altering gene regulatory sequences. We hypothesised that SNPs within the binding sites of the transcriptional repressor REST alter the degree of repression of target genes. Given that changes in the effective concentration of REST contribute to several pathologies—various cancers, Huntington's disease, cardiac hypertrophy, vascular smooth muscle proliferation—these SNPs should alter disease-susceptibility in carriers. We devised a strategy to identify SNPs that affect the recruitment of REST to target genes through the alteration of its DNA recognition element, the RE1. A multi-step screen combining genetic, genomic, and experimental filters yielded 56 polymorphic RE1 sequences with robust and statistically significant differences of affinity between alleles. These SNPs have a considerable effect on the the functional recruitment of REST to DNA in a range of in vitro, reporter gene, and in vivo analyses. Furthermore, we observe allele-specific biases in deeply sequenced chromatin immunoprecipitation data, consistent with predicted differenes in RE1 affinity. Amongst the targets of polymorphic RE1 elements are important disease genes including NPPA, PTPRT, and CDH4. Thus, considerable genetic variation exists in the DNA motifs that connect gene regulatory networks. Recently available ChIP–seq data allow the annotation of human genetic polymorphisms with regulatory information to generate prior hypotheses about their disease-causing mechanism.

Common human diseases such as cancer, heart disease, or epilepsy have a genetic component that predisposes particular individuals to suffer from them. Huge sums have been invested to map the regions of the human genome where small DNA variations, or SNPs (“single-nucleotide polymorphisms”), determine the probability of developing these diseases. A major problem with this approach, however, is that, once the culprit SNPs are discovered, we know very little about how they cause disease—which is critical if we are to use this information to develop drugs and therapies. In this study, we demonstrate a new approach, employing functional maps of the human genome that have recently been published. We begin with regions of the genome recognised by a gene repressor protein—REST—that is involved in a number of important human diseases. Using information on where REST binds in the human genome, we predict and validate common DNA variations that increase or decrease this binding. By affecting how much REST is recruited to important genes, these variations may predispose or protect individuals from a number of diseases. Studies like this show how we can use genomic information to gain a deeper understanding of the genetics behind human disease.

Widespread site-dependent buffering of human regulatory polymorphism

The average individual is expected to harbor thousands of variants within non-coding genomic regions involved in gene regulation. However, it is currently not possible to interpret reliably the functional consequences of genetic variation within any given transcription factor recognition sequence. To address this, we comprehensively analyzed heritable genome-wide binding patterns of a major sequence-specific regulator (CTCF) in relation to genetic variability in binding site sequences across a multi-generational pedigree. We localized and quantified CTCF occupancy by ChIP-seq in 12 related and unrelated individuals spanning three generations, followed by comprehensive targeted resequencing of the entire CTCF–binding landscape across all individuals. We identified hundreds of variants with reproducible quantitative effects on CTCF occupancy (both positive and negative). While these effects paralleled protein–DNA recognition energetics when averaged, they were extensively buffered by striking local context dependencies. In the significant majority of cases buffering was complete, resulting in silent variants spanning every position within the DNA recognition interface irrespective of level of binding energy or evolutionary constraint. The prevalence of complex partial or complete buffering effects severely constrained the ability to predict reliably the impact of variation within any given binding site instance. Surprisingly, 40% of variants that increased CTCF occupancy occurred at positions of human–chimp divergence, challenging the expectation that the vast majority of functional regulatory variants should be deleterious. Our results suggest that, even in the presence of “perfect” genetic information afforded by resequencing and parallel studies in multiple related individuals, genomic site-specific prediction of the consequences of individual variation in regulatory DNA will require systematic coupling with empirical functional genomic measurements.

A comprehensive understanding of the contribution of individual genome sequences to disease and quantitative traits will require the general ability to predict consequences of genetic variation in non-protein-coding regions, particularly those involved in gene regulation. Here we tested the power to predict such consequences when presented with “complete” information encompassing the genomic DNA binding site patterns of a well-studied regulatory protein across multiple related individuals, coupled with all individual genome sequences at the binding positions. We find that, while there is reasonable ability to predict the average effects of variation within the consensus recognition sequence of a transcriptional regulator, it is not possible to determine reliably the consequences of variation at any given genomic instance. This suggests that the interpretation of individual genome sequences will require comprehensive complementation with functional genomic studies.

p53 Isoforms: An Intracellular Microprocessor?

Normal function of the p53 pathway is ubiquitously lost in cancers either through mutation or inactivating interaction with viral or cellular proteins. However, it is difficult in clinical studies to link p53 mutation status to cancer treatment and clinical outcome, suggesting that the p53 pathway is not fully understood. We have recently reported that the human p53 gene expresses not only 1 but 12 different p53 proteins (isoforms) due to alternative splicing, alternative initiation of translation, and alternative promoter usage. p53 isoform proteins thus contain distinct protein domains. They are expressed in normal human tissues but are abnormally expressed in a wide range of cancer types. We have recently reported that p53 isoform expression is associated with breast cancer prognosis, suggesting that they play a role in carcinogenesis. Indeed, the cellular response to damages can be switched from cell cycle arrest to apoptosis by only manipulating p53 isoform expression. This may provide an explanation to the hitherto inconsistent relationship between p53 mutation, treatment response, and outcome in breast cancer. However, the molecular mechanism is still unknown. Recent reports suggest that it involves modulation of gene expression in a p53-dependent and -independent manner. In this review, we summarize our current knowledge about the biological activities of p53 isoforms and propose a molecular mechanism conciliating our current knowledge on p53 and integrating p63 and p73 isoforms in the p53 pathway.

The Toll-like receptor gene family is integrated into human DNA damage and p53 networks

In recent years the functions that the p53 tumor suppressor plays in human biology have been greatly extended beyond “guardian of the genome.” Our studies of promoter response element sequences targeted by the p53 master regulatory transcription factor suggest a general role for this DNA damage and stress-responsive regulator in the control of human Toll-like receptor (TLR) gene expression. The TLR gene family mediates innate immunity to a wide variety of pathogenic threats through recognition of conserved pathogen-associated molecular motifs. Using primary human immune cells, we have examined expression of the entire TLR gene family following exposure to anti-cancer agents that induce the p53 network. Expression of all TLR genes, TLR1 to TLR10, in blood lymphocytes and alveolar macrophages from healthy volunteers can be induced by DNA metabolic stressors. However, there is considerable inter-individual variability. Most of the TLR genes respond to p53 via canonical as well as noncanonical promoter binding sites. Importantly, the integration of the TLR gene family into the p53 network is unique to primates, a recurrent theme raised for other gene families in our previous studies. Furthermore, a polymorphism in a TLR8 response element provides the first human example of a p53 target sequence specifically responsible for endogenous gene induction. These findings—demonstrating that the human innate immune system, including downstream induction of cytokines, can be modulated by DNA metabolic stress—have many implications for health and disease, as well as for understanding the evolution of damage and p53 responsive networks.

Among the most prominently studied regulators of gene function is the p53 tumor suppressor, which has many roles in human biology. The transcriptional master regulator p53 directly targets expression of >200 genes. Previously, we sought to define the p53 network in terms of functionality, specifically the ability of target response element sequences (REs) to support p53 transactivation. Here we identify p53 target canonical and noncanonical REs in the family of Toll-like Receptor (TLR) innate immune response genes and establish p53 regulation of most TLR genes. We address p53 responsiveness in primary human lymphocytes and alveolar macrophages collected from healthy volunteers. Notably, all TLR genes show responses to DNA damage, and most are p53-mediated. However, there is considerable variability between individuals, suggesting that DNA and p53 metabolic stresses can markedly differ in impact on the innate immune system as well as downstream appearance of cytokines. Indeed, we report a SNP in a p53 RE within the TLR8 promoter that alters p53 responsiveness in primary human cells. Furthermore, the p53-mediated expression of TLRs is unique to primates. Overall, these findings identify a new, pivotal role for the well-known human tumor suppressor p53, namely, integration of DNA damage and innate immune responses.

Formation of stress-specific p53 binding patterns is influenced by chromatin but not by modulation of p53 binding affinity to response elements

The p53 protein is crucial for adapting programs of gene expression in response to stress. Recently, we revealed that this occurs partly through the formation of stress-specific p53 binding patterns. However, the mechanisms that generate these binding patterns remain largely unknown. It is not established whether the selective binding of p53 is achieved through modulation of its binding affinity to certain response elements (REs) or via a chromatin-dependent mechanism. To shed light on this issue, we used a microsphere assay for protein-DNA binding to measure p53 binding patterns on naked DNA. In parallel, we measured p53 binding patterns within chromatin using chromatin immunoprecipitation and DNase I coupled to ligation-mediated polymerase chain reaction footprinting. Through this experimental approach, we revealed that UVB and Nutlin-3 doses, which lead to different cellular outcomes, induce similar p53 binding patterns on naked DNA. Conversely, the same treatments lead to stress-specific p53 binding patterns on chromatin. We show further that altering chromatin remodeling using an histone acetyltransferase inhibitor reduces p53 binding to REs. Altogether, our results reveal that the formation of p53 binding patterns is not due to the modulation of sequence-specific p53 binding affinity. Rather, we propose that chromatin and chromatin remodeling are required in this process.

Human single-nucleotide polymorphisms alter p53 sequence-specific binding at gene regulatory elements

p53 coordinates the expression of an intricate network of genes in response to stress signals. Sequence-specific DNA binding is essential for p53-mediated tumor suppression. We evaluated the impact of single-nucleotide polymorphisms (SNPs) in p53 response elements (p53RE) on DNA binding and gene expression in response to DNA damage. Using a bioinformatics approach based on incorporating p53 binding strength into a position weight matrix, we selected 32 SNPs in putative and validated p53REs. The microsphere assay for protein–DNA binding (MAPD) and allele-specific expression analysis was employed to assess the impact of SNPs on p53-DNA binding and gene expression, respectively. Comparing activated p53 binding in nuclear extracts from doxorubicin- or ionizing radiation (IR)-treated human cells, we observed little difference in binding profiles. Significant p53 binding was observed for most polymorphic REs and several displayed binding comparable to the p21 RE. SNP alleles predicted to lower p53 binding indeed reduced binding in 25 of the 32 sequences. Chromatin immunoprecipitation-sequencing in lymphoblastoid cells confirmed p53 binding to seven polymorphic p53 REs in response to doxorubicin. In addition, five polymorphisms were associated with altered gene expression following doxorubicin treatment. Our findings demonstrate an effective strategy to identify and evaluate SNPs that may alter p53-mediated stress responses.

Lysine120 interactions with p53 response elements can allosterically direct p53 organization

p53 can serve as a paradigm in studies aiming to figure out how allosteric perturbations in transcription factors (TFs) triggered by small changes in DNA response element (RE) sequences, can spell selectivity in co-factor recruitment. p53-REs are 20-base pair (bp) DNA segments specifying diverse functions. They may be located near the transcription start sites or thousands of bps away in the genome. Their number has been estimated to be in the thousands, and they all share a common motif. A key question is then how does the p53 protein recognize a particular p53-RE sequence among all the similar ones? Here, representative p53-REs regulating diverse functions including cell cycle arrest, DNA repair, and apoptosis were simulated in explicit solvent. Among the major interactions between p53 and its REs involving Lys120, Arg280 and Arg248, the bps interacting with Lys120 vary while the interacting partners of other residues are less so. We observe that each p53-RE quarter site sequence has a unique pattern of interactions with p53 Lys120. The allosteric, DNA sequence-induced conformational and dynamic changes of the altered Lys120 interactions are amplified by the perturbation of other p53-DNA interactions. The combined subtle RE sequence-specific allosteric effects propagate in the p53 and in the DNA. The resulting amplified allosteric effects far away are reflected in changes in the overall p53 organization and in the p53 surface topology and residue fluctuations which play key roles in selective co-factor recruitment. As such, these observations suggest how similar p53-RE sequences can spell the preferred co-factor binding, which is the key to the selective gene transactivation and consequently different functional effects.

Single-nucleotide polymorphisms in the p53 signaling pathway

The p53 tumor suppressor pathway is central both in reducing cancer frequency in vertebrates and in mediating the response of commonly used cancer therapies. This article aims to summarize and discuss a large body of evidence suggesting that the p53 pathway harbors functional inherited single-nucleotide polymorphisms (SNPs) that affect p53 signaling in cells, resulting in differences in cancer risk and clinical outcome in humans. The insights gained through these studies into how the functional p53 pathway SNPs could help in the tailoring of cancer therapies to the individual are discussed. Moreover, recent work is discussed that suggests that many more functional p53 pathway SNPs are yet to be fully characterized and that a thorough analysis of the functional human genetics of this important tumor suppressor pathway is required.

Conservation of DNA-binding specificity and oligomerisation properties within the p53 family

Background

Transcription factors activate their target genes by binding to specific response elements. Many transcription factor families evolved from a common ancestor by gene duplication and subsequent divergent evolution. Members of the p53 family, which play key roles in cell-cycle control and development, share conserved DNA binding and oligomerisation domains but exhibit distinct functions. In this study, the molecular basis of the functional divergence of related transcription factors was investigated.

Results

We characterised the DNA-binding specificity and oligomerisation properties of human p53, p63 and p73, as well as p53 from other organisms using novel biophysical approaches. All p53 family members bound DNA cooperatively as tetramers with high affinity. Despite structural differences in the oligomerisation domain, the dissociation constants of the tetramers was in the low nanomolar range for all family members, indicating that the strength of tetramerisation was evolutionarily conserved. However, small differences in the oligomerisation properties were observed, which may play a regulatory role. Intriguingly, the DNA-binding specificity of p53 family members was highly conserved even for evolutionarily distant species. Additionally, DNA recognition was only weakly affected by CpG methylation. Prediction of p53/p63/p73 binding sites in the genome showed almost complete overlap between the different homologs.

Conclusion

Diversity of biological function of p53 family members is not reflected in differences in sequence-specific DNA binding. Hence, additional specificity factors must exist, which allowed the acquisition of novel functions during evolution while preserving original roles.

The expanding universe of p53 targets

The p53 tumour suppressor is modified through mutation or changes in expression in most cancers, leading to the altered regulation of hundreds of genes that are directly influenced by this sequence-specific transcription factor. Central to the p53 master regulatory network are the target response element (RE) sequences. The extent of p53 transactivation and transcriptional repression is influenced by many factors, including p53 levels, cofactors and the specific RE sequences, all of which contribute to the role that p53 has in the aetiology of cancer. This Review describes the identification and functionality of REs and highlights the inclusion of non-canonical REs that expand the universe of genes and regulation mechanisms in the p53 tumour suppressor network.