Crystal structure of a multidomain human p53 tetramer bound to the natural CDKN1A (p21) p53-response element

Germline variant affecting p53β isoforms predisposes to familial cancer

Germline and somatic TP53 variants play a crucial role during tumorigenesis. However, genetic variations that solely affect the alternatively spliced p53 isoforms, p53β and p53γ, are not fully considered in the molecular diagnosis of Li-Fraumeni syndrome and cancer. In our search for additional cancer predisposing variants, we identify a heterozygous stop-lost variant affecting the p53β isoforms (p.*342Serext*17) in four families suspected of an autosomal dominant cancer syndrome with colorectal, breast and papillary thyroid cancers. The stop-lost variant leads to the 17 amino-acid extension of the p53β isoforms, which increases oligomerization to canonical p53α and dysregulates the expression of p53's transcriptional targets. Our study reveals the capacity of p53β mutants to influence p53 signalling and contribute to the susceptibility of different cancer types. These findings underscore the significance of p53 isoforms and the necessity of comprehensive investigation into the entire TP53 gene in understanding cancer predisposition.

Distinct sets of molecular characteristics define tumor-rejecting neoantigens

Challenges in identifying tumor-rejecting neoantigens limit the efficacy of neoantigen vaccines to treat cancers, including cutaneous squamous cell carcinoma (cSCC). A minority of human cSCC tumors shared neoantigens, supporting the need for personalized vaccines. Using a UV-induced mouse cSCC model which recapitulated the mutational signature and driver mutations found in human disease, we found that CD8 T cells constrain cSCC. Two MHC class I neoantigens were identified that constrained cSCC growth. Compared to the wild-type peptides, one tumor-rejecting neoantigen exhibited improved MHC binding and the other had increased solvent accessibility of the mutated residue. Across known neoantigens that do not impact MHC binding, structural modeling of the peptide/MHC complexes indicated that increased solvent accessibility, which will facilitate TCR recognition of the neoantigen, distinguished tumor-rejecting from non-immunogenic neoantigens. This work reveals characteristics of tumor-rejecting neoantigens that may be of considerable importance in identifying optimal vaccine candidates in cSCC and other cancers.

Structural assessment of the full-length wild-type tumor suppressor protein p53 by mass spectrometry-guided computational modeling

The tetrameric tumor suppressor p53 represents a great challenge for 3D-structural analysis due to its high degree of intrinsic disorder (ca. 40%). We aim to shed light on the structural and functional roles of p53's C-terminal region in full-length, wild-type human p53 tetramer and their importance for DNA binding. For this, we employed complementary techniques of structural mass spectrometry (MS) in an integrated approach with computational modeling. Our results show no major conformational differences in p53 between DNA-bound and DNA-free states, but reveal a substantial compaction of p53's C-terminal region. This supports the proposed mechanism of unspecific DNA binding to the C-terminal region of p53 prior to transcription initiation by specific DNA binding to the core domain of p53. The synergies between complementary structural MS techniques and computational modeling as pursued in our integrative approach is envisioned to serve as general strategy for studying intrinsically disordered proteins (IDPs) and intrinsically disordered region (IDRs).

Mapping Interactions of the Intrinsically Disordered C-Terminal Regions of Tetrameric p53 by Segmental Isotope Labeling and NMR

The C-terminal region of the tumor suppressor protein p53 contains three domains, nuclear localization signal (NLS), tetramerization domain (TET), and C-terminal regulatory domain (CTD), which are essential for p53 function. Characterization of the structure and interactions of these domains within full-length p53 has been limited by the overall size and flexibility of the p53 tetramer. Using trans-intein splicing, we have generated full-length p53 constructs in which the C-terminal region is isotopically labeled with 15N for NMR analysis, allowing us to obtain atomic-level information on the C-terminal domains in the context of the full-length protein. Resonances of NLS and CTD residues have narrow linewidths, showing that these regions are largely solvent-exposed and dynamically disordered, whereas resonances from the folded TET are broadened beyond detection. Two regions of the CTD, spanning residues 369-374 and 381-388 and with high lysine content, make dynamic and sequence-independent interactions with DNA in regions that flank the p53 recognition element. The population of DNA-bound states increases as the length of the flanking regions is extended up to approximately 20 base pairs on either side of the recognition element. Acetylation of K372, K373, and K382, using a construct of the transcriptional coactivator CBP containing the TAZ2 and acetyltransferase domains, inhibits interaction of the CTD with DNA. This work provides high-resolution insights into the behavior of the intrinsically disordered C-terminal regions of p53 within the full-length tetramer and the molecular basis by which the CTD mediates DNA binding and specificity.

Insights into Allosteric Mechanisms of the Lung-Enriched p53 Mutants V157F and R158L

Lung cancer is a leading fatal malignancy in humans. p53 mutants exhibit not only loss of tumor suppressor capability but also oncogenic gain-of-function, contributing to lung cancer initiation, progression and therapeutic resistance. Research shows that p53 mutants V157F and R158L occur with high frequency in lung squamous cell carcinomas. Revealing their conformational dynamics is critical for developing novel lung therapies. Here, we used all-atom molecular dynamics (MD) simulations to investigate the effect of V157F and R158L substitutions on the structural properties of the p53 core domain (p53C). Compared to wild-type (WT) p53C, both V157F and R158L mutants display slightly lesser β-sheet structure, larger radius of gyration, larger volume and larger exposed surface area, showing aggregation-prone structural characteristics. The aggregation-prone fragments (residues 249–267 and 268–282) of two mutants are more exposed to water solution than that of WT p53C. V157F and R158L mutation sites can affect the conformation switch of loop 1 through long-range associations. Simulations also reveal that the local structure and conformation around the V157F and R158L mutation sites are in a dynamic equilibrium between the misfolded and properly folded conformations. These results provide molecular mechanistic insights into allosteric mechanisms of the lung-enriched p53 mutants.

Structural Basis of Mutation-Dependent p53 Tetramerization Deficiency

The formation of a tetrameric assembly is essential for the ability of the tumor suppressor protein p53 to act as a transcription factor. Such a quaternary conformation is driven by a specific tetramerization domain, separated from the central DNA-binding domain by a flexible linker. Despite the distance, functional crosstalk between the two domains has been reported. This phenomenon can explain the pathogenicity of some inherited or somatically acquired mutations in the tetramerization domain, including the widespread R337H missense mutation present in the population in south Brazil. In this work, we combined computational predictions through extended all-atom molecular dynamics simulations with functional assays in a genetically defined yeast-based model system to reveal structural features of p53 tetramerization domains and their transactivation capacity and specificity. In addition to the germline and cancer-associated R337H and R337C, other rationally designed missense mutations targeting a significant salt-bridge interaction that stabilizes the p53 tetramerization domain were studied (i.e., R337D, D352R, and the double-mutation R337D plus D352R). The simulations revealed a destabilizing effect of the pathogenic mutations within the p53 tetramerization domain and highlighted the importance of electrostatic interactions between residues 337 and 352. The transactivation assay, performed in yeast by tuning the expression of wild-type and mutant p53 proteins, revealed that p53 tetramerization mutations could decrease the transactivation potential and alter transactivation specificity, in particular by better tolerating negative features in weak DNA-binding sites. These results establish the effect of naturally occurring variations at positions 337 and 352 on p53’s conformational stability and function.

Structural diversity of p63 and p73 isoforms

Abstract

The p53 protein family is the most studied protein family of all. Sequence analysis and structure determination have revealed a high similarity of crucial domains between p53, p63 and p73. Functional studies, however, have shown a wide variety of different tasks in tumor suppression, quality control and development. Here we review the structure and organization of the individual domains of p63 and p73, the interaction of these domains in the context of full-length proteins and discuss the evolutionary origin of this protein family.

Facts

Open questions

Protein of a thousand faces: The tumor-suppressive and oncogenic responses of p53

The p53 protein is a pleiotropic regulator working as a tumor suppressor and as an oncogene. Depending on the cellular insult and the mutational status, p53 may trigger opposing activities such as cell death or survival, senescence and cell cycle arrest or proliferative signals, antioxidant or prooxidant activation, glycolysis, or oxidative phosphorylation, among others. By augmenting or repressing specific target genes or directly interacting with cellular partners, p53 accomplishes a particular set of activities. The mechanism in which p53 is activated depends on increased stability through post-translational modifications (PTMs) and the formation of higher-order structures (HOS). The intricate cell death and metabolic p53 response are reviewed in light of gaining stability via PTM and HOS formation in health and disease.

Association between Predicted Effects of TP53 Missense Variants on Protein Conformation and Their Phenotypic Presentation as Li-Fraumeni Syndrome or Hereditary Breast Cancer

Rare germline pathogenic TP53 missense variants often predispose to a wide spectrum of tumors characterized by Li-Fraumeni syndrome (LFS) but a subset of variants is also seen in families with exclusively hereditary breast cancer (HBC) outcomes. We have developed a logistic regression model with the aim of predicting LFS and HBC outcomes, based on the predicted effects of individual TP53 variants on aspects of protein conformation. A total of 48 missense variants either unique for LFS (n = 24) or exclusively reported in HBC (n = 24) were included. LFS-variants were over-represented in residues tending to be buried in the core of the tertiary structure of TP53 (p = 0.0014). The favored logistic regression model describes disease outcome in terms of explanatory variables related to the surface or buried status of residues as well as their propensity to contribute to protein compactness or protein-protein interactions. Reduced, internally validated models discriminated well between LFS and HBC (C-statistic = 0.78−0.84; equivalent to the area under the ROC (receiver operating characteristic) curve), had a low risk for over-fitting and were well calibrated in relation to the known outcome risk. In conclusion, this study presents a phenotypic prediction model of LFS and HBC risk for germline TP53 missense variants, in an attempt to provide a complementary tool for future decision making and clinical handling.

Structure of the p53/RNA polymerase II assembly

The tumor suppressor p53 protein activates expression of a vast gene network in response to stress stimuli for cellular integrity. The molecular mechanism underlying how p53 targets RNA polymerase II (Pol II) to regulate transcription remains unclear. To elucidate the p53/Pol II interaction, we have determined a 4.6 Å resolution structure of the human p53/Pol II assembly via single particle cryo-electron microscopy. Our structure reveals that p53's DNA binding domain targets the upstream DNA binding site within Pol II. This association introduces conformational changes of the Pol II clamp into a further-closed state. A cavity was identified between p53 and Pol II that could possibly host DNA. The transactivation domain of p53 binds the surface of Pol II's jaw that contacts downstream DNA. These findings suggest that p53's functional domains directly regulate DNA binding activity of Pol II to mediate transcription, thereby providing insights into p53-regulated gene expression.

An integrated view of p53 dynamics, function, and reactivation

The tumor suppressor p53 plays a vital role in responding to cell stressors such as DNA damage, hypoxia, and tumor formation by inducing cell-cycle arrest, senescence, or apoptosis. Expression level alterations and mutational frequency implicates p53 in most human cancers. In this review, we show how both computational and experimental methods have been used to provide an integrated view of p53 dynamics, function, and reactivation potential. We argue that p53 serves as an exceptional case study for developing methods in modeling intrinsically disordered proteins. We describe how these methods can be leveraged to improve p53 reactivation molecule design and other novel therapeutic modalities, such as PROteolysis TARgeting Chimeras (PROTACs).

Markov state models and NMR uncover an overlooked allosteric loop in p53

The tumor suppressor p53 is the most frequently mutated gene in human cancer, and thus reactivation of mutated p53 is a promising avenue for cancer therapy. Analysis of wildtype p53 and the Y220C cancer mutant long-timescale molecular dynamics simulations with Markov state models and validation by NMR relaxation studies has uncovered the involvement of loop L6 in the slowest motions of the protein. Due to its distant location from the DNA-binding surface, the conformational dynamics of this loop has so far remained largely unexplored. We observe mutation-induced stabilization of alternate L6 conformations, distinct from all experimentally-determined structures, in which the loop is both extended and located further away from the DNA-interacting surface. Additionally, the effect of the L6-adjacent Y220C mutation on the conformational landscape of the functionally-important loop L1 suggests an allosteric role to this dynamic loop and the inactivation mechanism of the mutation. Finally, the simulations reveal a novel Y220C cryptic pocket that can be targeted for p53 rescue efforts. Our approach exemplifies the power of the MSM methodology for uncovering intrinsic dynamic and kinetic differences among distinct protein ensembles, such as for the investigation of mutation effects on protein function.

Wildtype and Y220C L1 and L6 loops conformational landscape, with MSM-identified L6 states highlighted on the right.

Liquid-like droplet formation by tumor suppressor p53 induced by multivalent electrostatic interactions between two disordered domains

Early in vivo studies demonstrated the involvement of a tumor-suppressing transcription factor, p53, into cellular droplets such as Cajal and promyelocytic leukemia protein bodies, suggesting that the liquid-liquid phase separation (LLPS) might be involved in the cellular functions of p53. To examine this possibility, we conducted extensive investigations on the droplet formation of p53 in vitro. First, p53 itself was found to form liquid-like droplets at neutral and slightly acidic pH and at low salt concentrations. Truncated p53 mutants modulated droplet formation, suggesting the importance of multivalent electrostatic interactions among the N-terminal and C-terminal domains. Second, FRET efficiency measurements for the dimer mutants of p53 revealed that distances between the core domains and between the C-terminal domains were modulated in an opposite manner within the droplets. Third, the molecular crowding agents were found to promote droplet formation, whereas ssDNA, dsDNA, and ATP, to suppress it. Finally, the p53 mutant mimicking posttranslational phosphorylation did not form the droplets. We conclude that p53 itself has a potential to form droplets that can be controlled by cellular molecules and by posttranslational modifications, suggesting that LLPS might be involved in p53 function.

The R882H substitution in the human de novo DNA methyltransferase DNMT3A disrupts allosteric regulation by the tumor supressor p53

A myriad of protein partners modulate the activity of the human DNA methyltransferase 3A (DNMT3A), whose interactions with these other proteins are frequently altered during oncogenesis. We show here that the tumor suppressor p53 decreases DNMT3A activity by forming a heterotetramer complex with DNMT3A. Mutational and modeling experiments suggested that p53 interacts with the same region in DNMT3A as does the structurally characterized DNMT3L. We observed that the p53-mediated repression of DNMT3A activity is blocked by amino acid substitutions within this interface, but surprisingly, also by a distal DNMT3A residue, R882H. DNMT3A R882H occurs frequently in various cancers, including acute myeloid leukemia, and our results suggest that the effects of R882H and other DNMT3A mutations may go beyond changes in DNMT3A methylation activity. To further understand the dynamics of how protein-protein interactions modulate DNMT3A activity, we determined that p53 has a greater affinity for DNMT3A than for DNMT3L and that p53 readily displaces DNMT3L from the DNMT3A:DNMT3L heterotetramer. Interestingly, this occurred even when the preformed DNMT3A:DNMT3L complex was actively methylating DNA. The frequently identified p53 substitutions (R248W and R273H), whereas able to regulate DNMT3A function when forming the DNMT3A:p53 heterotetramer, no longer displaced DNMT3L from the DNMT3A:DNMT3L heterotetramer. The results of our work highlight the complex interplay between DNMT3A, p53, and DNMT3L and how these interactions are further modulated by clinically derived mutations in each of the interacting partners.

iASPP mediates p53 selectivity through a modular mechanism fine-tuning DNA recognition

The most frequently mutated protein in human cancer is p53, a transcription factor (TF) that regulates myriad genes instrumental in diverse cellular outcomes including growth arrest and cell death. Cell context-dependent p53 modulation is critical for this life-or-death balance, yet remains incompletely understood. Here we identify sequence signatures enriched in genomic p53-binding sites modulated by the transcription cofactor iASPP. Moreover, our p53-iASPP crystal structure reveals that iASPP displaces the p53 L1 loop-which mediates sequence-specific interactions with the signature-corresponding base-without perturbing other DNA-recognizing modules of the p53 DNA-binding domain. A TF commonly uses multiple structural modules to recognize its cognate DNA, and thus this mechanism of a cofactor fine-tuning TF-DNA interactions through targeting a particular module is likely widespread. Previously, all tumor suppressors and oncoproteins that associate with the p53 DNA-binding domain-except the oncogenic E6 from human papillomaviruses (HPVs)-structurally cluster at the DNA-binding site of p53, complicating drug design. By contrast, iASPP inhibits p53 through a distinct surface overlapping the E6 footprint, opening prospects for p53-targeting precision medicine to improve cancer therapy.

Glucocorticoid receptor quaternary structure drives chromatin occupancy and transcriptional outcome

Most transcription factors, including nuclear receptors, are widely modeled as binding regulatory elements as monomers, homodimers, or heterodimers. Recent findings in live cells show that the glucocorticoid receptor NR3C1 (also known as GR) forms tetramers on enhancers, owing to an allosteric alteration induced by DNA binding, and suggest that higher oligomerization states are important for the gene regulatory responses of GR. By using a variant (GRtetra) that mimics this allosteric transition, we performed genome-wide studies using a GR knockout cell line with reintroduced wild-type GR or reintroduced GRtetra. GRtetra acts as a super receptor by binding to response elements not accessible to the wild-type receptor and both induces and represses more genes than GRwt. These results argue that DNA binding induces a structural transition to the tetrameric state, forming a transient higher-order structure that drives both the activating and repressive actions of glucocorticoids.

Dynamics and Molecular Mechanisms of p53 Transcriptional Activation

The "guardian of the genome", p53, functions as a tumor suppressor that responds to cell stressors such as DNA damage, hypoxia, and tumor formation by inducing cell-cycle arrest, senescence, or apoptosis. Mutation of p53 disrupts its tumor suppressor function, leading to various types of human cancers. One particular mutant, R175H, is a structural mutant that inactivates the DNA damage response pathway and acquires oncogenic functions that promotes both cancer and drug resistance. Our current work aims to understand how p53 wild-type function is disrupted due to the R175H mutation. We use a series of atomistic integrative models built previously from crystal structures of the full-length p53 tetramer bound to DNA and model the R175H mutant using in silico site-directed mutagenesis. Explicitly solvated all-atom molecular dynamics (MD) simulations on wild-type and the R175H mutant p53 reveal insights into how wild-type p53 searches and recognizes DNA, and how this mechanism is disrupted as a result of the R175H mutation. Specifically, our work reveals the optimal quaternary DNA binding mode of the DNA binding domain and shows how this binding mode is altered via symmetry loss as a result of the R175H mutation, indicating a recognition mechanism that is reminiscent of the asymmetry seen in wild type p53 binding to nonspecific genomic elements. Altogether our work sheds new light into the hitherto unseen molecular mechanisms governing transcription factor, DNA recognition.

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.

Biophysical simulations and structure-based modeling of residue interaction networks in the tumor suppressor proteins reveal functional role of cancer mutation hotspots in molecular communication

In the current study, we have combined molecular simulations and energetic analysis with dynamics-based network modeling and perturbation response scanning to determine molecular signatures of mutational hotspot residues in the p53, PTEN, and SMAD4 tumor suppressor proteins. By examining structure, energetics and dynamics of these proteins, we have shown that inactivating mutations preferentially target a group of structurally stable residues that play a fundamental role in global propagation of dynamic fluctuations and mediating allosteric interaction networks. Through integration of long-range perturbation dynamics and network-based approaches, we have quantified allosteric potential of residues in the studied proteins. The results have revealed that mutational hotspot sites often correspond to high centrality mediating centers of the residue interaction networks that are responsible for coordination of global dynamic changes and allosteric signaling. Our findings have also suggested that structurally stable mutational hotpots can act as major effectors of allosteric interactions and mutations in these positions are typically associated with severe phenotype. Modeling of shortest inter-residue pathways has shown that mutational hotspot sites can also serve as key mediating bridges of allosteric communication in the p53 and PTEN protein structures. Multiple regression models have indicated that functional significance of mutational hotspots can be strongly associated with the network signatures serving as robust predictors of critical regulatory positions responsible for loss-of-function phenotype. The results of this computational investigation are compared with the experimental studies and reveal molecular signatures of mutational hotspots, providing a plausible rationale for explaining and localizing disease-causing mutations in tumor suppressor genes.

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.

A new era of studying p53-mediated transcription activation

To prevent tumorigenesis, p53 stimulates transcription by facilitating the recruitment of the transcription machinery on target gene promoters. Cryo-Electron Microscopy studies on p53-bound RNA Polymerase II (Pol II) reveal that p53 structurally regulates Pol II to affect its DNA binding and elongation, providing new insights into p53-mediated transcriptional regulation.

Structural visualization of the p53/RNA polymerase II assembly

Singh et al. dissected the human p53/Pol II interaction via single-particle cryo-electron microscopy, structural docking, and biochemical analyses. These findings indicate that p53 may structurally regulate DNA-binding functions of Pol II via the clamp domain, thereby providing insights into p53-regulated Pol II transcription.

The master tumor suppressor p53 activates transcription in response to various cellular stresses in part by facilitating recruitment of the transcription machinery to DNA. Recent studies have documented a direct yet poorly characterized interaction between p53 and RNA polymerase II (Pol II). Therefore, we dissected the human p53/Pol II interaction via single-particle cryo-electron microscopy, structural docking, and biochemical analyses. This study reveals that p53 binds Pol II via the Rpb1 and Rpb2 subunits, bridging the DNA-binding cleft of Pol II proximal to the upstream DNA entry site. In addition, the key DNA-binding surface of p53, frequently disrupted in various cancers, remains exposed within the assembly. Furthermore, the p53/Pol II cocomplex displays a closed conformation as defined by the position of the Pol II clamp domain. Notably, the interaction of p53 and Pol II leads to increased Pol II elongation activity. These findings indicate that p53 may structurally regulate DNA-binding functions of Pol II via the clamp domain, thereby providing insights into p53-regulated Pol II transcription.

Cluster Analysis of p53 Binding Site Sequences Reveals Subsets with Different Functions

p53 is an important regulator of cell cycle arrest, senescence, apoptosis and metabolism, and is frequently mutated in tumors. It functions as a tetramer, where each component dimer binds to a decameric DNA region known as a response element. We identify p53 binding site subtypes and examine the functional and evolutionary properties of these subtypes. We start with over 1700 known binding sites and, with no prior labeling, identify two sets of response elements by unsupervised clustering. When combined, they give rise to three types of p53 binding sites. We find that probabilistic and alignment-based assessments of cross-species conservation show no strong evidence of differential conservation between types of binding sites. In contrast, functional analysis of the genes most proximal to the binding sites provides strong bioinformatic evidence of functional differentiation between the three types of binding sites. Our results are consistent with recent structural data identifying two conformations of the L1 loop in the DNA binding domain, suggesting that they reflect biologically meaningful groups imposed by the p53 protein structure.

Full-length p53 tetramer bound to DNA and its quaternary dynamics

P53 is a major tumor suppressor that is mutated and inactivated in ~50% of all human cancers. Thus, reactivation of mutant p53 using small molecules has been a long sought-after anticancer therapeutic strategy. Full structural characterization of the full-length oligomeric p53 is challenging because of its complex architecture and multiple highly flexible regions. To explore p53 structural dynamics, here we developed a series of atomistic integrative models with available crystal structures of the full-length p53 (fl-p53) tetramer bound to three different DNA sequences: a p21 response element, a puma response element and a nonspecific DNA sequence. Explicitly solvated, all-atom molecular dynamics simulations of the three complexes (totaling nearly 1 μs of aggregate simulation time) yield final structures consistent with electron microscopy maps and, for the first time, show the direct interactions of the p53 C-terminal with DNA. Through a collective principal component analysis, we identify sequence-dependent differential quaternary binding modes of the p53 tetramer interfacing with DNA. Additionally, L1 loop dynamics of fl-p53 in the presence of DNA is revealed, and druggable pockets of p53 are identified via solvent mapping to aid future drug discovery studies.

Motif co-regulation and co-operativity are common mechanisms in transcriptional, post-transcriptional and post-translational regulation

A substantial portion of the regulatory interactions in the higher eukaryotic cell are mediated by simple sequence motifs in the regulatory segments of genes and (pre-)mRNAs, and in the intrinsically disordered regions of proteins. Although these regulatory modules are physicochemically distinct, they share an evolutionary plasticity that has facilitated a rapid growth of their use and resulted in their ubiquity in complex organisms. The ease of motif acquisition simplifies access to basal housekeeping functions, facilitates the co-regulation of multiple biomolecules allowing them to respond in a coordinated manner to changes in the cell state, and supports the integration of multiple signals for combinatorial decision-making. Consequently, motifs are indispensable for temporal, spatial, conditional and basal regulation at the transcriptional, post-transcriptional and post-translational level. In this review, we highlight that many of the key regulatory pathways of the cell are recruited by motifs and that the ease of motif acquisition has resulted in large networks of co-regulated biomolecules. We discuss how co-operativity allows simple static motifs to perform the conditional regulation that underlies decision-making in higher eukaryotic biological systems. We observe that each gene and its products have a unique set of DNA, RNA or protein motifs that encode a regulatory program to define the logical circuitry that guides the life cycle of these biomolecules, from transcription to degradation. Finally, we contrast the regulatory properties of protein motifs and the regulatory elements of DNA and (pre-)mRNAs, advocating that co-regulation, co-operativity, and motif-driven regulatory programs are common mechanisms that emerge from the use of simple, evolutionarily plastic regulatory modules.

Sequence Variation in the Response Element Determines Binding by the Transcription Factor p73

How the sequence of a response element affects the binding of a transcription factor and, ultimately, the differential rate of transcription of genes under its control is not well-understood. In the case of the p73 transcription factor, it binds to >200 response elements to trigger developmental, cell arrest, and apoptotic pathways. The p73 response elements match the 20 bp consensus sequence of the p53 response elements that are formed by two 10 bp half-sites, where each half-site is an inverted repeat of two 5 bp quarter-sites. Using sedimentation velocity and fluorescence anisotropy experiments, we studied how systematic variations in the sequence of a half-site response element modify the DNA binding affinity of the p73 DNA-binding domain. We observed that each nucleotide position in the response element has a different influence in determining the binding of the p73 DNA-binding domain. The cytosine in the fourth position of each quarter-site is the largest determinant of DNA binding, followed by the nucleotide in the fifth position, and last, the first three positions show a slight regulatory preference for purines. Together with previous structural and functional results, our data suggest a hierarchical model of binding in which some nucleotide positions in the response element are more important than others in determining the binding of the transcription factor.

The variation game: Cracking complex genetic disorders with NGS and omics data

Tremendous advances in Next Generation Sequencing (NGS) and high-throughput omics methods have brought us one step closer towards mechanistic understanding of the complex disease at the molecular level. In this review, we discuss four basic regulatory mechanisms implicated in complex genetic diseases, such as cancer, neurological disorders, heart disease, diabetes, and many others. The mechanisms, including genetic variations, copy-number variations, posttranscriptional variations, and epigenetic variations, can be detected using a variety of NGS methods. We propose that malfunctions detected in these mechanisms are not necessarily independent, since these malfunctions are often found associated with the same disease and targeting the same gene, group of genes, or functional pathway. As an example, we discuss possible rewiring effects of the cancer-associated genetic, structural, and posttranscriptional variations on the protein-protein interaction (PPI) network centered around P53 protein. The review highlights multi-layered complexity of common genetic disorders and suggests that integration of NGS and omics data is a critical step in developing new computational methods capable of deciphering this complexity.

An Integrated Spin-Labeling/Computational-Modeling Approach for Mapping Global Structures of Nucleic Acids

The technique of site-directed spin labeling (SDSL) provides unique information on biomolecules by monitoring the behavior of a stable radical tag (i.e., spin label) using electron paramagnetic resonance (EPR) spectroscopy. In this chapter, we describe an approach in which SDSL is integrated with computational modeling to map conformations of nucleic acids. This approach builds upon a SDSL tool kit previously developed and validated, which includes three components: (i) a nucleotide-independent nitroxide probe, designated as R5, which can be efficiently attached at defined sites within arbitrary nucleic acid sequences; (ii) inter-R5 distances in the nanometer range, measured via pulsed EPR; and (iii) an efficient program, called NASNOX, that computes inter-R5 distances on given nucleic acid structures. Following a general framework of data mining, our approach uses multiple sets of measured inter-R5 distances to retrieve "correct" all-atom models from a large ensemble of models. The pool of models can be generated independently without relying on the inter-R5 distances, thus allowing a large degree of flexibility in integrating the SDSL-measured distances with a modeling approach best suited for the specific system under investigation. As such, the integrative experimental/computational approach described here represents a hybrid method for determining all-atom models based on experimentally-derived distance measurements.

Mutant TP53 posttranslational modifications: challenges and opportunities

The wild-type (WT) human p53 (TP53) tumor suppressor can be posttranslationally modified at over 60 of its 393 residues. These modifications contribute to changes in TP53 stability and in its activity as a transcription factor in response to a wide variety of intrinsic and extrinsic stresses in part through regulation of protein-protein and protein-DNA interactions. The TP53 gene frequently is mutated in cancers, and in contrast to most other tumor suppressors, the mutations are mostly missense often resulting in the accumulation of mutant (MUT) protein, which may have novel or altered functions. Most MUT TP53s can be posttranslationally modified at the same residues as in WT TP53. Strikingly, however, codons for modified residues are rarely mutated in human tumors, suggesting that TP53 modifications are not essential for tumor suppression activity. Nevertheless, these modifications might alter MUT TP53 activity and contribute to a gain-of-function leading to increased metastasis and tumor progression. Furthermore, many of the signal transduction pathways that result in TP53 modifications are altered or disrupted in cancers. Understanding the signaling pathways that result in TP53 modification and the functions of these modifications in both WT TP53 and its many MUT forms may contribute to more effective cancer therapies.

Structural studies on mechanisms to activate mutant p53

The design of a broad-spectrum cancer drug would provide enormous clinical benefits to treat cancer patients. Most of cancerous cells have a mutation in the p53 gene that results in an inactive mutant p53 protein. For this reason, p53 is a prime target for the development of a broad-spectrum cancer drug. To provide the atomic information to rationally design a drug to recover p53 activity is the main goal of the structural studies on mutant p53. We review three mechanisms that influence p53 activity and provide information about how reactivation of mutant p53 can be achieved: stabilization of the active conformation of the DNA-binding domain of the protein, suppression of missense mutations in the DNA-binding domain by a second-site mutation, and increased transactivation.

Reversal of the DNA-binding-induced loop L1 conformational switch in an engineered human p53 protein

The gene encoding the p53 tumor suppressor protein, a sequence-specific DNA binding transcription factor, is the most frequently mutated gene in human cancer. Crystal structures of homo-oligomerizing p53 polypeptides with specific DNA suggest that DNA binding is associated with a conformational switch. Specifically, in the absence of DNA, loop L1 of the p53 DNA binding domain adopts an extended conformation, whereas two p53 subunits switch to a recessed loop L1 conformation when bound to DNA as a tetramer. We previously designed a p53 protein, p53FG, with amino substitutions S121F and V122G targeting loop L1. These two substitutions enhanced the affinity of p53 for specific DNA yet, counterintuitively, decreased the residency time of p53 on DNA. Here, we confirmed these DNA binding properties of p53FG using a different method. We also determined by crystallography the structure of p53FG in its free state and bound to DNA as a tetramer. In the free state, loop L1 adopted a recessed conformation, whereas upon DNA binding, two subunits switched to the extended loop L1 conformation, resulting in a final structure that was very similar to that of wild-type p53 bound to DNA. Thus, altering the apo structure of p53 changed its DNA binding properties, even though the DNA-bound structure was not altered.

Conformations of p53 response elements in solution deduced using site-directed spin labeling and Monte Carlo sampling

The tumor suppressor protein p53 regulates numerous signaling pathways by specifically recognizing diverse p53 response elements (REs). Understanding the mechanisms of p53-DNA interaction requires structural information on p53 REs. However, such information is limited as a 3D structure of any RE in the unbound form is not available yet. Here, site-directed spin labeling was used to probe the solution structures of REs involved in p53 regulation of the p21 and Bax genes. Multiple nanometer distances in the p21-RE and BAX-RE, measured using a nucleotide-independent nitroxide probe and double-electron-electron-resonance spectroscopy, were used to derive molecular models of unbound REs from pools of all-atom structures generated by Monte-Carlo simulations, thus enabling analyses to reveal sequence-dependent DNA shape features of unbound REs in solution. The data revealed distinct RE conformational changes on binding to the p53 core domain, and support the hypothesis that sequence-dependent properties encoded in REs are exploited by p53 to achieve the energetically most favorable mode of deformation, consequently enhancing binding specificity. This work reveals mechanisms of p53-DNA recognition, and establishes a new experimental/computational approach for studying DNA shape in solution that has far-reaching implications for studying protein-DNA interactions.

Mapping the structural and dynamical features of multiple p53 DNA binding domains: insights into loop 1 intrinsic dynamics

The transcription factor p53 regulates cellular integrity in response to stress. p53 is mutated in more than half of cancerous cells, with a majority of the mutations localized to the DNA binding domain (DBD). In order to map the structural and dynamical features of the DBD, we carried out multiple copy molecular dynamics simulations (totaling 0.8 μs). Simulations show the loop 1 to be the most dynamic element among the DNA-contacting loops (loops 1-3). Loop 1 occupies two major conformational states: extended and recessed; the former but not the latter displays correlations in atomic fluctuations with those of loop 2 (~24 Å apart). Since loop 1 binds to the major groove whereas loop 2 binds to the minor groove of DNA, our results begin to provide some insight into the possible mechanism underpinning the cooperative nature of DBD binding to DNA. We propose (1) a novel mechanism underlying the dynamics of loop 1 and the possible tread-milling of p53 on DNA and (2) possible mutations on loop 1 residues to restore the transcriptional activity of an oncogenic mutation at a distant site.

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.

Transactivation specificity is conserved among p53 family proteins and depends on a response element sequence code

Structural and biochemical studies have demonstrated that p73, p63 and p53 recognize DNA with identical amino acids and similar binding affinity. Here, measuring transactivation activity for a large number of response elements (REs) in yeast and human cell lines, we show that p53 family proteins also have overlapping transactivation profiles. We identified mutations at conserved amino acids of loops L1 and L3 in the DNA-binding domain that tune the transactivation potential nearly equally in p73, p63 and p53. For example, the mutant S139F in p73 has higher transactivation potential towards selected REs, enhanced DNA-binding cooperativity in vitro and a flexible loop L1 as seen in the crystal structure of the protein–DNA complex. By studying, how variations in the RE sequence affect transactivation specificity, we discovered a RE-transactivation code that predicts enhanced transactivation; this correlation is stronger for promoters of genes associated with apoptosis.

Structure of p53 binding to the BAX response element reveals DNA unwinding and compression to accommodate base-pair insertion

The p53 core domain binds to response elements (REs) that contain two continuous half-sites as a cooperative tetramer, but how p53 recognizes discontinuous REs is not well understood. Here we describe the crystal structure of the p53 core domain bound to a naturally occurring RE located at the promoter of the Bcl-2-associated X protein (BAX) gene, which contains a one base-pair insertion between the two half-sites. Surprisingly, p53 forms a tetramer on the BAX-RE that is nearly identical to what has been reported on other REs with a 0-bp spacer. Each p53 dimer of the tetramer binds in register to a half-site and maintains the same protein–DNA interactions as previously observed, and the two dimers retain all the protein–protein contacts without undergoing rotation or translation. To accommodate the additional base pair, the DNA is deformed and partially disordered around the spacer region, resulting in an apparent unwinding and compression, such that the interactions between the dimers are maintained. Furthermore, DNA deformation within the p53-bound BAX-RE is confirmed in solution by site-directed spin labeling measurements. Our results provide a structural insight into the mechanism by which p53 binds to discontinuous sites with one base-pair spacer.

Computational identification of a transiently open L1/S3 pocket for reactivation of mutant p53

The tumour suppressor p53 is the most frequently mutated gene in human cancer. Reactivation of mutant p53 by small molecules is an exciting potential cancer therapy. Although several compounds restore wild-type function to mutant p53, their binding sites and mechanisms of action are elusive. Here computational methods identify a transiently open binding pocket between loop L1 and sheet S3 of the p53 core domain. Mutation of residue Cys124, located at the centre of the pocket, abolishes p53 reactivation of mutant R175H by PRIMA-1, a known reactivation compound. Ensemble-based virtual screening against this newly revealed pocket selects stictic acid as a potential p53 reactivation compound. In human osteosarcoma cells, stictic acid exhibits dose-dependent reactivation of p21 expression for mutant R175H more strongly than does PRIMA-1. These results indicate the L1/S3 pocket as a target for pharmaceutical reactivation of p53 mutants.

About 40% of human cancers carry missense mutations in the tumour suppressor protein p53. Here the authors identify a transiently open pocket in the protein, and by targeting a small molecule to it, partially restore mutant p53 tumour suppressor activity.

Crystal structures of the DNA-binding domain tetramer of the p53 tumor suppressor family member p73 bound to different full-site response elements

How cells choose between developmental pathways remains a fundamental biological question. In the case of the p53 protein family, its three transcription factors (p73, p63, and p53) each trigger a gene expression pattern that leads to specific cellular pathways. At the same time, these transcription factors recognize the same response element (RE) consensus sequences, and their transactivation of target genes overlaps. We aimed to understand target gene selectivity at the molecular level by determining the crystal structures of the p73 DNA-binding domain (DBD) in complex with full-site REs that vary in sequence. We report two structures of the p73 DBD bound as a tetramer to 20-bp full-site REs based on two distinct quarter-sites: GAACA and GAACC. Our study confirms that the DNA-binding residues are conserved within the p53 family, whereas the dimerization and tetramerization interfaces diverge. Moreover, a conserved lysine residue in loop L1 of the DBD senses the presence of guanines in positions 2 and 3 of the quarter-site RE, whereas a conserved arginine in loop 3 adapts to changes in position 5. Sequence variations in the RE elicit a p73 conformational response that might explain target gene specificity.

Structure of p73 DNA-binding domain tetramer modulates p73 transactivation

The transcription factor p73 triggers developmental pathways and overlaps stress-induced p53 transcriptional pathways. How p53-family response elements determine and regulate transcriptional specificity remains an unsolved problem. In this work, we have determined the first crystal structures of p73 DNA-binding domain tetramer bound to response elements with spacers of different length. The structure and function of the adaptable tetramer are determined by the distance between two half-sites. The structures with zero and one base-pair spacers show compact p73 DNA-binding domain tetramers with large tetramerization interfaces; a two base-pair spacer results in DNA unwinding and a smaller tetramerization interface, whereas a four base-pair spacer hinders tetramerization. Functionally, p73 is more sensitive to spacer length than p53, with one base-pair spacer reducing 90% of transactivation activity and longer spacers reducing transactivation to basal levels. Our results establish the quaternary structure of the p73 DNA-binding domain required as a scaffold to promote transactivation.