DNA structure determines protein binding and transcriptional efficiency of the proenkephalin cAMP-responsive enhancer

The AP-1 transcriptional complex: Local switch or remote command?

The ubiquitous family of AP-1 dimeric transcription complexes is involved in virtually all cellular and physiological functions. It is paramount for cells to reprogram gene expression in response to cues of many sorts and is involved in many tumorigenic processes. How AP-1 controls gene transcription has largely remained elusive till recently. The advent of the "omics" technologies permitting genome-wide studies of transcription factors has however changed and improved our view of AP-1 mechanistical actions. If these studies confirm that AP-1 can sometimes act as a local transcriptional switch operating in the vicinity of transcription start sites (TSS), they strikingly indicate that AP-1 principally operates as a remote command binding to distal enhancers, placing chromatin architecture dynamics at the heart of its transcriptional actions. They also unveil novel constraints operating on AP-1, as well as novel mechanisms used to regulate gene expression via transcription-pioneering-, chromatin-remodeling- and chromatin accessibility maintenance effects.

Comparison of discriminative motif optimization using matrix and DNA shape-based models

Background

Transcription factor (TF) binding site specificity is commonly represented by some form of matrix model in which the positions in the binding site are assumed to contribute independently to the site’s activity. The independence assumption is known to be an approximation, often a good one but sometimes poor. Alternative approaches have been developed that use k-mers (DNA “words” of length k) to account for the non-independence, and more recently DNA structural parameters have been incorporated into the models. ChIP-seq data are often used to assess the discriminatory power of motifs and to compare different models. However, to measure the improvement due to using more complex models, one must compare to optimized matrix models.

Results

We describe a program “Discriminative Additive Model Optimization” (DAMO) that uses positive and negative examples, as in ChIP-seq data, and finds the additive position weight matrix (PWM) that maximizes the Area Under the Receiver Operating Characteristic Curve (AUROC). We compare to a recent study where structural parameters, serving as features in a gradient boosting classifier algorithm, are shown to improve the AUROC over JASPAR position frequency matrices (PFMs). In agreement with the previous results, we find that adding structural parameters gives the largest improvement, but most of the gain can be obtained by an optimized PWM and nearly all of the gain can be obtained with a di-nucleotide extension to the PWM.

Conclusion

To appropriately compare different models for TF bind sites, optimized models must be used. PWMs and their extensions are good representations of binding specificity for most TFs, and more complex models, including the incorporation of DNA shape features and gradient boosting classifiers, provide only moderate improvements for a few TFs.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2104-7) contains supplementary material, which is available to authorized users.

High Level of Fasting Plasma Proenkephalin-A Predicts Deterioration of Kidney Function and Incidence of CKD

High levels of proenkephalin-A (pro-ENK) have been associated with decreased eGFR in an acute setting. Here, we examined whether pro-ENK levels predict CKD and decline of renal function in a prospective cohort of 2568 participants without CKD (eGFR>60 ml/min per 1.73 m²) at baseline. During a mean follow-up of 16.6 years, 31.7% of participants developed CKD. Participants with baseline pro-ENK levels in the highest tertile had significantly greater yearly mean decline of eGFR (P_trend<0.001) and rise of cystatin C (P_trend=0.01) and creatinine (P_trend<0.001) levels. Furthermore, compared with participants in the lowest tertile, participants in the highest tertile of baseline pro-ENK concentration had increased CKD incidence (odds ratio, 1.51; 95% confidence interval, 1.18 to 1.94) when adjusted for multiple factors. Adding pro-ENK to a model of conventional risk factors in net reclassification improvement analysis resulted in reclassification of 14.14% of participants. Genome-wide association analysis in 4150 participants of the same cohort revealed the strongest association of pro-ENK levels with rs1012178 near the PENK gene, where the minor T-allele associated with a 0.057 pmol/L higher pro-ENK level per allele (P=4.67x10^-21). Furthermore, the T-allele associated with a 19% increased risk of CKD per allele (P=0.03) and a significant decrease in the instrumental variable estimator for eGFR (P<0.01) in a Mendelian randomization analysis. In conclusion, circulating plasma pro-ENK level predicts incident CKD and may aid in identifying subjects in need of primary preventive regimens. Additionally, the Mendelian randomization analysis suggests a causal relationship between pro-ENK level and deterioration of kidney function over time.

A genome-wide screen of CREB occupancy identifies the RhoA inhibitors Par6C and Rnd3 as regulators of BDNF-induced synaptogenesis

Neurotrophin-regulated gene expression is believed to play a key role in long-term changes in synaptic structure and the formation of dendritic spines. Brain-derived neurotrophic factor (BDNF) has been shown to induce increases in dendritic spine formation, and this process is thought to function in part by stimulating CREB-dependent transcriptional changes. To identify CREB-regulated genes linked to BDNF-induced synaptogenesis, we profiled transcriptional occupancy of CREB in hippocampal neurons. Interestingly, de novo motif analysis of hippocampal ChIP-Seq data identified a non-canonical CRE motif (TGGCG) that was enriched at CREB target regions and conferred CREB-responsiveness. Because cytoskeletal remodeling is an essential element of the formation of dendritic spines, within our screens we focused our attention on genes previously identified as inhibitors of RhoA GTPase. Bioinformatic analyses identified dozens of candidate CREB target genes known to regulate synaptic architecture and function. We showed that two of these, the RhoA inhibitors Par6C (Pard6A) and Rnd3 (RhoE), are BDNF-induced CREB-regulated genes. Interestingly, CREB occupied a cluster of non-canonical CRE motifs in the Rnd3 promoter region. Lastly, we show that BDNF-stimulated synaptogenesis requires the expression of Par6C and Rnd3, and that overexpression of either protein is sufficient to increase synaptogenesis. Thus, we propose that BDNF can regulate formation of functional synapses by increasing the expression of the RhoA inhibitors, Par6C and Rnd3. This study shows that genome-wide analyses of CREB target genes can facilitate the discovery of new regulators of synaptogenesis.

Energetic coupling along an allosteric communication channel drives the binding of Jun-Fos heterodimeric transcription factor to DNA

Although allostery plays a central role in driving protein-DNA interactions, the physical basis of such cooperative behavior remains poorly understood. In the present study, using isothermal titration calorimetry in conjunction with site-directed mutagenesis, we provide evidence that an intricate network of energetically-coupled residues within the basic regions of the Jun-Fos heterodimeric transcription factor accounts for its allosteric binding to DNA. Remarkably, energetic coupling is prevalent in residues that are both close in space, as well as residues distant in space, implicating the role of both short- and long-range cooperative interactions in driving the assembly of this key protein-DNA interaction. Unexpectedly, many of the energetically-coupled residues involved in orchestrating such a cooperative network of interactions are poorly conserved across other members of the basic zipper family, emphasizing the importance of basic residues in dictating the specificity of basic zipper-DNA interactions. Collectively, our thermodynamic analysis maps an allosteric communication channel driving a key protein-DNA interaction central to cellular functions in health and disease.

Structural recognition of DNA by poly(ADP-ribose)polymerase-like zinc finger families

PARP-like zinc fingers (zf-PARPs) are protein domains apt to the recognition of multiple DNA secondary structures. They were initially described as the DNA-binding, nick-sensor domains of poly(ADP-ribose)polymerases (PARPs). It now appears that zf-PARPs are evolutionary conserved in the eukaryotic lineage and associated with various enzymes implicated in nucleic acid transactions. In the present study, we discuss the functional and structural data of zf-PARPSs in the light of a comparative analysis of the protein family. Sequence and structural analyses allow the definition of the conserved features of the zf-PARP domain and the identification of five distinct phylogenetic groups. Differences among the groups accumulate on the putative DNA binding surface of the PARP zinc-finger fold. These observations suggest that different zf-PARP types have distinctive recognition properties for DNA secondary structures. A comparison of various functional studies confirms that the different finger types can accomplish a selective recognition of DNA structures.

Increases in preproenkephalin mRNA levels in the Syrian hamster: The influence of glucocorticoids is dependent on age and tissue

In adult hamsters, basal proenkephalin (Penk) gene expression in adrenals is independent of glucocorticoids and glucocorticoid receptor blockade, by RU 486, increases striatal preproenkephalin (PPenk) mRNA levels. However, glucocorticoids maintain both basal and induced Penk gene expression in rat adrenal (medulla) and striatum. This suggests species and tissue-specific differences in Penk gene regulation. Since studies show temporal coordination in Penk gene expression in developing hamster adrenal and striatum, we tested the hypothesis that increasing PPenk mRNA levels are dependent, while basal levels are independent of glucocorticoids in developing hamsters. To facilitate this study, we examined the influence of glucocorticoids on the temporal increases in developing hamster PPenk mRNA observed in adrenals between postnatal days 0 and 4 and in striatum between postnatal days 12 and 48. PPenk mRNA levels were determined in hamster pups after treatment with increasing doses of metyrapone (an 11beta hydroxylase inhibitor) or with the glucocorticoid receptor antagonist RU 486 +/- metyrapone between postnatal days 2 and 4. Levels were also determined 36 days after hypophysectomy at age 16-17 days. Although plasma glucocorticoid levels and/or the influence from glucocorticoids were reduced, only developmental increases in PPenk mRNA are influenced by glucocorticoids in hamster adrenals, while basal adrenal mRNA levels are unchanged. However, pituitary influence on striatal PPenk mRNA levels appears complex and may involve steroid and/or non-steroid factors. These results suggest that glucocorticoids regulate hamster Penk gene expression via a mechanism that varies with age and tissue and functions during the induction of the Penk gene and not to maintain basal gene expression. Possible mechanisms and species variation are discussed.

A simple motif for protein recognition in DNA secondary structures

DNA in a single-stranded form (ssDNA) exists transiently within the cell and comprises the telomeres of linear chromosomes and the genomes of some DNA viruses. As with RNA, in the single-stranded state, some DNA sequences are able to fold into complex secondary and tertiary structures that may be recognized by proteins and participate in gene regulation. To better understand how such DNA elements might fold and interact with proteins, and to compare recognition features to those of a structured RNA, we used in vitro selection to identify ssDNAs that bind an RNA-binding peptide from the HIV Rev protein with high affinity and specificity. The large majority of selected binders contain a non-Watson-Crick G.T base-pair and an adjacent C:G base-pair and both are essential for binding. This GT motif can be presented in different DNA contexts, including a nearly perfect duplex and a branched three-helix structure, and appears to be recognized in large part by arginine residues separated by one turn of an alpha-helix. Interestingly, a very similar GT motif is necessary also for protein binding and function of a well-characterized model ssDNA regulatory element from the proenkephalin promoter.

Nuclease-deficient FEN-1 blocks Rad51/BRCA1-mediated repair and causes trinucleotide repeat instability

Previous studies have shown that expansion-prone repeats form structures that inhibit human flap endonuclease (FEN-1). We report here that faulty processing by FEN-1 initiates repeat instability in mammalian cells. Disease-length CAG tracts in Huntington's disease mice heterozygous for FEN-1 display a tendency toward expansions over contractions during intergenerational inheritance compared to those in homozygous wild-type mice. Further, with regard to human cells expressing a nuclease-defective FEN-1, we provide direct evidence that an unprocessed FEN-1 substrate is a precursor to instability. In cells with no endogenous defects in DNA repair, exogenous nuclease-defective FEN-1 causes repeat instability and aberrant DNA repair. Inefficient flap processing blocks the formation of Rad51/BRCA1 complexes but invokes repair by other pathways.

Calmodulin kinase II attenuation of gene transcription by preventing cAMP response element-binding protein (CREB) dimerization and binding of the CREB-binding protein

Calmodulin Kinase II (CamKII) inhibits the transcription of many CRE-dependent genes, but the mechanism of dominant transcriptional inhibition is unknown. Here we show that phosphorylation of serine 142 in CREB by CamKII leads to dissociation of the CREB dimer without impeding DNA binding capacity. CamKII-modified CREB binds to DNA efficiently as a monomer; however, monomeric CREB is unable to recruit the CREB-binding protein (CBP) even when phosphorylated at serine 133. Thus, CamKII confers a dominant inhibitory effect on transcription by preventing dimerization of CREB, and this mechanism may account for the attenuation of gene expression.

In vitro and in vivo ligation-mediated polymerase chain reaction analysis of a polypurine/polypyrimidine sequence upstream of the mouse metallothionein-I gene

The mouse metallothionein-I homopurine/homopyrimidine (MT-I R/Y) sequence is a 128-base pair element located approximately 1.2 kilobase pairs upstream of the MT-I gene. Previous in vitro studies of this sequence in purified plasmids indicated the formation of a non-B DNA structure stabilized by acidic pH and negative supercoiling. We now present a detailed in vitro and in vivo analysis of the MT-I R/Y sequence using chemical probes of DNA structure and ligation-mediated polymerase chain reaction. In vivo analysis suggests neither profound base unpairing nor protein binding within the MT-I R/Y sequence before or after metal induction of MT-I. We conclude for this element that the propensity to adopt an unusual DNA structure in vitro does not imply the occurrence of such a structure in vivo. We were able to show both in purified genomic DNA and in vivo that only isolated thymines and the 3' terminal thymine in strings of consecutive thymines are modified significantly by KMnO(4), indicating an altered thymine accessibility pattern within the R/Y sequence. This KMnO(4) reactivity pattern is more consistent and predictable within the R/Y sequence when compared with flanking sequences. We propose a simple steric interference model to explain the observed pattern of KMnO(4) modification of thymines.

Keratin 10 gene expression during differentiation of mouse epidermis requires transcription factors C/EBP and AP-2

The epidermis forms a vital barrier composed of stratified keratinocytes and their differentiated products. One of these products, keratin K10, is critical to epidermal integrity, because mutations in k10 lead to abnormal blistering. For the normal expression of k10, differentiation-associated transcription factors C/EBPalpha, C/EBPbeta, and AP-2 are well positioned to play an important role. Here, regulation of the k10 gene is examined in keratinocytes in the skin of normal mice and in transgenic mice carrying targeted deletions of c/ebpbeta and ap-2alpha. In cultured cells, C/EBPalpha and C/EBPbeta are each capable of activating the k10 promoter via three binding sites, identified by site-directed mutagenesis. In a given epidermal cell in vivo, however, the selection of C/EBPalpha versus C/EBPbeta for k10 regulation is determined via a third transcription factor, AP-2. This novel regulatory scheme involves: (1) unique gradients of expression for each transcription factor, i.e., C/EBPbeta and AP-2 most abundant in the lower epidermis, C/EBPalpha in the upper; (2) C/EBP-binding sites in the ap-2alpha gene promoter, through which C/EBPbeta stimulates ap-2alpha; and (3) AP-2 binding sites in the c/ebpalpha promoter, through which AP-2 represses c/ebpalpha. Promoter-analysis and gene-expression data presented herein support a regulatory model in which C/EBPbeta activates and maintains AP-2 expression in basal keratinocytes, whereas AP-2 represses C/EBPalpha in those cells. In response to differentiation signals, loss of AP-2 expression leads to derepression of the c/ebpalpha promoter and activation of k10 as cells migrate upward.

Sequence and DNA structural determinants of N4 virion RNA polymerase-promoter recognition

Coliphage N4-coded, virion-encapsidated RNA polymerase (vRNAP) is able to bind to and transcribe promoter-containing double-stranded DNAs when the template is supercoiled and Escherichia coli single-stranded DNA-binding protein (Eco SSB) is present. We report that vRNAP-promoter recognition and activity on these templates require specific sequences and a hairpin structure on the template strand. Hairpin extrusion, induced by Mg(II) and physiological superhelical density, is essential to provide the correct DNA structure for polymerase recognition, as mutant promoters that do not form hairpins show reduced in vitro activity. Therefore, a supercoil-induced DNA structural transition regulates N4 vRNAP transcription. Eco SSB activates transcription at physiological superhelical densities by stabilizing the template-strand hairpin. Specific sequences at the promoters are conserved to provide proper contacts for vRNAP, to support hairpin extrusion, or both. We propose a model for in vivo utilization of the vRNAP promoters, and discuss the roles of DNA supercoiling and Eco SSB in promoter activation.

cAMP response element-binding protein monomers cooperatively assemble to form dimers on DNA

We have analyzed the properties of cAMP response element-binding protein (CREB) in solution with emphasis on dimerization and effects of phosphorylation. Using a purified CREB fusion protein, a novel dye-label technique, and sedimentation equilibrium analysis, we directly and conclusively demonstrate that, unlike Jun and Fos, CREB dimerization is DNA-dependent. CREB exists primarily as a monomer in solution and cooperatively assembles on DNA to form dimers. Sedimentation equilibrium analysis also indicates that dimerization is unaffected by cAMP-dependent protein kinase-phosphorylation or by the symmetry of the cAMP-responsive element binding site. Filter binding assays reveal that CREB binding is unaffected by phosphorylation regardless of the symmetry of the cAMP-responsive element binding site. Our results suggest that structurally similar members of the same bZIP superfamily may differ significantly in their regulation at the level of dimerization.

Influence of hairpins on template reannealing at trinucleotide repeat duplexes: a model for slipped DNA

Hairpin stabilization of polymerase slippage has been proposed as part of the mechanism for large-scale expansion of CG-rich (CNG, where N = A, T, G, or C) trinucleotide repeats. However, hairpin formation does not entirely account for why long repeats but not short repeats or palindromes expand. Using ultraviolet spectroscopic methods, we examine the thermodynamic and kinetic properties of repeating trinucleotides to evaluate their behavior at a slippage site. We find that CNG trinucleotide repeats associated with expansion form stable hairpins whether they are short (with as few as 10 repeats) or long. However, long repeating stretches exist as single strands up to 2 orders of magnitude longer than sequences with either short repeats or random DNA. Thus, long hairpins have long lifetimes even in the presence of their complementary strands and inhibit duplex reannealing at a slippage site. The kinetic properties explain why expansion occurs with high frequency at long repeats but not at short repeats or palindromes.

Cruciform-extruding regulatory element controls cell-specific activity of the tyrosine hydroxylase gene promoter

Tyrosine hydroxylase (TH) is expressed specifically in catecholaminergic cells. We have identified a novel regulatory sequence in the upstream region of the bovine TH gene promoter formed by a dyad symmetry element (DSE1;-352/-307 bp). DSE1 supports TH promoter activity in TH-expressing bovine adrenal medulla chromaffin (BAMC) cells and inhibits promoter activity in non-expressing TE671 cells. DNase I footprinting of relaxed TH promoter DNA showed weak binding of nuclear BAMC cell proteins to a short sequence in the right DSE1 arm. In BAMC cells, deletion of the right arm markedly reduced the expression of luciferase from the TH promoter. However, deletion of the left DSE1 arm or its reversed orientation (RevL) also inactivated the TH promoter. In supercoiled TH promoter, DSE1 assumes a cruciform-like conformation i.e., it binds cruciform-specific 2D3 antibody, and S1 nuclease-cleavage and OsO4-modification assays have identified an imperfect cruciform extruded by the DSE1. DNase I footprinting of supercoiled plasmid showed that cruciformed DSE1 is targeted by nuclear proteins more efficiently than the linear duplex isomer and that the protected site encompasses the left arm and center of DSE1. Our results suggest that the disruption of intrastrand base-pairing preventing cruciform formation and protein binding to DSE1 is responsible for its inactivation in DSE1 mutants. DSE1 cruciform may act as a target site for activator (BAMC cells) and repressor (TE671) proteins. Its extrusion emerges as a novel mechanism that controls cell-specific promoter activity.

Initiation of DNA replication at palindromic telomeres is mediated by a duplex-to-hairpin transition induced by the minute virus of mice nonstructural protein NS1

The linear single-stranded DNA genome of the minute virus of mice (MVM) is replicated via a double-stranded replicative form (RF) intermediate. Amplification of this RF is initiated by the folding-back of palindromic sequences serving as primers for strand-displacement synthesis and formation of dimeric RF DNA. Using an in vitro replication assay and a cloned MVM DNA template, we observed hairpin-primed DNA replication at both MVM DNA termini, with a bias toward right-end initiation. Initiation of DNA replication is favored by nuclear components of A9 cell extract and highly stimulated by the MVM nonstructural protein NS1. Hairpin-primed DNA replication is also observed in the presence of NS1 and the Klenow fragment of the Escherichia coli DNA polymerase I. Addition of ATPgammaS (adenosine 5'-O-(thiotriphosphate)) blocks the initiation of DNA replication but not the extension of pre-existing hairpin primers formed in the presence of NS1 only. The NS1-mediated unwinding of the right-end palindrome may account for the recently reported capacity of NS1 for driving dimer RF synthesis in vitro.

Switching of DNA secondary structure in proenkephalin transcriptional regulation

Proper transcriptional regulation of the proenkephalin gene requires a switch between distinct factor binding sites that cannot exist at the same time. Each of the sites is formed from a nearly palindromic region that contains two functionally defined cAMP response elements. The region can switch between cruciform and linear duplex. Formation of the cruciform creates an alternative binding site for mediators of second messenger-directed transcription and abolishes the site present in the native duplex form. Use of the cruciform site has been shown to correlate with activated transcription. Analysis of DNA structure, protein binding, and gene expression from plasmids with mutant enhancers shows, however, that both sites are required for regulation of transcription. The two distinct structures form within the same enhancer. Shifting the balance between the two alters transcriptional response.