A molecular mechanism for the "digital" response of p53 to stress.

The tumor suppressor protein p53 accumulates in response to cellular stress and consequently orchestrates the expression of multiple genes in a p53-level and time-dependent manner to overcome stress consequences, for which a molecular mechanism is currently unknown. Previously, we reported that DNA torsional flexibility distinguishes among p53 response elements (REs) and that transactivation at basal p53 levels is correlated with p53 REs flexibility. Here, we calculated the flexibility of ~200 p53 REs. By connecting functional outcomes of p53-target genes' activation to the calculated flexibility of their REs, we show that genes known to belong to pathways that are activated rapidly upon stress contain REs that are significantly more flexible relative to REs of genes known to be involved in pathways that are activated later in the response to stress. The global structural properties of several p53 REs belonging to different pathways were experimentally validated. Additionally, reporter-gene expression driven by flexible p53 REs occurred at lower p53 levels and with faster rates than expression from rigid REs. Furthermore, analysis of published endogenous mRNA levels of p53-target genes as a function of REs' flexibility showed that early versus late genes differ significantly in their flexibility properties of their REs and that highly flexible p53 REs enable high-activation level exclusively to early-response genes. Overall, we demonstrate that DNA flexibility of p53 REs contributes significantly to functional selectivity in the p53 system by facilitating the initial steps of p53-dependent target-genes expression, thereby contributing to survival versus death decisions in the p53 system.

Evolutionary history of the p53 family DNA-binding domain: insights from an Alvinella pompejana homolog.

The extremophile Alvinella pompejana, an annelid worm living on the edge of hydrothermal vents in the Pacific Ocean, is an excellent model system for studying factors that govern protein stability. Low intrinsic stability is a crucial factor for the susceptibility of the transcription factor p53 to inactivating mutations in human cancer. Understanding its molecular basis may facilitate the design of novel therapeutic strategies targeting mutant p53. By analyzing expressed sequence tag (EST) data, we discovered a p53 family gene in A. pompejana. Protein crystallography and biophysical studies showed that it has a p53/p63-like DNA-binding domain (DBD) that is more thermostable than all vertebrate p53 DBDs tested so far, but not as stable as that of human p63. We also identified features associated with its increased thermostability. In addition, the A. pompejana homolog shares DNA-binding properties with human p53 family DBDs, despite its evolutionary distance, consistent with a potential role in maintaining genome integrity. Through extensive structural and phylogenetic analyses, we could further trace key evolutionary events that shaped the structure, stability, and function of the p53 family DBD over time, leading to a potent but vulnerable tumor suppressor in humans.

The complex architecture of p53 binding sites.

Sequence-specific protein-DNA interactions are at the heart of the response of the tumor-suppressor p53 to numerous physiological and stress-related signals. Large variability has been previously reported in p53 binding to and transactivating from p53 response elements (REs) due, at least in part, to changes in direct (base) and indirect (shape) readouts of p53 REs. Here, we dissect p53 REs to decipher the mechanism by which p53 optimizes this highly regulated variable level of interaction with its DNA binding sites. We show that hemi-specific binding is more prevalent in p53 REs than previously envisioned. We reveal that sequences flanking the REs modulate p53 binding and activity and show that these effects extend to 4-5 bp from the REs. Moreover, we show here that the arrangement of p53 half-sites within its REs, relative to transcription direction, has been fine-tuned by selection pressure to optimize and regulate the response levels from p53 REs. This directionality in the REs arrangement is at least partly encoded in the structural properties of the REs. Furthermore, we show here that in the p21-5' RE the orientation of the half-sites is such that the effect of the flanking sequences is minimized and we discuss its advantages.

DNA mismatches reveal conformational penalties in protein-DNA recognition.

Transcription factors recognize specific genomic sequences to regulate complex gene-expression programs. Although it is well-established that transcription factors bind to specific DNA sequences using a combination of base readout and shape recognition, some fundamental aspects of protein-DNA binding remain poorly understood1,2. Many DNA-binding proteins induce changes in the structure of the DNA outside the intrinsic B-DNA envelope. However, how the energetic cost that is associated with distorting the DNA contributes to recognition has proven difficult to study, because the distorted DNA exists in low abundance in the unbound ensemble3-9. Here we use a high-throughput assay that we term SaMBA (saturation mismatch-binding assay) to investigate the role of DNA conformational penalties in transcription factor-DNA recognition. In SaMBA, mismatched base pairs are introduced to pre-induce structural distortions in the DNA that are much larger than those induced by changes in the Watson-Crick sequence. Notably, approximately 10% of mismatches increased transcription factor binding, and for each of the 22 transcription factors that were examined, at least one mismatch was found that increased the binding affinity. Mismatches also converted non-specific sites into high-affinity sites, and high-affinity sites into 'super sites' that exhibit stronger affinity than any known canonical binding site. Determination of high-resolution X-ray structures, combined with nuclear magnetic resonance measurements and structural analyses, showed that many of the DNA mismatches that increase binding induce distortions that are similar to those induced by protein binding-thus prepaying some of the energetic cost incurred from deforming the DNA. Our work indicates that conformational penalties are a major determinant of protein-DNA recognition, and reveals mechanisms by which mismatches can recruit transcription factors and thus modulate replication and repair activities in the cell10,11.

New Insights into the Role of DNA Shape on Its Recognition by p53 Proteins.

The tumor suppressor p53 acts as a transcription factor recognizing diverse DNA response elements (REs). Previous structural studies of p53-DNA complexes revealed non-canonical Hoogsteen geometry of A/T base pairs at conserved CATG motifs leading to changes in DNA shape and its interface with p53. To study the effects of DNA shape on binding characteristics, we designed REs with modified base pairs "locked" into either Hoogsteen or Watson-Crick form. Here we present crystal structures of these complexes and their thermodynamic and kinetic parameters, demonstrating that complexes with Hoogsteen base pairs are stabilized relative to those with all-Watson-Crick base pairs. CATG motifs are abundant in p53REs such as GADD45 and p53R2 related to cell-cycle arrest and DNA repair. The high-resolution structures of these complexes validate their propensity to adopt the unique Hoogsteen-induced structure, thus providing insights into the functional role of DNA shape and broadening the mechanisms that contribute to DNA recognition by proteins.

How mutations shape p53 interactions with the genome to promote tumorigenesis and drug resistance.

The tumor suppressive transcription factor p53 regulates a wide array of cellular processes that confer upon cells an essential protection against cancer development. Wild-type p53 regulates gene expression by directly binding to DNA in a sequence-specific manner. p53 missense mutations are the most common mutations in malignant cells and can be regarded as synonymous with anticancer drug resistance and poor prognosis. The current review provides an overview of how the extraordinary variety of more than 2000 different mutant p53 proteins, known as the p53 mutome, affect the interaction of p53 with DNA. We discuss how the classification of p53 mutations to loss of function (LOF), gain of function (GOF), and dominant-negative (DN) inhibition of a remaining wild-type allele, hides a complex p53 mutation spectrum that depends on the distinctive nature of each mutant protein, requiring different therapeutic strategies for each mutant p53 protein. We propose to regard the different mutant p53 categories as continuous variables, that may not be independent of each other. In particular, we suggest here to consider GOF mutations as a special subset of LOF mutations, especially when mutant p53 binds to DNA through cooperation with other transcription factors, and we present a model for GOF mechanism that consolidates many observations on the GOF phenomenon. We review how novel mutant p53 targeting approaches aim to restore a wild-type-like DNA interaction and to overcome resistance to cancer therapy.

Diverse p53/DNA binding modes expand the repertoire of p53 response elements.

The tumor suppressor protein p53 acts as a transcription factor, binding sequence-specifically to defined DNA sites, thereby activating the expression of genes leading to diverse cellular outcomes. Canonical p53 response elements (REs) are made of two decameric half-sites separated by a variable number of base pairs (spacers). Fifty percent of all validated p53 REs contain spacers between 1 and 18 bp; however, their functional significance is unclear at present. Here, we show that p53 forms two different tetrameric complexes with consensus or natural REs, both with long spacers: a fully specific complex where two p53 dimers bind to two specific half-sites, and a hemispecific complex where one dimer binds to a specific half-site and the second binds to an adjacent spacer sequence. The two types of complexes have comparable binding affinity and specificity, as judged from binding competition against bulk genomic DNA. Structural analysis of the p53 REs in solution shows that these sites are not bent in both their free and p53-bound states when the two half-sites are either abutting or separated by spacers. Cell-based assay supports the physiological relevance of our findings. We propose that p53 REs with long spacers comprise separate specific half-sites that can lead to several different tetrameric complexes. This finding expands the universe of p53 binding sites and demonstrates that even isolated p53 half-sites can form tetrameric complexes. Moreover, it explains the manner in which p53 binds to clusters of more than one canonical binding site, common in many natural REs.

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.

Sequence-dependent cooperative binding of p53 to DNA targets and its relationship to the structural properties of the DNA targets.

The prime mechanism by which p53 acts as a tumor suppressor is as a transcription factor regulating the expression of diverse downstream genes. The DNA-binding domain of p53 (p53DBD) interacts with defined DNA sites and is the main target for mutations in human primary tumors. Here, we show that the CWWG motif, found in the center of each consensus p53 half-site, is a key player in p53/DNA interactions. Gel-mobility-shift assays provide a unique opportunity to directly observe the various oligomeric complexes formed between p53DBD and its target sites. We demonstrate that p53DBD binds to p53 consensus sites containing CATG with relatively low cooperativity, as both dimers and tetramers, and with even lower cooperativity to such sites containing spacer sequences. p53DBD binds to sites containing CAAG and CTAG with measurable affinity only when imbedded in two contiguous p53 half-sites and only as tetramers (with very high cooperativity). There are three orders-of-magnitude difference in the cooperativity of interaction between sites differing in their non-contacted step, and further two orders-of-magnitude difference as a function of spacer sequences. By experimentally measuring the global structural properties of these sites, by cyclization kinetics of DNA minicircles, we correlate these differences with the torsional flexibility of the binding sites.

The coexistence of the nucleosome positioning code with the genetic code on eukaryotic genomes.

It is known that there are several codes residing simultaneously on the DNA double helix. The two best-characterized codes are the genetic code--the code for protein production, and the code for DNA packaging into nucleosomes. Since these codes have to coexist simultaneously on the same DNA region, both must be degenerate to allow this coexistence. A-tracts are homopolymeric stretches of several adjacent deoxyadenosines on one strand of the double helix, having unusual structural properties, which were shown to exclude nucleosomes and as such are instrumental in setting the translational positioning of DNA within nucleosomes. We observe, cross-kingdoms, a strong codon bias toward the avoidance of long A-tracts in exon regions, which enables the formation of high density of nucleosomes in these regions. Moreover, long A-tract avoidance is restricted exclusively to nucleosome-occupied exon regions. We show that this bias in codon usage is sufficient for enabling DNA organization within nucleosomes without constraints on the actual code for proteins. Thus, there is inter-dependency of the two major codes within DNA to allow their coexistence. Furthermore, we show that modulation of A-tract occurrences in exon versus non-exon regions may result in a unique alternation of the diameter of the '30-nm' fiber model.

The unique structure of A-tracts and intrinsic DNA bending.

Short runs of adenines are a ubiquitous DNA element in regulatory regions of many organisms. When runs of 4-6 adenine base pairs ('A-tracts') are repeated with the helical periodicity, they give rise to global curvature of the DNA double helix, which can be macroscopically characterized by anomalously slow migration on polyacrylamide gels. The molecular structure of these DNA tracts is unusual and distinct from that of canonical B-DNA. We review here our current knowledge about the molecular details of A-tract structure and its interaction with sequences flanking them of either side and with the environment. Various molecular models were proposed to describe A-tract structure and how it causes global deflection of the DNA helical axis. We review old and recent findings that enable us to amalgamate the various findings to one model that conforms to the experimental data. Sequences containing phased repeats of A-tracts have from the very beginning been synonymous with global intrinsic DNA bending. In this review, we show that very often it is the unique structure of A-tracts that is at the basis of their widespread occurrence in regulatory regions of many organisms. Thus, the biological importance of A-tracts may often be residing in their distinct structure rather than in the global curvature that they induce on sequences containing them.

Nearest-neighbor non-additivity versus long-range non-additivity in TATA-box structure and its implications for TBP-binding mechanism.

TBP recognizes its target sites, TATA boxes, by recognizing their sequence-dependent structure and flexibility. Studying this mode of TATA-box recognition, termed 'indirect readout', is important for elucidating the binding mechanism in this system, as well as for developing methods to locate new binding sites in genomic DNA. We determined the binding stability and TBP-induced TATA-box bending for consensus-like TATA boxes. In addition, we calculated the individual information score of all studied sequences. We show that various non-additive effects exist in TATA boxes, dependent on their structural properties. By several criterions, we divide TATA boxes to two main groups. The first group contains sequences with 3-4 consecutive adenines. Sequences in this group have a rigid context-independent cooperative structure, best described by a nearest-neighbor non-additive model. Sequences in the second group have a flexible, context-dependent conformation, which cannot be described by an additive model or by a nearest-neighbor non-additive model. Classifying TATA boxes by these and other structural rules clarifies the different recognition pathways and binding mechanisms used by TBP upon binding to different TATA boxes. We discuss the structural and evolutionary sources of the difficulties in predicting new binding sites by probabilistic weight-matrix methods for proteins in which indirect readout is dominant.

Structural basis of DNA recognition by p53 tetramers.

The tumor-suppressor protein p53 is among the most effective of the cell's natural defenses against cancer. In response to cellular stress, p53 binds as a tetramer to diverse DNA targets containing two decameric half-sites, thereby activating the expression of genes involved in cell-cycle arrest or apoptosis. Here we present high-resolution crystal structures of sequence-specific complexes between the core domain of human p53 and different DNA half-sites. In all structures, four p53 molecules self-assemble on two DNA half-sites to form a tetramer that is a dimer of dimers, stabilized by protein-protein and base-stacking interactions. The protein-DNA interface varies as a function of the specific base sequence in correlation with the measured binding affinities of the complexes. The new data establish a structural framework for understanding the mechanisms of specificity, affinity, and cooperativity of DNA binding by p53 and suggest a model for its regulation by regions outside the sequence-specific DNA binding domain.

TBP flanking sequences: asymmetry of binding, long-range effects and consensus sequences.

We carried out in vitro selection experiments to systematically probe the effects of TATA-box flanking sequences on its interaction with the TATA-box binding protein (TBP). This study validates our previous hypothesis that the effect of the flanking sequences on TBP/TATA-box interactions is much more significant when the TATA box has a context-dependent DNA structure. Several interesting observations, with implications for protein-DNA interactions in general, came out of this study. (i) Selected sequences are selection-method specific and TATA-box dependent. (ii) The variability in binding stability as a function of the flanking sequences for (T-A)4 boxes is as large as the variability in binding stability as a function of the core TATA box itself. Thus, for (T-A)4 boxes the flanking sequences completely dominate and determine the binding interaction. (iii) Binding stabilities of all but one of the individual selected sequences of the (T-A)4 form is significantly higher than that of their mononucleotide-based consensus sequence. (iv) Even though the (T-A)4 sequence is symmetric the flanking sequence pattern is asymmetric. We propose that the plasticity of (T-A)n sequences increases the number of conformationally distinct TATA boxes without the need to extent the TBP contact region beyond the eight-base-pair long TATA box.

Characteristics of migration patterns of DNA oligomers in gels and the relationship to the question of intrinsic DNA bending.

We have developed a methodology that is capable of quantitatively describing the electrophoretic mobility patterns of oligomeric B-DNA through polyacrylamide gels (PAG) in the presence of varying concentration of the organic solvent 2-methyl-2,4-pentanediol (MPD), used routinely to induce DNA crystallization. The model includes the ion atmosphere and its polarization, electrostatic excluded volume, hydrodynamic interactions, and fluctuation effects that characterize the overall size of the migrating polyion. Using this model, and by critically examining the mobility patterns of linear random-sequence B-DNA molecules in PAG as a function of MPD, we address the question of the discrepancy between current models used to explain the molecular origins of A-tract-induced DNA bending. Direct analysis of the mobility of B-DNA oligomers on PAG, and comparison to the mobility of A-tract-containing oligomers, shows a significant apparent effect of MPD on the mobility of generic B-DNA sequences, which is larger than the effect on A-tract-containing oligomers. The effect is chain-length dependent, especially at lower MPD concentration. Thus, the apparent reduction in gel mobility, as a function of MPD, is not unique to A-tract regions or A-tract-containing molecules. However, our analysis suggests that MPD molecules are probably excluded from the surface of both B-DNA and A-tract molecules. This is supported by circular dichroism studies on A-tract and B-DNA molecules in solutions containing various MPD concentrations.

A-tract polarity dominate the curvature in flanking sequences.

It is well-known, but little understood, that the nucleotide sequences between phased A(4-6)-tracts (at 10-11 bp intervals) have only a slight effect on overall curvature. To explore this phenomenon, we have examined the gel-migration properties of sequences containing both A-tracts as well as G-tracts (i.e., sequences of the form G(n)C(m) or C(n)G(m), n + m > 4) in various relative positioning. We show that the composite bend of these sequences depends on their relative arrangement. When G-tracts are placed between two A-tracts, such that both tracts are repeated in phase to themselves (e.g., G(5)A(6)G(5)A(5)), or adjacent to the 3'-side of A-tracts (e.g., A(6)G(5)N(10)), they have minimal influence on the extent of bending of the composite sequence. When G-tracts are placed one helical repeat away from A-tracts (e.g., G(5)N(5)A(6)N(6)), or are adjacent only to the 5'-side of A-tracts (e.g., G(5)A(6)N(10)) their influence on the composite bend is larger. The differential behavior of AG- versus GA-tracts means that A-tracts influence their flanking sequences in a polar manner. Whereas they suppress, or make constant, the intrinsic bending characteristics of any sequence placed immediately 3' to them (and hence by definition any sequence placed between two phased A-tracts), sequences adjoining them on their 5'-side are free to modulate the overall curvature. We interpret these results as evidence for the dominant nature of the unique and nonuniform structure adopted by tracts of four adenines or more. The effects of A-tracts extend at least five base pairs into the adjoining 3' region. This is further evidence for the complexity of DNA structure and the inadequacy of simple nearest-neighbor models to explain all its manifestations.

Dynamics of curved DNA molecules: prediction and experiment.

We have developed a quantitative predictive model capable of describing the dynamics of migration of intrinsically curved DNA fragments on polyacrylamide gels. The model takes into account structural features of DNA, end-to-end distance, screening of hydrodynamic interactions, ionic strength of buffer, electrostatic persistence length, structural fluctuations of the macromolecule, counter condensation, and variation of dielectric constant and viscosity of water with MPD. In doing so, we have also addressed a decade old issue on the effect of the organic solvent 2-methyl-2,4-pentanediol on gel migration of phased A-tracts. We show here that A-tract-solvent interactions are less favored compared with A-tract-A-tract and solvent-solvent interactions.

DNA binding and 3'-5' exonuclease activity in the murine alternatively-spliced p53 protein.

In this study we show that the naturally occurring C-terminally alternative spliced p53 (referred to as AS-p53) is active as a sequence-specific DNA binding protein as well as a 3'-5'-exonuclease in the presence of Mg2+ ions. The two activities are positively correlated as the sequence-specific DNA target is more efficiently degraded than a non-specific target. In contrast, a mutated AS-p53 protein that is deficient in DNA binding lacks exonuclease activity. The use of modified p53 binding sites, where the 3'-phosphate is replaced by a phosphorothioate group, enabled the inhibition of DNA degradation under the binding conditions. We demonstrate that AS-p53 interacts with its specific DNA target by two distinct binding modes: a high-affinity mode characterized by a low-mobility protein-DNA complex at the nanomolar range, and a low-affinity mode shown by a high-mobility complex at the micromolar range. Comparison of the data on the natural and the modified p53 binding sites suggests that the high-affinity mode is related to AS-p53 function as a transcription factor and that the low-affinity mode is associated with its exonuclease activity. The implications of these findings to a specific cellular role of AS-p53 are discussed.