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

Evaluating chemotherapeutic potential of soya-isoflavonoids against high penetrance genes in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is the most aggressive molecular subtype of breast cancer (BC) associated with a poor prognosis. Owing to the structural similarity with 17-β-estradiol, consumption of soya-isoflavonoids are associated with a reduced rate of hormone-receptive BC incidence, but their role in TNBC is not deciphered in detail. This present study thus aims to investigate the therapeutic binding dynamics of dietary soya-flavonoids with the six high penetrance (HP) receptors in TNBC, viz. BRCA1, BRCA2, PALB2, PTEN, STK11 and TP53. Out of the 14 soya-flavonoids screened based on ADMET descriptors and several other physicochemical, bioavailability, drug and lead-likeness properties, four hits were shortlisted (Daidzein, Genistein, Glycitein and Biochanin A). Docking and molecular dynamics (MD) simulation revealed Genistein as the most potential multi-target inhibitor of the six TNBC HP genes. Additionally, Genistein exhibited excellent binding specificity with PTEN, a potent mediator of the PI3K signaling pathway in TNBC. The binding interaction of PTEN and Genistein was further compared against a standardized FDA-approved chemotherapeutic inhibitor, Olaparib, computed through various MD trajectory analysis, principal component analysis and computation of free energy landscape. This study reveals a comparatively better binding dynamics of PTEN-Genistein than PTEN-Olaparib. With a significant global surge in biomarker-based precision therapeutics in oncology, the results of this exhaustive in-silico study thus encourage the prospect of validating PTEN as a druggable target of Genistein, a unique drug-receptor combination in the future.Communicated by Ramaswamy H. Sarma.

Investigating the Effects of Chelidonic Acid on Oxidative Stress-Induced Premature Cellular Senescence in Human Skin Fibroblast Cells

This study investigated the effects of chelidonic acid (CA) on hydrogen peroxide (H2O2) induced cellular senescence in human skin fibroblast cells (BJ). Cellular senescence is a critical mechanism that is linked to age-related diseases and chronic conditions. CA, a γ-pyrone compound known for its broad pharmacological activity, was assessed for its potential to mitigate oxidative stress and alter senescence markers. A stress-induced premature senescence (SIPS) model was designed in BJ fibroblast cells using the oxidative stress agent H2O2. After this treatment, cells were treated with CA, and the potential effect of CA on senescence was evaluated using senescence-related β-galactosidase, 4',6-diamino-2-phenylindole (DAPI), acridine-orange staining (AO), comet assay, molecular docking assays, gene expression, and protein analysis. These results demonstrate that CA effectively reduces senescence markers, including senescence-associated β-galactosidase activity, DNA damage, lysosomal activity, and oxidative stress indicators such as malondialdehyde. Molecular docking revealed CA's potential interactions with critical proteins involved in senescence signalling pathways, suggesting mechanisms by which CA may exert its effects. Gene expression and protein analyses corroborated the observed anti-senescent effects, with CA modulating p16, p21, and pRB1 expressions and reducing oxidative stress markers. In conclusion, CA appeared to have senolytic and senomorphic potential in vitro, which could mitigate and reverse SIPS markers in BJ fibroblasts.

Navigating the complexity of p53-DNA binding: implications for cancer therapy

Abstract

The tumor suppressor protein p53, a transcription factor playing a key role in cancer prevention, interacts with DNA as its primary means of determining cell fate in the event of DNA damage. When it becomes mutated, it opens damaged cells to the possibility of reproducing unchecked, which can lead to formation of cancerous tumors. Despite its critical role, therapies at the molecular level to restore p53 native function remain elusive, due to its complex nature. Nevertheless, considerable information has been amassed, and new means of investigating the problem have become available.

Objectives

We consider structural, biophysical, and bioinformatic insights and their implications for the role of direct and indirect readout and how they contribute to binding site recognition, particularly those of low consensus. We then pivot to consider advances in computational approaches to drug discovery.

Materials and methods

We have conducted a review of recent literature pertinent to the p53 protein.

Results

Considerable literature corroborates the idea that p53 is a complex allosteric protein that discriminates its binding sites not only via consensus sequence through direct H-bond contacts, but also a complex combination of factors involving the flexibility of the binding site. New computational methods have emerged capable of capturing such information, which can then be utilized as input to machine learning algorithms towards the goal of more intelligent and efficient de novo allosteric drug design.

Conclusions

Recent improvements in machine learning coupled with graph theory and sector analysis hold promise for advances to more intelligently design allosteric effectors that may be able to restore native p53-DNA binding activity to mutant proteins.

Clinical relevance

The ideas brought to light by this review constitute a significant advance that can be applied to ongoing biophysical studies of drugs for p53, paving the way for the continued development of new methodologies for allosteric drugs. Our discoveries hold promise to provide molecular therapeutics which restore p53 native activity, thereby offering new insights for cancer therapies.

Graphical Abstract

Structural representation of the p53 DBD (PDBID 1TUP). DNA consensus sequence is shown in gray, and the protein is shown in blue. Red beads indicate hotspot residue mutations, green beads represent DNA interacting residues, and yellow beads represent both.

Single molecule studies characterize the kinetic mechanism of tetrameric p53 binding to different native response elements

The transcriptional activator p53 is a tumor suppressor protein that controls cellular pathways important for cell fate decisions, including cell cycle arrest, senescence, and apoptosis. It functions as a tetramer by binding to specific DNA sequences known as response elements (REs) to control transcription via interactions with co-regulatory complexes. Despite its biological importance, the mechanism by which p53 binds REs remains unclear. To address this, we have used an in vitro single molecule fluorescence approach to quantify the dynamic binding of full-length human p53 to five native REs in real time under equilibrium conditions. Our approach enabled us to quantify the oligomeric state of DNA-bound p53. We found little evidence that dimer/DNA complexes form as intermediates en route to binding or dissociation of p53 tetramer/DNA complexes. Interestingly, however, at some REs dimers can rapidly exchange from tetramer/DNA complexes. Real time kinetic measurements enabled us to determine rate constants for association and dissociation at all five REs, which revealed two kinetically distinct populations of tetrameric p53/RE complexes. For the less stable population, the rate constants for dissociation were larger at REs closest to consensus, showing that the more favorable binding sequences form the least kinetically stable complexes. Together our single molecule measurements provide new insight into mechanisms by which tetrameric p53 forms complexes on different native REs.

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.

P63 and P73 Activation in Cancers with p53 Mutation

The members of the p53 family comprise p53, p63, and p73, and full-length isoforms of the p53 family have a tumor suppressor function. However, p53, but not p63 or p73, has a high mutation rate in cancers causing it to lose its tumor suppressor function. The top and second-most prevalent p53 mutations are missense and nonsense mutations, respectively. In this review, we discuss possible drug therapies for nonsense mutation and a missense mutation in p53. p63 and p73 activators may be able to replace mutant p53 and act as anti-cancer drugs. Herein, these p63 and p73 activators are summarized and how to improve these activator responses, particularly focusing on p53 gain-of-function mutants, is discussed.

Folding Intermediates, Heterogeneous Native Ensembles and Protein Function

Single domain proteins fold via diverse mechanisms emphasizing the intricate relationship between energetics and structure, which is a direct consequence of functional constraints and demands imposed at the level of sequence. On the other hand, elucidating the interplay between folding mechanisms and function is challenging in large proteins, given the inherent shortcomings in identifying metastable states experimentally and the sampling limitations associated with computational methods. Here, we show that free energy profiles and surfaces of large systems (>150 residues), as predicted by a statistical mechanical model, display a wide array of folding mechanisms with ubiquitous folding intermediates and heterogeneous native ensembles. Importantly, residues around the ligand binding or enzyme active site display a larger tendency to partially unfold and this manifests as intermediates or excited states along the folding coordinate in ligand binding domains, transcription repressors, and representative enzymes from all the six classes, including the SARS-CoV-2 receptor binding domain (RBD) of the spike protein and the protease M^pro. It thus appears that it is relatively easier to distill the imprints of function on the folding landscape of larger proteins as opposed to smaller systems. We discuss how an understanding of energetic-entropic features in ordered proteins can pinpoint specific avenues through which folding mechanisms, populations of partially structured states and function can be engineered.

Protein-DNA complex structure modeling based on structural template

DNA-binding is an important feature of proteins, and protein-DNA interaction involves in many life processes. Various computational methods have been developed to predict protein-DNA complex structures due to the difficulty of experimentally obtaining protein-DNA complex structures. However, prediction of protein-DNA complex is still a challenging problem compared with prediction of protein-RNA complex, this may be due to the large conformational changes between bound and unbound structure in both protein and DNA. We extend PRIME 2.0 to PRIME 2.0.1 to model protein-DNA complex structures. By comparing sequence and structure alignment methods, we found that structure-based methods can find more templates than sequence-based methods. The results of all-to-all structure alignments showed that DNA structure plays an important role in prediction of protein-DNA complex structure. By exploring the relationship of sequence and structure, we found that in protein-DNA interaction, numerous structures with dissimilar sequences have similar 3D structures and perform the similar function.

The Changes in the p53 Protein across the Animal Kingdom Point to Its Involvement in Longevity

Recently, the quest for the mythical fountain of youth has produced extensive research programs that aim to extend the healthy lifespan of humans. Despite advances in our understanding of the aging process, the surprisingly extended lifespan and cancer resistance of some animal species remain unexplained. The p53 protein plays a crucial role in tumor suppression, tissue homeostasis, and aging. Long-lived, cancer-free African elephants have 20 copies of the TP53 gene, including 19 retrogenes (38 alleles), which are partially active, whereas humans possess only one copy of TP53 and have an estimated cancer mortality rate of 11–25%. The mechanism through which p53 contributes to the resolution of the Peto’s paradox in Animalia remains vague. Thus, in this work, we took advantage of the available datasets and inspected the p53 amino acid sequence of phylogenetically related organisms that show variations in their lifespans. We discovered new correlations between specific amino acid deviations in p53 and the lifespans across different animal species. We found that species with extended lifespans have certain characteristic amino acid substitutions in the p53 DNA-binding domain that alter its function, as depicted from the Phenotypic Annotation of p53 Mutations, using the PROVEAN tool or SWISS-MODEL workflow. In addition, the loop 2 region of the human p53 DNA-binding domain was identified as the longest region that was associated with longevity. The 3D model revealed variations in the loop 2 structure in long-lived species when compared with human p53. Our findings show a direct association between specific amino acid residues in p53 protein, changes in p53 functionality, and the extended animal lifespan, and further highlight the importance of p53 protein in aging.

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.

Single molecule studies reveal that p53 tetramers dynamically bind response elements containing one or two half sites

The tumor suppressor protein p53 is critical for cell fate decisions, including apoptosis, senescence, and cell cycle arrest. p53 is a tetrameric transcription factor that binds DNA response elements to regulate transcription of target genes. p53 response elements consist of two decameric half-sites, and data suggest one p53 dimer in the tetramer binds to each half-site. Despite a broad literature describing p53 binding DNA, unanswered questions remain, due partly to the need for more quantitative and structural studies with full length protein. Here we describe a single molecule fluorescence system to visualize full length p53 tetramers binding DNA in real time. The data revealed a dynamic interaction in which tetrameric p53/DNA complexes assembled and disassembled without a dimer/DNA intermediate. On a wild type DNA containing two half sites, p53/DNA complexes existed in two kinetically distinct populations. p53 tetramers bound response elements containing only one half site to form a single population of complexes with reduced kinetic stability. Altering the spacing and helical phasing between two half sites affected both the population distribution of p53/DNA complexes and their kinetic stability. Our real time single molecule measurements of full length p53 tetramers binding DNA reveal the parameters that define the stability of p53/DNA complexes, and provide insight into the pathways by which those complexes assemble.

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.

Roles of p53 Family Structure and Function in Non-Canonical Response Element Binding and Activation

The p53 canonical consensus sequence is a 10-bp repeat of PuPuPuC(A/T)(A/T)GPyPyPy, separated by a spacer with up to 13 bases. C(A/T)(A/T)G is the core sequence and purine (Pu) and pyrimidine (Py) bases comprise the flanking sequence. However, in the p53 noncanonical sequences, there are many variations, such as length of consensus sequence, variance of core sequence or flanking sequence, and variance in number of bases making up the spacer or AT gap composition. In comparison to p53, the p53 family members p63 and p73 have been found to have more tolerance to bind and activate several of these noncanonical sequences. The p53 protein forms monomers, dimers, and tetramers, and its nonspecific binding domain is well-defined; however, those for p63 or p73 are still not fully understood. Study of p63 and p73 structure to determine the monomers, dimers or tetramers to bind and regulate noncanonical sequence is a new challenge which is crucial to obtaining a complete picture of structure and function in order to understand how p63 and p73 regulate genes differently from p53. In this review, we will summarize the rules of p53 family non-canonical sequences, especially focusing on the structure of p53 family members in the regulation of specific target genes. In addition, we will compare different software programs for prediction of p53 family responsive elements containing parameters with canonical or non-canonical sequences.

Quantifying the two-state facilitated diffusion model of protein-DNA interactions

The current report extends the facilitated diffusion model to account for conflict between the search and recognition binding modes adopted by DNA-binding proteins (DBPs) as they search DNA and subsequently recognize and bind to their specific binding site. The speed of the search dynamics is governed by the energetic ruggedness of the protein-DNA landscape, whereas the rate for the recognition process is mostly dictated by the free energy barrier for the transition between the DBP's search and recognition binding modes. We show that these two modes are negatively coupled, such that fast 1D sliding and rapid target site recognition probabilities are unlikely to coexist. Thus, a tradeoff occurs between optimizing the timescales for finding and binding the target site. We find that these two kinetic properties can be balanced to produce a fast timescale for the total target search and recognition process by optimizing frustration. Quantification of the facilitated diffusion model by including a frustration term enables it to explain several experimental observations concerning search and recognition speeds. The extended model captures experimental estimate of the energetic ruggedness of the protein-DNA landscape and predicts how various molecular properties of protein-DNA binding affect recognition kinetics. Particularly, point mutations may change the frustration and so affect protein association with DNA, thus providing a means to modulate protein-DNA affinity by manipulating the protein's association or dissociation reactions.

Mutantelec: An In Silico mutation simulation platform for comparative electrostatic potential profiling of proteins

The electrostatic potential plays a key role in many biological processes like determining the affinity of a ligand to a given protein target, and they are responsible for the catalytic activity of many enzymes. Understanding the effect that amino acid mutations will have on the electrostatic potential of a protein, will allow a thorough understanding of which residues are the most important in a protein. MutantElec, is a friendly web application for in silico generation of site-directed mutagenesis of proteins and the comparison of electrostatic potential between the wild type protein and the mutant(s), based on the three-dimensional structure of the protein. The effect of the mutation is evaluated using different approach to the traditional surface map. MutantElec provides a graphical display of the results that allows the visualization of changes occurring at close distance from the mutation and thus uncovers the local and global impact of a specific change. © 2017 Wiley Periodicals, Inc.

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.

Visualization and analysis of non-covalent contacts using the Protein Contacts Atlas

Visualizations of biomolecular structures empower us to gain insights into biological functions, generate testable hypotheses, and communicate biological concepts. Typical visualizations (such as ball and stick) primarily depict covalent bonds. In contrast, non-covalent contacts between atoms, which govern normal physiology, pathogenesis, and drug action, are seldom visualized. We present the Protein Contacts Atlas, an interactive resource of non-covalent contacts from over 100,000 PDB crystal structures. We developed multiple representations for visualization and analysis of non-covalent contacts at different scales of organization: atoms, residues, secondary structure, subunits, and entire complexes. The Protein Contacts Atlas enables researchers from different disciplines to investigate diverse questions in the framework of non-covalent contacts, including the interpretation of allostery, disease mutations and polymorphisms, by exploring individual subunits, interfaces, and protein-ligand contacts and by mapping external information. The Protein Contacts Atlas is available at http://www.mrc-lmb.cam.ac.uk/pca/ and also through PDBe.

p53 shades of Hippo

The three p53 family members, p53, p63 and p73, are structurally similar and share many biochemical activities. Yet, along with their common fundamental role in protecting genomic fidelity, each has acquired distinct functions related to diverse cell autonomous and non-autonomous processes. Similar to the p53 family, the Hippo signaling pathway impacts a multitude of cellular processes, spanning from cell cycle and metabolism to development and tumor suppression. The core Hippo module consists of the tumor-suppressive MST-LATS kinases and oncogenic transcriptional co-effectors YAP and TAZ. A wealth of accumulated data suggests a complex and delicate regulatory network connecting the p53 and Hippo pathways, in a highly context-specific manner. This generates multiple layers of interaction, ranging from interdependent and collaborative signaling to apparent antagonistic activity. Furthermore, genetic and epigenetic alterations can disrupt this homeostatic network, paving the way to genomic instability and cancer. This strengthens the need to better understand the nuances that control the molecular function of each component and the cross-talk between the different components. Here, we review interactions between the p53 and Hippo pathways within a subset of physiological contexts, focusing on normal stem cells and development, as well as regulation of apoptosis, senescence and metabolism in transformed cells.

A DNA Structural Alphabet Distinguishes Structural Features of DNA Bound to Regulatory Proteins and in the Nucleosome Core Particle

We analyzed the structural behavior of DNA complexed with regulatory proteins and the nucleosome core particle (NCP). The three-dimensional structures of almost 25 thousand dinucleotide steps from more than 500 sequentially non-redundant crystal structures were classified by using DNA structural alphabet CANA (Conformational Alphabet of Nucleic Acids) and associations between ten CANA letters and sixteen dinucleotide sequences were investigated. The associations showed features discriminating between specific and non-specific binding of DNA to proteins. Important is the specific role of two DNA structural forms, A-DNA, and BII-DNA, represented by the CANA letters AAA and BB2: AAA structures are avoided in non-specific NCP complexes, where the wrapping of the DNA duplex is explained by the periodic occurrence of BB2 every 10.3 steps. In both regulatory and NCP complexes, the extent of bending of the DNA local helical axis does not influence proportional representation of the CANA alphabet letters, namely the relative incidences of AAA and BB2 remain constant in bent and straight duplexes.

p53 Dynamically Directs TFIID Assembly on Target Gene Promoters

p53 is a central regulator that turns on vast gene networks to maintain cellular integrity in the presence of various stimuli. p53 activates transcription initiation in part by aiding recruitment of TFIID to the promoter. However, the precise means by which p53 dynamically interacts with TFIID to facilitate assembly on target gene promoters remains elusive. To address this key issue, we have undertaken an integrated approach involving single-molecule fluorescence microscopy, single-particle cryo-electron microscopy, and biochemistry. Our real-time single-molecule imaging data demonstrate that TFIID alone binds poorly to native p53 target promoters. p53 unlocks TFIID's ability to bind DNA by stabilizing TFIID contacts with both the core promoter and a region within p53's response element. Analysis of single-molecule dissociation kinetics reveals that TFIID interacts with promoters via transient and prolonged DNA binding modes that are each regulated by p53. Importantly, our structural work reveals that TFIID's conversion to a rearranged DNA binding conformation is enhanced in the presence of DNA and p53. Notably, TFIID's interaction with DNA induces p53 to rapidly dissociate, which likely leads to additional rounds of p53-mediated recruitment of other basal factors. Collectively, these findings indicate that p53 dynamically escorts and loads TFIID onto its target promoters.

Structural Evolution and Dynamics of the p53 Proteins

The family of the p53 tumor suppressive transcription factors includes p73 and p63 in addition to p53 itself. Given the high degree of amino-acid-sequence homology and structural organization shared by the p53 family members, they display some common features (i.e., induction of cell death, cell-cycle arrest, senescence, and metabolic regulation in response to cellular stress) as well as several distinct properties. Here, we describe the structural evolution of the family members with recent advances on the molecular dynamic studies of p53 itself. A crucial role of the carboxy-terminal domain in regulating the properties of the DNA-binding domain (DBD) supports an induced-fit mechanism, in which the binding of p53 on individual promoters is preferentially regulated by the K_OFF over K_ON.

p53 Proteoforms and Intrinsic Disorder: An Illustration of the Protein Structure-Function Continuum Concept

Although it is one of the most studied proteins, p53 continues to be an enigma. This protein has numerous biological functions, possesses intrinsically disordered regions crucial for its functionality, can form both homo-tetramers and isoform-based hetero-tetramers, and is able to interact with many binding partners. It contains numerous posttranslational modifications, has several isoforms generated by alternative splicing, alternative promoter usage or alternative initiation of translation, and is commonly mutated in different cancers. Therefore, p53 serves as an important illustration of the protein structure–function continuum concept, where the generation of multiple proteoforms by various mechanisms defines the ability of this protein to have a multitude of structurally and functionally different states. Considering p53 in the light of a proteoform-based structure–function continuum represents a non-canonical and conceptually new contemplation of structure, regulation, and functionality of this important protein.

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.

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.

Mitochondrial biology, targets, and drug delivery

In recent years, mitochondrial medicine has emerged as a new discipline resting at the intersection of mitochondrial biology, pathology, and pharmaceutics. The central role of mitochondria in critical cellular processes such as metabolism and apoptosis has placed mitochondria at the forefront of cell science. Advances in mitochondrial biology have revealed that these organelles continually undergo fusion and fission while functioning independently and in complex cellular networks, establishing direct membrane contacts with each other and with other organelles. Understanding the diverse cellular functions of mitochondria has contributed to understanding mitochondrial dysfunction in disease states. Polyplasmy and heteroplasmy contribute to mitochondrial phenotypes and associated dysfunction. Residing at the center of cell biology, cellular functions, and disease pathology and being laden with receptors and targets, mitochondria are beacons for pharmaceutical modification. This review presents the current state of mitochondrial medicine with a focus on mitochondrial function, dysfunction, and common disease; mitochondrial receptors, targets, and substrates; and mitochondrial drug design and drug delivery with a focus on the application of nanotechnology to mitochondrial medicine. Mitochondrial medicine is at the precipice of clinical translation; the objective of this review is to aid in the advancement of mitochondrial medicine from infancy to application.