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

JinLiDa granules alleviates cardiac hypertrophy and inflammation in diabetic cardiomyopathy by regulating TP53

BACKGROUND

JinLiDa granules (JLD) is a traditional Chinese medicine (TCM) used to treat type 2 diabetes mellitus with Qi and Yin deficiency. Clinical evidence has shown that JLD can alleviate diabetic cardiomyopathy, but the exact mechanism is not yet clear.

PURPOSE

The purpose of this study was to examine the potential role and mechanism of JLD in the treatment of diabetic cardiomyopathy through network pharmacological analysis and basic experiments.

METHODS

The targets of JLD associated with diabetic cardiomyopathy were examined by network pharmacology. Protein interaction analysis was performed on the targets, and the associated pathways were searched by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Diabetic mice were treated with low or high doses of JLD by gavage, and AC16 and H9C2 cardiomyocytes exposed to high-glucose conditions were treated with JLD. The analysis results were verified by various experimental techniques to examine molecular mechanisms.

RESULTS

Network pharmacological analysis revealed that JLD acted on the tumor suppressor p53 (TP53) during inflammation and fibrosis associated with diabetic cardiomyopathy. The results of basic experiments showed that after JLD treatment, ventricular wall thickening in diabetic mouse hearts was attenuated, cardiac hypertrophy and myocardial inflammation were alleviated, and the expression of cardiac hypertrophy- and inflammation-related factors in cardiomyocytes exposed to a high-glucose environment was decreased. Cardiomyocyte morphology also improved after JLD treatment. TP53 expression and the tumor necrosis factor (TNF) and transforming growth factor beta-1 (TGFβ1) signaling pathways were significantly altered, and inhibiting TP53 expression effectively alleviated the activation of the TNF and TGFβ1 signaling pathways under high glucose conditions. Overexpression of TP53 activated these signaling pathways.

CONCLUSIONS

JLD acted on TP53 to regulate the TNF and TGFβ1 signaling pathways, effectively alleviating cardiomyocyte hypertrophy and inflammation in high glucose and diabetic conditions. Our study provides a solid foundation for the future treatment of diabetic cardiomyopathy with JLD.

Gastric Inhibitory Polypeptide Receptor (GIPR) Overexpression Reduces the Tumorigenic Potential of Retinoblastoma Cells

Retinoblastoma (RB) is the most common malignant intraocular tumor in early childhood. Gene expression profiling revealed that the gastric inhibitory polypeptide receptor (GIPR) is upregulated following trefoil factor family peptide 1 (TFF1) overexpression in RB cells. In the study presented, we found this G protein-coupled transmembrane receptor to be co-expressed with TFF1, a new diagnostic and prognostic RB biomarker for advanced subtype 2 RBs. Functional analyses in two RB cell lines revealed a significant reduction in cell viability and growth and a concomitant increase in apoptosis following stable, lentiviral GIPR overexpression, matching the effects seen after TFF1 overexpression. In chicken chorioallantoic membrane (CAM) assays, GIPR-overexpressing RB cells developed significantly smaller CAM tumors. The effect of GIPR overexpression in RB cells was reversed by the GIPR inhibitor MK0893. The administration of recombinant TFF1 did not augment GIPR overexpression effects, suggesting that GIPR does not serve as a TFF1 receptor. Investigations of potential GIPR up- and downstream mediators suggest the involvement of miR-542-5p and p53 in GIPR signaling. Our results indicate a tumor suppressor role of GIPR in RB, suggesting its pathway as a new potential target for future retinoblastoma therapy.

Orexins in apoptosis: a dual regulatory role

The orexins, also referred to as hypocretins, are neuropeptides that originate from the lateral hypothalamus (LH) region of the brain. They are composed of two small peptides, orexin-A, and orexin-B, which are broadly distributed throughout the central and peripheral nervous systems. Orexins are recognized to regulate diverse functions, involving energy homeostasis, the sleep-wake cycle, stress responses, and reward-seeking behaviors. Additionally, it is suggested that orexin-A deficiency is linked to sleepiness and narcolepsy. The orexins bind to their respective receptors, the orexin receptor type 1 (OX1R) and type 2 (OX2R), and activate different signaling pathways, which results in the mediation of various physiological functions. Orexin receptors are widely expressed in different parts of the body, including the skin, muscles, lungs, and bone marrow. The expression levels of orexins and their receptors play a crucial role in apoptosis, which makes them a potential target for clinical treatment of various disorders. This article delves into the significance of orexins and orexin receptors in the process of apoptosis, highlighting their expression levels and their potential contributions to different diseases. The article offers an overview of the existing understanding of the orexin/receptor system and how it influences the regulation of apoptosis.

Brewing COFFEE: A Sequence-Specific Coarse-Grained Energy Function for Simulations of DNA-Protein Complexes

DNA-protein interactions are pervasive in a number of biophysical processes ranging from transcription and gene expression to chromosome folding. To describe the structural and dynamic properties underlying these processes accurately, it is important to create transferable computational models. Toward this end, we introduce Coarse-grained Force Field for Energy Estimation, COFFEE, a robust framework for simulating DNA-protein complexes. To brew COFFEE, we integrated the energy function in the self-organized polymer model with side-chains for proteins and the three interaction site model for DNA in a modular fashion, without recalibrating any of the parameters in the original force-fields. A unique feature of COFFEE is that it describes sequence-specific DNA-protein interactions using a statistical potential (SP) derived from a data set of high-resolution crystal structures. The only parameter in COFFEE is the strength (λDNAPRO) of the DNA-protein contact potential. For an optimal choice of λDNAPRO, the crystallographic B-factors for DNA-protein complexes with varying sizes and topologies are quantitatively reproduced. Without any further readjustments to the force-field parameters, COFFEE predicts scattering profiles that are in quantitative agreement with small-angle X-ray scattering experiments, as well as chemical shifts that are consistent with NMR. We also show that COFFEE accurately describes the salt-induced unraveling of nucleosomes. Strikingly, our nucleosome simulations explain the destabilization effect of ARG to LYS mutations, which do not alter the balance of electrostatic interactions but affect chemical interactions in subtle ways. The range of applications attests to the transferability of COFFEE, and we anticipate that it would be a promising framework for simulating DNA-protein complexes at the molecular length-scale.

Expanding Roles of the E2F-RB-p53 Pathway in Tumor Suppression

The transcription factor E2F links the RB pathway to the p53 pathway upon loss of function of pRB, thereby playing a pivotal role in the suppression of tumorigenesis. E2F fulfills a major role in cell proliferation by controlling a variety of growth-associated genes. The activity of E2F is controlled by the tumor suppressor pRB, which binds to E2F and actively suppresses target gene expression, thereby restraining cell proliferation. Signaling pathways originating from growth stimulative and growth suppressive signals converge on pRB (the RB pathway) to regulate E2F activity. In most cancers, the function of pRB is compromised by oncogenic mutations, and E2F activity is enhanced, thereby facilitating cell proliferation to promote tumorigenesis. Upon such events, E2F activates the Arf tumor suppressor gene, leading to activation of the tumor suppressor p53 to protect cells from tumorigenesis. ARF inactivates MDM2, which facilitates degradation of p53 through proteasome by ubiquitination (the p53 pathway). P53 suppresses tumorigenesis by inducing cellular senescence or apoptosis. Hence, in almost all cancers, the p53 pathway is also disabled. Here we will introduce the canonical functions of the RB-E2F-p53 pathway first and then the non-classical functions of each component, which may be relevant to cancer biology.

PDBminer to Find and Annotate Protein Structures for Computational Analysis

Computational methods relying on protein structure strongly depend on the structure selected for investigation. Typical sources of protein structures include experimental structures available at the Protein Data Bank (PDB) and high-quality in silico model structures, such as those available at the AlphaFold Protein Structure Database. Either option has significant advantages and drawbacks, and exploring the wealth of available structures to identify the most suitable ones for specific applications can be a daunting task. We provide an open-source software package, PDBminer, with the purpose of making structure identification and selection easier, faster, and less error prone. PDBminer searches the AlphaFold Database and the PDB for available structures of interest and provides an up-to-date, quality-ranked table of structures applicable for further use. PDBminer provides an overview of the available protein structures to one or more input proteins, parallelizing the runs if multiple cores are specified. The output table reports the coverage of the protein structures aligned to the UniProt sequence, overcoming numbering differences in PDB structures and providing information regarding model quality, protein complexes, ligands, and nucleic acid chain binding. The PDBminer2coverage and PDBminer2network tools assist in visualizing the results. PDBminer can be applied to overcome the tedious task of choosing a PDB structure without losing the wealth of additional information available in the PDB. Here, we showcase the main functionalities of the package on the p53 tumor suppressor protein. The package is available at http://github.com/ELELAB/PDBminer.

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

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

Brewing COFFEE: A sequence-specific coarse-grained energy function for simulations of DNA-protein complexes

DNA-protein interactions are pervasive in a number of biophysical processes ranging from transcription, gene expression, to chromosome folding. To describe the structural and dynamic properties underlying these processes accurately, it is important to create transferable computational models. Toward this end, we introduce Coarse grained force field for energy estimation, COFFEE, a robust framework for simulating DNA-protein complexes. To brew COFFEE, we integrated the energy function in the Self-Organized Polymer model with Side Chains for proteins and the Three Interaction Site model for DNA in a modular fashion, without re-calibrating any of the parameters in the original force-fields. A unique feature of COFFEE is that it describes sequence-specific DNA-protein interactions using a statistical potential (SP) derived from a dataset of high-resolution crystal structures. The only parameter in COFFEE is the strength (λ D N A P R O ) of the DNA-protein contact potential. For an optimal choice of λ D N A P R O , the crystallographic B-factors for DNA-protein complexes, with varying sizes and topologies, are quantitatively reproduced. Without any further readjustments to the force-field parameters, COFFEE predicts the scattering profiles that are in quantitative agreement with SAXS experiments as well as chemical shifts that are consistent with NMR. We also show that COFFEE accurately describes the salt-induced unraveling of nucleosomes. Strikingly, our nucleosome simulations explain the destabilization effect of ARG to LYS mutations, which does not alter the balance of electrostatic interactions, but affects chemical interactions in subtle ways. The range of applications attests to the transferability of COFFEE, and we anticipate that it would be a promising framework for simulating DNA-protein complexes at the molecular length-scale.

Targeting p53 pathways: mechanisms, structures, and advances in therapy

The TP53 tumor suppressor is the most frequently altered gene in human cancers, and has been a major focus of oncology research. The p53 protein is a transcription factor that can activate the expression of multiple target genes and plays critical roles in regulating cell cycle, apoptosis, and genomic stability, and is widely regarded as the "guardian of the genome". Accumulating evidence has shown that p53 also regulates cell metabolism, ferroptosis, tumor microenvironment, autophagy and so on, all of which contribute to tumor suppression. Mutations in TP53 not only impair its tumor suppressor function, but also confer oncogenic properties to p53 mutants. Since p53 is mutated and inactivated in most malignant tumors, it has been a very attractive target for developing new anti-cancer drugs. However, until recently, p53 was considered an "undruggable" target and little progress has been made with p53-targeted therapies. Here, we provide a systematic review of the diverse molecular mechanisms of the p53 signaling pathway and how TP53 mutations impact tumor progression. We also discuss key structural features of the p53 protein and its inactivation by oncogenic mutations. In addition, we review the efforts that have been made in p53-targeted therapies, and discuss the challenges that have been encountered in clinical development.

Markov State Modeling Analysis Captures Changes in the Temperature-Sensitive N-Terminal and β-Turn Regions of the p53 DNA-Binding Domain

The transcription factor p53 is one of the most widely studied cancer proteins. Its temperature-sensitive nature suggests reduction in functionality at physiological temperatures. Temperature-induced conformational variations and their impact on its functional ability still remain unexplored. A total of 20.8 μs molecular dynamics simulations of wildtype p53 in the apo and the DNA-bound states have been performed at 300 K and 310 K. Further, Markov State Modeling (MSM) analyses were performed, considering C_α-C_α distances as reaction coordinates. Filtering of these distances based on correlation with the time-independent components (tICs) resulted in 16 and 32 distances for apo and DNA-bound systems, respectively. Individual MSM analyses using these filtered distances were performed for both p53 systems. These C_α-C_α residue pairs belonged to the N-terminal, S6/7 β-turn, loop L2, loop L3, and hydrophobic core residues. At physiological temperatures, apo-p53 exhibits exposure of its hydrophobic core, where the temperature-sensitive hotspot residues were also located. This exposure was the result of the S6/7 β-turn and N-terminal moving apart. In the DNA-bound p53 system, loop L1 attains an open conformation at physiological temperatures, which weakens the DNA binding. It is already known that p53 mutants that lack DNA binding also tend to show similar conformational variations. The S6/7 β-turn along with the already known functionally important loop L2 may pose as regions to be targeted to overcome the loss in binding of temperature-sensitive wildtype p53. Rescue strategies directed toward these temperature-sensitive regions may be useful to recuperate its strong binding at physiological temperatures.

Group Theory of Syntactical Freedom in DNA Transcription and Genome Decoding

Transcription factors (TFs) are proteins that recognize specific DNA fragments in order to decode the genome and ensure its optimal functioning. TFs work at the local and global scales by specifying cell type, cell growth and death, cell migration, organization and timely tasks. We investigate the structure of DNA-binding motifs with the theory of finitely generated groups. The DNA ‘word’ in the binding domain—the motif—may be seen as the generator of a finitely generated group Fdna on four letters, the bases A, T, G and C. It is shown that, most of the time, the DNA-binding motifs have subgroup structures close to free groups of rank three or less, a property that we call ‘syntactical freedom’. Such a property is associated with the aperiodicity of the motif when it is seen as a substitution sequence. Examples are provided for the major families of TFs, such as leucine zipper factors, zinc finger factors, homeo-domain factors, etc. We also discuss the exceptions to the existence of such DNA syntactical rules and their functional roles. This includes the TATA box in the promoter region of some genes, the single-nucleotide markers (SNP) and the motifs of some genes of ubiquitous roles in transcription and regulation.

Most Probable Druggable Pockets in Mutant p53-Arg175His Clusters Extracted from Gaussian Accelerated Molecular Dynamics Simulations

p53, a tumor suppressor protein, is essential for preventing cancer development. Enhancing our understanding of the human p53 function and its modifications in carcinogenesis will aid in developing more highly effective strategies for cancer prevention and treatment. In this study, we have modeled five human p53 forms, namely, inactive, distal-active, proximal-active, distal-Arg175His mutant, and proximal-Arg175His mutant forms. These forms have been investigated using Gaussian accelerated molecular dynamics (GaMD) simulations in OPC water model at physiological temperature and pH. Our observations, obtained throughout [Formula: see text] of production run, are in good agreement with the relevant results in the classical molecular dynamics (MD) studies. Therefore, GaMD method is more economic and efficient method than the classical MD method for studying biomolecular systems. The featured dynamics of the five human p53-DBD forms include noticeable conformational changes of L1 and [Formula: see text]-[Formula: see text] loops as well as [Formula: see text]-[Formula: see text] and [Formula: see text]-[Formula: see text] turns. We have identified two clusters that represent two distinct conformational states in each p53-DBD form. The free-energy profiles of these clusters demonstrate the flexibility of the protein to undergo a conformational transition between the two clusters. We have predicted two out of seven possible druggability pockets on the clusters of the Arg175His forms. These two druggability pockets are near the mutation site and are expected to be actual pockets, which will be helpful for the compound clinical progression studies.

Structural insight into the molecular mechanism of p53-mediated mitochondrial apoptosis

The tumor suppressor p53 is mutated in approximately half of all human cancers. p53 can induce apoptosis through mitochondrial membrane permeabilization by interacting with and antagonizing the anti-apoptotic proteins BCL-xL and BCL-2. However, the mechanisms by which p53 induces mitochondrial apoptosis remain elusive. Here, we report a 2.5 Å crystal structure of human p53/BCL-xL complex. In this structure, two p53 molecules interact as a homodimer, and bind one BCL-xL molecule to form a ternary complex with a 2:1 stoichiometry. Mutations at the p53 dimer interface or p53/BCL-xL interface disrupt p53/BCL-xL interaction and p53-mediated apoptosis. Overall, our current findings of the bona fide structure of p53/BCL-xL complex reveal the molecular basis of the interaction between p53 and BCL-xL, and provide insight into p53-mediated mitochondrial apoptosis.

The structure of human tumor suppressor p53 in complex with the antiapoptotic protein BCL-xL reveals the basis of the p53–BCL-xL interaction and provides insight into the mechanisms of p53-mediated mitochondrial apoptosis.

The complex architecture of p53 binding sites

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

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.

Reactive Oxygen Species Induced p53 Activation: DNA Damage, Redox Signaling, or Both?

Significance: The p53 tumor suppressor has been dubbed the "guardian of genome" because of its various roles in the response to DNA damage such as DNA damage repair, cell cycle arrest, senescence, and apoptosis, all of which are in place to prevent mutations from being passed on down the lineage. Recent Advances: Reactive oxygen species (ROS), for instance hydrogen peroxide derived from mitochondrial respiration, have long been regarded mainly as a major source of cellular damage to DNA and other macromolecules. Critical Issues: More recently, ROS have been shown to also play important physiological roles as second messengers in so-called redox signaling. It is, therefore, not clear whether the observed activation of p53 by ROS is mediated through the DNA damage response, redox signaling, or both. In this review, we will discuss the similarities and differences between p53 activation in response to DNA damage and redox signaling in terms of upstream signaling and downstream transcriptional program activation. Future Directions: Understanding whether and how DNA damage and redox signaling-dependent p53 activation can be dissected could be useful to develop anti-cancer therapeutic p53-reactivation strategies that do not depend on the induction of DNA damage and the resulting additional mutational load.

The roles of TET family proteins in development and stem cells

Ten-eleven translocation (TET) methylcytosine dioxygenases are enzymes that catalyze the demethylation of 5-methylcytosine on DNA. Through global and site-specific demethylation, they regulate cell fate decisions during development and in embryonic stem cells by maintaining pluripotency or by regulating differentiation. In this Primer, we provide an updated overview of TET functions in development and stem cells. We discuss the catalytic and non-catalytic activities of TETs, and their roles as epigenetic regulators of both DNA and RNA hydroxymethylation, highlighting how TET proteins function in regulating gene expression at both the transcriptional and post-transcriptional levels.

Cell proliferation activity delineated by molecular docking of four new compounds isolated from the aerial parts of Suaeda monoica Forssk. ex. J.F. Gmel

Using different chromatographic methods, four new compounds were isolated from the aerial parts of Suaeda monoica (Chenopodiaceae) along with 2-hydroxy-1-naphthoic acid (SCM-3). The structures of the new compounds were established as 6'-hydroxy-10'-geranilanyl naphtha-1-oate (SMC-1), 4,4,8β,10β-Tetramethyl-9β-isobutanyl decalin-13-ol-13-O-β-D-xylopyranoside (SCM-2), 6'-(2-hydroxynaphthalen-3-yl) hexanoic acid (SCM-4) and 1'-(2-Methoxy-3-naphthyl)-4'-(2''-methylbenzoyl)-n-butane (SMC-5) by IR, EIMS and NMR (1 & 2D) analyses. All compounds (50 μg/mL) were tested for cell proliferative potential on cultured human liver cell HepG2 cells by MTT assay. The results revealed a marked cell proliferative potential of all compounds (1.42-1.48 fold) as compared to untreated control. The results of molecular docking and binding with specific proteins such as PTEN (Phosphatase and Tensin homolog) and p53 also justify the cell proliferative potential of the isolated compounds. Glide program with Schrodinger suit 2018 was used to evaluate the binding between SMC compounds and proteins (PTEN and p53). The binding affinity of all compounds was in order of 10⁴-10⁵ M^-1 towards both PTEN and p53. All the SMC compounds have been found to bind at the active site of PTEN, thereby may prevent the binding of phosphatidylinositiol 3,4,5-triphosphate (PI3P). In the locked position, PTEN would not be able to hydrolyze PI3P and hence the PI3P regulated signaling pathway remains active. Similarly, SMC molecules were found to interact with the amino acid residues (Ser99, Thr170, Gly199, and Asp224) which are critically involved in the formation of tetrameric p53. The blockage of p53 to attain its active conformation thus may prevent the recruitment of p53 on DNA and hence may promote cell proliferation.

Mutant p53 drives cancer chemotherapy resistance due to loss of function on activating transcription of PUMA

P53 is a critical tumor suppressor gene, activating p53 and its downstream targets to induce apoptosis is a promising way for cancer therapy. However, more than 50% of cancer patients have p53 mutations, which may cause cancer therapy resistance, and the underline mechanism is poorly understood. Here, we found that cell viability decrease and apoptosis induced by p53-dependent traditional drugs in colon cancer cells were eliminated in p53 mutant cells. Mutant p53 did not up-regulate the expression of its direct downstream targets PUMA and p21, due to the inhibition of PUMA transcription. Furthermore, mutant p53 could not bind to the promoter of PUMA to activate its transcription like WT p53 did, while overexpressed WT p53 rescued PUMA-induced subsequent apoptosis. In conclusion, our findings demonstrate mutant p53 may cause chemo-resistance of tumor because of inactivating PUMA transcription, which prompts some new insights for clinical therapy of cancer patients with mutant p53.Abbreviations: CRC: Colorectal cancer; CDKs: Cyclin-dependent kinases; PUMA: p53 up-regulated modulator of apoptosis; PDGF: the platelet-derived growth factor; WT p53: wild-type p53 protein; mutp53: mutant p53 proteins; BAX: Bcl-2-associated X protein; NOXA: Phorbol-12-myristate-13-acetate-induced protein 1.

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.

GenProBiS: web server for mapping of sequence variants to protein binding sites

Discovery of potentially deleterious sequence variants is important and has wide implications for research and generation of new hypotheses in human and veterinary medicine, and drug discovery. The GenProBiS web server maps sequence variants to protein structures from the Protein Data Bank (PDB), and further to protein–protein, protein–nucleic acid, protein–compound, and protein–metal ion binding sites. The concept of a protein–compound binding site is understood in the broadest sense, which includes glycosylation and other post-translational modification sites. Binding sites were defined by local structural comparisons of whole protein structures using the Protein Binding Sites (ProBiS) algorithm and transposition of ligands from the similar binding sites found to the query protein using the ProBiS-ligands approach with new improvements introduced in GenProBiS. Binding site surfaces were generated as three-dimensional grids encompassing the space occupied by predicted ligands. The server allows intuitive visual exploration of comprehensively mapped variants, such as human somatic mis-sense mutations related to cancer and non-synonymous single nucleotide polymorphisms from 21 species, within the predicted binding sites regions for about 80 000 PDB protein structures using fast WebGL graphics. The GenProBiS web server is open and free to all users at http://genprobis.insilab.org.

Building Transcription Factor Binding Site Models to Understand Gene Regulation in Plants

Transcription factors (TFs) are key cellular components that control gene expression. They recognize specific DNA sequences, the TF binding sites (TFBSs), and thus are targeted to specific regions of the genome where they can recruit transcriptional co-factors and/or chromatin regulators to fine-tune spatiotemporal gene regulation. Therefore, the identification of TFBSs in genomic sequences and their subsequent quantitative modeling is of crucial importance for understanding and predicting gene expression. Here, we review how TFBSs can be determined experimentally, how the TFBS models can be constructed in silico, and how they can be optimized by taking into account features such as position interdependence within TFBSs, DNA shape, and/or by introducing state-of-the-art computational algorithms such as deep learning methods. In addition, we discuss the integration of context variables into the TFBS modeling, including nucleosome positioning, chromatin states, methylation patterns, 3D genome architectures, and TF cooperative binding, in order to better predict TF binding under cellular contexts. Finally, we explore the possibilities of combining the optimized TFBS model with technological advances, such as targeted TFBS perturbation by CRISPR, to better understand gene regulation, evolution, and plant diversity.

The molecular mechanism of action of methylene quinuclidinone and its effects on the structure of p53 mutants

One of the most important tumor suppressor proteins in eukaryotic cells is the transcription factor called p53. The importance of this protein in cells comes from the fact that it regulates a wide variety of cellular processes including the cell cycle, metabolism, DNA repair, senescence and apoptosis. In cancer cells, p53 is a major target as the most mutated protein, which has led to the search for potential activators of the mutant protein. Currently, the only mutated-p53 activator in clinical trials is a small molecule called APR-246. There is evidence that the active metabolite of APR-246 binds covalently to mutant p53 and restores its wild-type (wt) activity. In this work, we created atomistic in silico models of the wt, mutant and drugged mutant p53 proteins each in complex with DNA. Using molecular dynamics simulations we generated equilibrated models of the complexes. Detailed analysis revealed that the binding of the APR-246 active metabolite to the mutant proteins alters their interaction with DNA. In particular, the binding of the molecule at loop L1 of the protein allows the loop to anchor the protein to DNA similarly to wt p53. Several important p53-DNA interactions lost due to mutation were also restored in the drugged mutants. These findings, not only provide a possible mechanism of action of this drug, but also criteria to use in virtual screening campaigns for other p53 activators.

Virtual screening using covalent docking to find activators for G245S mutant p53

TP53 is the most mutated gene in all cancers. The mutant protein also accumulates in cells. The high frequency of p53 mutations makes the protein a promising target for anti-cancer therapy. Only a few molecules have been found, using in vitro screening, to reactivate the mutant protein. APR-246 is currently the most successful mutant p53 activator, which reactivates the transcriptional activity of p53 by covalently binding to C124 of the protein. We have recently created in silico models of G245S-mp53 in its apo and DNA-bound forms. In this paper we further report on our in silico screening for potential activators of G245S-mp53. We filtered the ZINC15 database (13 million compounds) to only include drug-like molecules with moderate to standard reactivity. Our filtered database of 130,000 compounds was screened using the DOCKTITE protocol in the Molecular Operating Environment software. We performed covalent docking at C124 of G245S-mp53 to identify potential activators of the mutant protein. The docked compounds were ranked using a consensus scoring approach. We also used ADMET Predictor™ to predict pharmacokinetics and the possible toxicities of the compounds. Our screening procedure has identified compounds, mostly thiosemicarbazones and halo-carbonyls, with the best potential as G245S-mp53 activators, which are described in this work. Based on its binding scores and ADMET risk score, compound 2 is likely to have the best potential as a G245S-mp53 activator compared to the other top hits.

DelPhiForce web server: electrostatic forces and energy calculations and visualization

Summary

Electrostatic force is an essential component of the total force acting between atoms and macromolecules. Therefore, accurate calculations of electrostatic forces are crucial for revealing the mechanisms of many biological processes. We developed a DelPhiForce web server to calculate and visualize the electrostatic forces at molecular level. DelPhiForce web server enables modeling of electrostatic forces on individual atoms, residues, domains and molecules, and generates an output that can be visualized by VMD software. Here we demonstrate the usage of the server for various biological problems including protein-cofactor, domain-domain, protein-protein, protein-DNA and protein-RNA interactions.

Availability and implementation

The DelPhiForce web server is available at: http://compbio.clemson.edu/delphi-force.

Contact

delphi@clemson.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

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

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

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

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

Expanding the repertoire of DNA shape features for genome-scale studies of transcription factor binding

Uncovering the mechanisms that affect the binding specificity of transcription factors (TFs) is critical for understanding the principles of gene regulation. Although sequence-based models have been used successfully to predict TF binding specificities, we found that including DNA shape information in these models improved their accuracy and interpretability. Previously, we developed a method for modeling DNA binding specificities based on DNA shape features extracted from Monte Carlo (MC) simulations. Prediction accuracies of our models, however, have not yet been compared to accuracies of models incorporating DNA shape information extracted from X-ray crystallography (XRC) data or Molecular Dynamics (MD) simulations. Here, we integrated DNA shape information extracted from MC or MD simulations and XRC data into predictive models of TF binding and compared their performance. Models that incorporated structural information consistently showed improved performance over sequence-based models regardless of data source. Furthermore, we derived and validated nine additional DNA shape features beyond our original set of four features. The expanded repertoire of 13 distinct DNA shape features, including six intra-base pair and six inter-base pair parameters and minor groove width, is available in our R/Bioconductor package DNAshapeR and enables a comprehensive structural description of the double helix on a genome-wide scale.

Bonding of Butylparaben, Bis(2-ethylhexyl)-phthalate, and Perfluorooctanesulfonic Acid to DNA: Comparison with Benzo[a]pyrene Shows Low Probability for Strong Noncovalent DNA Intercalation

Parabens, phthalates, and perfluorinated compounds are pollutant compounds used in cosmetics, plastics, and fire-fighting foams. All three compounds have been studied over several years for toxicity mechanism; however, a clear view of their ability to bind to DNA has not been supplied empirically. In this work, a simulation study is done to reveal the interaction of three of these pollutants, bis(2-ethylhexyl)-phthalate (DEHP), butylparaben (BPRB), and the protonated form of perfluorooctanesulfonic acid (PFOS(H)), with DNA. The results show that the DEHP, PFOS(H), and BPRB bind with a probability of 1/5 to DNA, with respective bonding energies -23.96 kJ/mol (PFOS(H)), -94.92 kJ/mol (BPRB), and -216.52 kJ/mol (DEHP). The positive control, benzo[a]pyrene diol epoxide (BAP), which is known for its notorious DNA intercalation, binds at a rate of 3/5 simulations, with bonding energies of -6544.52, -7034.66, and -7578.67 kJ/mol. The results are compared to empirical studies and conclusively show that all these pollutants can interfere with transcription and DNA related mechanisms by forming noncovalent interactions with DNA. The results show also that these pollutants are unlikely to undergo strong noncovalent intercalation to DNA, such as BAP, and do not possess the frontier orbital profiles to undergo adduct formation. After many years of research and several unanswered questions on the action of these pollutants on DNA, a calculation on their properties hence to the DNA confirms that there is a low probability for these to undergo a strong intercalation with DNA. Literature shows however that the pollutants are strongly interfering with the protein machinery and receptors on the cell surface and are therefore still priority pollutants for ecotoxicity research.

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.

Cancer-associated arginine-to-histidine mutations confer a gain in pH sensing to mutant proteins

The intracellular pH (pHi) of most cancers is constitutively higher than that of normal cells and enhances proliferation and cell survival. We found that increased pHi enabled the tumorigenic behaviors caused by somatic arginine-to-histidine mutations, which are frequent in cancer and confer pH sensing not seen with wild-type proteins. Experimentally raising the pHi increased the activity of R776H mutant epidermal growth factor receptor (EGFR-R776H), thereby increasing proliferation and causing transformation in fibroblasts. An Arg-to-Gly mutation did not confer these effects. Molecular dynamics simulations of EGFR suggested that decreased protonation of His⁷⁷⁶ at high pH causes conformational changes in the αC helix that may stabilize the active form of the kinase. An Arg-to-His, but not Arg-to-Lys, mutation in the transcription factor p53 (p53-R273H) decreased its transcriptional activity and attenuated the DNA damage response in fibroblasts and breast cancer cells with high pHi. Lowering pHi attenuated the tumorigenic effects of both EGFR-R776H and p53-R273H. Our data suggest that some somatic mutations may confer a fitness advantage to the higher pHi of cancer cells.

Insights into the Effect of the G245S Single Point Mutation on the Structure of p53 and the Binding of the Protein to DNA

The transcription factor p53 is a potent tumor suppressor dubbed as the "guardian of the genome" because of its ability to orchestrate protective biological outputs in response to a variety of oncogenic stresses. Mutation and thus inactivation of p53 can be found in 50% of human tumors. The majority are missense mutations located in the DNA binding region. Among them, G245S is known to be a structural hotspot mutation. To understand the behaviors and differences between the wild-type and mutant, both a dimer of the wild type p53 (wt-p53) and its G245S mutant (G245S-mp53), complexed with DNA, were simulated using molecular dynamics for more than 1 μs. wt-p53 and G245S-mp53 apo monomers were simulated for 1 μs as well. Conformational analyses and binding energy evaluations performed underline important differences and therefore provide insights to understand the G245S-mp53 loss of function. Our results indicate that the G245S mutation destabilizes several structural regions in the protein that are crucial for DNA binding when found in its apo form and highlight differences in the mutant-DNA complex structure compared to the wt protein. These findings not only provide means that can be applied to other p53 mutants but also serve as structural basis for further studies aimed at the development of cancer therapies based on restoring the function of p53.

Proteins Recognizing DNA: Structural Uniqueness and Versatility of DNA-Binding Domains in Stem Cell Transcription Factors

Proteins in the form of transcription factors (TFs) bind to specific DNA sites that regulate cell growth, differentiation, and cell development. The interactions between proteins and DNA are important toward maintaining and expressing genetic information. Without knowing TFs structures and DNA-binding properties, it is difficult to completely understand the mechanisms by which genetic information is transferred between DNA and proteins. The increasing availability of structural data on protein-DNA complexes and recognition mechanisms provides deeper insights into the nature of protein-DNA interactions and therefore, allows their manipulation. TFs utilize different mechanisms to recognize their cognate DNA (direct and indirect readouts). In this review, we focus on these recognition mechanisms as well as on the analysis of the DNA-binding domains of stem cell TFs, discussing the relative role of various amino acids toward facilitating such interactions. Unveiling such mechanisms will improve our understanding of the molecular pathways through which TFs are involved in repressing and activating gene expression.

Two different, mutually exclusively distributed, TP53 mutations in ovarian and peritoneal tumor tissues of a serous ovarian cancer patient: indicative for tumor origin?

High-grade serous ovarian cancer (HGSOC) is characterized by a TP53 mutation rate of up to 96.7% and associated with a more aggressive tumor biology. The origin of HGSOC is thought to arise either from fallopian tube secretory cells or the ovarian surface epithelium/inclusion cysts, the former with more evidence. Peritoneal tumor spread is heterogeneous, either excessive in the peritoneum (with miliary appearance) or more confined to the ovaries with only few (bigger and exophytically growing) peritoneal implants. Using RNA sequencing and DNA digital droplet polymerase chain reaction (PCR), we identified two different functional TP53 mutations in one HGSOC patient: one exclusively in the ovarian tumor mass and the other exclusively in ascites tumor cells, peritoneal tumor masses, and a lymph node metastasis. In blood, both mutations could be detected, the one from the peritoneal tumors with much higher frequency, presumably because of the higher tumor load. We conclude that this mutually exclusive distribution of two different TP53 mutations in different tumor tissues indicates the development of two independent carcinomas in the peritoneal cavity, probably one originating from a precancerous lesion in the fallopian tube and the other from the ovaries. In addition, in the patient's ascites CD45 and EpCAM, double-positive cells were found-proliferating but testing negative for the above-mentioned TP53 mutations. This mutually exclusive distribution of two TP53 mutations is probably further evidence that HGSOC can originate either from the fallopian tube or (more seldom) the ovaries, the former more prone for excessive peritoneal tumor spread.

Protein-induced DNA linking number change by sequence-specific DNA binding proteins and its biological effects

Sequence-specific DNA-binding proteins play essential roles in many fundamental biological events such as DNA replication, recombination, and transcription. One common feature of sequence-specific DNA-binding proteins is to introduce structural changes to their DNA recognition sites including DNA-bending and DNA linking number change (ΔLk). In this article, I review recent progress in studying protein-induced ΔLk by several sequence-specific DNA-binding proteins, such as E. coli cAMP receptor protein (CRP) and lactose repressor (LacI). It was demonstrated recently that protein-induced ΔLk is an intrinsic property for sequence-specific DNA-binding proteins and does not correlate to protein-induced other structural changes, such as DNA bending. For instance, although CRP bends its DNA recognition site by 90°, it was not able to introduce a ΔLk to it. However, LacI was able to simultaneously bend and introduce a ΔLk to its DNA binding sites. Intriguingly, LacI also constrained superhelicity within LacI-lac O1 complexes if (-) supercoiled DNA templates were provided. I also discuss how protein-induced ΔLk help sequence-specific DNA-binding proteins regulate their biological functions. For example, it was shown recently that LacI utilizes the constrained superhelicity (ΔLk) in LacI-lac O1 complexes and serves as a topological barrier to constrain free, unconstrained (-) supercoils within the 401-bp DNA loop. These constrained (-) supercoils enhance LacI's binding affinity and therefore the repression of the lac promoter. Other biological functions include how DNA replication initiators λ O and DnaA use the induced ΔLk to open/melt bacterial DNA replication origins.

Design of a molecular support for cryo-EM structure determination

Despite the recent rapid progress in cryo-electron microscopy (cryo-EM), there still exist ample opportunities for improvement in sample preparation. Macromolecular complexes may disassociate or adopt nonrandom orientations against the extended air-water interface that exists for a short time before the sample is frozen. We designed a hollow support structure using 3D DNA origami to protect complexes from the detrimental effects of cryo-EM sample preparation. For a first proof-of-principle, we concentrated on the transcription factor p53, which binds to specific DNA sequences on double-stranded DNA. The support structures spontaneously form monolayers of preoriented particles in a thin film of water, and offer advantages in particle picking and sorting. By controlling the position of the binding sequence on a single helix that spans the hollow support structure, we also sought to control the orientation of individual p53 complexes. Although the latter did not yet yield the desired results, the support structures did provide partial information about the relative orientations of individual p53 complexes. We used this information to calculate a tomographic 3D reconstruction, and refined this structure to a final resolution of ∼15 Å. This structure settles an ongoing debate about the symmetry of the p53 tetramer bound to DNA.

Protein-induced DNA linking number change by sequence-specific DNA binding proteins and its biological effects

Sequence-specific DNA-binding proteins play essential roles in many fundamental biological events such as DNA replication, recombination, and transcription. One common feature of sequence-specific DNA-binding proteins is to introduce structural changes to their DNA recognition sites including DNA-bending and DNA linking number change (ΔLk). In this article, I review recent progress in studying protein-induced ΔLk by several sequence-specific DNA-binding proteins, such as E. coli cAMP receptor protein (CRP) and lactose repressor (LacI). It was demonstrated recently that protein-induced ΔLk is an intrinsic property for sequence-specific DNA-binding proteins and does not correlate to protein-induced other structural changes, such as DNA bending. For instance, although CRP bends its DNA recognition site by 90°, it was not able to introduce a ΔLk to it. However, LacI was able to simultaneously bend and introduce a ΔLk to its DNA binding sites. Intriguingly, LacI also constrained superhelicity within LacI-lac O1 complexes if (-) supercoiled DNA templates were provided. I also discuss how protein-induced ΔLk help sequence-specific DNA-binding proteins regulate their biological functions. For example, it was shown recently that LacI utilizes the constrained superhelicity (ΔLk) in LacI-lac O1 complexes and serves as a topological barrier to constrain free, unconstrained (-) supercoils within the 401-bp DNA loop. These constrained (-) supercoils enhance LacI's binding affinity and therefore the repression of the lac promoter. Other biological functions include how DNA replication initiators λ O and DnaA use the induced ΔLk to open/melt bacterial DNA replication origins.

Functional Diversification after Gene Duplication: Paralog Specific Regions of Structural Disorder and Phosphorylation in p53, p63, and p73

Conformational and functional flexibility promote protein evolvability. High evolvability allows related proteins to functionally diverge and perhaps to neostructuralize. p53 is a multifunctional protein frequently referred to as the Guardian of the Genome–a hub for e.g. incoming and outgoing signals in apoptosis and DNA repair. p53 has been found to be structurally disordered, an extreme form of conformational flexibility. Here, p53, and its paralogs p63 and p73, were studied for further insights into the evolutionary dynamics of structural disorder, secondary structure, and phosphorylation. This study is focused on the post gene duplication phase for the p53 family in vertebrates, but also visits the origin of the protein family and the early domain loss and gain events. Functional divergence, measured by rapid evolutionary dynamics of protein domains, structural properties, and phosphorylation propensity, is inferred across vertebrate p53 proteins, in p63 and p73 from fish, and between the three paralogs. In particular, structurally disordered regions are redistributed among paralogs, but within clades redistribution of structural disorder also appears to be an ongoing process. Despite its deemed importance as the Guardian of the Genome, p53 is indeed a protein with high evolvability as seen not only in rearranged structural disorder, but also in fluctuating domain sequence signatures among lineages.

Binding of transcription factor GabR to DNA requires recognition of DNA shape at a location distinct from its cognate binding site

Mechanisms for transcription factor recognition of specific DNA base sequences are well characterized and recent studies demonstrate that the shape of these cognate binding sites is also important. Here, we uncover a new mechanism where the transcription factor GabR simultaneously recognizes two cognate binding sites and the shape of a 29 bp DNA sequence that bridges these sites. Small-angle X-ray scattering and multi-angle laser light scattering are consistent with a model where the DNA undergoes a conformational change to bend around GabR during binding. In silico predictions suggest that the bridging DNA sequence is likely to be bendable in one direction and kinetic analysis of mutant DNA sequences with biolayer interferometry, allowed the independent quantification of the relative contribution of DNA base and shape recognition in the GabR-DNA interaction. These indicate that the two cognate binding sites as well as the bendability of the DNA sequence in between these sites are required to form a stable complex. The mechanism of GabR-DNA interaction provides an example where the correct shape of DNA, at a clearly distinct location from the cognate binding site, is required for transcription factor binding and has implications for bioinformatics searches for novel binding sites.

Conformations of Human Telomeric G-Quadruplex Studied Using a Nucleotide-Independent Nitroxide Label

Guanine-rich oligonucleotides can form a unique G-quadruplex (GQ) structure with stacking units of four guanine bases organized in a plane through Hoogsteen bonding. GQ structures have been detected in vivo and shown to exert their roles in maintaining genome integrity and regulating gene expression. Understanding GQ conformation is important for understanding its inherent biological role and for devising strategies to control and manipulate functions based on targeting GQ. Although a number of biophysical methods have been used to investigate structure and dynamics of GQs, our understanding is far from complete. As such, this work explores the use of the site-directed spin labeling technique, complemented by molecular dynamics simulations, for investigating GQ conformations. A nucleotide-independent nitroxide label (R5), which has been previously applied for probing conformations of noncoding RNA and DNA duplexes, is attached to multiple sites in a 22-nucleotide DNA strand derived from the human telomeric sequence (hTel-22) that is known to form GQ. The R5 labels are shown to minimally impact GQ folding, and inter-R5 distances measured using double electron-electron resonance spectroscopy are shown to adequately distinguish the different topological conformations of hTel-22 and report variations in their occupancies in response to changes of the environment variables such as salt, crowding agent, and small molecule ligand. The work demonstrates that the R5 label is able to probe GQ conformation and establishes the base for using R5 to study more complex sequences, such as those that may potentially form multimeric GQs in long telomeric repeats.

p53 dynamics upon response element recognition explored by molecular simulations

p53 is a representative transcription factor that activates multiple target genes. To realize stimulus-dependent specificities, p53 has to recognize targets with structural variety, of which molecular mechanisms are largely unknown. Here, we conducted a series of long-time scale (totally more than 100-ms) coarse-grained molecular dynamics simulations, uncovering structure and dynamics of full-length p53 tetramer that recognizes its response element (RE). We obtained structures of a full-length p53 tetramer that binds to the RE, which is strikingly different from the structure of p53 at search. These structures are not only consistent with previous low-resolution or partial structural information, but also give access to previously unreachable detail, such as the preferential distribution of intrinsically disordered regions, the contacts between core domains, the DNA bending, and the connectivity of linker regions. We also explored how the RE variation affects the structure of the p53-RE complex. Further analysis of simulation trajectories revealed how p53 finds out the RE and how post-translational modifications affect the search mechanism.

Significant Differences in the Development of Acquired Resistance to the MDM2 Inhibitor SAR405838 between In Vitro and In Vivo Drug Treatment

SAR405838 is a potent and specific MDM2 inhibitor currently being evaluated in Phase I clinical trials for the treatment of human cancer. Using the SJSA-1 osteosarcoma cell line which harbors an amplified MDM2 gene and wild-type p53, we have investigated the acquired resistance mechanisms both in vitro and in vivo to SAR405838. Treatment of SJSA-1 cells with SAR405838 in vitro leads to dose-dependent cell growth inhibition, cell cycle arrest and robust apoptosis. However, prolonged treatment of SJSA-1 cells in vitro with SAR405838 results in profound acquired resistance to the drug. Analysis of in vitro-derived resistant cell lines showed that p53 is mutated in the DNA binding domain and can no longer be activated by SAR405838. Treatment of the parental SJSA-1 xenograft tumors with SAR405838 in mice yields rapid tumor regression but the tumors eventually regrow. Culturing the regrown tumors established a number of sublines, which showed only modest (3–5 times) loss of sensitivity to SAR405838 in vitro. Sequencing of the p53 showed that it retains its wild-type status in these in vivo sublines, with the exception of one subline, which harbors a single heterozygous C176F p53 mutation. Using xenograft models of two in vivo derived sublines, which has either wild-type p53 or p53 containing a single heterozygous C176F mutation, we showed that while SAR405838 effectively achieves partial tumor regression in these models, it no longer induces complete tumor regression and tumors resume growth once the treatment is stopped. Harvesting and culturing tumors obtained from a prolonged treatment with SAR405838 in mice established additional in vivo sublines, which all contain a single heterozygous C176F mutation with no additional p53 mutation detected. Interestingly, SAR405838 can still effectively activate p53 in all sublines containing a single heterozygous C176F mutation, with a moderately reduced potency as compared to that in the parental cell line. Consistently, SAR405838 is 3–5 times less effective in all the in vivo derived sublines containing a single heterozygous C176F p53 mutation than in the SJSA-1 parental cell line in assays of cell growth and apoptosis. Computational modeling suggested that a p53 tetramer containing two wild-type p53 molecules and two C176F mutated molecules can maintain the structural stability and interactions with DNA by formation of additional hydrophobic and cation-π interactions which compensate for the loss of sulphur-zinc coordination. Our data thus show that SJSA-1 tumor cells acquire very different levels of resistance in vitro and in vivo to the MDM2 inhibitor SAR405838. Our present study may have a significant implication for the investigation of resistant mechanisms for other classes of anticancer drugs.

FUSE Binding Protein 1 Facilitates Persistent Hepatitis C Virus Replication in Hepatoma Cells by Regulating Tumor Suppressor p53

Hepatitis C virus (HCV) is a leading cause of chronic hepatitis C (CHC), liver cirrhosis, and hepatocellular carcinoma (HCC). Immunohistochemistry of archived HCC tumors showed abundant FBP1 expression in HCC tumors with the CHC background. Oncomine data analysis of normal versus HCC tumors with the CHC background indicated a 4-fold increase in FBP1 expression with a concomitant 2.5-fold decrease in the expression of p53. We found that FBP1 promotes HCV replication by inhibiting p53 and regulating BCCIP and TCTP, which are positive and negative regulators of p53, respectively. The severe inhibition of HCV replication in FBP1-knockdown Huh7.5 cells was restored to a normal level by downregulation of either p53 or BCCIP. Although p53 in Huh7.5 cells is transcriptionally inactive as a result of Y220C mutation, we found that the activation and DNA binding ability of Y220C p53 were strongly suppressed by FBP1 but significantly activated upon knockdown of FBP1. Transient expression of FBP1 in FBP1 knockdown cells fully restored the control phenotype in which the DNA binding ability of p53 was strongly suppressed. Using electrophoretic mobility shift assay (EMSA) and isothermal titration calorimetry (ITC), we found no significant difference in in vitro target DNA binding affinity of recombinant wild-type p53 and its Y220C mutant p53. However, in the presence of recombinant FBP1, the DNA binding ability of p53 is strongly inhibited. We confirmed that FBP1 downregulates BCCIP, p21, and p53 and upregulates TCTP under radiation-induced stress. Since FBP1 is overexpressed in most HCC tumors with an HCV background, it may have a role in promoting persistent virus infection and tumorigenesis.

IMPORTANCE

It is our novel finding that FUSE binding protein 1 (FBP1) strongly inhibits the function of tumor suppressor p53 and is an essential host cell factor required for HCV replication. Oncomine data analysis of a large number of samples has revealed that overexpression of FBP1 in most HCC tumors with chronic hepatitis C is significantly linked with the decreased expression level of p53. The most significant finding is that FBP1 not only physically interacts with p53 and interferes with its binding to the target DNA but also functions as a negative regulator of p53 under cellular stress. FBP1 is barely detectable in normal differentiated cells; its overexpression in HCC tumors with the CHC background suggests that FBP1 has an important role in promoting HCV infection and HCC tumors by suppressing p53.

High-resolution profiling of Drosophila replication start sites reveals a DNA shape and chromatin signature of metazoan origins

At every cell cycle, faithful inheritance of metazoan genomes requires the concerted activation of thousands of DNA replication origins. However, the genetic and chromatin features defining metazoan replication start sites remain largely unknown. Here, we delineate the origin repertoire of the Drosophila genome at high resolution. We address the role of origin-proximal G-quadruplexes and suggest that they transiently stall replication forks in vivo. We dissect the chromatin configuration of replication origins and identify a rich spatial organization of chromatin features at initiation sites. DNA shape and chromatin configurations, not strict sequence motifs, mark and predict origins in higher eukaryotes. We further examine the link between transcription and origin firing and reveal that modulation of origin activity across cell types is intimately linked to cell-type-specific transcriptional programs. Our study unravels conserved origin features and provides unique insights into the relationship among DNA topology, chromatin, transcription, and replication initiation across metazoa.

Quantitative modeling of transcription factor binding specificities using DNA shape

DNA binding specificities of transcription factors (TFs) are a key component of gene regulatory processes. Underlying mechanisms that explain the highly specific binding of TFs to their genomic target sites are poorly understood. A better understanding of TF-DNA binding requires the ability to quantitatively model TF binding to accessible DNA as its basic step, before additional in vivo components can be considered. Traditionally, these models were built based on nucleotide sequence. Here, we integrated 3D DNA shape information derived with a high-throughput approach into the modeling of TF binding specificities. Using support vector regression, we trained quantitative models of TF binding specificity based on protein binding microarray (PBM) data for 68 mammalian TFs. The evaluation of our models included cross-validation on specific PBM array designs, testing across different PBM array designs, and using PBM-trained models to predict relative binding affinities derived from in vitro selection combined with deep sequencing (SELEX-seq). Our results showed that shape-augmented models compared favorably to sequence-based models. Although both k-mer and DNA shape features can encode interdependencies between nucleotide positions of the binding site, using DNA shape features reduced the dimensionality of the feature space. In addition, analyzing the feature weights of DNA shape-augmented models uncovered TF family-specific structural readout mechanisms that were not revealed by the DNA sequence. As such, this work combines knowledge from structural biology and genomics, and suggests a new path toward understanding TF binding and genome function.

Oxidation of p53 through DNA charge transport involves a network of disulfides within the DNA-binding domain

Transcription factor p53 plays a critical role in the cellular response to stress stimuli. We have seen that p53 dissociates selectively from various promoter sites as a result of oxidation at long-range through DNA-mediated charge transport (CT). Here, we examine this chemical oxidation and determine the residues in p53 that are essential for oxidative dissociation, focusing on the network of cysteine residues adjacent to the DNA-binding site. Of the eight mutants studied, only the C275S mutation shows decreased affinity for the Gadd45 promoter site. However, both mutations C275S and C277S result in substantial attenuation of oxidative dissociation, with C275S causing the most severe attenuation. Differential thiol labeling was used to determine the oxidation states of cysteine residues within p53 after DNA-mediated oxidation. Reduced cysteines were iodoacetamide-labeled, whereas oxidized cysteines participating in disulfide bonds were (13)C2D2-iodoacetamide-labeled. Intensities of respective iodoacetamide-modified peptide fragments were analyzed by mass spectrometry. A distinct shift in peptide labeling toward (13)C2D2-iodoacetamide-labeled cysteines is observed in oxidized samples, confirming that chemical oxidation of p53 occurs at long range. All observable cysteine residues trend toward the heavy label under conditions of DNA CT, indicating the formation of multiple disulfide bonds among the cysteine network. On the basis of these data, it is proposed that disulfide formation involving C275 is critical for inducing oxidative dissociation of p53 from DNA.

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.

DNA binding by FOXP3 domain-swapped dimer suggests mechanisms of long-range chromosomal interactions

FOXP3 is a lineage-specific transcription factor that is required for regulatory T cell development and function. In this study, we determined the crystal structure of the FOXP3 forkhead domain bound to DNA. The structure reveals that FOXP3 can form a stable domain-swapped dimer to bridge DNA in the absence of cofactors, suggesting that FOXP3 may play a role in long-range gene interactions. To test this hypothesis, we used circular chromosome conformation capture coupled with high throughput sequencing (4C-seq) to analyze FOXP3-dependent genomic contacts around a known FOXP3-bound locus, Ptpn22. Our studies reveal that FOXP3 induces significant changes in the chromatin contacts between the Ptpn22 locus and other Foxp3-regulated genes, reflecting a mechanism by which FOXP3 reorganizes the genome architecture to coordinate the expression of its target genes. Our results suggest that FOXP3 mediates long-range chromatin interactions as part of its mechanisms to regulate specific gene expression in regulatory T cells.