Using protein-binding microarrays to study transcription factor specificity: homologs, isoforms and complexes

Experimental approaches to investigate biophysical interactions between homeodomain transcription factors and DNA

Homeodomain transcription factors (TFs) bind to specific DNA sequences to regulate the expression of target genes. Structural work has provided insight into molecular identities and aided in unraveling structural features of these TFs. However, the detailed affinity and specificity by which these TFs bind to DNA sequences is still largely unknown. Qualitative methods, such as DNA footprinting, Electrophoretic Mobility Shift Assays (EMSAs), Systematic Evolution of Ligands by Exponential Enrichment (SELEX), Bacterial One Hybrid (B1H) systems, Surface Plasmon Resonance (SPR), and Protein Binding Microarrays (PBMs) have been widely used to investigate the biochemical characteristics of TF-DNA binding events. In addition to these qualitative methods, bioinformatic approaches have also assisted in TF binding site discovery. Here we discuss the advantages and limitations of these different approaches, as well as the benefits of utilizing more quantitative approaches, such as Mechanically Induced Trapping of Molecular Interactions (MITOMI), Microscale Thermophoresis (MST) and Isothermal Titration Calorimetry (ITC), in determining the biophysical basis of binding specificity of TF-DNA complexes and improving upon existing computational approaches aimed at affinity predictions.

Protein target search diffusion-association/dissociation free energy landscape around DNA binding site with flanking sequences

In this work we present a minimal structure-based model of protein diffusional search along local DNA amid protein binding and unbinding events on the DNA, taking into account protein-DNA electrostatic interactions and hydrogen-bonding (HB) interactions or contacts at the interface. We accordingly constructed the protein diffusion-association/dissociation free energy surface and mapped it to 1D as the protein slides along DNA, maintaining the protein-DNA interfacial HB contacts that presumably dictate the DNA sequence information detection. Upon DNA helical path correction, the protein 1D diffusion rates along local DNA can be physically derived to be consistent with experimental measurements. We also show that the sequence-dependent protein sliding or stepping patterns along DNA are regulated by collective interfacial HB dynamics, which also determines the ruggedness of the protein diffusion free energy landscape on the local DNA. In comparison, protein association or binding with DNA are generically dictated by the protein-DNA electrostatic interactions, with an interaction zone of nanometers around DNA. Extra degrees of freedom (DOFs) of the protein such as rotations and conformational fluctuations can be well accommodated within the protein-DNA electrostatic interaction zone. As such we demonstrate that the protein binding or association free energy profiling along DNA smoothens over the 1D diffusion free energy landscape, which leads to population variations for an order of magnitude upon a marginal free energetic smoothening around the specific or consensus sites. We further show that the protein unbinding or dissociation from a comparatively high-binding affinity DNA site is dominated by lateral diffusion to the flanking low-affinity sites. The results predict that experimental characterizations on the relative protein-DNA binding affinities or population profiling on the DNA are systematically and physically impacted by the extra DOFs of protein motions aside from 1D translation or helical tracking, as well as from flanking DNA sequences due to protein 1D diffusion and nonspecific binding/unbinding.

Using Synthetic DNA Libraries to Investigate Chromatin and Gene Regulation

Despite the recent explosion in genome-wide studies in chromatin and gene regulation, we are still far from extracting a set of genetic rules that can predict the function of the regulatory genome. One major reason for this deficiency is that gene regulation is a multi-layered process that involves an enormous variable space, which cannot be fully explored using native genomes. This problem can be partially solved by introducing synthetic DNA libraries into cells, a method that can test the regulatory roles of thousands to millions of sequences with limited variables. Here, we review recent applications of this method to study transcription factor (TF) binding, nucleosome positioning, and transcriptional activity. We discuss the design principles, experimental procedures, and major findings from these studies and compare the pros and cons of different approaches.

Prediction and Experimental Validation of a New Salinity-Responsive Cis-Regulatory Element (CRE) in a Tilapia Cell Line

Transcriptional regulation is a major mechanism by which organisms integrate gene x environment interactions. It can be achieved by coordinated interplay between cis-regulatory elements (CREs) and transcription factors (TFs). Euryhaline tilapia (Oreochromis mossambicus) tolerate a wide range of salinity and thus are an appropriate model to examine transcriptional regulatory mechanisms during salinity stress in fish. Quantitative proteomics in combination with the transcription inhibitor actinomycin D revealed 19 proteins that are transcriptionally upregulated by hyperosmolality in tilapia brain (OmB) cells. We searched the extended proximal promoter up to intron1 of each corresponding gene for common motifs using motif discovery tools. The top-ranked motif identified (STREME1) represents a binding site for the Forkhead box TF L1 (FoxL1). STREME1 function during hyperosmolality was experimentally validated by choosing two of the 19 genes, chloride intracellular channel 2 (clic2) and uridine phosphorylase 1 (upp1), that are enriched in STREME1 in their extended promoters. Transcriptional induction of these genes during hyperosmolality requires STREME1, as evidenced by motif mutagenesis. We conclude that STREME1 represents a new functional CRE that contributes to gene x environment interactions during salinity stress in tilapia. Moreover, our results indicate that FoxL1 family TFs are contribute to hyperosmotic induction of genes in euryhaline fish.

Transcriptional control of Arabidopsis seed development

The entire process of embryo development is under the tight control of various transcription factors. Together with other proteins, they act in a combinatorial manner and control distinct events during embryo development. Seed development is a complex process that proceeds through sequences of events regulated by the interplay of various genes, prominent among them being the transcription factors (TFs). The members of WOX, HD-ZIP III, ARF, and CUC families have a preferential role in embryonic patterning. While WOX TFs are required for initiating body axis, HD-ZIP III TFs and CUCs establish bilateral symmetry and SAM. And ARF5 performs a major role during embryonic root, ground tissue, and vasculature development. TFs such as LEC1, ABI3, FUS3, and LEC2 (LAFL) are considered the master regulators of seed maturation. Furthermore, several new TFs involved in seed storage reserves and dormancy have been identified in the last few years. Their association with those master regulators has been established in the model plant Arabidopsis. Also, using chromatin immunoprecipitation (ChIP) assay coupled with transcriptomics, genome-wide target genes of these master regulators have recently been proposed. Many seed-specific genes, including those encoding oleosins and albumins, have appeared as the direct target of LAFL. Also, several other TFs act downstream of LAFL TFs and perform their function during maturation. In this review, the function of different TFs in different phases of early embryogenesis and maturation is discussed in detail, including information about their genetic and molecular interactors and target genes. Such knowledge can further be leveraged to understand and manipulate the regulatory mechanisms involved in seed development. In addition, the genomics approaches and their utilization to identify TFs aiming to study embryo development are discussed.

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.

Gene Ontology representation for transcription factor functions

Transcription plays a central role in defining the identity and functionalities of cells, as well as in their responses to changes in the cellular environment. The Gene Ontology (GO) provides a rigorously defined set of concepts that describe the functions of gene products. A GO annotation is a statement about the function of a particular gene product, represented as an association between a gene product and the biological concept a GO term defines. Critically, each GO annotation is based on traceable scientific evidence. Here, we describe the different GO terms that are associated with proteins involved in transcription and its regulation, focusing on the standard of evidence required to support these associations. This article is intended to help users of GO annotations understand how to interpret the annotations and can contribute to the consistency of GO annotations. We distinguish between three classes of activities involved in transcription or directly regulating it - general transcription factors, DNA-binding transcription factors, and transcription co-regulators.

High-Throughput Analysis of the Cell and DNA Site-Specific Binding of Native NF-κB Dimers Using Nuclear Extract Protein-Binding Microarrays (NextPBMs)

Nuclear factor-kappa B (NF-κB) transcription factors coordinate gene expression in response to a broad array of cellular signals. In vertebrates, there are five NF-κB proteins (c-Rel, RelA/p65, RelB, p50, and p52) that can form various dimeric combinations exhibiting both common and dimer-specific DNA-binding specificity. In this chapter, we describe the use of the nuclear extract protein-binding microarray (nextPBM), a high-throughput method to characterize the DNA binding of transcription factors present in cell nuclear extracts. NextPBMs allow for sensitive analysis of the DNA binding of NF-κB dimers and their interactions with cell-specific cofactors.

Custom DNA Microarrays Reveal Diverse Binding Preferences of Proteins and Small Molecules to Thousands of G-Quadruplexes

Single-stranded DNA (ssDNA) containing four guanine repeats can form G-quadruplex (G4) structures. While cellular proteins and small molecules can bind G4s, it has been difficult to broadly assess their DNA-binding specificity. Here, we use custom DNA microarrays to examine the binding specificities of proteins, small molecules, and antibodies across ∼15,000 potential G4 structures. Molecules used include fluorescently labeled pyridostatin (Cy5-PDS, a small molecule), BG4 (Cy5-BG4, a G4-specific antibody), and eight proteins (GST-tagged nucleolin, IGF2, CNBP, FANCJ, PIF1, BLM, DHX36, and WRN). Cy5-PDS and Cy5-BG4 selectively bind sequences known to form G4s, confirming their formation on the microarrays. Cy5-PDS binding decreased when G4 formation was inhibited using lithium or when ssDNA features on the microarray were made double-stranded. Similar conditions inhibited the binding of all other molecules except for CNBP and PIF1. We report that proteins have different G4-binding preferences suggesting unique cellular functions. Finally, competition experiments are used to assess the binding specificity of an unlabeled small molecule, revealing the structural features in the G4 required to achieve selectivity. These data demonstrate that the microarray platform can be used to assess the binding preferences of molecules to G4s on a broad scale, helping to understand the properties that govern molecular recognition.

The physical basis and practical consequences of biological promiscuity

Proteins interact with metabolites, nucleic acids, and other proteins to orchestrate the myriad catalytic, structural and regulatory functions that support life from the simplest microbes to the most complex multicellular organisms. These molecular interactions are often exquisitely specific, but never perfectly so. Adventitious "promiscuous" interactions are ubiquitous due to the thousands of macromolecules and small molecules crowded together in cells. Such interactions may perturb protein function at the molecular level, but as long as they do not compromise organismal fitness, they will not be removed by natural selection. Although promiscuous interactions are physiologically irrelevant, they are important because they can provide a vast reservoir of potential functions that can provide the starting point for evolution of new functions, both in nature and in the laboratory.

NextPBM: a platform to study cell-specific transcription factor binding and cooperativity

High-throughput (HT) in vitro methods for measuring protein-DNA binding have become invaluable for characterizing transcription factor (TF) complexes and modeling gene regulation. However, current methods do not utilize endogenous proteins and, therefore, do not quantify the impact of cell-specific post-translational modifications (PTMs) and cooperative cofactors. We introduce the HT nextPBM (nuclear extract protein-binding microarray) approach to study DNA binding of native cellular TFs that accounts for PTMs and cell-specific cofactors. We integrate immune-depletion and phosphatase treatment steps into our nextPBM pipeline to characterize the impact of cofactors and phosphorylation on TF binding. We analyze binding of PU.1/SPI1 and IRF8 from human monocytes, delineate DNA-sequence determinants for their cooperativity, and show how PU.1 affinity correlates with enhancer status and the presence of cooperative and collaborative cofactors. We describe how nextPBMs, and our accompanying computational framework, can be used to discover cell-specific cofactors, screen for synthetic cooperative DNA elements, and characterize TF cooperativity.

DNA-binding landscape of IRF3, IRF5 and IRF7 dimers: implications for dimer-specific gene regulation

Transcription factors IRF3, IRF5 and IRF7 (IRF3/5/7) have overlapping, yet distinct, roles in the mammalian response to pathogens. To examine the role that DNA-binding specificity plays in delineating IRF3/5/7-specific gene regulation we used protein-binding microarrays (PBMs) to characterize the DNA binding of IRF3/5/7 homodimers. We identified both common and dimer-specific DNA binding sites, and show that DNA-binding differences can translate into dimer-specific gene regulation. Central to the antiviral response, IRF3/5/7 regulate type I interferon (IFN) genes. We show that IRF3 and IRF7 bind to many interferon-stimulated response element (ISRE)-type sites in the virus-response elements (VREs) of IFN promoters. However, strikingly, IRF5 does not bind the VREs, suggesting evolutionary selection against IRF5 homodimer binding. Mutational analysis reveals a critical specificity-determining residue that inhibits IRF5 binding to the ISRE-variants present in the IFN gene promoters. Integrating PBM and reporter gene data we find that both DNA-binding affinity and affinity-independent mechanisms determine the function of DNA-bound IRF dimers, suggesting that DNA-based allostery plays a role in IRF binding site function. Our results provide new insights into the role and limitations of DNA-binding affinity in delineating IRF3/5/7-specific gene expression.

Insights from resolving protein-DNA interactions at near base-pair resolution

One of the central goals in molecular biology is to understand how cell-type-specific expression patterns arise through selective recruitment of RNA polymerase II (Pol II) to a subset of gene promoters. Pol II needs to be recruited to a precise genomic position at the proper time to produce messenger RNA from a DNA template. Ostensibly, transcription is a relatively simple cellular process; yet, experimentally measuring and then understanding the combinatorial possibilities of transcriptional regulators remain a daunting task. Since its introduction in 1985, chromatin immunoprecipitation (ChIP) has remained a key tool for investigating protein-DNA contacts in vivo. Over 30 years of intensive research using ChIP have provided numerous insights into mechanisms of gene regulation. As functional genomic technologies improve, they present new opportunities to address key biological questions. ChIP-exo is a refined version of ChIP-seq that significantly reduces background signal, while providing near base-pair mapping resolution for protein-DNA interactions. This review discusses the evolution of the ChIP assay over the years; the methodological differences between ChIP-seq, ChIP-exo and ChIP-nexus; and highlight new insights into epigenetic and transcriptional mechanisms that were uniquely enabled with the near base-pair resolution of ChIP-exo.

Neurospora crassa developmental control mediated by the FLB-3 transcription factor

Here, we report that the Neurospora crassa FLB-3 protein, the ortholog of the Aspergillus nidulans FlbC transcription factor, is required for developmental control. Deletion of flb-3 leads to changes in hyphae morphology and affects sexual and asexual development. We identified, as putative FLB-3 targets, the N. crassa aba-1, wet-1 and vos-1 genes, orthologs of the ones involved in A. nidulans asexual development and that work downstream of FlbC (abaA, wetA and vosA). In N. crassa, these three genes require FLB-3 for proper expression; however, they appear not to be required for normal development, as demonstrated by gene expression analyses during vegetative growth and asexual development. Moreover, mutant strains in the three genes conidiate well and produce viable conidia. We also determined FLB-3 DNA-binding preferences via protein-binding microarrays (PBMs) and demonstrated by chromatin immunoprecipitation (ChIP) that FLB-3 binds the aba-1, wet-1 and vos-1 promoters. Our data support an important role for FLB-3 in N. crassa development and highlight differences between the regulatory pathways controlled by this transcription factor in different fungal species.

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.

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.