Defining the sequence-recognition profile of DNA-binding molecules

Dominant Malignant Clones Leverage Lineage Restricted Epigenomic Programs to Drive Ependymoma Development

ZFTA-RELA is the most recurrent genetic alteration seen in pediatric supratentorial ependymoma (EPN) and is sufficient to initiate tumors in mice. Despite ZFTA-RELA's potent oncogenic potential, ZFTA-RELA gene fusions are observed exclusively in childhood EPN, with tumors located distinctly in the supratentorial region of the central nervous system (CNS). We hypothesized that specific chromatin modules accessible during brain development would render distinct cell lineage programs at direct risk of transformation by ZFTA-RELA. To this end, we performed combined single cell ATAC and RNA-seq analysis (scMultiome) of the developing mouse forebrain as compared to ZR-driven mouse and human EPN. We demonstrate that specific developmental lineage programs present in radial glial cells and regulated by Plagl family transcription factors are at risk of neoplastic transformation. Binding of this chromatin network by ZFTA-RELA or other PLAGL family motif targeting fusion proteins leads to persistent chromatin accessibility at oncogenic loci and oncogene expression. Cross-species analysis of mouse and human EPN reveals significant cell type heterogeneity mirroring incomplete neurogenic and gliogenic differentiation, with a small percentage of cycling intermediate progenitor-like cells that establish a putative tumor cell hierarchy. In vivo lineage tracing studies reveal single neoplastic clones that aggressively dominate tumor growth and establish the entire EPN cellular hierarchy. These findings unravel developmental epigenomic states critical for fusion oncoprotein driven transformation and elucidate how these states continue to shape tumor progression.

HIGHLIGHTS

1. Specific chromatin modules accessible during brain development render distinct cell lineage programs at risk of transformation by pediatric fusion oncoproteins.2. Cross-species single cell ATAC and RNA (scMultiome) of mouse and human ependymoma (EPN) reveals diverse patterns of lineage differentiation programs that restrain oncogenic transformation.3. Early intermediate progenitor-like EPN cells establish a tumor cell hierarchy that mirrors neural differentiation programs.4. ZFTA-RELA transformation is compatible with distinct developmental epigenetic states requiring precise 'goldilocks' levels of fusion oncoprotein expression.5. Dominant tumor clones establish the entire EPN cellular hierarchy that reflects normal gliogenic and neurogenic differentiation programs.

Accelerated, high-quality photolithographic synthesis of RNA microarrays in situ

Nucleic acid photolithography is the only microarray fabrication process that has demonstrated chemical versatility accommodating any type of nucleic acid. The current approach to RNA microarray synthesis requires long coupling and photolysis times and suffers from unavoidable degradation postsynthesis. In this study, we developed a series of RNA phosphoramidites with improved chemical and photochemical protection of the 2'- and 5'-OH functions. In so doing, we reduced the coupling time by more than half and the photolysis time by a factor of 4. Sequence libraries that would otherwise take over 6 hours to synthesize can now be prepared in half the time. Degradation is substantially lowered, and concomitantly, hybridization signals can reach over seven times those of the previous state of the art. Under those conditions, high-density RNA microarrays and RNA libraries can now be synthesized at greatly accelerated rates. We also synthesized fluorogenic RNA Mango aptamers on microarrays and investigated the effect of sequence mutations on their fluorogenic properties.

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.

Nucleic acid-based small molecules as targeted transcription therapeutics for immunoregulation

Transcription therapy is an emerging approach that centers on identifying the factors associated with the malfunctioning gene transcription machinery that causes diseases and controlling them with designer agents. Until now, the primary research focus in therapeutic gene modulation has been on small-molecule drugs that target epigenetic enzymes and critical signaling pathways. However, nucleic acid-based small molecules have gained popularity in recent years for their amenability to be pre-designed and realize operative control over the dynamic transcription machinery that governs how the immune system responds to diseases. Pyrrole-imidazole polyamides (PIPs) are well-established DNA-based small-molecule gene regulators that overcome the limitations of their conventional counterparts owing to their sequence-targeted specificity, versatile regulatory efficiency, and biocompatibility. Here, we emphasize the rational design of PIPs, their functional mechanisms, and their potential as targeted transcription therapeutics for disease treatment by regulating the immune response. Furthermore, we also discuss the challenges and foresight of this approach in personalized immunotherapy in precision medicine.

Minding the Gap: Exploring Neuroinflammatory and Microglial Sex Differences in Alzheimer's Disease

Research into Alzheimer's Disease (AD) describes a link between AD and the resident immune cells of the brain, the microglia. Further, this suspected link is thought to have underlying sex effects, although the mechanisms of these effects are only just beginning to be understood. Many of these insights are the result of policies put in place by funding agencies such as the National Institutes of Health (NIH) to consider sex as a biological variable (SABV) and the move towards precision medicine due to continued lackluster therapeutic options. The purpose of this review is to provide an updated assessment of the current research that summarizes sex differences and the research pertaining to microglia and their varied responses in AD.

A rugged yet easily navigable fitness landscape

Fitness landscape theory predicts that rugged landscapes with multiple peaks impair Darwinian evolution, but experimental evidence is limited. In this study, we used genome editing to map the fitness of >260,000 genotypes of the key metabolic enzyme dihydrofolate reductase in the presence of the antibiotic trimethoprim, which targets this enzyme. The resulting landscape is highly rugged and harbors 514 fitness peaks. However, its highest peaks are accessible to evolving populations via abundant fitness-increasing paths. Different peaks share large basins of attraction that render the outcome of adaptive evolution highly contingent on chance events. Our work shows that ruggedness need not be an obstacle to Darwinian evolution but can reduce its predictability. If true in general, the complexity of optimization problems on realistic landscapes may require reappraisal.

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.

Dipodal Silanes Greatly Stabilize Glass Surface Functionalization for DNA Microarray Synthesis and High-Throughput Biological Assays

Glass is by far the most common substrate for biomolecular arrays, including high-throughput sequencing flow cells and microarrays. The native glass hydroxyl surface is modified by using silane chemistry to provide appropriate functional groups and reactivities for either in situ synthesis or surface immobilization of biologically or chemically synthesized biomolecules. These arrays, typically of oligonucleotides or peptides, are then subjected to long incubation times in warm aqueous buffers prior to fluorescence readout. Under these conditions, the siloxy bonds to the glass are susceptible to hydrolysis, resulting in significant loss of biomolecules and concomitant loss of signal from the assay. Here, we demonstrate that functionalization of glass surfaces with dipodal silanes results in greatly improved stability compared to equivalent functionalization with standard monopodal silanes. Using photolithographic in situ synthesis of DNA, we show that dipodal silanes are compatible with phosphoramidite chemistry and that hybridization performed on the resulting arrays provides greatly improved signal and signal-to-noise ratios compared with surfaces functionalized with monopodal silanes.

Genome-wide association and multi-omics studies identify MGMT as a novel risk gene for Alzheimer's disease among women

INTRODUCTION

Variants in the tau gene (MAPT) region are associated with breast cancer in women and Alzheimer's disease (AD) among persons lacking apolipoprotein E ε4 (ε4-).

METHODS

To identify novel genes associated with tau-related pathology, we conducted two genome-wide association studies (GWAS) for AD, one among 10,340 ε4- women in the Alzheimer's Disease Genetics Consortium (ADGC) and another in 31 members (22 women) of a consanguineous Hutterite kindred.

RESULTS

We identified novel associations of AD with MGMT variants in the ADGC (rs12775171, odds ratio [OR] = 1.4, P = 4.9 × 10-8) and Hutterite (rs12256016 and rs2803456, OR = 2.0, P = 1.9 × 10-14) datasets. Multi-omics analyses showed that the most significant and largest number of associations among the single nucleotide polymorphisms (SNPs), DNA-methylated CpGs, MGMT expression, and AD-related neuropathological traits were observed among women. Furthermore, promoter capture Hi-C analyses revealed long-range interactions of the MGMT promoter with MGMT SNPs and CpG sites.

DISCUSSION

These findings suggest that epigenetically regulated MGMT expression is involved in AD pathogenesis, especially in women.

Physicochemical models of protein-DNA binding with standard and modified base pairs

DNA-binding proteins play important roles in various cellular processes, but the mechanisms by which proteins recognize genomic target sites remain incompletely understood. Functional groups at the edges of the base pairs (bp) exposed in the DNA grooves represent physicochemical signatures. As these signatures enable proteins to form specific contacts between protein residues and bp, their study can provide mechanistic insights into protein-DNA binding. Existing experimental methods, such as X-ray crystallography, can reveal such mechanisms based on physicochemical interactions between proteins and their DNA target sites. However, the low throughput of structural biology methods limits mechanistic insights for selection of many genomic sites. High-throughput binding assays enable prediction of potential target sites by determining relative binding affinities of a protein to massive numbers of DNA sequences. Many currently available computational methods are based on the sequence of standard Watson-Crick bp. They assume that the contribution of overall binding affinity is independent for each base pair, or alternatively include dinucleotides or short k-mers. These methods cannot directly expand to physicochemical contacts, and they are not suitable to apply to DNA modifications or non-Watson-Crick bp. These variations include DNA methylation, and synthetic or mismatched bp. The proposed method, DeepRec, can predict relative binding affinities as function of physicochemical signatures and the effect of DNA methylation or other chemical modifications on binding. Sequence-based modeling methods are in comparison a coarse-grain description and cannot achieve such insights. Our chemistry-based modeling framework provides a path towards understanding genome function at a mechanistic level.

Flexibility of flanking DNA is a key determinant of transcription factor affinity for the core motif

Selective gene regulation is mediated by recognition of specific DNA sequences by transcription factors (TFs). The extremely challenging task of searching out specific cognate DNA binding sites among several million putative sites within the eukaryotic genome is achieved by complex molecular recognition mechanisms. Elements of this recognition code include the core binding sequence, the flanking sequence context, and the shape and conformational flexibility of the composite binding site. To unravel the extent to which DNA flexibility modulates TF binding, in this study, we employed experimentally guided molecular dynamics simulations of ternary complex of closely related Hox heterodimers Exd-Ubx and Exd-Scr with DNA. Results demonstrate that flexibility signatures embedded in the flanking sequences impact TF binding at the cognate binding site. A DNA sequence has intrinsic shape and flexibility features. While shape features are localized, our analyses reveal that flexibility features of the flanking sequences percolate several basepairs and allosterically modulate TF binding at the core. We also show that lack of flexibility in the motif context can render the cognate site resistant to protein-induced shape changes and subsequently lower TF binding affinity. Overall, this study suggests that flexibility-guided DNA shape, and not merely the static shape, is a key unexplored component of the complex DNA-TF recognition code.

The Effect of Gut Microbe Dysbiosis on the Pathogenesis of Alzheimer's Disease (AD) and related conditions

It has been hypothesized that the shift in gut microbiota composition, known as gut microbe dysbiosis, may be correlated with the onset of Alzheimer's disease (AD), which is the most common cause of dementia characterized by a gradual deterioration in cognitive function associated with the development of amyloid-beta (Aβ) plaques. The gut microbiota dysbiosis induces the release of significant amounts of amyloids, lipopolysaccharides, and neurotoxins, which might play a role in modulating signaling pathways and immune activation, leading to the production of proinflammatory cytokines related to the pathogenesis of AD. The dysbiosis of gut microbe is associated with various diseases such as type 2 diabetes, obesity, hypertension, and some neuropsychiatric disorders like depression, anxiety, and stress. It is conceivable that these diseases trigger the onset of AD. Thus, modifying the gut microbiota composition with probiotic and prebiotic supplementation can reduce depression and anxiety symptoms, lower stress reactivity, and improve memory. This narrative review aimed to examine the possible role of gut microbe dysbiosis in AD's pathogenesis.

CASCADE: high-throughput characterization of regulatory complex binding altered by non-coding variants

Non-coding DNA variants (NCVs) impact gene expression by altering binding sites for regulatory complexes. New high-throughput methods are needed to characterize the impact of NCVs on regulatory complexes. We developed CASCADE (Customizable Approach to Survey Complex Assembly at DNA Elements), an array-based high-throughput method to profile cofactor (COF) recruitment. CASCADE identifies DNA-bound transcription factor-cofactor (TF-COF) complexes in nuclear extracts and quantifies the impact of NCVs on their binding. We demonstrate CASCADE sensitivity in characterizing condition-specific recruitment of COFs p300 and RBBP5 (MLL subunit) to the CXCL10 promoter in lipopolysaccharide (LPS)-stimulated human macrophages and quantify the impact of all possible NCVs. To demonstrate applicability to NCV screens, we profile TF-COF binding to ~1,700 single-nucleotide polymorphism quantitative trait loci (SNP-QTLs) in human macrophages and identify perturbed ETS domain-containing complexes. CASCADE will facilitate high-throughput testing of molecular mechanisms of NCVs for diverse biological applications.

Bray et al. develop CASCADE, a method to profile transcription factor (TF)-cofactor (COF) complexes binding to DNA. They demonstrate the approach by profiling complex binding across the CXCL10 cytokine promoter and to ~1,700 single-nucleotide polymorphisms (SNPs). They anticipate that CASCADE can be applied to diverse biological systems to examine regulatory complex binding to DNA.

Sequence-dependent quenching of fluorescein fluorescence on single-stranded and double-stranded DNA

Fluorescein is commonly used to label macromolecules, particularly proteins and nucleic acids, but its fluorescence is known to be strongly dependent on its direct chemical environment.

Inferring primase-DNA specific recognition using a data driven approach

DNA–protein interactions play essential roles in all living cells. Understanding of how features embedded in the DNA sequence affect specific interactions with proteins is both challenging and important, since it may contribute to finding the means to regulate metabolic pathways involving DNA–protein interactions. Using a massive experimental benchmark dataset of binding scores for DNA sequences and a machine learning workflow, we describe the binding to DNA of T7 primase, as a model system for specific DNA–protein interactions. Effective binding of T7 primase to its specific DNA recognition sequences triggers the formation of RNA primers that serve as Okazaki fragment start sites during DNA replication.

Targeting a noncanonical, hairpin-containing G-quadruplex structure from the MYCN gene

The MYCN gene encodes the transcription factor N-Myc, a driver of neuroblastoma (NB). Targeting G-quadruplexes (G4s) with small molecules is attractive strategy to control the expression of undruggable proteins such as N-Myc. However, selective binders to G4s are challenging to identify due to the structural similarity of many G4s. Here, we report the discovery of a small molecule ligand (4) that targets the noncanonical, hairpin containing G4 structure found in the MYCN gene using small molecule microarrays (SMMs). Unlike many G4 binders, the compound was found to bind to a pocket at the base of the hairpin region of the MYCN G4. This compound stabilizes the G4 and has affinity of 3.5 ± 1.6 μM. Moreover, an improved analog, MY-8, suppressed levels of both MYCN and MYCNOS (a lncRNA embedded within the MYCN gene) in NBEB neuroblastoma cells. This work indicates that the approach of targeting complex, hybrid G4 structures that exist throughout the human genome may be an applicable strategy to achieve selectivity for targeting disease-relevant genes including protein coding (MYCN) as well as non-coding (MYCNOS) gene products.

Single position substitution of hairpin pyrrole-imidazole polyamides imparts distinct DNA-binding profiles across the human genome

Pyrrole-imidazole (Py-Im) polyamides are synthetic molecules that can be rationally designed to target specific DNA sequences to both disrupt and recruit transcriptional machinery. While in vitro binding has been extensively studied, in vivo effects are often difficult to predict using current models of DNA binding. Determining the impact of genomic architecture and the local chromatin landscape on polyamide-DNA sequence specificity remains an unresolved question that impedes their effective deployment in vivo. In this report we identified polyamide-DNA interaction sites across the entire genome, by covalently crosslinking and capturing these events in the nuclei of human LNCaP cells. This technique confirms the ability of two eight ring hairpin-polyamides, with similar architectures but differing at a single ring position (Py to Im), to retain in vitro specificities and display distinct genome-wide binding profiles.

Chemical Approaches to the Development of Artificial Transcription Factors Based on Pyrrole-Imidazole Polyamides

To maintain the functions of living organisms, cells have developed complex gene regulatory networks. Transcription factors have a central role in spatiotemporal control of gene expression and this has motivated us to develop artificial transcription factors that mimic their function. We found that three functions could be mimicked by applying our chemical approaches: i) efficient delivery into organelles that contain target DNA, ii) specific DNA binding to the target genomic region, and iii) regulation of gene expression by interaction with other transcription coregulators. We chose pyrrole-imidazole polyamides (PIPs), sequence-selective DNA binding molecules, as DNA binding domains, and have achieved each of the required functions by introducing other functional moieties. The developed artificial transcription factors have potential as chemical tools that can be used to artificially modulate gene expression to enable cell fate control and to correct abnormal gene regulation for therapeutic purposes.

DNA mismatches reveal conformational penalties in protein-DNA recognition

Transcription factors recognize specific genomic sequences to regulate complex gene-expression programs. Although it is well-established that transcription factors bind to specific DNA sequences using a combination of base readout and shape recognition, some fundamental aspects of protein-DNA binding remain poorly understood1,2. Many DNA-binding proteins induce changes in the structure of the DNA outside the intrinsic B-DNA envelope. However, how the energetic cost that is associated with distorting the DNA contributes to recognition has proven difficult to study, because the distorted DNA exists in low abundance in the unbound ensemble3-9. Here we use a high-throughput assay that we term SaMBA (saturation mismatch-binding assay) to investigate the role of DNA conformational penalties in transcription factor-DNA recognition. In SaMBA, mismatched base pairs are introduced to pre-induce structural distortions in the DNA that are much larger than those induced by changes in the Watson-Crick sequence. Notably, approximately 10% of mismatches increased transcription factor binding, and for each of the 22 transcription factors that were examined, at least one mismatch was found that increased the binding affinity. Mismatches also converted non-specific sites into high-affinity sites, and high-affinity sites into 'super sites' that exhibit stronger affinity than any known canonical binding site. Determination of high-resolution X-ray structures, combined with nuclear magnetic resonance measurements and structural analyses, showed that many of the DNA mismatches that increase binding induce distortions that are similar to those induced by protein binding-thus prepaying some of the energetic cost incurred from deforming the DNA. Our work indicates that conformational penalties are a major determinant of protein-DNA recognition, and reveals mechanisms by which mismatches can recruit transcription factors and thus modulate replication and repair activities in the cell10,11.

Quantitative analysis of transcription factor binding and expression using calling cards reporter arrays

We report a tool, Calling Cards Reporter Arrays (CCRA), that measures transcription factor (TF) binding and the consequences on gene expression for hundreds of synthetic promoters in yeast. Using Cbf1p and MAX, we demonstrate that the CCRA method is able to detect small changes in binding free energy with a sensitivity comparable to in vitro methods, enabling the measurement of energy landscapes in vivo. We then demonstrate the quantitative analysis of cooperative interactions by measuring Cbf1p binding at synthetic promoters with multiple sites. We find that the cooperativity between Cbf1p dimers varies sinusoidally with a period of 10.65 bp and energetic cost of 1.37 K_BT for sites that are positioned ‘out of phase’. Finally, we characterize the binding and expression of a group of TFs, Tye7p, Gcr1p and Gcr2p, that act together as a ‘TF collective’, an important but poorly characterized model of TF cooperativity. We demonstrate that Tye7p often binds promoters without its recognition site because it is recruited by other collective members, whereas these other members require their recognition sites, suggesting a hierarchy where these factors recruit Tye7p but not vice versa. Our experiments establish CCRA as a useful tool for quantitative investigations into TF binding and function.

From Interaction to Function in DNA-Templated Supramolecular Self-Assemblies

DNA-templated self-assembly represents a rich and growing subset of supramolecular chemistry where functional self-assemblies are programmed in a versatile manner using nucleic acids as readily-available and readily-tunable templates. In this review, we summarize the different DNA recognition modes and the basic supramolecular interactions at play in this context. We discuss the recent results that report the DNA-templated self-assembly of small molecules into complex yet precise nanoarrays, going from 1D to 3D architectures. Finally, we show their emerging functions as photonic/electronic nanowires, sensors, gene delivery vectors, and supramolecular catalysts, and their growing applications in a wide range of area from materials to biological sciences.

Large-Scale Photolithographic Synthesis of Chimeric DNA/RNA Hairpin Microarrays To Explore Sequence Specificity Landscapes of RNase HII Cleavage

Ribonuclease HII (RNase HII) is an essential endoribonuclease that binds to double-stranded DNA with RNA nucleotide incorporations and cleaves 5' of the ribonucleotide at RNA-DNA junctions. Thought to be present in all domains of life, RNase HII protects genomic integrity by initiating excision repair pathways that protect the encoded information from rapid degradation. There is sparse evidence that the enzyme cleaves some substrates better than others, but a large-scale study is missing. Such large-scale studies can be carried out on microarrays, and we employ chemical photolithography to synthesize very large combinatorial libraries of fluorescently labeled DNA/RNA chimeric sequences that self-anneal to form hairpin structures that are substrates for Escherichia coli RNase HII. The relative activity is determined by the loss of fluorescence upon cleavage. Each substrate includes a double-stranded 5 bp variable region with one to five consecutive ribonucleotide substitutions. We also examined the effect of all possible single and double mismatches, for a total of >9500 unique structures. Differences in cleavage efficiency indicate some level of substrate preference, and we identified the 5'-dC/rC-rA-dX-3' motif in well-cleaved substrates. The results significantly extend known patterns of RNase HII sequence specificity and serve as a template using large-scale photolithographic synthesis to comprehensively map landscapes of substrate specificity of nucleic acid-processing enzymes.

Recent insights into the genotype-phenotype relationship from massively parallel genetic assays

With the molecular revolution in Biology, a mechanistic understanding of the genotype-phenotype relationship became possible. Recently, advances in DNA synthesis and sequencing have enabled the development of deep mutational scanning assays, capable of scoring comprehensive libraries of genotypes for fitness and a variety of phenotypes in massively parallel fashion. The resulting empirical genotype-fitness maps pave the way to predictive models, potentially accelerating our ability to anticipate the behaviour of pathogen and cancerous cell populations from sequencing data. Besides from cellular fitness, phenotypes of direct application in industry (e.g. enzyme activity) and medicine (e.g. antibody binding) can be quantified and even selected directly by these assays. This review discusses the technological basis of and recent developments in massively parallel genetics, along with the trends it is uncovering in the genotype-phenotype relationship (distribution of mutation effects, epistasis), their possible mechanistic bases and future directions for advancing towards the goal of predictive genetics.

BindSpace decodes transcription factor binding signals by large-scale sequence embedding

Decoding transcription factor (TF) binding signals in genomic DNA is a fundamental problem. Here we present a prediction model called BindSpace that learns to embed DNA sequences and TF class/family labels into the same space. By training on binding data for hundreds of TFs and embedding over 1M DNA sequences, BindSpace achieves state-of-the-art multiclass binding prediction performance, in vitro and in vivo, and can distinguish signals of closely related TFs.

Low-Affinity Binding Sites and the Transcription Factor Specificity Paradox in Eukaryotes

Eukaryotic transcription factors (TFs) from the same structural family tend to bind similar DNA sequences, despite the ability of these TFs to execute distinct functions in vivo. The cell partly resolves this specificity paradox through combinatorial strategies and the use of low-affinity binding sites, which are better able to distinguish between similar TFs. However, because these sites have low affinity, it is challenging to understand how TFs recognize them in vivo. Here, we summarize recent findings and technological advancements that allow for the quantification and mechanistic interpretation of TF recognition across a wide range of affinities. We propose a model that integrates insights from the fields of genetics and cell biology to provide further conceptual understanding of TF binding specificity. We argue that in eukaryotes, target specificity is driven by an inhomogeneous 3D nuclear distribution of TFs and by variation in DNA binding affinity such that locally elevated TF concentration allows low-affinity binding sites to be functional.

Duplexed aptamers: history, design, theory, and application to biosensing

Nucleic acid aptamers are single stranded DNA or RNA sequences that specifically bind a cognate ligand. In addition to their widespread use as stand-alone affinity binding reagents in analytical chemistry, aptamers have been engineered into a variety of ligand-specific biosensors, termed aptasensors. One of the most common aptasensor formats is the duplexed aptamer (DA). As defined herein, DAs are aptasensors containing two nucleic acid elements coupled via Watson-Crick base pairing: (i) an aptamer sequence, which serves as a ligand-specific receptor, and (ii) an aptamer-complementary element (ACE), such as a short DNA oligonucleotide, which is designed to hybridize to the aptamer. The ACE competes with ligand binding, such that DAs generate a signal upon ligand-dependent ACE-aptamer dehybridization. DAs possess intrinsic advantages over other aptasensor designs. For example, DA biosensing designs generalize across DNA and RNA aptamers, DAs are compatible with many readout methods, and DAs are inherently tunable on the basis of nucleic acid hybridization. However, despite their utility and popularity, DAs have not been well defined in the literature, leading to confusion over the differences between DAs and other aptasensor formats. In this review, we introduce a framework for DAs based on ACEs, and use this framework to distinguish DAs from other aptasensor formats and to categorize cis- and trans-DA designs. We then explore the ligand binding dynamics and chemical properties that underpin DA systems, which fall under conformational selection and induced fit models, and which mirror classical SN1 and SN2 models of nucleophilic substitution reactions. We further review a variety of in vitro and in vivo applications of DAs in the chemical and biological sciences, including riboswitches and riboregulators. Finally, we present future directions of DAs as ligand-responsive nucleic acids. Owing to their tractability, versatility and ease of engineering, DA biosensors bear a great potential for the development of new applications and technologies in fields ranging from analytical chemistry and mechanistic modeling to medicine and synthetic biology.

The Phosphorylated Estrogen Receptor α (ER) Cistrome Identifies a Subset of Active Enhancers Enriched for Direct ER-DNA Binding and the Transcription Factor GRHL2

Posttranslational modifications are key regulators of protein function, providing cues that can alter protein interactions and cellular location. Phosphorylation of estrogen receptor α (ER) at serine 118 (pS118-ER) occurs in response to multiple stimuli and is involved in modulating ER-dependent gene transcription. While the cistrome of ER is well established, surprisingly little is understood about how phosphorylation impacts ER-DNA binding activity. To define the pS118-ER cistrome, chromatin immunoprecipitation sequencing was performed on pS118-ER and ER in MCF-7 cells treated with estrogen. pS118-ER occupied a subset of ER binding sites which were associated with an active enhancer mark, acetylated H3K27. Unlike ER, pS118-ER sites were enriched in GRHL2 DNA binding motifs, and estrogen treatment increased GRHL2 recruitment to sites occupied by pS118-ER. Additionally, pS118-ER occupancy sites showed greater enrichment of full-length estrogen response elements relative to ER sites. In an in vitro DNA binding array of genomic binding sites, pS118-ER was more commonly associated with direct DNA binding events than indirect binding events. These results indicate that phosphorylation of ER at serine 118 promotes direct DNA binding at active enhancers and is a distinguishing mark for associated transcription factor complexes on chromatin.

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.

High-Density RNA Microarrays Synthesized In Situ by Photolithography

While high-density DNA microarrays have been available for over three decades, the synthesis of equivalent RNA microarrays has proven intractable until now. Herein we describe the first in situ synthesis of mixed-based, high-density RNA microarrays using photolithography and light-sensitive RNA phosphoramidites. With coupling efficiencies comparable to those of DNA monomers, RNA oligonucleotides at least 30 nucleotides long can now efficiently be prepared using modified phosphoramidite chemistry. A two-step deprotection route unmasks the phosphodiester, the exocyclic amines and the 2' hydroxyl. Hybridization and enzymatic assays validate the quality and the identity of the surface-bound RNA. We show that high-density is feasible by synthesizing a complex RNA permutation library with 262144 unique sequences. We also introduce DNA/RNA chimeric microarrays and explore their applications by mapping the sequence specificity of RNase HII.

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.

Accurate and sensitive quantification of protein-DNA binding affinity

Transcription factors (TFs) control gene expression by binding to genomic DNA in a sequence-specific manner. Mutations in TF binding sites are increasingly found to be associated with human disease, yet we currently lack robust methods to predict these sites. Here, we developed a versatile maximum likelihood framework named No Read Left Behind (NRLB) that infers a biophysical model of protein-DNA recognition across the full affinity range from a library of in vitro selected DNA binding sites. NRLB predicts human Max homodimer binding in near-perfect agreement with existing low-throughput measurements. It can capture the specificity of the p53 tetramer and distinguish multiple binding modes within a single sample. Additionally, we confirm that newly identified low-affinity enhancer binding sites are functional in vivo, and that their contribution to gene expression matches their predicted affinity. Our results establish a powerful paradigm for identifying protein binding sites and interpreting gene regulatory sequences in eukaryotic genomes.

Sequence-specific DNA binding Pyrrole-imidazole polyamides and their applications

Pyrrole-imidazole polyamides (Py-Im polyamides) are cell-permeable compounds that bind to the minor groove of double-stranded DNA in a sequence-specific manner without causing denaturation of the DNA. These compounds can be used to control gene expression and to stain specific sequences in cells. Here, we review the history, structural variations, and functional investigations of Py-Im polyamides.

Identification of Chemical Inducers of the Phosphate-Starvation Signaling Pathway in A. thaliana Using Chemical Genetics

In spite of its importance for agriculture and 30 years of genetic studies, the phosphate-starvation signaling pathway, that allows plants to detect, respond, and adapt to changes in the phosphate concentration of the rhizosphere, remains poorly known. Chemical genetics has been increasingly and successfully used as a complementary approach to genetics for the dissection of signaling pathways in diverse organisms. Screens can be designed to identify chemicals interfering specifically with a pathway of interest. We designed a screen that led to the discovery of the first chemical capable to induce specifically the phosphate-starvation signaling pathway in Arabidopsis thaliana. This procedure, described here, can be adapted for the discovery of inducers or repressors of other pathways.

Mapping the affinity landscape of Thrombin-binding aptamers on 2΄F-ANA/DNA chimeric G-Quadruplex microarrays

In situ fabricated nucleic acids microarrays are versatile and very high-throughput platforms for aptamer optimization and discovery, but the chemical space that can be probed against a given target has largely been confined to DNA, while RNA and non-natural nucleic acid microarrays are still an essentially uncharted territory. 2΄-Fluoroarabinonucleic acid (2΄F-ANA) is a prime candidate for such use in microarrays. Indeed, 2΄F-ANA chemistry is readily amenable to photolithographic microarray synthesis and its potential in high affinity aptamers has been recently discovered. We thus synthesized the first microarrays containing 2΄F-ANA and 2΄F-ANA/DNA chimeric sequences to fully map the binding affinity landscape of the TBA1 thrombin-binding G-quadruplex aptamer containing all 32 768 possible DNA-to-2΄F-ANA mutations. The resulting microarray was screened against thrombin to identify a series of promising 2΄F-ANA-modified aptamer candidates with K_ds significantly lower than that of the unmodified control and which were found to adopt highly stable, antiparallel-folded G-quadruplex structures. The solution structure of the TBA1 aptamer modified with 2΄F-ANA at position T3 shows that fluorine substitution preorganizes the dinucleotide loop into the proper conformation for interaction with thrombin. Overall, our work strengthens the potential of 2΄F-ANA in aptamer research and further expands non-genomic applications of nucleic acids microarrays.

Network-based approaches that exploit inferred transcription factor activity to analyze the impact of genetic variation on gene expression

Over the past decade, a number of methods have emerged for inferring protein-level transcription factor activities in individual samples based on prior information about the structure of the gene regulatory network. We discuss how this has enabled new methods for dissecting trans-acting mechanisms that underpin genetic variation in gene expression.

Combinatorial bZIP dimers display complex DNA-binding specificity landscapes

How transcription factor dimerization impacts DNA-binding specificity is poorly understood. Guided by protein dimerization properties, we examined DNA binding specificities of 270 human bZIP pairs. DNA interactomes of 80 heterodimers and 22 homodimers revealed that 72% of heterodimer motifs correspond to conjoined half-sites preferred by partnering monomers. Remarkably, the remaining motifs are composed of variably-spaced half-sites (12%) or 'emergent' sites (16%) that cannot be readily inferred from half-site preferences of partnering monomers. These binding sites were biochemically validated by EMSA-FRET analysis and validated in vivo by ChIP-seq data from human cell lines. Focusing on ATF3, we observed distinct cognate site preferences conferred by different bZIP partners, and demonstrated that genome-wide binding of ATF3 is best explained by considering many dimers in which it participates. Importantly, our compendium of bZIP-DNA interactomes predicted bZIP binding to 156 disease associated SNPs, of which only 20 were previously annotated with known bZIP motifs.

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.

Tailor-made transcriptional biosensors for optimizing microbial cell factories

Monitoring cellular behavior and eventually properly adapting cellular processes is key to handle the enormous complexity of today's metabolic engineering questions. Hence, transcriptional biosensors bear the potential to augment and accelerate current metabolic engineering strategies, catalyzing vital advances in industrial biotechnology. The development of such transcriptional biosensors typically starts with exploring nature's richness. Hence, in a first part, the transcriptional biosensor architecture and the various modi operandi are briefly discussed, as well as experimental and computational methods and relevant ontologies to search for natural transcription factors and their corresponding binding sites. In the second part of this review, various engineering approaches are reviewed to tune the main characteristics of these (natural) transcriptional biosensors, i.e., the response curve and ligand specificity, in view of specific industrial biotechnology applications, which is illustrated using success stories of transcriptional biosensor engineering.

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.

Dry and wet approaches for genome-wide functional annotation of conventional and unconventional transcriptional activators

Transcription factors (TFs) are master gene products that regulate gene expression in response to a variety of stimuli. They interact with DNA in a sequence-specific manner using a variety of DNA-binding domain (DBD) modules. This allows to properly position their second domain, called "effector domain", to directly or indirectly recruit positively or negatively acting co-regulators including chromatin modifiers, thus modulating preinitiation complex formation as well as transcription elongation. At variance with the DBDs, which are comprised of well-defined and easily recognizable DNA binding motifs, effector domains are usually much less conserved and thus considerably more difficult to predict. Also not so easy to identify are the DNA-binding sites of TFs, especially on a genome-wide basis and in the case of overlapping binding regions. Another emerging issue, with many potential regulatory implications, is that of so-called "moonlighting" transcription factors, i.e., proteins with an annotated function unrelated to transcription and lacking any recognizable DBD or effector domain, that play a role in gene regulation as their second job. Starting from bioinformatic and experimental high-throughput tools for an unbiased, genome-wide identification and functional characterization of TFs (especially transcriptional activators), we describe both established (and usually well affordable) as well as newly developed platforms for DNA-binding site identification. Selected combinations of these search tools, some of which rely on next-generation sequencing approaches, allow delineating the entire repertoire of TFs and unconventional regulators encoded by the any sequenced genome.

Sequence-specific DNA binding by long hairpin pyrrole-imidazole polyamides containing an 8-amino-3,6-dioxaoctanoic acid unit

With the aim of improving aqueous solubility, we designed and synthesized five N-methylpyrrole (Py)-N-methylimidazole (Im) polyamides capable of recognizing 9-bp sequences. Their DNA-binding affinities and sequence specificities were evaluated by SPR and Bind-n-Seq analyses. The design of polyamide 1 was based on a conventional model, with three consecutive Py or Im rings separated by a β-alanine to match the curvature and twist of long DNA helices. Polyamides 2 and 3 contained an 8-amino-3,6-dioxaoctanoic acid (AO) unit, which has previously only been used as a linker within linear Py-Im polyamides or between Py-Im hairpin motifs for tandem hairpin. It is demonstrated herein that AO also functions as a linker element that can extend to 2-bp in hairpin motifs. Notably, although the AO-containing unit can fail to bind the expected sequence, polyamide 4, which has two AO units facing each other in a hairpin form, successfully showed the expected motif and a KD value of 16nM was recorded. Polyamide 5, containing a β-alanine-β-alanine unit instead of the AO of polyamide 2, was synthesized for comparison. The aqueous solubilities and nuclear localization of three of the polyamides were also examined. The results suggest the possibility of applying the AO unit in the core of Py-Im polyamide compounds.

How Good Are Statistical Models at Approximating Complex Fitness Landscapes?

Fitness landscapes determine the course of adaptation by constraining and shaping evolutionary trajectories. Knowledge of the structure of a fitness landscape can thus predict evolutionary outcomes. Empirical fitness landscapes, however, have so far only offered limited insight into real-world questions, as the high dimensionality of sequence spaces makes it impossible to exhaustively measure the fitness of all variants of biologically meaningful sequences. We must therefore revert to statistical descriptions of fitness landscapes that are based on a sparse sample of fitness measurements. It remains unclear, however, how much data are required for such statistical descriptions to be useful. Here, we assess the ability of regression models accounting for single and pairwise mutations to correctly approximate a complex quasi-empirical fitness landscape. We compare approximations based on various sampling regimes of an RNA landscape and find that the sampling regime strongly influences the quality of the regression. On the one hand it is generally impossible to generate sufficient samples to achieve a good approximation of the complete fitness landscape, and on the other hand systematic sampling schemes can only provide a good description of the immediate neighborhood of a sequence of interest. Nevertheless, we obtain a remarkably good and unbiased fit to the local landscape when using sequences from a population that has evolved under strong selection. Thus, current statistical methods can provide a good approximation to the landscape of naturally evolving populations.

Deciphering the genomic targets of alkylating polyamide conjugates using high-throughput sequencing

Chemically engineered small molecules targeting specific genomic sequences play an important role in drug development research. Pyrrole-imidazole polyamides (PIPs) are a group of molecules that can bind to the DNA minor-groove and can be engineered to target specific sequences. Their biological effects rely primarily on their selective DNA binding. However, the binding mechanism of PIPs at the chromatinized genome level is poorly understood. Herein, we report a method using high-throughput sequencing to identify the DNA-alkylating sites of PIP-indole-seco-CBI conjugates. High-throughput sequencing analysis of conjugate 2: showed highly similar DNA-alkylating sites on synthetic oligos (histone-free DNA) and on human genomes (chromatinized DNA context). To our knowledge, this is the first report identifying alkylation sites across genomic DNA by alkylating PIP conjugates using high-throughput sequencing.

Express photolithographic DNA microarray synthesis with optimized chemistry and high-efficiency photolabile groups

Background

DNA microarrays are a core element of modern genomics research and medical diagnostics, allowing the simple and simultaneous determination of the relative abundances of hundreds of thousands to millions of genomic DNA or RNA sequences in a sample. Photolithographic in situ synthesis, using light projection from a digitally-controlled array of micromirrors, has been successful at both commercial and laboratory scales. The advantages of this synthesis method are its ability to reliably produce high-quality custom microarrays with a very high spatial density of DNA features using a compact device with few moving parts. The phosphoramidite chemistry used in photolithographic synthesis is similar to that used in conventional solid-phase synthesis of oligonucleotides, but some unique differences require an independent optimization of the synthesis chemistry to achieve fast and low-cost synthesis without compromising microarray quality.

Results

High microarray quality could be maintained while reducing coupling time to a few seconds using DCI activator. Five coupling activators were compared, which resulted in microarray hybridization signals following the order ETT > Activator 42 > DCI ≫ BTT ≫ pyridinium chloride, but only the use of DCI led to both high signal and highly uniform feature intensities. The photodeprotection time was also reduced to a few seconds by replacing the NPPOC photolabile group with the new thiophenyl-NPPOC group. Other chemical parameters, such as oxidation and washing steps were also optimized.

Conclusions

Highly optimized and microarray-specific phosphoramidite chemistry, along with the use of the very photosensitive thiophenyl-NPPOC protecting group allow for the synthesis of high-complexity DNA arrays using coupling times of 15 s and deprotection times of 9 s. The resulting overall cycle time (coupling to coupling) of about 50 s, results in a three-fold reduction in synthesis time.

Sequence-Dependent Fluorescence of Cy3- and Cy5-Labeled Double-Stranded DNA

The fluorescent intensity of Cy3 and Cy5 dyes is strongly dependent on the nucleobase sequence of the labeled oligonucleotides. Sequence-dependent fluorescence may significantly influence the data obtained from many common experimental methods based on fluorescence detection of nucleic acids, such as sequencing, PCR, FRET, and FISH. To quantify sequence dependent fluorescence, we have measured the fluorescence intensity of Cy3 and Cy5 bound to the 5' end of all 1024 possible double-stranded DNA 5mers. The fluorescence intensity was also determined for these dyes bound to the 5' end of fixed-sequence double-stranded DNA with a variable sequence 3' overhang adjacent to the dye. The labeled DNA oligonucleotides were made using light-directed, in situ microarray synthesis. The results indicate that the fluorescence intensity of both dyes is sensitive to all five bases or base pairs, that the sequence dependence is stronger for double- (vs single-) stranded DNA, and that the dyes are sensitive to both the adjacent dsDNA sequence and the 3'-ssDNA overhang. Purine-rich sequences result in higher fluorescence. The results can be used to estimate measurement error in experiments with fluorescent-labeled DNA, as well as to optimize the fluorescent signal by considering the nucleobase environment of the labeling cyanine dye.