Specificity landscapes of DNA binding molecules elucidate biological function

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.

Site-Specific Arrangement and Structure Determination of Minor Groove Binding Molecules in Self-Assembled Three-Dimensional DNA Crystals

The structural analysis of guest molecules in rationally designed and self-assembling DNA crystals has proven an elusive goal since its conception. Oligonucleotide frameworks provide an especially attractive route toward studying DNA-binding molecules by using three-dimensional lattices with defined sequence and structure. In this work, we site-specifically position a suite of minor groove binding molecules, and solve their structures via X-ray crystallography as a proof-of-principle toward scaffolding larger guest species. Two crystal motifs were used to precisely immobilize the molecules DAPI, Hoechst, and netropsin at defined positions in the lattice, allowing us to control occupancy within the crystal. We also solved the structure of a three-ring imidazole-pyrrole-pyrrole polyamide molecule, which sequence-specifically packs in an antiparallel dimeric arrangement within the minor groove. Finally, we engineered a crystal designed to position both netropsin and the polyamide at two distinct locations within the same lattice. Our work elucidates the design principles for the spatial arrangement of functional guests within lattices and opens new potential opportunities for the use of DNA crystals to display and structurally characterize small molecules, peptides, and ultimately proteins of unknown structure.

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.

Site-specific arrangement and structure determination of minor groove binding molecules in self-assembled three-dimensional DNA crystals

The structural analysis of guest molecules in rationally designed and self-assembling DNA crystals has proven elusive since its conception. Oligonucleotide frameworks provide an especially attractive route towards studying DNA-binding molecules by using three-dimensional lattices with defined sequence and structure. In this work, we site-specifically position a suite of minor groove binding molecules, and solve their structures via x-ray crystallography, as a proof-of-principle towards scaffolding larger guest species. Two crystal motifs were used to precisely immobilize the molecules DAPI, Hoechst, and netropsin at defined positions in the lattice, allowing us to control occupancy within the crystal. We also solved the structure of a three-ring imidazole-pyrrole-pyrrole polyamide molecule, which sequence-specifically packs in an anti-parallel dimeric arrangement within the minor groove. Finally, we engineered a crystal designed to position both netropsin and the polyamide at two distinct locations within the same lattice. Our work elucidates the design principles for the spatial arrangement of functional guests within lattices and opens new potential opportunities for the use of DNA crystals to display and structurally characterize small molecules, peptides, and ultimately proteins of unknown structure.

Hidden modes of DNA binding by human nuclear receptors

Human nuclear receptors (NRs) are a superfamily of ligand-responsive transcription factors that have central roles in cellular function. Their malfunction is linked to numerous diseases, and the ability to modulate their activity with synthetic ligands has yielded 16% of all FDA-approved drugs. NRs regulate distinct gene networks, however they often function from genomic sites that lack known binding motifs. Here, to annotate genomic binding sites of known and unexamined NRs more accurately, we use high-throughput SELEX to comprehensively map DNA binding site preferences of all full-length human NRs, in complex with their ligands. Furthermore, to identify non-obvious binding sites buried in DNA-protein interactomes, we develop MinSeq Find, a search algorithm based on the MinTerm concept from electrical engineering and digital systems design. The resulting MinTerm sequence set (MinSeqs) reveal a constellation of binding sites that more effectively annotate NR-binding profiles in cells. MinSeqs also unmask binding sites created or disrupted by 52,106 single-nucleotide polymorphisms associated with human diseases. By implicating druggable NRs as hidden drivers of multiple human diseases, our results not only reveal new biological roles of NRs, but they also provide a resource for drug-repurposing and precision medicine.

Disease-associated non-coding variants alter NKX2-5 DNA-binding affinity

Genome-wide association studies (GWAS) have mapped over 90 % of disease- or trait-associated variants within the non-coding genome, like cis-regulatory elements (CREs). Non-coding single nucleotide polymorphisms (SNPs) are genomic variants that can change how DNA-binding regulatory proteins, like transcription factors (TFs), interact with the genome and regulate gene expression. NKX2-5 is a TF essential for proper heart development, and mutations affecting its function have been associated with congenital heart diseases (CHDs). However, establishing a causal mechanism between non-coding genomic variants and human disease remains challenging. To address this challenge, we identified 8475 SNPs predicted to alter NKX2-5 DNA-binding using a position weight matrix (PWM)-based predictive model. Five variants were prioritized for in vitro validation; four of them are associated with traits and diseases that impact cardiovascular health. The impact of these variants on NKX2-5 binding was evaluated with electrophoretic mobility shift assay (EMSA) using purified recombinant NKX2-5 homeodomain. Binding curves were constructed to determine changes in binding between variant and reference alleles. Variants rs7350789, rs7719885, rs747334, and rs3892630 increased binding affinity, whereas rs61216514 decreased binding by NKX2-5 when compared to the reference genome. Our findings suggest that differential TF-DNA binding affinity can be key in establishing a causal mechanism of pathogenic variants.

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.

Simple synthesis of massively parallel RNA microarrays via enzymatic conversion from DNA microarrays

RNA catalytic and binding interactions with proteins and small molecules are fundamental elements of cellular life processes as well as the basis for RNA therapeutics and molecular engineering. In the absence of quantitative predictive capacity for such bioaffinity interactions, high throughput experimental approaches are needed to sufficiently sample RNA sequence space. Here we report on a simple and highly accessible approach to convert commercially available customized DNA microarrays of any complexity and density to RNA microarrays via a T7 RNA polymerase-mediated extension of photocrosslinked methyl RNA primers and subsequent degradation of the DNA templates.

RNA microarrays have many potential applications, but are difficult to produce. Here, the AUs present a method for converting commercial, customizable DNA microarrays into RNA microarrays using an accessible three-step process involving primer photocrosslinking, extension, and template degradation.

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.

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.

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.

The Road Not Taken with Pyrrole-Imidazole Polyamides: Off-Target Effects and Genomic Binding

The high sequence specificity of minor groove-binding N-methylpyrrole-N-methylimidazole polyamides have made significant advances in cancer and disease biology, yet there have been few comprehensive reports on their off-target effects, most likely as a consequence of the lack of available tools in evaluating genomic binding, an essential aspect that has gone seriously underexplored. Compared to other N-heterocycles, the off-target effects of these polyamides and their specificity for the DNA minor groove and primary base pair recognition require the development of new analytical methods, which are missing in the field today. This review aims to highlight the current progress in deciphering the off-target effects of these N-heterocyclic molecules and suggests new ways that next-generating sequencing can be used in addressing off-target effects.

Getting a handle on chemical probes of chomatin readers

The dynamic nature of histone post-translational modifications such as methylation or acetylation makes possible the alteration of disease associated epigenetic states through the manipulation of the associated epigenetic machinery. One approach is through small molecule perturbation. Chemical probes of epigenetic reader domains have been critical in improving our understanding of the biological consequences of modulating their targets, while also enabling the development of novel probe-based reagents. By appending a functional handle to a reader domain probe, a chemical toolbox of reagents can be created to facilitate chemiprecipitation of epigenetic complexes, evaluate probe selectivity, develop in vitro screening assays, visualize cellular target localization, enable target degradation and recruit epigenetic machinery to a site within the genome in a highly controlled fashion.

Targeted suppression of metastasis regulatory transcription factor SOX2 in various cancer cell lines using a sequence-specific designer pyrrole-imidazole polyamide

Metastasis, a deadly feature of cancer, compromises the prognosis and accounts for mortality in the majority of cancer patients. SOX2, a well-known pluripotency transcription factor, plays a central role in cell fate determination and has an overlapping role as a regulatory factor in tumorigenesis and metastasis. The demand is increasing for clinically useful strategies for artificial control of SOX2 expression and its complex transcription machinery in cancer cells. N-Methylpyrrole (Py) and N-methylimidazole (Im) polyamides are small programmable designer ligands that can be pre-programmed to selectively recognize DNA sequence and control endogenous gene expression. Herein, we evaluated the anticancer activity of a designer ligand (SOX2i). SOX2i remarkably altered the expression of SOX2 at the mRNA and protein level in human cancer cell lines such as SW620 (colorectal adenocarcinoma), MKN45 (gastric adenocarcinoma), MCF7 (breast carcinoma), U2OS (osteosarcoma) and other cancer cell lines of different origin and type. Genome-wide transcriptome analysis and cell-based assays showed SOX2 to be a downregulated upstream regulator that alters cell proliferation, cell cycle progression, metabolism and apoptotic pathway. Studies in the mouse model confirmed the anti-metastatic property of SOX2i. SOX2i inhibited the expression of genes associated with EMT and stemness. Moreover, Wnt-canonical signaling was found to be downregulated in the SOX2i-treated group. Our proof-of-concept study supports the potential of DNA-based programmable small molecules for controlling the key regulatory factors associated with tumorigenesis and metastasis.

Nkx2-5 Second Heart Field Target Gene Ccdc117 Regulates DNA Metabolism and Proliferation

The cardiac transcription factor Nkx2-5 is essential for normal outflow tract (OFT) and right ventricle (RV) development. Nkx2-5-/- null mouse embryos display severe OFT and RV hypoplasia and a single ventricle phenotype due to decreased proliferation of Second Heart Field (SHF) cells, a pool of cardiac progenitors present in anterior pharyngeal arch mesoderm at mid-gestation. However, definition of the precise role of Nkx2-5 in facilitating SHF expansion is incomplete. We have found that Nkx2-5 positively and directly regulates a novel target gene, Ccdc117, in cells of the SHF at these stages. The nuclear/mitotic spindle associated protein Ccdc117 interacts with the MIP18/MMS19 cytoplasmic iron-sulfur (FeS) cluster assembly (CIA) complex, which transfers critical FeS clusters to several key enzymes with functions in DNA repair and replication. Loss of cellular Ccdc117 expression results in reduced proliferation rates associated with a delay at the G1-S transition, decreased rates of DNA synthesis, and unresolved DNA damage. These results implicate a novel role for Nkx2-5 in the regulation of cell cycle events in the developing heart, through Ccdc117's interaction with elements of the CIA pathway and the facilitation of DNA replication during SHF expansion.

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.

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.

Targeting RNA in mammalian systems with small molecules

The recognition of RNA functions beyond canonical protein synthesis has challenged the central dogma of molecular biology. Indeed, RNA is now known to directly regulate many important cellular processes, including transcription, splicing, translation, and epigenetic modifications. The misregulation of these processes in disease has led to an appreciation of RNA as a therapeutic target. This potential was first recognized in bacteria and viruses, but discoveries of new RNA classes following the sequencing of the human genome have invigorated exploration of its disease-related functions in mammals. As stable structure formation is evolving as a hallmark of mammalian RNAs, the prospect of utilizing small molecules to specifically probe the function of RNA structural domains and their interactions is gaining increased recognition. To date, researchers have discovered bioactive small molecules that modulate phenotypes by binding to expanded repeats, microRNAs, G-quadruplex structures, and RNA splice sites in neurological disorders, cancers, and other diseases. The lessons learned from achieving these successes both call for additional studies and encourage exploration of the plethora of mammalian RNAs whose precise mechanisms of action remain to be elucidated. Efforts toward understanding fundamental principles of small molecule-RNA recognition combined with advances in methodology development should pave the way toward targeting emerging RNA classes such as long noncoding RNAs. Together, these endeavors can unlock the full potential of small molecule-based probing of RNA-regulated processes and enable us to discover new biology and underexplored avenues for therapeutic intervention in human disease. This article is categorized under: RNA Methods > RNA Analyses In Vitro and In Silico RNA Interactions with Proteins and Other Molecules > Small Molecule-RNA Interactions RNA in Disease and Development > RNA in Disease.

Global pairwise RNA interaction landscapes reveal core features of protein recognition

RNA–protein interactions permeate biology. Transcription, translation, and splicing all hinge on the recognition of structured RNA elements by RNA-binding proteins. Models of RNA–protein interactions are generally limited to short linear motifs and structures because of the vast sequence sampling required to access longer elements. Here, we develop an integrated approach that calculates global pairwise interaction scores from in vitro selection and high-throughput sequencing. We examine four RNA-binding proteins of phage, viral, and human origin. Our approach reveals regulatory motifs, discriminates between regulated and non-regulated RNAs within their native genomic context, and correctly predicts the consequence of mutational events on binding activity. We design binding elements that improve binding activity in cells and infer mutational pathways that reveal permissive versus disruptive evolutionary trajectories between regulated motifs. These coupling landscapes are broadly applicable for the discovery and characterization of protein–RNA recognition at single nucleotide resolution.

RNA–protein interactions often depend on the recognition of extended RNA elements but the identification of these motifs is challenging. Here, the authors present a global integrated approach to analyze RNA–protein binding landscapes, mapping extended RNA interaction motifs for four RNA-binding proteins.

Comprehensive, high-resolution binding energy landscapes reveal context dependencies of transcription factor binding

Transcription factors (TFs) are primary regulators of gene expression in cells, where they bind specific genomic target sites to control transcription. Quantitative measurements of TF-DNA binding energies can improve the accuracy of predictions of TF occupancy and downstream gene expression in vivo and shed light on how transcriptional networks are rewired throughout evolution. Here, we present a sequencing-based TF binding assay and analysis pipeline (BET-seq, for Binding Energy Topography by sequencing) capable of providing quantitative estimates of binding energies for more than one million DNA sequences in parallel at high energetic resolution. Using this platform, we measured the binding energies associated with all possible combinations of 10 nucleotides flanking the known consensus DNA target interacting with two model yeast TFs, Pho4 and Cbf1. A large fraction of these flanking mutations change overall binding energies by an amount equal to or greater than consensus site mutations, suggesting that current definitions of TF binding sites may be too restrictive. By systematically comparing estimates of binding energies output by deep neural networks (NNs) and biophysical models trained on these data, we establish that dinucleotide (DN) specificities are sufficient to explain essentially all variance in observed binding behavior, with Cbf1 binding exhibiting significantly more nonadditivity than Pho4. NN-derived binding energies agree with orthogonal biochemical measurements and reveal that dynamically occupied sites in vivo are both energetically and mutationally distant from the highest affinity sites.

Reprogramming cell fate with artificial transcription factors

Transcription factors (TFs) reprogram cell states by exerting control over gene regulatory networks and the epigenetic landscape of a cell. Artificial transcription factors (ATFs) are designer regulatory proteins comprised of modular units that can be customized to overcome challenges faced by natural TFs in establishing and maintaining desired cell states. Decades of research on DNA-binding proteins and synthetic molecules has provided a molecular toolkit for ATF design and the construction of genome-scale libraries of ATFs capable of phenotypic manipulation and reprogramming of cell states. Here, we compare the unique strengths and limitations of different ATF platforms, highlight the advantages of cooperative assembly, and present the potential of ATF libraries in revealing gene regulatory networks that govern cell fate choices.

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.

Conservation of Specificity in Two Low-Specificity Proteins

Many regulatory proteins bind peptide regions of target proteins and modulate their activity. Such regulatory proteins can often interact with highly diverse target peptides. In many instances, it is not known if the peptide-binding interface discriminates targets in a biological context, or whether biological specificity is achieved exclusively through external factors such as subcellular localization. We used an evolutionary biochemical approach to distinguish these possibilities for two such low-specificity proteins: S100A5 and S100A6. We used isothermal titration calorimetry to study the binding of peptides with diverse sequence and biochemistry to human S100A5 and S100A6. These proteins bound distinct, but overlapping, sets of peptide targets. We then studied the peptide binding properties of orthologs sampled from across five amniote species. Binding specificity was conserved along all lineages, for the last 320 million years, despite the low specificity of each protein. We used ancestral sequence reconstruction to determine the binding specificity of the last common ancestor of the paralogs. The ancestor bound the entire set of peptides bound by modern S100A5 and S100A6 proteins, suggesting that paralog specificity evolved via subfunctionalization. To rule out the possibility that specificity is conserved because it is difficult to modify, we identified a single historical mutation that, when reverted in human S100A5, gave it the ability to bind an S100A6-specific peptide. These results reveal strong evolutionary constraints on peptide binding specificity. Despite being able to bind a large number of targets, the specificity of S100 peptide interfaces is likely important for the biology of these proteins.

Synthetic transcription elongation factors license transcription across repressive chromatin

The release of paused RNA polymerase II into productive elongation is highly regulated, especially at genes that affect human development and disease. To exert control over this rate-limiting step, we designed sequence-specific synthetic transcription elongation factors (Syn-TEFs). These molecules are composed of programmable DNA-binding ligands flexibly tethered to a small molecule that engages the transcription elongation machinery. By limiting activity to targeted loci, Syn-TEFs convert constituent modules from broad-spectrum inhibitors of transcription into gene-specific stimulators. Here we present Syn-TEF1, a molecule that actively enables transcription across repressive GAA repeats that silence frataxin expression in Friedreich's ataxia, a terminal neurodegenerative disease with no effective therapy. The modular design of Syn-TEF1 defines a general framework for developing a class of molecules that license transcription elongation at targeted genomic loci.

A synthetic DNA-binding inhibitor of SOX2 guides human induced pluripotent stem cells to differentiate into mesoderm

Targeted differentiation of human induced pluripotent stem cells (hiPSCs) using only chemicals would have value-added clinical potential in the regeneration of complex cell types including cardiomyocytes. Despite the availability of several chemical inhibitors targeting proteins involved in signaling pathways, no bioactive synthetic DNA-binding inhibitors, targeting key cell fate-controlling genes such as SOX2, are yet available. Here, we demonstrate a novel DNA-based chemical approach to guide the differentiation of hiPSCs using pyrrole-imidazole polyamides (PIPs), which are sequence-selective DNA-binding synthetic molecules. Harnessing knowledge about key transcriptional changes during the induction of cardiomyocyte, we developed a DNA-binding inhibitor termed PIP-S2, targeting the 5'-CTTTGTT-3' and demonstrated that inhibition of SOX2-DNA interaction by PIP-S2 triggers the mesoderm induction in hiPSCs. Genome-wide gene expression analyses revealed that PIP-S2 induced mesoderm by targeted alterations in SOX2-associated gene regulatory networks. Also, employment of PIP-S2 along with a Wnt/β-catenin inhibitor successfully generated spontaneously contracting cardiomyocytes, validating our concept that DNA-binding inhibitors could drive the directed differentiation of hiPSCs. Because PIPs can be fine-tuned to target specific DNA sequences, our DNA-based approach could be expanded to target and regulate key transcription factors specifically associated with desired cell types.

Massively Parallel Biophysical Analysis of CRISPR-Cas Complexes on Next Generation Sequencing Chips

CRISPR-Cas nucleoproteins target foreign DNA via base pairing with a crRNA. However, a quantitative description of protein binding and nuclease activation at off-target DNA sequences remains elusive. Here, we describe a chip-hybridized association-mapping platform (CHAMP) that repurposes next-generation sequencing chips to simultaneously measure the interactions between proteins and ∼10⁷ unique DNA sequences. Using CHAMP, we provide the first comprehensive survey of DNA recognition by a type I-E CRISPR-Cas (Cascade) complex and Cas3 nuclease. Analysis of mutated target sequences and human genomic DNA reveal that Cascade recognizes an extended protospacer adjacent motif (PAM). Cascade recognizes DNA with a surprising 3-nt periodicity. The identity of the PAM and the PAM-proximal nucleotides control Cas3 recruitment by releasing the Cse1 subunit. These findings are used to develop a model for the biophysical constraints governing off-target DNA binding. CHAMP provides a framework for high-throughput, quantitative analysis of protein-DNA interactions on synthetic and genomic DNA. PAPERCLIP.

Modeling protein-DNA binding via high-throughput in vitro technologies

Protein-DNA binding plays a central role in gene regulation and by that in all processes in the living cell. Novel experimental and computational approaches facilitate better understanding of protein-DNA binding preferences via high-throughput measurement of protein binding to a large number of DNA sequences and inference of binding models from them. Here we review the state of the art in measuring protein-DNA binding in vitro, emphasizing the advantages and limitations of different technologies. In addition, we describe models for representing protein-DNA binding preferences and key computational approaches to learn those from high-throughput data. Using large experimental data sets, we test the performance of different models based on different measuring techniques. We conclude with pertinent open problems.

Knockdown of Broad-Complex Gene Expression of Bombyx mori by Oligopyrrole Carboxamides Enhances Silk Production

Bombyx mori (B. mori) is important due to its major role in the silk production. Though DNA binding ligands often influence gene expression, no attempt has been made to exploit their use in sericulture. The telomeric heterochromatin of B. mori is enriched with 5'-TTAGG-3' sequences. These sequences were also found to be present in several genes in the euchromatic regions. We examined three synthetic oligopyrrole carboxamides that target 5'-TTAGG-3' sequences in controlling the gene expression in B. mori. The ligands did not show any defect or feeding difference in the larval stage, crucial for silk production. The ligands caused silencing of various isoforms of the broad-complex transcription factor and cuticle proteins which resulted in late pupal developmental defects. Furthermore, treatment with such drugs resulted in statistically enhanced cocoon weight, shell weight, and silk yield. This study shows for the first time use of oligopyrrole carboxamide drugs in controlling gene expression in B. mori and their long term use in enhancing silk production.

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.

A thousand empirical adaptive landscapes and their navigability

The adaptive landscape is an iconic metaphor that pervades evolutionary biology. It was mostly applied in theoretical models until recent years, when empirical data began to allow partial landscape reconstructions. Here, we exhaustively analyse 1,137 complete landscapes from 129 eukaryotic species, each describing the binding affinity of a transcription factor to all possible short DNA sequences. We find that the navigability of these landscapes through single mutations is intermediate to that of additive and shuffled null models, suggesting that binding affinity-and thereby gene expression-is readily fine-tuned via mutations in transcription factor binding sites. The landscapes have few peaks that vary in their accessibility and in the number of sequences they contain. Binding sites in the mouse genome are enriched in sequences found in the peaks of especially navigable landscapes and the genetic diversity of binding sites in yeast increases with the number of sequences in a peak. Our findings suggest that landscape navigability may have contributed to the enormous success of transcriptional regulation as a source of evolutionary adaptations and innovations.

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.

Sliding on DNA: From Peptides to Small Molecules

Many DNA binding proteins utilize one-dimensional (1D) diffusion along DNA to accelerate their DNA target recognition. Although 1D diffusion of proteins along DNA has been studied for decades, a quantitative understanding is only beginning to emerge and few chemical tools are available to apply 1D diffusion as a design principle. Recently, we discovered that peptides can bind and slide along DNA-even transporting cargo along DNA. Such molecules are known as molecular sleds. Here, to advance our understanding of structure-function relationships governing sequence nonspecific DNA interaction of natural molecular sleds and to explore the potential for controlling sliding activity, we test the DNA binding and sliding activities of chemically modified peptides and analogs, and show that synthetic small molecules can slide on DNA. We found new ways to control molecular sled activity, novel small-molecule synthetic sleds, and molecular sled activity in N-methylpyrrole/N-methylimidazole polyamides that helps explain how these molecules locate rare target sites.

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.

Integrated analysis of RNA-binding protein complexes using in vitro selection and high-throughput sequencing and sequence specificity landscapes (SEQRS)

RNA-binding proteins (RBPs) collaborate to control virtually every aspect of RNA function. Tremendous progress has been made in the area of global assessment of RBP specificity using next-generation sequencing approaches both in vivo and in vitro. Understanding how protein-protein interactions enable precise combinatorial regulation of RNA remains a significant problem. Addressing this challenge requires tools that can quantitatively determine the specificities of both individual proteins and multimeric complexes in an unbiased and comprehensive way. One approach utilizes in vitro selection, high-throughput sequencing, and sequence-specificity landscapes (SEQRS). We outline a SEQRS experiment focused on obtaining the specificity of a multi-protein complex between Drosophila RBPs Pumilio (Pum) and Nanos (Nos). We discuss the necessary controls in this type of experiment and examine how the resulting data can be complemented with structural and cell-based reporter assays. Additionally, SEQRS data can be integrated with functional genomics data to uncover biological function. Finally, we propose extensions of the technique that will enhance our understanding of multi-protein regulatory complexes assembled onto RNA.

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.

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.

Zinc finger protein binding to DNA: an energy perspective using molecular dynamics simulation and free energy calculations on mutants of both zinc finger domains and their specific DNA bases

Energy calculations based on MM-GBSA were employed to study various zinc finger protein (ZF) motifs binding to DNA. Mutants of both the DNA bound to their specific amino acids were studied. Calculated energies gave evidence for a relationship between binding energy and affinity of ZF motifs to their sites on DNA. ΔG values were -15.82(12), -3.66(12), and -12.14(11.6) kcal/mol for finger one, finger two, and finger three, respectively. The mutations in the DNA bases reduced the value of the negative energies of binding (maximum value for ΔΔG = 42Kcal/mol for F1 when GCG mutated to GGG, and ΔΔG = 22 kcal/mol for F2, the loss in total energy of binding originated in the loss in electrostatic energies upon mutation (r = .98). The mutations in key amino acids in the ZF motif in positions-1, 2, 3, and 6 showed reduced binding energies to DNA with correlation coefficients between total free energy and electrostatic was .99 and with Van der Waal was .93. Results agree with experimentally found selectivity which showed that Arginine in position-1 is specific to G, while Aspartic acid (D) in position 2 plays a complicated role in binding. There is a correlation between the MD calculated free energies of binding and those obtained experimentally for prepared ZF motifs bound to triplet bases in other reports (), our results may help in the design of ZF motifs based on the established recognition codes based on energies and contributing energies to the total energy.

Genome-wide Mapping of Drug-DNA Interactions in Cells with COSMIC (Crosslinking of Small Molecules to Isolate Chromatin)

The genome is the target of some of the most effective chemotherapeutics, but most of these drugs lack DNA sequence specificity, which leads to dose-limiting toxicity and many adverse side effects. Targeting the genome with sequence-specific small molecules may enable molecules with increased therapeutic index and fewer off-target effects. N-methylpyrrole/N-methylimidazole polyamides are molecules that can be rationally designed to target specific DNA sequences with exquisite precision. And unlike most natural transcription factors, polyamides can bind to methylated and chromatinized DNA without a loss in affinity. The sequence specificity of polyamides has been extensively studied in vitro with cognate site identification (CSI) and with traditional biochemical and biophysical approaches, but the study of polyamide binding to genomic targets in cells remains elusive. Here we report a method, the crosslinking of small molecules to isolate chromatin (COSMIC), that identifies polyamide binding sites across the genome. COSMIC is similar to chromatin immunoprecipitation (ChIP), but differs in two important ways: (1) a photocrosslinker is employed to enable selective, temporally-controlled capture of polyamide binding events, and (2) the biotin affinity handle is used to purify polyamide-DNA conjugates under semi-denaturing conditions to decrease DNA that is non-covalently bound. COSMIC is a general strategy that can be used to reveal the genome-wide binding events of polyamides and other genome-targeting chemotherapeutic agents.

Peptidylprolyl Isomerase Pin1 Directly Enhances the DNA Binding Functions of Estrogen Receptor α

The transcriptional activity of estrogen receptor α (ERα), the key driver of breast cancer proliferation, is enhanced by multiple cellular interactions, including phosphorylation-dependent interaction with Pin1, a proline isomerase, which mediates cis-trans isomerization of the N-terminal Ser(P)(118)-Pro(119) in the intrinsically disordered AF1 (activation function 1) domain of ERα. Because both ERα and Pin1 have multiple cellular partners, it is unclear how Pin1 assists in the regulation of ERα transactivation mechanisms and whether the functional effects of Pin1 on ERα signaling are direct or indirect. Here, we tested the specific action of Pin1 on an essential step in ERα transactivation, binding to specific DNA sites. DNA binding analysis demonstrates that stable overexpression of Pin1 increases endogenous ERα DNA binding activity when activated by estrogen but not by tamoxifen or EGF. Increased DNA binding affinity is a direct effect of Pin1 on ERα because it is observed in solution-based assays with purified components. Further, our data indicate that isomerization is required for Pin1-modulation of ERα-DNA interactions. In an unbiased in vitro DNA binding microarray with hundreds of thousands of permutations of ERα-binding elements, Pin1 selectively enhances the binding affinity of ERα to consensus DNA elements. These studies reveal that Pin1 isomerization of phosphorylated ERα can directly regulate the function of the adjacent DNA binding domain, and this interaction is further modulated by ligand binding in the ligand-binding domain, providing evidence for Pin1-dependent allosteric regulation of ERα function.

Targets of the Sex Inducer homeodomain proteins are required for fungal development and virulence in Cryptococcus neoformans

In the yeast Saccharomyces cerevisiae, the regulation of cell types by homeodomain transcription factors is a key paradigm; however, many questions remain regarding this class of developmental regulators in other fungi. In the human fungal pathogen Cryptococcus neoformans, the homeodomain transcription factors Sxi1α and Sxi2a are required for sexual development that produces infectious spores, but the molecular mechanisms by which they drive this process are unknown. To better understand homeodomain control of fungal development, we determined the targets of the Sxi2a-Sxi1α heterodimer using whole genome expression analyses paired with in silico and in vitro binding site identification methods. We identified Sxi-regulated genes that contained a site bound directly by the Sxi proteins that is required for full regulation in vivo. Among the targets of the Sxi2a-Sxi1α complex were many genes known to be involved in sexual reproduction, as well as several well-studied virulence genes. Our findings suggest that genes involved in sexual development are also important in mammalian disease. Our work advances the understanding of how homeodomain transcription factors control complex developmental events and suggests an intimate link between fungal development and virulence.

Multivalent DNA recognition by self-assembled clusters: deciphering structural effects by fragments screening and evaluation as siRNA vectors

Fragment self-assembly was used for producing clusters with a variety of scaffolds and ligands, and an effective siRNA vector was identified.

Identification of a small molecule that turns ON the pluripotency gene circuitry in human fibroblasts

A nontransgenic approach to reprogram mouse somatic cells into induced pluripotent stem cells using only small molecules got achieved to propose a potential clinical-friendly cellular reprogramming strategy. Consequently, the screening and identification of small molecules capable of inducing pluripotency genes in human cells are increasingly a focus of research. Because cellular reprogramming is multifactorial in nature, there is a need for versatile small molecules capable of modulating the complicated gene networks associated with pluripotency. We have developed a targeting small molecule called SAHA-PIP comprising the histone deacetylase inhibitor SAHA and the sequence-specific DNA binding pyrrole-imidazole polyamides for modulating distinct gene networks. Here, we report the identification of a SAHA-PIP termed Ì that could trigger genome-wide epigenetic reprogramming and turn ON the typically conserved core pluripotency gene network. Through independent lines of evidence, we report for the first time a synthetic small molecule inducer that target and activate the OCT-3/4 regulated pluripotency genes in human dermal fibroblasts.

Placental Nkx2-5 and target gene expression in early-onset and severe preeclampsia

OBJECTIVE

Preeclampsia (PE) affects 2-8% of pregnancies worldwide and is a significant source of maternal and neonatal morbidity and mortality. However, the mechanisms underlying PE are poorly understood and major questions regarding etiology and risk factors remain to be addressed. Our objective was to examine whether abnormal expression of the cardiovascular developmental transcription factor, Nkx2-5, was associated with early onset and severe preeclampsia (EOSPE).

METHODS

Using qPCR and immunohistochemical assay, we examined expression of Nkx2-5 and target gene expression in EOSPE and control placental tissue. We tested resulting mechanistic hypotheses in cultured cells using shRNA knockdown, qPCR, and western blot.

RESULTS

Nkx2-5 is highly expressed in racially disparate fashion (Caucasians > African Americans) in a subset of early EOSPE placentae. Nkx2-5 mRNA expression is highly correlated (Caucasians > African Americans) to mRNA expression of the preeclampsia marker sFlt-1, and of the Nkx2-5 target and RNA splicing factor, Sam68. Knockdown of Sam68 expression in cultured cells significantly impacts sFlt-1 mRNA isoform generation in vitro, supporting a mechanistic hypothesis that Nkx2-5 impacts EOSPE severity in a subset of patients via upregulation of Sam68 to increase sFlt-1 expression. Expression of additional Nkx2-5 targets potentially regulating metabolic stress response is also elevated in a racially disparate fashion in EOSPE.

CONCLUSIONS

Expression of Nkx2-5 and its target genes may directly influence the genesis and racially disparate severity, and define a mechanistically distinct subclass of EOSPE.

Mapping polyamide-DNA interactions in human cells reveals a new design strategy for effective targeting of genomic sites

Targeting the genome with sequence-specific synthetic molecules is a major goal at the interface of chemistry, biology, and personalized medicine. Pyrrole/imidazole-based polyamides can be rationally designed to target specific DNA sequences with exquisite precision in vitro; yet, the biological outcomes are often difficult to interpret using current models of binding energetics. To directly identify the binding sites of polyamides across the genome, we designed, synthesized, and tested polyamide derivatives that enabled covalent crosslinking and localization of polyamide-DNA interaction sites in live human cells. Bioinformatic analysis of the data reveals that clustered binding sites, spanning a broad range of affinities, best predict occupancy in cells. In contrast to the prevailing paradigm of targeting single high-affinity sites, our results point to a new design principle to deploy polyamides and perhaps other synthetic molecules to effectively target desired genomic sites in vivo.