p53 binding to nucleosomal DNA depends on the rotational positioning of DNA response element

Determinants of p53 DNA binding, gene regulation, and cell fate decisions

The extent to which transcription factors read and respond to specific information content within short DNA sequences remains an important question that the tumor suppressor p53 is helping us answer. We discuss recent insights into how local information content at p53 binding sites might control modes of p53 target gene activation and cell fate decisions. Significant prior work has yielded data supporting two potential models of how p53 determines cell fate through its target genes: a selective target gene binding and activation model and a p53 level threshold model. Both of these models largely revolve around an analogy of whether p53 is acting in a "smart" or "dumb" manner. Here, we synthesize recent and past studies on p53 decoding of DNA sequence, chromatin context, and cellular signaling cascades to elicit variable cell fates critical in human development, homeostasis, and disease.

Spectral analysis of DNA superhelical dynamics from molecular minicircle simulations

Torsional and bending deformations of DNA molecules often occur in vivo and are important for biological functions. DNA "under stress" is a conformational state, which is by far the most frequent state during DNA-protein and gene regulation. In DNA minicircles of length <100 base pairs (bp), the combined effect of torsional and bending stresses can cause local unusual conformations, with certain base pair steps often absorbing most of the stress, leaving other steps close to their relaxed conformation. To better understand the superhelical dynamics of DNA under stress, molecular simulations of 94 bp minicircles with different torsional linking numbers were interpreted using Fourier analyses and principal component analyses. Sharp localized bends of nearly 90° in the helical axis were observed, which in turn decreased fluctuations of the rotational register and helped redistribute the torsional stress into writhe, i.e., superhelical turn up to 360°. In these kinked minicircles, only two-thirds of the DNA molecule bends and writhes and the remaining segment stays close to straight and preserves a conformational flexibility typical of canonical B-DNA (bending of 39° ± 17° distributed parsimoniously across 36 bp), which was confirmed and visualized by principal component analysis. These results confirm that stressed DNA molecules are highly heterogeneous along their sequence, with segments designed to locally store and release stress so that nearby segments can stay relaxed.

Symphony of the DNA flexibility and sequence environment orchestrates p53 binding to its responsive elements

The p53 tumor suppressor protein maintains the genome fidelity and integrity by modulating several cellular activities. It regulates these events by interacting with a heterogeneous set of response elements (REs) of regulatory genes in the background of chromatin configuration. At the p53-RE interface, both the base readout and torsional-flexibility of DNA account for high-affinity binding. However, DNA structure is an entanglement of a multitude of physicochemical features, both local and global structure should be considered for dealing with DNA-protein interactions. The goal of current research work is to conceptualize and abstract basic principles of p53-RE binding affinity as a function of structural alterations in DNA such as bending, twisting, and stretching flexibility and shape. For this purpose, we have exploited high throughput in-vitro relative affinity information of responsive elements and genome binding events of p53 from HT-Selex and ChIP-Seq experiments respectively. Our results confirm the role of torsional flexibility in p53 binding, and further, we reveal that DNA axial bending, stretching stiffness, propeller twist, and wedge angles are intimately linked to p53 binding affinity when compared to homeodomain, bZIP, and bHLH proteins. Besides, a similar DNA structural environment is observed in the distal sequences encompassing the actual binding sites of p53 cistrome genes. Additionally, we revealed that p53 cistrome target genes have unique promoter architecture, and the DNA flexibility of genomic sequences around REs in cancer and normal cell types display major differences. Altogether, our work provides a keynote on DNA structural features of REs that shape up the in-vitro and in-vivo high-affinity binding of the p53 transcription factor.

Pioneer factor Foxa2 enables ligand-dependent activation of type II nuclear receptors FXR and LXRα

Objective

Type II nuclear hormone receptors, including farnesoid X receptors (FXR), liver X receptors (LXR), and peroxisome proliferator-activated receptors (PPAR), which serve as drug targets for metabolic diseases, are permanently positioned in the nucleus and thought to be bound to DNA regardless of the ligand status. However, recent genome-wide location analysis showed that LXRα and PPARα binding in the liver is largely ligand-dependent. We hypothesized that pioneer factor Foxa2 evicts nucleosomes to enable ligand-dependent binding of type II nuclear receptors and performed genome-wide studies to test this hypothesis.

Methods

ATAC-Seq was used to profile chromatin accessibility; ChIP-Seq was performed to assess transcription factors (Foxa2, FXR, LXRα, and PPARα) binding; and RNA-Seq analysis determined differentially expressed genes in wildtype and Foxa2 mutants treated with a ligand (GW4064 for FXR, GW3965, and T09 for LXRα).

Results

We reveal that chromatin accessibility, FXR binding, LXRα occupancy, and ligand-responsive activation of gene expression by FXR and LXRα require Foxa2. Unexpectedly, Foxa2 occupancy is drastically increased when either receptor, FXR or LXRα, is bound by an agonist. In addition, co-immunoprecipitation experiments demonstrate that Foxa2 interacts with either receptor in a ligand-dependent manner, suggesting that Foxa2 and the receptor, bind DNA as an interdependent complex during ligand activation. Furthermore, PPARα binding is induced in Foxa2 mutants treated with FXR and LXR ligands, leading to the activation of PPARα targets.

Conclusions

Our model requires pioneering activity for ligand activation that challenges the existing ligand-independent binding mechanism. We also demonstrate that Foxa2 is required to achieve activation of the proper receptor – one that binds the added ligand – by repressing the activity of a competing receptor.

•

Foxa2 opens chromatin for FXR and LXRα binding during acute ligand activation.

•

Ligand-dependent activation of FXR & LXR-dependent gene expression requires Foxa2.

•

Foxa2 interacts with FXR and LXRα in a ligand-dependent manner.

•

Foxa2 restricts binding of competing receptor PPARα to ensure proper ligand-dependent activation of FXR and LXRα.

p63 and p53: Collaborative Partners or Dueling Rivals?

The tumor suppressor p53 and its oncogenic sibling p63 (ΔNp63) direct opposing fates in tumor development. These paralog proteins are transcription factors that elicit their tumor suppressive and oncogenic capacity through the regulation of both shared and unique target genes. Both proteins predominantly function as activators of transcription, leading to a paradigm shift away from ΔNp63 as a dominant negative to p53 activity. The discovery of p53 and p63 as pioneer transcription factors regulating chromatin structure revealed new insights into how these paralogs can both positively and negatively influence each other to direct cell fate. The previous view of a strict rivalry between the siblings needs to be revisited, as p53 and p63 can also work together toward a common goal.

Machine learning predicts nucleosome binding modes of transcription factors

Background

Most transcription factors (TFs) compete with nucleosomes to gain access to their cognate binding sites. Recent studies have identified several TF-nucleosome interaction modes including end binding (EB), oriented binding, periodic binding, dyad binding, groove binding, and gyre spanning. However, there are substantial experimental challenges in measuring nucleosome binding modes for thousands of TFs in different species.

Results

We present a computational prediction of the binding modes based on TF protein sequences. With a nested cross-validation procedure, our model outperforms several fine-tuned off-the-shelf machine learning (ML) methods in the multi-label classification task. Our binary classifier for the EB mode performs better than these ML methods with the area under precision-recall curve achieving 75%. The end preference of most TFs is consistent with low nucleosome occupancy around their binding site in GM12878 cells. The nucleosome occupancy data is used as an alternative dataset to confirm the superiority of our EB classifier.

Conclusions

We develop the first ML-based approach for efficient and comprehensive analysis of nucleosome binding modes of TFs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04093-9.

Temporal gene regulation by p53 is associated with the rotational setting of its binding sites in nucleosomes

The tumor suppressor protein p53 is a DNA-binding transcription factor (TF) that, once activated, coordinates the expression of thousands of target genes. Increased p53 binding to gene promoters occurs shortly after p53 activation. Intriguingly, gene transcription exhibits differential kinetics with some genes being induced early (early genes) and others being induced late (late genes). To understand pre-binding factors contributing to the temporal gene regulation by p53, we performed time-course RNA sequencing experiments in human colon cancer cell line HCT116 treated with fluorouracil to identify early and late genes. Published p53 ChIP fragments co-localized with the early or late genes were used to uncover p53 binding sites (BS). We demonstrate that the BS associated with early genes are clustered around gene starts with decreased nucleosome occupancy. DNA analysis shows that these BS are likely exposed on nucleosomal surface if wrapped into nucleosomes, thereby facilitating stable interactions with and fast induction by p53. By contrast, p53 BS associated with late genes are distributed uniformly across the genes with increased nucleosome occupancy. Predicted rotational settings of these BS show limited accessibility. We therefore propose a hypothetical model in which the BS are fully, partially or not accessible to p53 in the nucleosomal context. The partial accessibility of the BS allows subunits of a p53 tetramer to bind, but the resulting p53-DNA complex may not be stable enough to recruit cofactors, which leads to delayed induction. Our work highlights the importance of DNA conformations of p53 BS in gene expression dynamics.

Tumor suppressor p53: from engaging DNA to target gene regulation

The p53 transcription factor confers its potent tumor suppressor functions primarily through the regulation of a large network of target genes. The recent explosion of next generation sequencing protocols has enabled the study of the p53 gene regulatory network (GRN) and underlying mechanisms at an unprecedented depth and scale, helping us to understand precisely how p53 controls gene regulation. Here, we discuss our current understanding of where and how p53 binds to DNA and chromatin, its pioneer-like role, and how this affects gene regulation. We provide an overview of the p53 GRN and the direct and indirect mechanisms through which p53 affects gene regulation. In particular, we focus on delineating the ubiquitous and cell type-specific network of regulatory elements that p53 engages; reviewing our understanding of how, where, and when p53 binds to DNA and the mechanisms through which these events regulate transcription. Finally, we discuss the evolution of the p53 GRN and how recent work has revealed remarkable differences between vertebrates, which are of particular importance to cancer researchers using mouse models.

STARR-seq identifies active, chromatin-masked, and dormant enhancers in pluripotent mouse embryonic stem cells

BACKGROUND

Enhancers are distal regulators of gene expression that shape cell identity and control cell fate transitions. In mouse embryonic stem cells (mESCs), the pluripotency network is maintained by the function of a complex network of enhancers, that are drastically altered upon differentiation. Genome-wide chromatin accessibility and histone modification assays are commonly used as a proxy for identifying putative enhancers and for describing their activity levels and dynamics.

RESULTS

Here, we applied STARR-seq, a genome-wide plasmid-based assay, as a read-out for the enhancer landscape in "ground-state" (2i+LIF; 2iL) and "metastable" (serum+LIF; SL) mESCs. This analysis reveals that active STARR-seq loci show modest overlap with enhancer locations derived from peak calling of ChIP-seq libraries for common enhancer marks. We unveil ZIC3-bound loci with significant STARR-seq activity in SL-ESCs. Knock-out of Zic3 removes STARR-seq activity only in SL-ESCs and increases their propensity to differentiate towards the endodermal fate. STARR-seq also reveals enhancers that are not accessible, masked by a repressive chromatin signature. We describe a class of dormant, p53 bound enhancers that gain H3K27ac under specific conditions, such as after treatment with Nocodazol, or transiently during reprogramming from fibroblasts to pluripotency.

CONCLUSIONS

In conclusion, loci identified as active by STARR-seq often overlap with those identified by chromatin accessibility and active epigenetic marking, yet a significant fraction is epigenetically repressed or display condition-specific enhancer activity.

Nucleosome allostery in pioneer transcription factor binding

While recent experiments revealed that some pioneer transcription factors (TFs) can bind to their target DNA sequences inside a nucleosome, the binding dynamics of their target recognitions are poorly understood. Here we used the latest coarse-grained models and molecular dynamics simulations to study the nucleosome-binding procedure of the two pioneer TFs, Sox2 and Oct4. In the simulations for a strongly positioning nucleosome, Sox2 selected its target DNA sequence only when the target was exposed. Otherwise, Sox2 entropically bound to the dyad region nonspecifically. In contrast, Oct4 plastically bound on the nucleosome mainly in two ways. First, the two POU domains of Oct4 separately bound to the two parallel gyres of the nucleosomal DNA, supporting the previous experimental results of the partial motif recognition. Second, the POU_S domain of Oct4 favored binding on the acidic patch of histones. Then, simulating the TFs binding to a genomic nucleosome, the LIN28B nucleosome, we found that the recognition of a pseudo motif by Sox2 induced the local DNA bending and shifted the population of the rotational position of the nucleosomal DNA. The redistributed DNA phase, in turn, changed the accessibility of a distant TF binding site, which consequently affected the binding probability of a second Sox2 or Oct4. These results revealed a nucleosomal DNA-mediated allosteric mechanism, through which one TF binding event can change the global conformation, and effectively regulate the binding of another TF at distant sites. Our simulations provide insights into the binding mechanism of single and multiple TFs on the nucleosome.

Pioneer factors and their in vitro identification methods

Pioneer transcription factors are a special group of transcription factors that can interact with nucleosomal DNA and initiate regulatory events. Their binding to regulatory regions is the first event in gene activation and can occur in silent or heterochromatin regions. Several research groups have endeavored to define pioneer factors and study their binding characteristics using various techniques. In this review, we describe the in vitro methods used to define and characterize pioneer factors, paying particular attention to differences in methodologies and how these differences can affect results.

Nucleosomes as allosteric scaffolds for genetic regulation

Nucleosomes are stable yet highly dynamic complexes exhibiting diverse types of motions, such as sliding, DNA unwrapping, and disassembly, encoding a landscape with a large number of metastable states. In this review, describing recent studies on these nucleosome structure changes, we propose that the nucleosome can be viewed as an ideal allosteric scaffold: regulated by effector molecules such as transcription factors and chromatin remodelers, the nucleosome controls the downstream gene activity. Binding of transcription factors to the nucleosome can enhance DNA unwrapping or slide the DNA, altering either the binding or the unbinding of other transcription factors to nearby sites. ATP-dependent chromatin remodelers induce a series of DNA deformations, which allosterically propagate throughout the nucleosome to induce DNA sliding or histone exchange.

Picking a nucleosome lock: Sequence- and structure-specific recognition of the nucleosome

The nucleosome presents a formidable barrier to DNA-templated transcription by the RNA polymerase II machinery. Overcoming this transcriptional barrier in a locus-specific manner requires sequence-specific recognition of nucleosomal DNA by 'pioneer' transcription factors (TFs). Cell fate decisions, in turn, depend on the coordinated action of pioneer TFs at cell lineage-specific gene regulatory elements. Although it is already appreciated that pioneer factors play a critical role in cell differentiation, our understanding of the structural and biochemical mechanisms by which they act is still rapidly expanding. Recent research has revealed novel insight into modes of nucleosome-TF binding and uncovered kinetic principles by which nucleosomal DNA compaction affects both TF binding and residence time. Here, we review progress and argue that these structural and kinetic studies suggest new models of gene regulation by pioneer TFs.

The nucleosome position-encoding WW/SS sequence pattern is depleted in mammalian genes relative to other eukaryotes

Nucleosomal DNA sequences generally follow a well-known pattern with ∼10-bp periodic WW (where W is A or T) dinucleotides that oscillate in phase with each other and out of phase with SS (where S is G or C) dinucleotides. However, nucleosomes with other DNA patterns have not been systematically analyzed. Here, we focus on an opposite pattern, namely anti-WW/SS pattern, in which WW dinucleotides preferentially occur at DNA sites that bend into major grooves and SS (where S is G or C) dinucleotides are often found at sites that bend into minor grooves. Nucleosomes with the anti-WW/SS pattern are widespread and exhibit a species- and context-specific distribution in eukaryotic genomes. Unlike non-mammals (yeast, nematode and fly), there is a positive correlation between the enrichment of anti-WW/SS nucleosomes and RNA Pol II transcriptional levels in mammals (mouse and human). Interestingly, such enrichment is not due to underlying DNA sequence. In addition, chromatin remodeling complexes have an impact on the abundance but not on the distribution of anti-WW/SS nucleosomes in yeast. Our data reveal distinct roles of cis- and trans-acting factors in the rotational positioning of nucleosomes between non-mammals and mammals. Implications of the anti-WW/SS sequence pattern for RNA Pol II transcription are discussed.

RBPJ binds to consensus and methylated cis elements within phased nucleosomes and controls gene expression in human aortic smooth muscle cells in cooperation with SRF

Given our previous demonstration that RBPJ binds a methylated repressor element and regulates smooth muscle cell (SMC)-specific gene expression, we used genome-wide approaches to identify RBPJ binding regions in human aortic SMC and to assess RBPJ's effects on chromatin structure and gene expression. RBPJ bound to consensus cis elements, but also to TCmGGGA sequences within Alu repeats that were less transcriptionally active as assessed by DNAse hypersensitivity, H3K9 acetylation, and Notch3 and RNA Pol II binding. Interestingly, RBPJ binding was frequently detected at the borders of open chromatin, and a large fraction of genes induced or repressed by RBPJ depletion were associated with this cluster of RBPJ binding sites. RBPJ binding dramatically co-localized with serum response factor (SRF) and RNA seq experiments in RBPJ- and SRF-depleted SMC demonstrated that these factors interact functionally to regulate the contraction and inflammatory gene programs that help define SMC phenotype. Finally, we showed that RBPJ bound preferentially to phased nucleosomes independent of active chromatin marks and to cis elements positioned at the beginning and middle of the nucleosome dyad. These novel findings add important insight into RBPJ's role in chromatin structure and gene expression in SMC.

Nucleosome Crowding in Chromatin Slows the Diffusion but Can Promote Target Search of Proteins

Dynamics of nuclear proteins in crowded chromatin has only been poorly understood. Here, we address the diffusion, target search, and structural dynamics of three proteins in a model chromatin using coarse-grained molecular simulations run on the K computer. We prepared two structures of chromatin made of 20 nucleosomes with different nucleosome densities and investigated dynamics of two transcription factors, HMGB1 and p53, and one signaling protein, ERK, embedded in the chromatin. We found fast and normal diffusion of the nuclear proteins in the low-density chromatins and slow and subdiffusional movements in the high-density chromatin. The diffusion of the largest transcription factor, p53, is slowed by high-density chromatin most markedly. The on rates and off rates for DNA binding are increased and decreased, respectively, in the high-density chromatin. To our surprise, the DNA sequence search was faster in chromatin with high nucleosome density, though the diffusion is slower. We also found that the three nuclear proteins preferred to bind on the linker DNA and the entry and exit regions of nucleosomal DNA. In addition to these regions, HMGB1 and p53 also bound to the dyad.

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.

The mechanistic basis for chromatin regulation by pioneer transcription factors

Pioneer transcription factors play a primary role in establishing competence for gene expression and initiating cellular programming and reprogramming, and their dysregulation causes severe effects on human health, such as promoting tumorigenesis. Although more than 200 transcription factors are expressed in each cell type, only a small number of transcription factors are necessary to elicit dramatic cell‐fate changes in embryonic development and cell‐fate conversion. Among these key transcription factors, a subset called “pioneer transcription factors” have a remarkable ability to target nucleosomal DNA, or closed chromatin, early in development, often leading to the local opening of chromatin, thereby establishing competence for gene expression. Although more key transcription factors have been identified as pioneer transcription factors, the molecular mechanisms behind their special properties are only beginning to be revealed. Understanding the pioneering mechanisms will enhance our ability to precisely control cell fate at will for research and therapeutic purposes.

This article is categorized under:

Biological Mechanisms > Cell Fates

Biological Mechanisms > Regulatory Biology

Developmental Biology > Lineages

Mechanisms of transcriptional regulation by p53

p53 is a transcription factor that suppresses tumor growth through regulation of dozens of target genes with diverse biological functions. The activity of this master transcription factor is inactivated in nearly all tumors, either by mutations in the TP53 locus or by oncogenic events that decrease the activity of the wild-type protein, such as overexpression of the p53 repressor MDM2. However, despite decades of intensive research, our collective understanding of the p53 signaling cascade remains incomplete. In this review, we focus on recent advances in our understanding of mechanisms of p53-dependent transcriptional control as they relate to five key areas: (1) the functionally distinct N-terminal transactivation domains, (2) the diverse regulatory roles of its C-terminal domain, (3) evidence that p53 is solely a direct transcriptional activator, not a direct repressor, (4) the ability of p53 to recognize many of its enhancers across diverse chromatin environments, and (5) mechanisms that modify the p53-dependent transcriptional program in a context-dependent manner.

p53 binding sites in normal and cancer cells are characterized by distinct chromatin context

The tumor suppressor protein p53 interacts with DNA in a sequence-dependent manner. Thousands of p53 binding sites have been mapped genome-wide in normal and cancer cells. However, the way p53 selectively binds its cognate sites in different types of cells is not fully understood. Here, we performed a comprehensive analysis of 25 published p53 cistromes and identified 3,551 and 6,039 ‘high-confidence’ binding sites in normal and cancer cells, respectively. Our analysis revealed 2 distinct epigenetic features underlying p53-DNA interactions in vivo. First, p53 binding sites are associated with transcriptionally active histone marks (H3K4me3 and H3K36me3) in normal-cell chromatin, but with repressive histone marks (H3K27me3) in cancer-cell chromatin. Second, p53 binding sites in cancer cells are characterized by a lower level of DNA methylation than their counterparts in normal cells, probably related to global hypomethylation in cancers. Intriguingly, regardless of the cell type, p53 sites are highly enriched in the endogenous retroviral elements of the ERV1 family, highlighting the importance of this repeat family in shaping the transcriptional network of p53. Moreover, the p53 sites exhibit an unusual combination of chromatin patterns: high nucleosome occupancy and, at the same time, high sensitivity to DNase I. Our results suggest that p53 can access its target sites in a chromatin environment that is non-permissive to most DNA-binding transcription factors, which may allow p53 to act as a pioneer transcription factor in the context of chromatin.

One-Dimensional Search Dynamics of Tumor Suppressor p53 Regulated by a Disordered C-Terminal Domain

Tumor suppressor p53 slides along DNA and finds its target sequence in drastically different and changing cellular conditions. To elucidate how p53 maintains efficient target search at different concentrations of divalent cations such as Ca²⁺ and Mg²⁺, we prepared two mutants of p53, each possessing one of its two DNA-binding domains, the CoreTet mutant having the structured core domain plus the tetramerization (Tet) domain, and the TetCT mutant having Tet plus the disordered C-terminal domain. We investigated their equilibrium and kinetic dissociation from DNA and search dynamics along DNA at various [Mg²⁺]. Although binding of CoreTet to DNA becomes markedly weaker at higher [Mg²⁺], binding of TetCT depends slightly on [Mg²⁺]. Single-molecule fluorescence measurements revealed that the one-dimensional diffusion of CoreTet along DNA consists of fast and slow search modes, the ratio of which depends strongly on [Mg²⁺]. In contrast, diffusion of TetCT consisted of only the fast mode. The disordered C-terminal domain can associate with DNA irrespective of [Mg²⁺], and can maintain an equilibrium balance of the two search modes and the p53 search distance. These results suggest that p53 modulates the quaternary structure of the complex between p53 and DNA under different [Mg²⁺] and that it maintains the target search along DNA.

Transcriptional Regulation by Wild-Type and Cancer-Related Mutant Forms of p53

TP53 missense mutations produce a mutant p53 protein that cannot activate the p53 tumor suppressive transcriptional response, which is the primary selective pressure for TP53 mutation. Specific codons of TP53, termed hotspot mutants, are mutated at elevated frequency. Hotspot forms of mutant p53 possess oncogenic properties in addition to being deficient in tumor suppression. Such p53 mutants accumulate to high levels in the cells they inhabit, causing transcriptional alterations that produce pro-oncogenic activities, such as increased pro-growth signaling, invasiveness, and metastases. These forms of mutant p53 very likely use features of wild-type p53, such as interactions with the transcriptional machinery, to produce oncogenic effects. In this review, we discuss commonalities between wild-type and mutant p53 proteins with an emphasis on transcriptional processes.

Regulation of mammalian transcription and splicing by Nuclear RNAi

RNA interference (RNAi) is well known as a mechanism for controlling mammalian mRNA translation in the cytoplasm, but what would be the consequences if it also functions in cell nuclei? Although RNAi has also been found in nuclei of plants, yeast, and other organisms, there has been relatively little progress towards understanding the potential involvement of mammalian RNAi factors in nuclear processes including transcription and splicing. This review summarizes evidence for mammalian RNAi factors in cell nuclei and mechanisms that might contribute to the control of gene expression. When RNAi factors bind small RNAs, they form ribonucleoprotein complexes that can be selective for target sequences within different classes of nuclear RNA substrates. The versatility of nuclear RNAi may supply a previously underappreciated layer of regulation to transcription, splicing, and other nuclear processes.

TP53 engagement with the genome occurs in distinct local chromatin environments via pioneer factor activity

Despite overwhelming evidence that transcriptional activation by TP53 is critical for its tumor suppressive activity, the mechanisms by which TP53 engages the genome in the context of chromatin to activate transcription are not well understood. Using a compendium of novel and existing genome-wide data sets, we examined the relationship between TP53 binding and the dynamics of the local chromatin environment. Our analysis revealed three distinct categories of TP53 binding events that differ based on the dynamics of the local chromatin environment. The first class of TP53 binding events occurs near transcriptional start sites (TSS) and is defined by previously characterized promoter-associated chromatin modifications. The second class comprises a large cohort of preestablished, promoter-distal enhancer elements that demonstrates dynamic histone acetylation and transcription upon TP53 binding. The third class of TP53 binding sites is devoid of classic chromatin modifications and, remarkably, fall within regions of inaccessible chromatin, suggesting that TP53 has intrinsic pioneer factor activity and binds within structurally inaccessible regions of chromatin. Intriguingly, these inaccessible TP53 binding sites feature several enhancer-like properties in cell types within the epithelial lineage, indicating that TP53 binding events include a group of “proto-enhancers” that become active enhancers given the appropriate cellular context. These data indicate that TP53, along with TP63, may act as pioneer factors to specify epithelial enhancers. Further, these findings suggest that rather than following a global cell-type invariant stress response program, TP53 may tune its response based on the lineage-specific epigenomic landscape.

DNA topology confers sequence specificity to nonspecific architectural proteins

Topological constraints placed on short fragments of DNA change the disorder found in chain molecules randomly decorated by nonspecific, architectural proteins into tightly organized 3D structures. The bacterial heat-unstable (HU) protein builds up, counter to expectations, in greater quantities and at particular sites along simulated DNA minicircles and loops. Moreover, the placement of HU along loops with the "wild-type" spacing found in the Escherichia coli lactose (lac) and galactose (gal) operons precludes access to key recognition elements on DNA. The HU protein introduces a unique spatial pathway in the DNA upon closure. The many ways in which the protein induces nearly the same closed circular configuration point to the statistical advantage of its nonspecificity. The rotational settings imposed on DNA by the repressor proteins, by contrast, introduce sequential specificity in HU placement, with the nonspecific protein accumulating at particular loci on the constrained duplex. Thus, an architectural protein with no discernible DNA sequence-recognizing features becomes site-specific and potentially assumes a functional role upon loop formation. The locations of HU on the closed DNA reflect long-range mechanical correlations. The protein responds to DNA shape and deformability—the stiff, naturally straight double-helical structure—rather than to the unique features of the constituent base pairs. The structures of the simulated loops suggest that HU architecture, like nucleosomal architecture, which modulates the ability of regulatory proteins to recognize their binding sites in the context of chromatin, may influence repressor-operator interactions in the context of the bacterial nucleoid.

RBPJ, the major transcriptional effector of Notch signaling, remains associated with chromatin throughout mitosis, suggesting a role in mitotic bookmarking

Mechanisms that maintain transcriptional memory through cell division are important to maintain cell identity, and sequence-specific transcription factors that remain associated with mitotic chromatin are emerging as key players in transcriptional memory propagation. Here, we show that the major transcriptional effector of Notch signaling, RBPJ, is retained on mitotic chromatin, and that this mitotic chromatin association is mediated through the direct association of RBPJ with DNA. We further demonstrate that RBPJ binds directly to nucleosomal DNA in vitro, with a preference for sites close to the entry/exit position of the nucleosomal DNA. Genome-wide analysis in the murine embryonal-carcinoma cell line F9 revealed that roughly 60% of the sites occupied by RBPJ in asynchronous cells were also occupied in mitotic cells. Among them, we found that a fraction of RBPJ occupancy sites shifted between interphase and mitosis, suggesting that RBPJ can be retained on mitotic chromatin by sliding on DNA rather than disengaging from chromatin during mitotic chromatin condensation. We propose that RBPJ can function as a mitotic bookmark, marking genes for efficient transcriptional activation or repression upon mitotic exit. Strikingly, we found that sites of RBPJ occupancy were enriched for CTCF-binding motifs in addition to RBPJ-binding motifs, and that RBPJ and CTCF interact. Given that CTCF regulates transcription and bridges long-range chromatin interactions, our results raise the intriguing hypothesis that by collaborating with CTCF, RBPJ may participate in establishing chromatin domains and/or long-range chromatin interactions that could be propagated through cell division to maintain gene expression programs.

How does a cell remember what it should be after cell division? One mechanism that is beginning to emerge is the retention of a few key regulatory proteins on the highly condensed mitotic chromatin during cell division. These proteins are called mitotic bookmarks, as they are believed to offer critical information as to how genetic information should be read immediately after mitosis. We have found that a protein called RBPJ, which plays pivotal roles in regulating cell-fate choices, is retained on mitotic chromatin. RBPJ transmits to DNA signals elicited by the Notch pathway: a pathway that conveys information resulting from the communication between two adjacent cells. Unlike many other factors, we found that RBPJ can bind to nucleosomes, which are the basic unit of packaged DNA consisting of DNA wrapped around eight histone proteins. We also found that RBPJ interacts with and binds to DNA sites regulated by the CTCF protein, which plays important roles in regulating long-range DNA interactions. Together, our results suggest that RBPJ can function as a mitotic bookmarking factor, to help maintain genetic programs, higher-order structural information and consequently the memory of cell identity through cell division.

Rotational positioning of nucleosomes facilitates selective binding of p53 to response elements associated with cell cycle arrest

The tumor suppressor protein p53 exhibits high affinity to the response elements regulating cell cycle arrest genes (CCA-sites), but relatively low affinity to the sites associated with apoptosis (Apo-sites). This in vivo tendency cannot be explained solely by the p53-DNA binding constants measured in vitro. Since p53 can bind nucleosomal DNA, we sought to understand if the two groups of p53 sites differ in their accessibility when embedded in nucleosomes. To this aim, we analyzed the sequence-dependent bending anisotropy of human genomic DNA containing p53 sites. For the 20 CCA-sites, we calculated rotational positioning patterns predicting that most of the sites are exposed on the nucleosomal surface. This is consistent with experimentally observed positioning of human nucleosomes. Remarkably, the sequence-dependent DNA anisotropy of both the p53 sites and flanking DNA work in concert producing strong positioning signals. By contrast, both the predicted and observed rotational settings of the 38 Apo-sites in nucleosomes suggest that many of these sites are buried inside, thus preventing immediate p53 recognition and delaying gene induction. The distinct chromatin organization of the CCA response elements appears to be one of the key factors facilitating p53-DNA binding and subsequent activation of genes associated with cell cycle arrest.

Contribution of nucleosome binding preferences and co-occurring DNA sequences to transcription factor binding

BACKGROUND

Chromatin plays a critical role in regulating transcription factors (TFs) binding to their canonical transcription factor binding sites (TFBS). Recent studies in vertebrates show that many TFs preferentially bind to genomic regions that are well bound by nucleosomes in vitro. Co-occurring secondary motifs sometimes correlated with functional TFBS.

RESULTS

We used a logistic regression to evaluate how well the propensity for nucleosome binding and co-occurrence of a secondary motif identify which canonical motifs are bound in vivo. We used ChIP-seq data for three transcription factors binding to their canonical motifs: c-Jun binding the AP-1 motif (TGA(C)/(G)TCA), GR (glucocorticoid receptor) binding the GR motif (G-ACA---(T)/(C)GT-C), and Hoxa2 (homeobox a2) binding the Pbx (Pre-B-cell leukemia homeobox) motif (TGATTGAT). For all canonical TFBS in the mouse genome, we calculated intrinsic nucleosome occupancy scores (INOS) for its surrounding 150-bps DNA and examined the relationship with in vivo TF binding. In mouse mammary 3134 cells, c-Jun and GR proteins preferentially bound regions calculated to be well-bound by nucleosomes in vitro with the canonical AP-1 and GR motifs themselves contributing to the high INOS. Functional GR motifs are enriched for AP-1 motifs if they are within a nucleosome-sized 150-bps region. GR and Hoxa2 also bind motifs with low INOS, perhaps indicating a different mechanism of action.

CONCLUSION

Our analysis quantified the contribution of INOS and co-occurring sequence to the identification of functional canonical motifs in the genome. This analysis revealed an inherent competition between some TFs and nucleosomes for binding canonical TFBS. GR and c-Jun cooperate if they are within 150-bps. Binding of Hoxa2 and a fraction of GR to motifs with low INOS values suggesting they are not in competition with nucleosomes and may function using different mechanisms.

Combined micrococcal nuclease and exonuclease III digestion reveals precise positions of the nucleosome core/linker junctions: implications for high-resolution nucleosome mapping

Micrococcal nuclease (MNase) is extensively used in genome-wide mapping of nucleosomes but its preference for AT-rich DNA leads to errors in establishing precise positions of nucleosomes. Here, we show that the MNase digestion of nucleosomes assembled on a strong nucleosome positioning sequence, Widom's clone 601, releases nucleosome cores whose sizes are strongly affected by the linker DNA sequence. Our experiments produced nucleosomal DNA sizes varying between 147 and 155 bp, with positions of the MNase cuts reflecting positions of the A⋅T pairs rather than the nucleosome core/linker junctions determined by X-ray crystallography. Extent of chromatosomal DNA protection by linker histone H1 also depends on the linker DNA sequence. Remarkably, we found that a combined treatment with MNase and exonuclease III (exoIII) overcomes MNase sequence preference producing nucleosomal DNA trimmed symmetrically and precisely at the core/linker junctions regardless of the underlying DNA sequence. We propose that combined MNase/exoIII digestion can be applied to in situ chromatin for unbiased genome-wide mapping of nucleosome positions that is not influenced by DNA sequences at the core/linker junctions. The same approach can be also used for the precise mapping of the extent of linker DNA protection by H1 and other protein factors associated with nucleosome linkers.

Radioprobing the conformation of DNA in a p53-DNA complex

PURPOSE

The frequency of DNA strand breaks produced by the decay of Auger electron-emitting radionuclides is inversely proportional to the distance of DNA nucleotides from the decay site; and thus is very sensitive to changes in the local conformation of the DNA. Analysis of the frequency of DNA breaks, or radioprobing, gives valuable information about the local DNA structure. More than 10 years ago, we demonstrated the feasibility of radioprobing using a DNA-repressor complex with a known structure. Herein, we used radioprobing to study the conformation of DNA in complex with the tumor suppressor protein 53 (p53). Several structures of p53-DNA complexes have been solved by X-ray crystallography. These structures, obtained with the p53 DNA binding domain, a truncated form, laid the groundwork for understanding p53-DNA interactions and their relation to p53 functions. However, whether all observed stereochemical details are relevant to the native p53-DNA complex remains unclear. A common theme of the crystallographic structures is the lack of significant bending in the central part of the DNA response element. In contrast, gel electrophoresis and electron microscopy data showed strong DNA bending and overtwisting upon binding to the native p53 tetramer.

METHODS

To analyze DNA in complex with p53, we incorporated (125)I-dCTP in two different positions of synthetic duplexes containing the consensus p53-binding site.

RESULTS

The most significant changes in the break frequency distributions were detected close to the center of the binding site, which is consistent with an increase in DNA twisting in this region and local DNA bending and sliding.

CONCLUSIONS

Our data confirm the main results of the studies made in solution and lay a foundation for systematic examination of interactions between DNA and native p53 using (125)I radioprobing.

Binding site specificity and factor redundancy in activator protein-1-driven human papillomavirus chromatin-dependent transcription

Activator protein-1 (AP-1) regulates diverse gene responses triggered by environmental cues and virus-induced cellular stress. Although many signaling events leading to AP-1 activation have been described, the fundamental features underlying binding site selection and factor recruitment of dimeric AP-1 complexes to their target genes remain mostly uncharacterized. Using recombinant full-length human AP-1 dimers formed between c-Jun and Fos family members (c-Fos, FosB, Fra-1, Fra-2) for DNA binding and transcriptional analysis, we found that each of these AP-1 complex exhibits differential activity for distinct non-consensus AP-1 sites present in human papillomavirus (HPV), and each AP-1 complex is capable of activating transcription from in vitro-reconstituted HPV chromatin in a p300- and acetyl-CoA-dependent manner. Transcription from HPV chromatin requires AP-1-dependent and contact-driven recruitment of p300. Acetylation of dimeric AP-1 complexes by p300 enhances AP-1 binding to DNA. Using a human C-33A cervical cancer-derived cell line harboring the episomal HPV type 11 genome, we illustrate binding site selectivity recognized by c-Jun, JunB, JunD, and various Fos family members in a combinatorial and unique pattern, highlighting the diversity and importance of non-canonical binding site recognition by various AP-1 family proteins.

p53 binding to nucleosomes within the p21 promoter in vivo leads to nucleosome loss and transcriptional activation

It is well established that p53 contacts DNA in a sequence-dependent manner in order to transactivate its myriad target genes. Yet little is known about how p53 interacts with its binding site/response element (RE) within such genes in vivo in the context of nucleosomal DNA. In this study we demonstrate that both distal (5') and proximal (3') p53 REs within the promoter of the p21 gene in unstressed HCT116 colon carcinoma cells are localized within a region of relatively high nucleosome occupancy. In the absence of cellular stress, p53 is prebound to both p21 REs within nucleosomal DNA in these cells. Treatment of cells with the DNA-damaging drug doxorubicin or the p53 stabilizing agent Nutlin-3, however, is accompanied by p53-dependent subsequent loss of nucleosomes associated with such p53 REs. We show that in vitro p53 can bind to mononucleosomal DNA containing the distal p21 RE, provided the binding site is not close to the diad center of the nucleosome. In line with this, our data indicate that the p53 distal RE within the p21 gene is located close to the end of the nucleosome. Thus, low- and high-resolution mapping of nucleosome boundaries around p53 REs within the p21 promoter have provided insight into the mechanism of p53 binding to its sites in cells and the consequent changes in nucleosome occupancy at such sites.

Working the kinks out of nucleosomal DNA

Condensation of DNA in the nucleosome takes advantage of its double-helical architecture. The DNA deforms at sites where the base pairs face the histone octamer. The largest so-called kink-and-slide deformations occur in the vicinity of arginines that penetrate the minor groove. Nucleosome structures formed from the 601 positioning sequence differ subtly from those incorporating an AT-rich human α-satellite DNA. Restraints imposed by the histone arginines on the displacement of base pairs can modulate the sequence-dependent deformability of DNA and potentially contribute to the unique features of the different nucleosomes. Steric barriers mimicking constraints found in the nucleosome induce the simulated large-scale rearrangement of canonical B DNA to kink-and-slide states. The pathway to these states shows nonharmonic behavior consistent with bending profiles inferred from AFM measurements.

Impact of Alu repeats on the evolution of human p53 binding sites

BACKGROUND

The p53 tumor suppressor protein is involved in a complicated regulatory network, mediating expression of ~1000 human genes. Recent studies have shown that many p53 in vivo binding sites (BSs) reside in transposable repeats. The relationship between these BSs and functional p53 response elements (REs) remains unknown, however. We sought to understand whether the p53 REs also reside in transposable elements and particularly in the most-abundant Alu repeats.

RESULTS

We have analyzed ~160 functional p53 REs identified so far and found that 24 of them occur in repeats. More than half of these repeat-associated REs reside in Alu elements. In addition, using a position weight matrix approach, we found ~400,000 potential p53 BSs in Alu elements genome-wide. Importantly, these putative BSs are located in the same regions of Alu repeats as the functional p53 REs - namely, in the vicinity of Boxes A/A' and B of the internal RNA polymerase III promoter. Earlier nucleosome-mapping experiments showed that the Boxes A/A' and B have a different chromatin environment, which is critical for the binding of p53 to DNA. Here, we compare the Alu-residing p53 sites with the corresponding Alu consensus sequences and conclude that the p53 sites likely evolved through two different mechanisms - the sites overlapping with the Boxes A/A' were generated by CG → TG mutations; the other sites apparently pre-existed in the progenitors of several Alu subfamilies, such as AluJo and AluSq. The binding affinity of p53 to the Alu-residing sites generally correlates with the age of Alu subfamilies, so that the strongest sites are embedded in the 'relatively young' Alu repeats.

CONCLUSIONS

The primate-specific Alu repeats play an important role in shaping the p53 regulatory network in the context of chromatin. One of the selective factors responsible for the frequent occurrence of Alu repeats in introns may be related to the p53-mediated regulation of Alu transcription, which, in turn, influences expression of the host genes.

Sequence-dependent Kink-and-Slide deformations of nucleosomal DNA facilitated by histone arginines bound in the minor groove

In addition to bending and twisting deformabilities, the lateral displacements of the DNA axis (Kink-and-Slide) play an important role in DNA wrapping around the histone core (M. Y. Tolstorukov, A. V. Colasanti, D. M. McCandlish, W. K. Olson, V. B. Zhurkin, J. Mol. Biol. 371, 725-738 (2007)). Here, we show that these Kink-and-Slide deformations are likely to be stabilized by the arginine residues of histones interacting with the minor groove of DNA. The arginines are positioned asymmetrically in the minor groove, being closer to one strand. The asymmetric arginine-DNA interactions facilitate lateral displacement of base pairs across the DNA grooves, thus leading to a stepwise accumulation of the superhelical pitch of nucleosomal DNA. To understand the sequence dependence of such Kink-and-Slide deformations, we performed all-atom calculations of DNA hexamers with the YR and RY steps in the center. We found that when the unrestrained DNA deformations are allowed, the YR steps tend to bend into the major groove, and RY steps bend into the minor groove. However, when the nucleosomal Kink-and-Slide deformation is considered, the YR steps prove to be more favorable for bending into the minor groove. Overall, the Kink-and-Slide deformation energy of DNA increases in the order TA < CA < CG < GC < AC < AT. We propose a simple stereochemical model accounting for this sequence dependence. Our results agree with experimental data indicating that the TA step most frequently occurs in the minor-groove kink positions in the most stable nucleosomes. Our computations demonstrate that the Kink-and-Slide distortion is accompanied by the BI to BII transition. This fact, together with irregularities in the two-dimensional (Roll, Slide) energy contour maps, suggest that the Kink-and-Slide deformations represent a nonharmonic behavior of the duplex. This explains the difference between the two estimates of the DNA deformation energy in nucleosome - the earlier one made using knowledge-based elastic energy functions, and the current one based on all-atom calculations. Our findings are useful for refining the score functions for the prediction of nucleosome positioning. In addition, the reverse bending behavior of the YR and RY steps revealed under the Kink-and-Slide constraint is important for understanding the molecular mechanisms of binding transcription factors (such as p53) to DNA exposed on the surface of nucleosome.