S1-END-seq reveals DNA secondary structures in human cells.

DNA becomes single stranded (ssDNA) during replication, transcription, and repair. Transiently formed ssDNA segments can adopt alternative conformations, including cruciforms, triplexes, and quadruplexes. To determine whether there are stable regions of ssDNA in the human genome, we utilized S1-END-seq to convert ssDNA regions to DNA double-strand breaks, which were then processed for high-throughput sequencing. This approach revealed two predominant non-B DNA structures: cruciform DNA formed by expanded (TA)n repeats that accumulate in microsatellite unstable human cancer cell lines and DNA triplexes (H-DNA) formed by homopurine/homopyrimidine mirror repeats common across a variety of cell lines. We show that H-DNA is enriched during replication, that its genomic location is highly conserved, and that H-DNA formed by (GAA)n repeats can be disrupted by treatment with a (GAA)n-binding polyamide. Finally, we show that triplex-forming repeats are hotspots for mutagenesis. Our results identify dynamic DNA secondary structures in vivo that contribute to elevated genome instability.

A partner evokes latent differences between Hox proteins.

Hox transcription factors bind highly related DNA sequences in vitro, yet they regulate different genes and play distinct roles in anterior-posterior patterning in animals. Slattery et al. report that a common cofactor, Exd, accentuates latent sequence specificities of all eight Hox proteins and directs binding to relevant sites across the genome.

Using synthetic genome readers/regulators to interrogate chromatin processes: A brief review.

Aberrant gene expression underlies numerous human ailments. Hence, developing small molecules to target and remedy dysfunctional gene regulation has been a long-standing goal at the interface of chemistry and medicine. A major challenge for designing small molecule therapeutics aimed at targeting desired genomic loci is the minimization of widescale disruption of genomic functions. To address this challenge, we rationally design polyamide-based multi-functional molecules, i.e., Synthetic Genome Readers/Regulators (SynGRs), which, by design, target distinct sequences in the genome. Herein, we briefly review how SynGRs access chromatin-bound and chromatin-free genomic sites, then highlight the methods for the study of chromatin processes using SynGRs on positioned nucleosomes in vitro or disease-causing repressive genomic loci in vivo.

Engineered Covalent Inactivation of TFIIH-Kinase Reveals an Elongation Checkpoint and Results in Widespread mRNA Stabilization.

During transcription initiation, the TFIIH-kinase Kin28/Cdk7 marks RNA polymerase II (Pol II) by phosphorylating the C-terminal domain (CTD) of its largest subunit. Here we describe a structure-guided chemical approach to covalently and specifically inactivate Kin28 kinase activity in vivo. This method of irreversible inactivation recapitulates both the lethal phenotype and the key molecular signatures that result from genetically disrupting Kin28 function in vivo. Inactivating Kin28 impacts promoter release to differing degrees and reveals a "checkpoint" during the transition to productive elongation. While promoter-proximal pausing is not observed in budding yeast, inhibition of Kin28 attenuates elongation-licensing signals, resulting in Pol II accumulation at the +2 nucleosome and reduced transition to productive elongation. Furthermore, upon inhibition, global stabilization of mRNA masks different degrees of reduction in nascent transcription. This study resolves long-standing controversies on the role of Kin28 in transcription and provides a rational approach to irreversibly inhibit other kinases in vivo.

Domain-Selective BET Ligands Yield Next-Generation Synthetic Genome Readers/Regulators with Nonidentical Cellular Functions.

SynTEF1, a prototype synthetic genome reader/regulator (SynGR), was designed to target GAA triplet repeats and restore the expression of frataxin (FXN) in Friedreich's ataxia patients. It achieves this complex task by recruiting BRD4, via a pan-BET ligand (JQ1), to the GAA repeats by using a sequence-selective DNA-binding polyamide. When bound to specific genomic loci in this way, JQ1 functions as a chemical prosthetic for acetyl-lysine residues that are natural targets of the two tandem bromodomains (BD1 and BD2) in bromo- and extra-terminal domain (BET) proteins. As next-generation BET ligands were disclosed, we tested a select set with improved physicochemical, pharmacological, and bromodomain-selective properties as substitutes for JQ1 in the SynGR design. Here, we report two unexpected findings: (1) SynGRs bearing pan-BET or BD2-selective ligands license transcription at the FXN locus, whereas those bearing BD1-selective ligands do not, and (2) rather than being neutral or inhibitory, an untethered BD1-selective ligand (GSK778) substantively enhances the activity of all active SynGRs. The failure of BD1-selective SynGRs to recruit BRD4/BET proteins suggests that rather than functioning as "epigenetic/chromatin mimics," active SynGRs mimic the functions of natural transcription factors in engaging BET proteins through BD2 binding. Moreover, the enhanced activity of SynGRs upon cotreatment with the BD1-selective ligand suggests that natural transcription factors compete for a limited pool of nonchromatin-bound BET proteins, and blocking BD1 directs pan-BET ligands to more effectively engage BD2. Taken together, SynGRs as chemical probes provide unique insights into the molecular recognition principles utilized by natural factors to precisely regulate gene expression, and they guide the design of more sophisticated synthetic gene regulators with greater therapeutic potential.

Noncanonical CTD kinases regulate RNA polymerase II in a gene-class-specific manner.

Phosphorylation of the carboxyl-terminal domain (CTD) of the largest subunit of RNA polymerase II (Pol II) governs stage-specific interactions with different cellular machines. The CTD consists of Y1S2P3T4S5P6S7 heptad repeats and sequential phosphorylations of Ser7, Ser5 and Ser2 occur universally at Pol II-transcribed genes. Phosphorylation of Thr4, however, appears to selectively modulate transcription of specific classes of genes. Here, we identify ten new Thr4 kinases from different kinase structural groups. Irreversible chemical inhibition of the most active Thr4 kinase, Hrr25, reveals a novel role for this kinase in transcription termination of specific class of noncoding snoRNA genes. Genome-wide profiles of Hrr25 reveal a selective enrichment at 3' regions of noncoding genes that display termination defects. Importantly, phospho-Thr4 marks placed by Hrr25 are recognized by Rtt103, a key component of the termination machinery. Our results suggest that these uncommon CTD kinases place phospho-Thr4 marks to regulate expression of targeted genes.

De novo design of programmable inducible promoters.

Ligand-responsive allosteric transcription factors (aTF) play a vital role in genetic circuits and high-throughput screening because they transduce biochemical signals into gene expression changes. Programmable control of gene expression from aTF-regulated promoter is important because different downstream effector genes function optimally at different expression levels. However, tuning gene expression of native promoters is difficult due to complex layers of homeostatic regulation encoded within them. We engineered synthetic promoters de novo by embedding operator sites with varying affinities and radically reshaped binding preferences within a minimal, constitutive Escherichia coli promoter. Multiplexed cell-based screening of promoters for three TetR-like aTFs generated with this approach gave rich diversity of gene expression levels, dynamic ranges and ligand sensitivities and were 50- to 100-fold more active over their respective native promoters. Machine learning on our dataset revealed that relative position of the core motif and bases flanking the core motif play an important role in modulating induction response. Our generalized approach yields customizable and programmable aTF-regulated promoters for engineering cellular pathways and enables the discovery of new small molecule biosensors.

Specificity landscapes unmask submaximal binding site preferences of transcription factors.

We have developed Differential Specificity and Energy Landscape (DiSEL) analysis to comprehensively compare DNA-protein interactomes (DPIs) obtained by high-throughput experimental platforms and cutting edge computational methods. While high-affinity DNA binding sites are identified by most methods, DiSEL uncovered nuanced sequence preferences displayed by homologous transcription factors. Pairwise analysis of 726 DPIs uncovered homolog-specific differences at moderate- to low-affinity binding sites (submaximal sites). DiSEL analysis of variants of 41 transcription factors revealed that many disease-causing mutations result in allele-specific changes in binding site preferences. We focused on a set of highly homologous factors that have different biological roles but "read" DNA using identical amino acid side chains. Rather than direct readout, our results indicate that DNA noncontacting side chains allosterically contribute to sculpt distinct sequence preferences among closely related members of transcription factor families.

A chemoprobe tracks its target.

Small-molecule inhibitors of histone-modifying enzymes have significant clinical utility for managing diseases such as cancer. These inhibitors are usually identified and monitored through their effects on the gain or loss of specific histone marks. In cells, multiple related enzymes can place or remove a specific mark; therefore, relying on an indirect measure of inhibitor engagement can be misleading. Mascaró et al. describe a luminescence-based ELISA approach that directly monitors binding of inhibitors to the histone lysine demethylase KDM1A.

Flexibility and structure of flanking DNA impact transcription factor affinity for its core motif.

Spatial and temporal expression of genes is essential for maintaining phenotype integrity. Transcription factors (TFs) modulate expression patterns by binding to specific DNA sequences in the genome. Along with the core binding motif, the flanking sequence context can play a role in DNA-TF recognition. Here, we employ high-throughput in vitro and in silico analyses to understand the influence of sequences flanking the cognate sites in binding of three most prevalent eukaryotic TF families (zinc finger, homeodomain and bZIP). In vitro binding preferences of each TF toward the entire DNA sequence space were correlated with a wide range of DNA structural parameters, including DNA flexibility. Results demonstrate that conformational plasticity of flanking regions modulates binding affinity of certain TF families. DNA duplex stability and minor groove width also play an important role in DNA-TF recognition but differ in how exactly they influence the binding in each specific case. Our analyses further reveal that the structural features of preferred flanking sequences are not universal, as similar DNA-binding folds can employ distinct DNA recognition modes.

Different phosphoisoforms of RNA polymerase II engage the Rtt103 termination factor in a structurally analogous manner.

The carboxyl-terminal domain (CTD) of the largest subunit of RNA polymerase II (Pol II) orchestrates dynamic recruitment of specific cellular machines during different stages of transcription. Signature phosphorylation patterns of Y1S2P3T4S5P6S7 heptapeptide repeats of the CTD engage specific "readers." Whereas phospho-Ser5 and phospho-Ser2 marks are ubiquitous, phospho-Thr4 is reported to only impact specific genes. Here, we identify a role for phospho-Thr4 in transcription termination at noncoding small nucleolar RNA (snoRNA) genes. Quantitative proteomics reveals an interactome of known readers as well as protein complexes that were not known to rely on Thr4 for association with Pol II. The data indicate a key role for Thr4 in engaging the machinery used for transcription elongation and termination. We focus on Rtt103, a protein that binds phospho-Ser2 and phospho-Thr4 marks and facilitates transcription termination at protein-coding genes. To elucidate how Rtt103 engages two distinct CTD modifications that are differentially enriched at noncoding genes, we relied on NMR analysis of Rtt103 in complex with phospho-Thr4- or phospho-Ser2-bearing CTD peptides. The structural data reveal that Rtt103 interacts with phospho-Thr4 in a manner analogous to its interaction with phospho-Ser2-modified CTD. The same set of hydrogen bonds involving either the oxygen on phospho-Thr4 and the hydroxyl on Ser2, or the phosphate on Ser2 and the Thr4 hydroxyl, can be formed by rotation of an arginine side chain, leaving the intermolecular interface otherwise unperturbed. This economy of design enables Rtt103 to engage Pol II at distinct sets of genes with differentially enriched CTD marks.

Reprogramming cell fate with a genome-scale library of artificial transcription factors.

Artificial transcription factors (ATFs) are precision-tailored molecules designed to bind DNA and regulate transcription in a preprogrammed manner. Libraries of ATFs enable the high-throughput screening of gene networks that trigger cell fate decisions or phenotypic changes. We developed a genome-scale library of ATFs that display an engineered interaction domain (ID) to enable cooperative assembly and synergistic gene expression at targeted sites. We used this ATF library to screen for key regulators of the pluripotency network and discovered three combinations of ATFs capable of inducing pluripotency without exogenous expression of Oct4 (POU domain, class 5, TF 1). Cognate site identification, global transcriptional profiling, and identification of ATF binding sites reveal that the ATFs do not directly target Oct4; instead, they target distinct nodes that converge to stimulate the endogenous pluripotency network. This forward genetic approach enables cell type conversions without a priori knowledge of potential key regulators and reveals unanticipated gene network dynamics that drive cell fate choices.

Synthetic genome readers target clustered binding sites across diverse chromatin states.

Targeting the genome with sequence-specific DNA-binding molecules is a major goal at the interface of chemistry, biology, and precision medicine. Polyamides, composed of N-methylpyrrole and N-methylimidazole monomers, are a class of synthetic molecules that can be rationally designed to "read" specific DNA sequences. However, the impact of different chromatin states on polyamide binding in live cells remains an unresolved question that impedes their deployment in vivo. Here, we use cross-linking of small molecules to isolate chromatin coupled to sequencing to map the binding of two bioactive and structurally distinct polyamides to genomes directly within live H1 human embryonic stem cells. This genome-wide view from live cells reveals that polyamide-based synthetic genome readers bind cognate sites that span a range of binding affinities. Polyamides can access cognate sites within repressive heterochromatin. The occupancy patterns suggest that polyamides could be harnessed to target loci within regions of the genome that are inaccessible to other DNA-targeting molecules.

The conserved bromo-adjacent homology domain of yeast Orc1 functions in the selection of DNA replication origins within chromatin.

The origin recognition complex (ORC) binds to the specific positions on chromosomes that serve as DNA replication origins. Although ORC is conserved from yeast to humans, the DNA sequence elements that specify ORC binding are not. In particular, metazoan ORC shows no obvious DNA sequence specificity, whereas yeast ORC binds to a specific DNA sequence within all yeast origins. Thus, whereas chromatin must play an important role in metazoan ORC's ability to recognize origins, it is unclear whether chromatin plays a role in yeast ORC's recognition of origins. This study focused on the role of the conserved N-terminal bromo-adjacent homology domain of yeast Orc1 (Orc1BAH). Recent studies indicate that BAH domains are chromatin-binding modules. We show that the Orc1BAH domain was necessary for ORC's stable association with yeast chromosomes, and was physiologically relevant to DNA replication in vivo. This replication role was separable from the Orc1BAH domain's previously defined role in transcriptional silencing. Genome-wide analyses of ORC binding in ORC1 and orc1bahDelta cells revealed that the Orc1BAH domain contributed to ORC's association with most yeast origins, including a class of origins highly dependent on the Orc1BAH domain for ORC association (orc1bahDelta-sensitive origins). Orc1bahDelta-sensitive origins required the Orc1BAH domain for normal activity on chromosomes and plasmids, and were associated with a distinct local nucleosome structure. These data provide molecular insights into how the Orc1BAH domain contributes to ORC's selection of replication origins, as well as new tools for examining conserved mechanisms governing ORC's selection of origins within eukaryotic chromosomes.

Regulator trafficking on bacterial transcription units in vivo.

The trafficking patterns of the bacterial regulators of transcript elongation sigma(70), rho, NusA, and NusG on genes in vivo and the explanation for promoter-proximal peaks of RNA polymerase (RNAP) are unknown. Genome-wide, E. coli ChIP-chip revealed distinct association patterns of regulators as RNAP transcribes away from promoters (rho first, then NusA, then NusG). However, the interactions of elongating complexes with these regulators did not differ significantly among most transcription units. A modest variation of NusG signal among genes reflected increased NusG interaction as transcription progresses, rather than functional specialization of elongating complexes. Promoter-proximal RNAP peaks were offset from sigma(70) peaks in the direction of transcription and co-occurred with NusA and rho peaks, suggesting that the RNAP peaks reflected elongating, rather than initiating, complexes. However, inhibition of rho did not increase RNAP levels within genes downstream from the RNAP peaks, suggesting the peaks are caused by a mechanism other than rho-dependent attenuation.

TFIIH kinase places bivalent marks on the carboxy-terminal domain of RNA polymerase II.

Posttranslational modifications of the carboxy-terminal domain (CTD) of the largest subunit of RNA polymerase II (Pol II) specify a molecular recognition code that is deciphered by proteins involved in RNA biogenesis. The CTD is comprised of a repeating heptapeptide (Y(1)S(2)P(3)T(4)S(5)P(6)S(7)). Recently, phosphorylation of serine 7 was shown to be important for cotranscriptional processing of two snRNAs in mammalian cells. Here we report that Kin28/Cdk7, a subunit of the evolutionarily conserved TFIIH complex, is a Ser7 kinase. The ability of Kin28/Cdk7 to phosphorylate Ser7 is particularly surprising because this kinase functions at promoters of protein-coding genes, rather than being restricted to promoter-distal regions of snRNA genes. Kin28/Cdk7 is also known to phosphorylate Ser5 residues of the CTD at gene promoters. Taken together, our results implicate the TFIIH kinase in placing bivalent Ser5 and Ser7 marks early in gene transcription. These bivalent CTD marks, in concert with cues within nascent transcripts, specify the cotranscriptional engagement of the relevant RNA processing machinery.

A library of yeast transcription factor motifs reveals a widespread function for Rsc3 in targeting nucleosome exclusion at promoters.

The sequence specificity of DNA-binding proteins is the primary mechanism by which the cell recognizes genomic features. Here, we describe systematic determination of yeast transcription factor DNA-binding specificities. We obtained binding specificities for 112 DNA-binding proteins representing 19 distinct structural classes. One-third of the binding specificities have not been previously reported. Several binding sequences have striking genomic distributions relative to transcription start sites, supporting their biological relevance and suggesting a role in promoter architecture. Among these are Rsc3 binding sequences, containing the core CGCG, which are found preferentially approximately 100 bp upstream of transcription start sites. Mutation of RSC3 results in a dramatic increase in nucleosome occupancy in hundreds of proximal promoters containing a Rsc3 binding element, but has little impact on promoters lacking Rsc3 binding sequences, indicating that Rsc3 plays a broad role in targeting nucleosome exclusion at yeast promoters.

The bacterial response regulator ArcA uses a diverse binding site architecture to regulate carbon oxidation globally.

Despite the importance of maintaining redox homeostasis for cellular viability, how cells control redox balance globally is poorly understood. Here we provide new mechanistic insight into how the balance between reduced and oxidized electron carriers is regulated at the level of gene expression by mapping the regulon of the response regulator ArcA from Escherichia coli, which responds to the quinone/quinol redox couple via its membrane-bound sensor kinase, ArcB. Our genome-wide analysis reveals that ArcA reprograms metabolism under anaerobic conditions such that carbon oxidation pathways that recycle redox carriers via respiration are transcriptionally repressed by ArcA. We propose that this strategy favors use of catabolic pathways that recycle redox carriers via fermentation akin to lactate production in mammalian cells. Unexpectedly, bioinformatic analysis of the sequences bound by ArcA in ChIP-seq revealed that most ArcA binding sites contain additional direct repeat elements beyond the two required for binding an ArcA dimer. DNase I footprinting assays suggest that non-canonical arrangements of cis-regulatory modules dictate both the length and concentration-sensitive occupancy of DNA sites. We propose that this plasticity in ArcA binding site architecture provides both an efficient means of encoding binding sites for ArcA, σ(70)-RNAP and perhaps other transcription factors within the same narrow sequence space and an effective mechanism for global control of carbon metabolism to maintain redox homeostasis.