Lysine120 interactions with p53 response elements can allosterically direct p53 organization

Raf/MEK/ERK Signaling Pathway Is Involved in the Inhibition of Glioma Cell Proliferation and Invasion in the Ketogenic Microenvironment

OBJECTIVE

A high-fat, low-carbohydrate ketogenic diet has been used to treat malignant glioma, in which the Raf/MEK/ERK signaling pathway is overactivated. However, whether the Raf/MEK/ERK signaling pathway is involved in the therapeutic effect of ketone bodies remains unknown. In this study, we investigated the effects of a major ketone body, 3-hydroxybutyric acid (3-HBA), on the proliferation and metastasis of malignant glioblastoma cells and the underlying mechanism.

METHODS

Two human malignant glioblastoma cell lines (U87 and U251) were treated with different concentrations of 3-HBA with or without the Raf inhibitor PAF C-16 for 24 h. Cell proliferation, cell cycle, cell invasion, and phospholipase D1 (PLD1) activity were determined. Protein and gene expression levels of Raf/MEK/ERK signaling pathway members were examined.

RESULTS

3-HBA significantly decreased cell proliferation, invasion, and intracellular PLD1 activity in both U87 and U251 glioblastoma cell lines. 3-HBA treatment significantly increased the proportion of cells in the G1 phase and decreased the proportion of cells in S phase in U87 cells. In the U251 line, the proportion of treated cells in S phase was increased and proportion of cells in G2 was decreased. 3-HBA treatment also significantly decreased the protein expression levels of Raf, MEK, p-MEK, ERK, p-ERK, and PLD1 while increasing p53 expression; an effect that was similar to treatment with the Raf inhibitor. Co-treatment of 3-HBA with the Raf inhibitor further enhanced the effects of the 3-HBA in both cell lines.

CONCLUSION

We confirmed that a ketogenic microenvironment can inhibit glioma cell proliferation and invasion by downregulating the expression of PLD1 through the Raf/MEK/ERK signaling pathway.

Distinct mechanisms control genome recognition by p53 at its target genes linked to different cell fates

The tumor suppressor p53 integrates stress response pathways by selectively engaging one of several potential transcriptomes, thereby triggering cell fate decisions (e.g., cell cycle arrest, apoptosis). Foundational to this process is the binding of tetrameric p53 to 20-bp response elements (REs) in the genome (RRRCWWGYYYN_0-13RRRCWWGYYY). In general, REs at cell cycle arrest targets (e.g. p21) are of higher affinity than those at apoptosis targets (e.g., BAX). However, the RE sequence code underlying selectivity remains undeciphered. Here, we identify molecular mechanisms mediating p53 binding to high- and low-affinity REs by showing that key determinants of the code are embedded in the DNA shape. We further demonstrate that differences in minor/major groove widths, encoded by G/C or A/T bp content at positions 3, 8, 13, and 18 in the RE, determine distinct p53 DNA-binding modes by inducing different Arg248 and Lys120 conformations and interactions. The predictive capacity of this code was confirmed in vivo using genome editing at the BAX RE to interconvert the DNA-binding modes, transcription pattern, and cell fate outcome.

Investigating Conformational Dynamics and Allostery in the p53 DNA-Binding Domain Using Molecular Simulations

The p53 tumor suppressor is a multifaceted context-dependent protein, which is involved in multiple cellular pathways, with the ability to either keep the cells alive or to kill them through mechanisms such as apoptosis. To complicate this picture, cancer cells that express mutant p53 becomes addicted to the mutant activity, so that the mutant variant features a myriad of gain-of-function activities, opening different venues for therapy. This makes essential to think outside the box and apply new approaches to the study of p53 structure-(mis)function relationship to find new critical components of its pathway or to understand how known parts are interconnected, compete, or cooperate. In this context, I will here illustrate how to integrate different computational methods to the identification of possible allosteric effects transmitted from the DNA binding interface of p53 to regions for cofactor recruitment. The protocol can be extended to any other cases of study. Indeed, it does not necessarily apply only to the study of DNA-induced effects, but more broadly to the investigation of long-range effects induced by a biological partner that binds to a biomolecule of interest.

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.

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.

Predicting conformational ensembles and genome-wide transcription factor binding sites from DNA sequences

DNA shape is emerging as an important determinant of transcription factor binding beyond just the DNA sequence. The only tool for large scale DNA shape estimates, DNAshape was derived from Monte-Carlo simulations and predicts four broad and static DNA shape features, Propeller twist, Helical twist, Minor groove width and Roll. The contributions of other shape features e.g. Shift, Slide and Opening cannot be evaluated using DNAshape. Here, we report a novel method DynaSeq, which predicts molecular dynamics-derived ensembles of a more exhaustive set of DNA shape features. We compared the DNAshape and DynaSeq predictions for the common features and applied both to predict the genome-wide binding sites of 1312 TFs available from protein interaction quantification (PIQ) data. The results indicate a good agreement between the two methods for the common shape features and point to advantages in using DynaSeq. Predictive models employing ensembles from individual conformational parameters revealed that base-pair opening - known to be important in strand separation - was the best predictor of transcription factor-binding sites (TFBS) followed by features employed by DNAshape. Of note, TFBS could be predicted not only from the features at the target motif sites, but also from those as far as 200 nucleotides away from the motif.

DNA-binding protects p53 from interactions with cofactors involved in transcription-independent functions

Binding-induced conformational changes of a protein at regions distant from the binding site may play crucial roles in protein function and regulation. The p53 tumour suppressor is an example of such an allosterically regulated protein. Little is known, however, about how DNA binding can affect distal sites for transcription factors. Furthermore, the molecular details of how a local perturbation is transmitted through a protein structure are generally elusive and occur on timescales hard to explore by simulations. Thus, we employed state-of-the-art enhanced sampling atomistic simulations to unveil DNA-induced effects on p53 structure and dynamics that modulate the recruitment of cofactors and the impact of phosphorylation at Ser215. We show that DNA interaction promotes a conformational change in a region 3 nm away from the DNA binding site. Specifically, binding to DNA increases the population of an occluded minor state at this distal site by more than 4-fold, whereas phosphorylation traps the protein in its major state. In the minor conformation, the interface of p53 that binds biological partners related to p53 transcription-independent functions is not accessible. Significantly, our study reveals a mechanism of DNA-mediated protection of p53 from interactions with partners involved in the p53 transcription-independent signalling. This also suggests that conformational dynamics is tightly related to p53 signalling.

The p53 tetramer shows an induced-fit interaction of the C-terminal domain with the DNA-binding domain

The Trp53 gene is the most frequently mutated gene in all human cancers. Its protein product p53 is a very powerful transcription factor that can activate different biochemical pathways and affect the regulation of metabolism, senescence, DNA damage response, cell cycle and cell death. The understanding of its function at the molecular level could be of pivotal relevance for therapy. Investigation of long-range intra- and interdomain communications in the p53 tetramer–DNA complex was performed by means of an atomistic model that included the tetramerization helices in the C-terminal domain, the DNA-binding domains and a consensus DNA-binding site of 18 base pairs. Nonsymmetric dynamics are illustrated in the four DNA-binding domains, with loop L1 switching from inward to outward conformations with respect to the DNA major groove. Direct intra- and intermonomeric long-range communications between the tetramerization and DNA-binding domains are noted. These long-distance conformational changes link the C terminus with the DNA-binding domain and provide a biophysical rationale for the reported functional regulation of the p53 C-terminal region. A fine characterization of the DNA deformation caused by p53 binding is obtained, with ‘static' deformations always present and measured by the slide parameter in the central thymine–adenine base pairs; we also detect ‘dynamic' deformations switched on and off by particular p53 tetrameric conformations and measured by the roll and twist parameters in the same base pairs. These different conformations can indeed modulate the electrostatic potential isosurfaces of the whole p53–DNA complex. These results provide a molecular/biophysical understanding of the evident role of the C terminus in post-translational modification that regulates the transcriptional function of p53. Furthermore, the unstructured C terminus is able to facilitate contacts between the core DNA-binding domains of the tetramer.

Structural and sequential context of p53: A review of experimental and theoretical evidence

Approximately 27 million people are suffering from cancer that contains either an inactivating missense mutation of TP53 gene or partially abrogated p53 signaling pathway. Concerted action of folded and intrinsically disordered domains accounts for multi-faceted role of p53. The intricacy of dynamic p53 structure is believed to shed light on its cellular activity for developing new cancer therapies. In this review, insights into structural details of p53, diverse single point mutations affecting its core domain, thermodynamic understanding and therapeutic strategies for pharmacological rescue of p53 function has been illustrated. An effort has been made here to bridge the structural and sequential evidence of p53 from experimental to computational studies. First, we focused on the individual domains and the crucial protein-protein or DNA-protein contacts that determine conformation and dynamic behavior of p53. Next, the oncogenic mutations associated with cancer and its contribution to thermodynamic fluctuation has been discussed. Thus the emerging anti-cancer strategies include targeting of destabilized cancer mutants with selective inhibition of its negative regulators. Recent advances in development of small molecule inhibitors and peptides exploiting p53-MDM2 interaction has been included. In a nutshell, this review attempts to describe structural biology of p53 which provide new openings for structure-guided rescue.

PyInteraph: a framework for the analysis of interaction networks in structural ensembles of proteins

In the last years, a growing interest has been gathering around the ability of Molecular Dynamics (MD) to provide insight into the paths of long-range structural communication in biomolecules. The knowledge of the mechanisms related to structural communication helps in the rationalization in atomistic details of the effects induced by mutations, ligand binding, and the intrinsic dynamics of proteins. We here present PyInteraph, a tool for the analysis of structural ensembles inspired by graph theory. PyInteraph is a software suite designed to analyze MD and structural ensembles with attention to binary interactions between residues, such as hydrogen bonds, salt bridges, and hydrophobic interactions. PyInteraph also allows the different classes of intra- and intermolecular interactions to be represented, combined or alone, in the form of interaction graphs, along with performing network analysis on the resulting interaction graphs. The program also integrates the network description with a knowledge-based force field to estimate the interaction energies between side chains in the protein. It can be used alone or together with the recently developed xPyder PyMOL plugin through an xPyder-compatible format. The software capabilities and associated protocols are here illustrated by biologically relevant cases of study. The program is available free of charge as Open Source software via the GPL v3 license at http://linux.btbs.unimib.it/pyinteraph/.

Regulation of the cell fate by DNA damage and hypoxia

In order to provide the means for the design of novel rational anti-cancer drug therapies research efforts are concentrated on unravelling the molecular circuits which induce programmed cell death and block proliferation of cancer cells. Modern therapeutic strategies are based on the understanding of the complexity of physiological functions such as differentiation, development, immune responses, cell-cycle arrest, DNA damage repair, apoptosis, autophagy, energy metabolism, and senescence. It has become evident that this knowledge will provide the means to target the components of the pathways involved in these processes in a specific and selective manner thus paving the way for the development of effective and personalised anti-cancer therapies. Transcription is a crucial cellular process that regulates a multitude of physiological functions, which are essential in disease progression and cellular response to therapy. Transcription factors such as the p53 tumor suppressor and the hypoxia-inducible factor-α (HIF-α) are key players in carcinogenesis and cellular response to cancer therapies. Both of these transcription factors regulate gene expression of genes involved in cell death and proliferation, in some cases cooperating towards producing the same outcome and in some others mediating opposing effects. It is thus apparent that fine tuning of the activity of these transcription factors is essential to determine the cellular response to therapeutic regimens, in other words whether tumor cells will commit to apoptosis or evade engagement with the anti-proliferative effects of drugs leading to drug resistance. Our observations support the notion that the functional crosstalk between HIF-1α and p53 pathways and thus the fine tuning of their transcriptional activity is mediated by cofactors shared between the two transcription factors such as components of the p300 co-activator multiprotein complex. In particular, there is evidence to suggest that differential composition of the co-modulatory protein complexes associated with p53 and HIF-1α under diverse types of stress conditions differentially regulate the expression of distinct subsets of p53 and HIF-1α target genes involved in processes such as cell cycle arrest, apoptosis, chronic inflammation, and cellular energy metabolism thereby determining the cellular fate under particular types of micro-environmental stress.

Molecular dynamics of the full-length p53 monomer

The p53 protein is frequently mutated in a very large proportion of human tumors, where it seems to acquire gain-of-function activity that facilitates tumor onset and progression. A possible mechanism is the ability of mutant p53 proteins to physically interact with other proteins, including members of the same family, namely p63 and p73, inactivating their function. Assuming that this interaction might occurs at the level of the monomer, to investigate the molecular basis for this interaction, here, we sample the structural flexibility of the wild-type p53 monomeric protein. The results show a strong stability up to 850 ns in the DNA binding domain, with major flexibility in the N-terminal transactivations domains (TAD1 and TAD2) as well as in the C-terminal region (tetramerization domain). Several stable hydrogen bonds have been detected between N-terminal or C-terminal and DNA binding domain, and also between N-terminal and C-terminal. Essential dynamics analysis highlights strongly correlated movements involving TAD1 and the proline-rich region in the N-terminal domain, the tetramerization region in the C-terminal domain; Lys120 in the DNA binding region. The herein presented model is a starting point for further investigation of the whole protein tetramer as well as of its mutants.

The HDAC inhibitor depsipeptide transactivates the p53/p21 pathway by inducing DNA damage

Histone deacetylase (HDAC) inhibitors have been proven to be effective therapeutic agents to kill cancer cells through inhibiting HDAC activity or altering the structure of chromatin. As a potent HDAC inhibitor, depsipeptide not only modulates histone deacetylation but also activates non-histone protein p53 to inhibit cancer cell growth. However, the mechanism of depsipeptide-induced p53 transactivity remains unknown. Here, we show that depsipeptide causes DNA damage through induction of reactive oxygen species (ROS) generation, as demonstrated by a comet assay and by detection of the phosphorylation of H2AX. Depsipeptide induced oxidative stress was confirmed to relate to a disturbance in reduction-oxidation (redox) reactions through inhibition of the transactivation of thioredoxin reductase (TrxR) in human cancer cells. Upon treatment with depsipeptide, p53 phosphorylation at threonine 18 (Thr18) was specifically induced. Furthermore, we also demonstrated that phosphorylation of p53 at Thr18 is required for p53 acetylation at lysine 373/382 and for p21 expression in response to depsipeptide treatment. Our results demonstrate that depsipeptide plays an anti-neoplastic role by generating ROS to elicit p53/p21 pathway activation.

The role of response elements organization in transcription factor selectivity: the IFN-β enhanceosome example

What is the mechanism through which transcription factors (TFs) assemble specifically along the enhancer DNA? The IFN-β enhanceosome provides a good model system: it is small; its components' crystal structures are available; and there are biochemical and cellular data. In the IFN-β enhanceosome, there are few protein-protein interactions even though consecutive DNA response elements (REs) overlap. Our molecular dynamics (MD) simulations on different motif combinations from the enhanceosome illustrate that cooperativity is achieved via unique organization of the REs: specific binding of one TF can enhance the binding of another TF to a neighboring RE and restrict others, through overlap of REs; the order of the REs can determine which complexes will form; and the alternation of consensus and non-consensus REs can regulate binding specificity by optimizing the interactions among partners. Our observations offer an explanation of how specificity and cooperativity can be attained despite the limited interactions between neighboring TFs on the enhancer DNA. To date, when addressing selective TF binding, attention has largely focused on RE sequences. Yet, the order of the REs on the DNA and the length of the spacers between them can be a key factor in specific combinatorial assembly of the TFs on the enhancer and thus in function. Our results emphasize cooperativity via RE binding sites organization.

Identification of a novel functional domain of ricin responsible for its potent toxicity

Ribosome-inactivating proteins (RIPs) are toxic N-glycosidases that depurinate the universally conserved α-sarcin loop of large rRNAs. They have received attention in biological and biomedical research because of their unique biological activities toward animals and human cells as cell-killing agents. A better understanding of the depurination mechanism of RIPs could allow us to develop potent neutralizing antibodies and to design efficient immunotoxins for clinical use. Among these RIPs, ricin exhibited remarkable efficacy in depurination activity and highly conserved tertiary structure with other RIPs. It can be considered as a prototype to investigate the depurination mechanism of RIPs. In the present study, we successfully identified a novel functional domain responsible for controlling the depurination activity of ricin, which is located far from the enzymatic active site reported previously. Our study indicated that ricin A-chain mAbs binding to this domain (an α-helix comprising the residues 99-106) exhibited an unusual potent neutralizing ability against ricin in vivo. To further investigate the potential role of the α-helix in regulating the catalytic activity of ricin, ricin A-chain variants with different flexibility of the α-helix were rationally designed. Our data clearly demonstrated that the flexibility of the α-helix is responsible for controlling the depurination activity of ricin and determining the extent of protein synthesis inhibition, suggesting that the conserved α-helix might be considered as a potential target for the prevention and treatment of RIP poisoning.