Structure and stability insights into tumour suppressor p53 evolutionary related proteins

Disease-related p63 DBD mutations impair DNA binding by distinct mechanisms and varying degree

The transcription factor p63 shares a high sequence identity with the tumour suppressor p53 which manifests itself in high structural similarity and preference for DNA sequences. Mutations in the DNA binding domain (DBD) of p53 have been studied in great detail, enabling a general mechanism-based classification. In this study we provide a detailed investigation of all currently known mutations in the p63 DBD, which are associated with developmental syndromes, by measuring their impact on transcriptional activity, DNA binding affinity, zinc binding capacity and thermodynamic stability. Some of the mutations we have further characterized with respect to their ability to convert human dermal fibroblasts into induced keratinocytes. Here we propose a classification of the p63 DBD mutations based on the four different mechanisms of DNA binding impairment which we identified: direct DNA contact, zinc finger region, H2 region, and dimer interface mutations. The data also demonstrate that, in contrast to p53 cancer mutations, no p63 mutation induces global unfolding and subsequent aggregation of the domain. The dimer interface mutations that affect the DNA binding affinity by disturbing the interaction between the individual DBDs retain partial DNA binding capacity which correlates with a milder patient phenotype.

Water Leakage Pathway Leads to Internal Hydration of the p53 Core Domain

The gene encoding the p53 tumor suppressor protein is the most frequently mutated oncogene in cancer patients; yet, generalized strategies for rescuing the function of different p53 mutants remain elusive. This work investigates factors that may contribute to the low inherent stability of the wild-type p53 core domain (p53C) and structurally compromised Y220C mutant. Pressure-induced unfolding of p53C was compared to p63C, the p53 family member with the highest stability, the engineered superstable p53C hexamutant (p53C HM), and lower stability p53C Y220C cancer-associated mutant. The following pressure unfolding values (P50% bar) were obtained: p53C 3346, p53C Y220C 2217, p53C HM 3943, and p63C 4326. Molecular dynamics (MD) simulations revealed that p53C Y220C was most prone to water infiltration, followed by p53C, whereas the interiors of p53C HM and p63C remained comparably dry. A strong correlation (r² = 0.92) between P50% and extent of interior hydration was observed. The pathways of individual water molecule entry and exit were mapped and analyzed, revealing a common route preserved across the p53 family involving a previously reported pocket, along with a novel surface cleft, both of which appear to be targetable by small molecules. Potential determinants of propensity to water incursion were assessed, including backbone hydrogen bond protection and combined sequence and structure similarity. Collectively, our results indicate that p53C has an intrinsic susceptibility to water leakage, which is exacerbated in a structural class mutant, suggesting that there may be a common avenue for rescuing p53 function.

Conformational Heterogeneity and Frustration of the Tumor Suppressor p53 as Tuned by Punctual Mutations

The conformational heterogeneity of the p53 tumor suppressor, the wild-type (p53wt) and mutated forms, was investigated by a computational approach, including the modeling and all atoms of the molecular dynamics (MD) simulations. Four different punctual mutations (p53R175H, p53R248Q, p53R273H, and p53R282W) which are known to affect the DNA binding and belong to the most frequent hot-spot mutations in human cancers, were taken into consideration. The MD trajectories of the wild-type and mutated p53 forms were analyzed by essential dynamics to extract the relevant collective motions and by the frustration method to evaluate the degeneracy of the energy landscape. We found that p53 is characterized by wide collective motions and its energy landscape exhibits a rather high frustration level, especially in the regions involved in the binding to physiological ligands. Punctual mutations give rise to a modulation of both the collective motions and the frustration of p53, with different effects depending on the mutation. The regions of p53wt and of the mutated forms characterized by a high frustration level are also largely involved in the collective motions. Such a correlation is discussed also in connection with the intrinsic disordered character of p53 and with its central functional role.

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

Hidden electrostatic energy contributions define dynamic allosteric communications within p53 during molecular recognition

Molecular recognition is fundamental to transcription regulation. As a transcription factor, the tumor suppressor p53 has to recognize either specific DNA sequences or repressor protein partners. However, the molecular mechanism underlying the p53 conformational switch from the DNA-bound to repressor-bound states is not fully characterized. The highly charged nature of these interacting molecules prompted us to explore the nonbonded energy contributions behind molecular recognition of either a DNA or the repressor protein iASPP by p53 DNA binding domain (p53DBD), using molecular dynamics simulation followed by rigorous analyses of energy terms. Our results illuminate the allosteric pathway by which iASPP binding to p53 diminishes binding affinity between p53 and DNA. Even though the p53DBD uses a common framework of residues for recognizing both DNA and iASPP, a comparison of the electrostatics in the two p53DBD complexes revealed significant differences in residue-wise contributions to the electrostatic energy. We found that an electrostatic allosteric communication path exists in the presence of both substrates. It consists of evolutionarily conserved residues, from residue K120 of the binding loop L1 to a distal residue R213 of p53DBD. K120 is near the DNA in the p53DBD-DNA complex, whereas iASPP binding moves it away from its DNA binding position in the p53DBD-iASPP complex. The "energy hubs" (the residues show a higher degree of connectivity with other residues in the electrostatic networks) determined from the electrostatic network analysis established that this conformational change in K120 completely rewires the electrostatic network from K120 to R213, thereby impeding DNA binding. Furthermore, we found shifting populations of hydrogen bonds and salt bridges reduce pairwise electrostatic energies within p53DBD in its DNA-bound state.

Temperature Modulation of the DBDp53 Structure as Monitored by Static and Time-Resolved Fluorescence Combined with Molecular Dynamics Simulations

Trp146 of the p53 DNA-binding domain (DBD) was investigated by static and time-resolved fluorescence combined with molecular dynamics (MD) simulations at different temperatures (25, 30, 37, and 45 °C). Static emission spectra exhibit an intensity maximum at 30 °C without any substantial peak shift, while the time-resolved fluorescence displays a peculiar stretched exponential decay, indicative of a structural disorder, at all of the investigated temperatures. The stretched exponential parameter was found to increase at 37 °C. An analysis of the MD simulation trajectories evidenced the occurrence of jumps in the temporal evolution of the distances between Trp146 and residues Arg110, Asp228, Cys229, and Gln144, which are mainly responsible for Trp146 fluorescence quenching. The times that these quenchers spend close to or far from Trp146 can provide an explanation for the static fluorescence behavior. Further essential dynamics analysis of the MD trajectories indicates a significant restriction of protein global motions above 37 °C. These results are consistent with a decrease in the structural heterogeneity of DBD as the temperature increases. The results are also discussed in view of understanding how temperature can modulate the p53 capability to binding partners, including DNA.

Molecular insights into the differential structure-dynamics-stability features of interleukin-8 orthologs: Implications to functional specificity

Chemokines are a sub-group of chemotactic cytokines that regulate the leukocyte migration by binding to G-protein coupled receptors (GPCRs) and cell surface glycosaminoglycans (GAGs). Interleukin-8 (CXCL8/IL8) is one of the most essential CXC chemokine that has been reported to be involved in various pathophysiological conditions. Structure-function relationships of human IL8 have been studied extensively. However, no such detailed information is available on IL8 orthologs, although they exhibit significant functional divergence. In order to unravel the differential structure-dynamics-stability-function relationship of IL8 orthologs, comparative molecular analysis was performed on canine (laurasians) and human (primates) IL8 proteins using in-silico molecular evolutionary analysis and solution NMR spectroscopy methods. The residue level NMR studies suggested that, although the overall structural architecture of canine IL8 is similar to that of human IL8, systematic differences were observed in their backbone dynamics and low-energy excited states due to amino acid substitutions. Further, these substitutions also resulted in attenuation of stability and heparin binding affinity in the canine IL8 as compared to its human counterpart. Indeed, structural and sequence analysis evidenced for specificity of molecular interactions with cognate receptor (CXCR1) and glycosaminoglycan (heparin), thus providing evidence for a noticeable functional specificity and divergence between the two IL8 orthologs.

TRIM3 inhibits P53 signaling in breast cancer cells

BACKGROUND

Beast cancer is the most common women cancer worldwide, while two third of them are ER alpha positive breast cancer. Among the ER alpha positive breast cancer, about 80% are P53 wild type, indicating the potential tumor suppression role in ER alpha positive breast cancer. Since P53 is an important safeguard to inhibit cell malignant transformation, reactivating P53 signaling could a plausible approach to treat breast cancer.

METHODS

TRIM3 protein levels were measured by western blot, while the P53 classical target genes were measured by real-time PCR. WST1 assay were used to measure cell proliferation, while cleaved caspase-3 was used to detect cell apoptosis. Protein stability and ubiquitin assay were used to detect the P53 protein ubiquitin and stability. The immuno-precipitation assays were used to detect the protein interactions. Immuno-staining was used to detect the protein localization of P53 and TRIM3, while the ubiquitin-based immuno-precipitation assays were used to detect the specific ubiquitination manner of P53.

RESULTS

In our study, we identified TRIM3 as an endogenous inhibitor for P53 signaling. TRIM3 depletion inhibited breast cancer cell proliferation and promoted apoptosis. In addition, TRIM3 depletion increased P53 protein level in breast cancer cell. Further investigation showed that TRIM3 could associate with P53 and promote P53 K48-linked ubiquitination and degradation.

CONCLUSION

Our study identified a novel post-translational modification mechanism between TRIM3 and P53. TRIM3 depletion or blockage could be a promising strategy to rescue P53 signaling and inhibit breast cancer progression.

Time-Resolved Fluorescence and Essential Dynamics Study on the Structural Heterogeneity of p53DBD Bound to the Anticancer p28 Peptide

Time-resolved fluorescence emission was combined with molecular dynamics (MD) simulations to investigate the DNA-binding domain (DBD) of the tumor suppressor p53 alone and its complex with the anticancer peptide p28 (DBD/p28). The fluorescence emission decay of the lone Trp residue, from both DBD and DBD/p28, was well-described by a stretched exponential function. Such a behavior was ascribed to heterogeneity in the Trp relaxation behavior, likely due to the coexistence of different conformational states. The increase of the stretching parameter, on passing from DBD to DBD/p28, indicates a reduced heterogeneity in the Trp146 environment for DBD/p28. Moreover, the effects of p28 on the global dynamics of DBD were analyzed by the essential dynamics method on 30 ns long MD trajectories of both DBD and DBD/p28. We found the establishment of wide-amplitude anharmonic modes throughout the DBD molecule, with a particularly high amplitude being detected in the DNA-binding region. These modes are significantly reduced when DBD is bound to p28, consistently with a structure stabilization. In summary, the results indicate that p28 binding has a strong effect on both the local and global heterogeneity of DBD, thus providing some hints to the understanding of its anticancer activity.

Roles of computational modelling in understanding p53 structure, biology, and its therapeutic targeting

The transcription factor p53 plays pivotal roles in numerous biological processes, including the suppression of tumours. The rich availability of biophysical data aimed at understanding its structure–function relationships since the 1990s has enabled the application of a variety of computational modelling techniques towards the establishment of mechanistic models. Together they have provided deep insights into the structure, mechanics, energetics, and dynamics of p53. In parallel, the observation that mutations in p53 or changes in its associated pathways characterize several human cancers has resulted in a race to develop therapeutic modulators of p53, some of which have entered clinical trials. This review describes how computational modelling has played key roles in understanding structural-dynamic aspects of p53, formulating hypotheses about domains that are beyond current experimental investigations, and the development of therapeutic molecules that target the p53 pathway.

Simulations of mutant p53 DNA binding domains reveal a novel druggable pocket

The DNA binding domain (DBD) of the tumor suppressor p53 is the site of several oncogenic mutations. A subset of these mutations lowers the unfolding temperature of the DBD. Unfolding leads to the exposure of a hydrophobic β-strand and nucleates aggregation which results in pathologies through loss of function and dominant negative/gain of function effects. Inspired by the hypothesis that structural changes that are associated with events initiating unfolding in DBD are likely to present opportunities for inhibition, we investigate the dynamics of the wild type (WT) and some aggregating mutants through extensive all atom explicit solvent MD simulations. Simulations reveal differential conformational sampling between the WT and the mutants of a turn region (S6-S7) that is contiguous to a known aggregation-prone region (APR). The conformational properties of the S6-S7 turn appear to be modulated by a network of interacting residues. We speculate that changes that take place in this network as a result of the mutational stress result in the events that destabilize the DBD and initiate unfolding. These perturbations also result in the emergence of a novel pocket that appears to have druggable characteristics. FDA approved drugs are computationally screened against this pocket.

Raman Evidence of p53-DBD Disorder Decrease upon Interaction with the Anticancer Protein Azurin

Raman spectroscopy, which is a suitable tool to elucidate the structural properties of intrinsically disordered proteins, was applied to investigate the changes in both the structure and the conformational heterogeneity of the DNA-binding domain (DBD) belonging to the intrinsically disordered protein p53 upon its binding to Azurin, an electron-transfer anticancer protein from Pseudomonas aeruginosa. The Raman spectra of the DBD and Azurin, isolated in solution or forming a complex, were analyzed by a combined analysis based on peak inspection, band convolution, and principal component analysis (PCA). In particular, our attention was focused on the Raman peaks of Tyrosine and Tryptophan residues, which are diagnostic markers of protein side chain environment, and on the Amide I band, of which the deconvolution allows us to extract information about α-helix, β-sheet, and random coil contents. The results show an increase of the secondary structure content of DBD concomitantly with a decrease of its conformational heterogeneity upon its binding to Azurin. These findings suggest an Azurin-induced conformational change of DBD structure with possible implications for p53 functionality.

Analysis of the TP53 Deleterious Single Nucleotide Polymorphisms Impact on Estrogen Receptor Alpha-p53 Interaction: A Machine Learning Approach

Breast cancer is a leading cancer type and one of the major health issues faced by women around the world. Some of its major risk factors include body mass index, hormone replacement therapy, family history and germline mutations. Of these risk factors, estrogen levels play a crucial role. Among the estrogen receptors, estrogen receptor alpha (ERα) is known to interact with tumor suppressor protein p53 directly thereby repressing its function. Previously, we have studied the impact of deleterious breast cancer-associated non-synonymous single nucleotide polymorphisms (nsnps) rs11540654 (R110P), rs17849781 (P278A) and rs28934874 (P151T) in TP53 gene on the p53 DNA-binding core domain. In the present study, we aimed to analyze the impact of these mutations on p53–ERα interaction. To this end, we, have modelled the full-length structure of human p53 and validated its quality using PROCHECK and subjected it to energy minimization using NOMAD-Ref web server. Three-dimensional structure of ERα activation function-2 (AF-2) domain was downloaded from the protein data bank. Interactions between the modelled native and mutant (R110P, P278A, P151T) p53 with ERα was studied using ZDOCK. Machine learning predictions on the interactions were performed using Weka software. Results from the protein–protein docking showed that the atoms, residues and solvent accessibility surface area (SASA) at the interface was increased in both p53 and ERα for R110P mutation compared to the native complexes indicating that the mutation R110P has more impact on the p53–ERα interaction compared to the other two mutants. Mutations P151T and P278A, on the other hand, showed a large deviation from the native p53-ERα complex in atoms and residues at the surface. Further, results from artificial neural network analysis showed that these structural features are important for predicting the impact of these three mutations on p53–ERα interaction. Overall, these three mutations showed a large deviation in total SASA in both p53 and ERα. In conclusion, results from our study will be crucial in making the decisions for hormone-based therapies against breast cancer.

Aggregation tendencies in the p53 family are modulated by backbone hydrogen bonds

The p53 family of proteins is comprised of p53, p63 and p73. Because the p53 DNA binding domain (DBD) is naturally unstable and possesses an amyloidogenic sequence, it is prone to form amyloid fibrils, causing loss of functions. To develop p53 therapies, it is necessary to understand the molecular basis of p53 instability and aggregation. Light scattering, thioflavin T (ThT) and high hydrostatic pressure (HHP) assays showed that p53 DBD aggregates faster and to a greater extent than p63 and p73 DBDs, and was more susceptible to denaturation. The aggregation tendencies of p53, p63, and p73 DBDs were strongly correlated with their thermal stabilities. Molecular Dynamics (MD) simulations indicated specific regions of structural heterogeneity unique to p53, which may be promoted by elevated incidence of exposed backbone hydrogen bonds (BHBs). The results indicate regions of structural vulnerability in the p53 DBD, suggesting new targetable sites for modulating p53 stability and aggregation, a potential approach to cancer therapy.

Identification of p53-target genes in Danio rerio

To orchestrate the genomic response to cellular stress signals, p53 recognizes and binds to DNA containing specific and well-characterized p53-responsive elements (REs). Differences in RE sequences can strongly affect the p53 transactivation capacity and occur even between closely related species. Therefore, the identification and characterization of a species-specific p53 Binding sistes (BS) consensus sequence and of the associated target genes may help to provide new insights into the evolution of the p53 regulatory networks across different species. Although p53 functions were studied in a wide range of species, little is known about the p53-mediated transcriptional signature in Danio rerio. Here, we designed and biochemically validated a computational approach to identify novel p53 target genes in Danio rerio genome. Screening all the Danio rerio genome by pattern-matching-based analysis, we found p53 RE-like patterns proximal to 979 annotated Danio rerio genes. Prioritization analysis identified a subset of 134 candidate pattern-related genes, 31 of which have been investigated in further biochemical assays. Our study identified runx1, axin1, traf4a, hspa8, col4a5, necab2, and dnajc9 genes as novel direct p53 targets and 12 additional p53-controlled genes in Danio rerio genome. The proposed combinatorial approach resulted to be highly sensitive and robust for identifying new p53 target genes also in additional animal species.

Evolution of p53 transactivation specificity through the lens of a yeast-based functional assay

Co-evolution of transcription factors (TFs) with their respective cis-regulatory network enhances functional diversity in the course of evolution. We present a new approach to investigate transactivation capacity of sequence-specific TFs in evolutionary studies. Saccharomyces cerevisiae was used as an in vivo test tube and p53 proteins derived from human and five commonly used animal models were chosen as proof of concept. p53 is a highly conserved master regulator of environmental stress responses. Previous reports indicated conserved p53 DNA binding specificity in vitro, even for evolutionary distant species. We used isogenic yeast strains where p53-dependent transactivation was measured towards chromosomally integrated p53 response elements (REs). Ten REs were chosen to sample a wide range of DNA binding affinity and transactivation capacity for human p53 and proteins were expressed at two levels using an inducible expression system. We showed that the assay is amenable to study thermo-sensitivity of frog p53, and that chimeric constructs containing an ectopic transactivation domain could be rapidly developed to enhance the activity of proteins, such as fruit fly p53, that are poorly effective in engaging the yeast transcriptional machinery. Changes in the profile of relative transactivation towards the ten REs were measured for each p53 protein and compared to the profile obtained with human p53. These results, which are largely independent from relative p53 protein levels, revealed widespread evolutionary divergence of p53 transactivation specificity, even between human and mouse p53. Fruit fly and human p53 exhibited the largest discrimination among REs while zebrafish p53 was the least selective.

Bioinformatics study of cancer-related mutations within p53 phosphorylation site motifs

p53 protein has about thirty phosphorylation sites located at the N- and C-termini and in the core domain. The phosphorylation sites are relatively less mutated than other residues in p53. To understand why and how p53 phosphorylation sites are rarely mutated in human cancer, using a bioinformatics approaches, we examined the phosphorylation site and its nearby flanking residues, focusing on the consensus phosphorylation motif pattern, amino-acid correlations within the phosphorylation motifs, the propensity of structural disorder of the phosphorylation motifs, and cancer mutations observed within the phosphorylation motifs. Many p53 phosphorylation sites are targets for several kinases. The phosphorylation sites match 17 consensus sequence motifs out of the 29 classified. In addition to proline, which is common in kinase specificity-determining sites, we found high propensity of acidic residues to be adjacent to phosphorylation sites. Analysis of human cancer mutations in the phosphorylation motifs revealed that motifs with adjacent acidic residues generally have fewer mutations, in contrast to phosphorylation sites near proline residues. p53 phosphorylation motifs are mostly disordered. However, human cancer mutations within phosphorylation motifs tend to decrease the disorder propensity. Our results suggest that combination of acidic residues Asp and Glu with phosphorylation sites provide charge redundancy which may safe guard against loss-of-function mutations, and that the natively disordered nature of p53 phosphorylation motifs may help reduce mutational damage. Our results further suggest that engineering acidic amino acids adjacent to potential phosphorylation sites could be a p53 gene therapy strategy.